Search NASASearch

SEARCH · Search NASA

Results for “EM simulation”

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 145 records · Page 8

Relativistic gas accretion onto supermassive black hole binaries from inspiral through merger

Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and electromagnetic (EM) radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to informing and guiding upcoming efforts aimed at detecting gravitational waves produced by these binaries. To this end, accurate numerical modeling is required to describe both the spacetime and the magnetized gas around the black holes. In this paper, we present two key advances in this field of research. First, we have developed a novel 3D general relativistic magnetohydrodynamics (GRMHD) framework that combines multiple numerical codes to simulate the inspiral and merger of supermassive black hole binaries starting from realistic initial data and running all the way through merger. Throughout the evolution, we adopt a simple but functional prescription to account for gas cooling through photon emission. Next, we have applied our new computational method to follow the time evolution of a circular, equal-mass, nonspinning black hole binary for ∼200 orbits, starting from a separation of 20⁢𝑟 𝑔 and reaching the postmerger evolutionary stage of the system. We have shown how mass continues to flow toward the binary even after the binary “decouples” from its surrounding disk, but the accretion rate onto the black holes diminishes. We have identified how the minidisks orbiting each black hole are slowly drained and eventually dissolve as the binary compresses. We confirm previous findings that the system’s luminosity decreases by a factor of a few during inspiral; however, we observe an abrupt increase by ∼50% in this quantity at the time of merger, likely accompanied by an equally abrupt change in spectrum. Lastly, we have demonstrated that during the inspiral, fluid ram pressure regulates the fraction of the magnetic flux transported to the binary that attaches to the black holes’ horizons.

Accretion disk & black-hole plasma

Development of a Robust Reference Electrode in Aggressive Chemical and Radiation Environments in the Hanford Waste Tanks

The Hanford site stores more than 200 million liters of radioactive and chemically hazardous wastes from the production of weapons materials. The wastes are stored in 177 underground carbon-steel storage tanks, separated between 149 single shell tanks (SSTs) and 28 double shell tanks (DSTs). The DSTs provide critical retrieval and interim storage before the waste is vitrified in the Waste Treatment and Isolation Plant (WTP). The tanks will need to remain in-service far beyond the initial 40-year design life, and effective corrosion control practices must remain in force to extend the tanks’ lifespans. This effort includes direct measurements of corrosion rate (e.g., ultrasonic measurements and corrosion coupons) and electrochemical processes (e.g., linear polarization measurements and open circuit potential measurements). The Hanford site began monitoring the corrosion potential in select DSTs in 2008. Of the 45 reference electrodes that have been installed, 29 have failed and 6 others provided unreliable results. DOE-EM is supporting a 3-year program to develop a chemical and radiation resistant reference electrode for application in the Hanford tanks. The first year of the program focused on understanding the failure mechanism for the reference electrodes and identification of candidate construction materials that would mitigate degradation of the electrodes in the waste environment. During the second year of the program, the objectives were to: 1) test candidate materials under simulated waste conditions, 2) design components that will extend the service life of the electrode, 3) fabricate materials for prototype reference electrodes, and 4) assemble prototype reference electrodes for accelerated testing. The reference electrode is constructed of four principal parts: 1) junction, 2) casing, 3) inner chamber backfill materials, and 4) the sensing wire. Principally, improvements of the junction, casing, and inner chamber backfill materials are being pursued. The junction material at the interface between the waste and the inner chamber of the reference electrode was identified as a critical component in the failure of the reference electrodes. Nine candidate replacement junction materials were tested under simulated waste conditions to evaluate permeation rate. These materials included a variety of polymeric and ceramic materials, some of which were 3-D printed. Thus far, porous polyvinylidene fluoride materials have performed satisfactorily and are being considered for prototype development. The commercial electrode casing materials in general have performed well. Additionally, 3-D printing of chemically and mechanically stable materials is being investigated as a means for further improvement in fabrication consistency. SRNL has also investigated altering the reference electrode design to extend the service life. The new design of the interior of the reference electrode casing creates a longer, more tortuous path between the junction material and the electrode sensing wire. A finite element model was used to optimize the design without adversely impacting the circuit resistance of the electrode during the measurements, thus preserving the measurement accuracy while enhancing the service life. The inner chamber back fill materials are also critical to the performance of the reference electrode. Materials that are resistant to intruding tank waste and provide a conductive path to the sensing wire were investigated. Gel and powder materials that are interspersed with a conductive chloride bearing material were tested for their influence on diffusion and electrode resistance. All the investigated materials and components will be assembled, with collaboration from commercial vendors, to fabricate the initial prototypes. Accelerated testing of the prototypes will be initiated in Year 2 of the program and will be completed in Year 3. A recommendation on the materials of construction and the design of the new robust reference electrode will be presented to the Hanford tank farm facility.

Sykes, Kiana [Savannah River National Laboratory (

Heat Pipes and Heat Rejection Component Testing at NASA Glenn Research Center

Titanium-water heat pipes are being evaluated for use in the heat rejection system for space fission power systems. The heat rejection syst em currently comprises heat pipes with a graphite saddle and a composite fin. The heat input is a pumped water loop from the cooling of the power conversion system. The National Aeronautics and Space Administration has been life testing titanium-water heat pipes as well as eval uating several heat pipe radiator designs. The testing includes thermal modeling and verification of model, material compatibility, frozen startup of heat pipe radiators, and simulating low-gravity environments. Future thermal testing of titanium-water heat pipes includes low-g ravity testing of thermosyphons, radiation testing of heat pipes and fin materials, water pump performance testing, as well as Small Busine ss Innovation Research funded deliverable prototype radiator panels.

Sanzi, James L.

CHANGO: A Software Tool for Boost Stage Guidance of the Space Launch System Exploration Mission 1

The Space Launch System (SLS) Exploration Mission 1 (EM-1) test flight will use open-loop guidance for Boost Stage (BS) flight. A table of attitude commands as a function of altitude, called the chi table, will be loaded onto the flight computers. The chi table will be generated using the measured winds on launch day by the Chi Angle Optimizer (CHANGO) software tool. Details of CHANGO’s design are given, including a Three Degrees-of-Freedom (3-DOF) simulation and a numerical minimization routine. CHANGO’s use in launch day operations is also described.

Ahmad, Naeem

Parameter Design and Optimal Control of an Open Core Flywheel Energy Storage System

In low earth orbit (LEO) satellite applications spacecraft power is provided by photovoltaic cells and batteries. To overcome battery shortcomings the University of Maryland, working in cooperation with NASA/GSFC and NASA/LeRC, has developed a magnetically suspended flywheel for energy storage applications. The system is referred to as an Open Core Composite Flywheel (OCCF) energy storage system. Successful application of flywheel energy storage requires integration of several technologies, viz. bearings, rotor design, motor/generator, power conditioning, and system control. In this paper we present a parameter design method which has been developed for analyzing the linear SISO model of the magnetic bearing controller for the OCCF. The objective of this continued research is to principally analyze the magnetic bearing system for nonlinear effects in order to increase the region of stability, as determined by high speed and large air gap control. This is achieved by four tasks: (1) physical modeling, design, prototyping, and testing of an improved magnetically suspended flywheel energy storage system, (2) identification of problems that limit performance and their corresponding solutions, (3) development of a design methodology for magnetic bearings, and (4) design of an optimal controller for future high speed applications. Both nonlinear SISO and MIMO models of the magnetic system were built to study limit cycle oscillations and power amplifier saturation phenomenon observed in experiments. The nonlinear models include the inductance of EM coils, the power amplifier saturation, and the physical limitation of the flywheel movement as discussed earlier. The control program EASY5 is used to study the nonlinear SISO and MIMO models. Our results have shown that the characteristics and frequency responses of the magnetic bearing system obtained from modeling are comparable to those obtained experimentally. Although magnetic saturation is shown in the bearings, there are good correlations between the theoretical model and experimental data. Both simulation and experiment confirm large variations of the magnetic bearing characteristics due to air gap growth. Therefore, the gap growth effect should be considered in the magnetic bearing system design. Additionally, the magnetic bearing control system will be compared to other design methods using not only parameter design but H-infinity optimal control and mu synthesis.

Pang, D.

Overview of Recent EAST Testing, Modeling & Analysis

Shock Layer Radiation at NASA Ames - Background: Complex aerothermal and thermochemical phenomena of planetary entry define convective and radiative heating. A spacecraft's TPS (Thermal Protective System) mitigates heat transfer to substructure. Successful TPS design relies on verifiable characterization of these phenomena in the anticipated flight environment. - Approach: EAST (Electric Arc Shock Tube) [facility at Ames] simulates high-enthalpy, real-gas phenomena encountered by hypersonic vehicles entering planetary atmospheres by spectrally imaging the flow behind a moving shock wave. - Goal: Validate aerothermal models (DPLR (Data Parallel Line Relaxation Code) and NEQAIR (Nonequilibrium Radiative Transport and Spectra Program)), inform model improvements, reduce uncertainty and quantify design uncertainties. - Recent Relevant Projects: MSL (Mars Science Laboratory) & Mars 2020, InSight, OSIRIS-REx (Origins Spectral Interpretation Resource Identification Security Regolith Explorer), Orion EFT-1 (Exploration Flight-1) and EM-1 (Exploration Mission-1) and New Frontiers.

Brandis, Aaron

Estimating Sparse Direct Effects in Multivariate Regression With the Spike-and-Slab LASSO

The multivariate regression interpretation of the Gaussian chain graph model simultaneously parametrizes (i) the direct effects of p predictors on q outcomes and (ii) the residual partial covariances between pairs of outcomes. We introduce a new method for fitting sparse versions of these models with spike-and-slab LASSO (SSL) priors. We develop an Expectation Conditional Maximization algorithm to obtain sparse estimates of the p × q matrix of direct effects and the q × q residual precision matrix. Our algorithm iteratively solves a sequence of penalized maximum likelihood problems with self-adaptive penalties that gradually filter out negligible regression coefficients and partial covariances. Because it adaptively penalizes individual model parameters, our method is seen to outperform fixed-penalty competitors on simulated data. We establish the posterior contraction rate for our model, buttressing our method’s excellent empirical performance with strong theoretical guarantees. Using our method, we estimated the direct effects of diet and residence type on the composition of the gut microbiome of elderly adults.

EM algorithm

The Transition of High-Resolution NASA MODIS Sea Surface Temperatures into the WRF Environmental Modeling System

The NASA Short-term Prediction Research and Transition (SPoRT) Center has developed a Moderate Resolution Imaging Spectroradiometer (MODIS) sea surface temperature (SST) composite at 2-km resolution that has been implemented in version 3 of the National Weather Service (NWS) Weather Research and Forecasting (WRF) Environmental Modeling System (EMS). The WRF EMS is a complete, full physics numerical weather prediction package that incorporates dynamical cores from both the Advanced Research WRF (ARW) and the Non-hydrostatic Mesoscale Model (NMM). The installation, configuration, and execution of either the ARW or NMM models is greatly simplified by the WRF EMS to encourage its use by NWS Weather Forecast Offices (WFOs) and the university community. The WRF EMS is easy to run on most Linux workstations and clusters without the need for compilers. Version 3 of the WRF EMS contains the most recent public release of the WRF-NMM and ARW modeling system (version 3 of the ARW is described in Skamarock et al. 2008), the WRF Pre-processing System (WPS) utilities, and the WRF Post-Processing program. The system is developed and maintained by the NWS National Science Operations Officer Science and Training Resource Coordinator. To initialize the WRF EMS with high-resolution MODIS SSTs, SPoRT developed the composite product consisting of MODIS SSTs over oceans and large lakes with the NCEP Real-Time Global (RTG) filling data over land points. Filling the land points is required due to minor inconsistencies between the WRF land-sea mask and that used to generate the MODIS SST composites. This methodology ensures a continuous field that adequately initializes all appropriate arrays in WRF. MODIS composites covering the Gulf of Mexico, western Atlantic Ocean and the Caribbean are generated daily at 0400, 0700, 1600, and 1900 UTC corresponding to overpass times of the NASA Aqua and Terra polar orbiting satellites. The MODIS SST product is output in gridded binary-1 (GRIB-1) data format for a seamless incorporation into WRF via the WPS utilities. The full-resolution, 1-km MODIS product is sub-sampled to 2-km grid spacing due to limitations in handling very large dimensions in the GRIB-1 data format. The GRIB-1 files are posted online at ftp://ftp.nsstc.org/sstcomp/WRF/, which is directly accessed by the WRF EMS scripts. The MODIS SST composites are also downloaded to the EMS data server, which is accessible by the WRF EMS users and NWS WFOs. The SPoRT MODIS SST composite provides the model with superior detail of the ocean gradients around Florida and surrounding waters, whereas the operational RTG SST typically depicts a relatively smooth field and is not able to capture sharp horizontal gradients in SST. Differences of 2-3 C are common over small horizontal distances, leading to enhanced SST gradients on either side of the Gulf Stream and along the edges of the cooler shelf waters. These sharper gradients can in turn produce atmospheric responses in simulated temperature and wind fields as depicted in LaCasse et al. Differences in atmospheric verification statistics over a several month study were generally small in the vicinity of south Florida; however, the validation of SSTs at specific buoy locations revealed important improvements in the biases and RMS errors, especially in the vicinity of the cooler shelf waters off the east-central Florida coast. A current weakness in the MODIS SST product is the occurrence of occasional discontinuities caused by high latency in SST coverage due to persistent cloud cover. An enhanced method developed by Jedlovec et al. (2009, GHRSST User Symposium) reduces the occurrence of these problems by adding Advanced Microwave Scanning Radiometer -- EOS (AMSR-E) SST data to the compositing process. Enhanced SST composites are produced over the ocean regions surrounding the Continental U.S. at four times each day corresponding to Terra and Aqua equator crossing times. For a given day and overpass time, both MODInd AMSR-E data from the previous seven days form a collection used in the compositing. At each MODIS pixel, cloud-free SST values from the collection are used to form a weighted average based on their latency (number of days from the current day). In this way, recent SST data are given more weight than older data. One of the primary issues involved in incorporating the AMSR-E microwave data in the composites is the tradeoff between the decreased spatial resolution of the AMSR-E data (25 km) and the increased coverage due to its near all-weather capability. Currently, the AMSR-E is given a weight of 20% compared to MODIS data, thereby preserving the spatial structure observed in the MODIS data. Day-time (night-time) AMSR-E SST data from Aqua are used with both Terra and Aqua MODIS day-time (night-time) SST data sets.

Case, Jonathan L.

Discrete-Event Model of WIPP Operations

The Waste Isolation Pilot Plant (WIPP) is the critical component of the Department of Energy's (DOE) Transuranic Radioactive Waste (TRU) disposition infrastructure. Quantifying the operations ofWIPP with a discrete-event model demonstrates the capability to assess efficiency, identity bottlenecks, and improve future operations. Such a model has been developed with the simulation software ExtendSim. This report outlines the structure of that model, summarizes the model's successful reproduction of the annual amount of emplaced waste, and demonstrates that WIPP is successfully receiving and emplacing waste at a rate consistent with the rate at which the waste arrives. The model serves as a first step toward future production enhancements at WIPP. Those enhancements will rely on close collaboration with the Carlsbad Field Office (CBFO), accurate interpretation and incorporation of the model results, and effective planning with other DOE Environmental Management (DOE-EM) entities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

Test Report: Direct and Indirect Lightning Effects on Composite Materials

Lightning tests were performed on composite materials as a part of an investigation of electromagnetic effects on the materials. Samples were subjected to direct and remote simulated lightning strikes. Samples included various thicknesses of graphite filament reinforced plastic (GFRP), material enhanced by expanded aluminum foil layers, and material with an aluminum honeycomb core. Shielding properties of the material and damage to the sample surfaces and joints were investigated. Adding expanded aluminum foil layers and increasing the thickness of GFRP improves the shielding effectiveness against lightning induced fields and the ability to withstand lightning strikes. A report describing the lightning strike tests performed by the U.S. Army Redstone Technical Test Center, Redstone Arsenal, AL, STERT-TE-E-EM, is included as an appendix.

Composite materials

Impact of Lake Okeechobee Sea Surface Temperatures on Numerical Predictions of Summertime Convective Systems over South Florida

The NASA Short-term Prediction Research and Transition (SPoRT) Center, the Florida Institute of Technology, and the NOAA/NWS Weather Forecast Office at Miami, FL (MFL) are collaborating on a project to investigate the impact of using high-resolution, 2-km Moderate Resolution Imaging Spectroradiometer (MODIS) sea surface temperature (SST) composites within the Weather Research and Forecasting (WRF) prediction system. The NWS MFL is currently running WRF in real-time to support daily forecast operations, using the National Centers for Environmental Prediction Nonhydrostatic Mesoscale Model dynamical core within the NWS Science and Training Resource Center's Environmental Modeling System (EMS) software. Twenty-seven hour forecasts are run daily initialized at 0300, 0900, 1500, and 2100 UTC on a domain with 4-km grid spacing covering the southern half of Florida and adjacent waters of the Gulf of Mexico and Atlantic Ocean. The SSTs are initialized with the NCEP Real-Time Global (RTG) analyses at 1/12deg resolution. The project objective is to determine whether more accurate specification of the lower-boundary forcing over water using the MODIS SST composites within the 4-km WRF runs will result in improved sea fluxes and hence, more accurate e\~olutiono f coastal mesoscale circulations and the associated sensible weather elements. SPoRT conducted parallel WRF EMS runs from February to August 2007 identical to the operational runs at NWS MFL except for the use of MODIS SST composites in place of the RTG product as the initial and boundary conditions over water. During the course of this evaluation, an intriguing case was examined from 6 May 2007, in which lake breezes and convection around Lake Okeechobee evolved quite differently when using the high-resolution SPoRT MODIS SST composites versus the lower-resolution RTG SSTs. This paper will analyze the differences in the 6 May simulations, as well as examine other cases from the summer 2007 in which the WRF-simulated Lake Okeechobee breezes evolved differently due to the SST initialization. The effects on wind fields and precipitation systems will be emphasized, including validation against surface mesonet observations and Stage IV precipitation grids.

Case, Jonathan L.

Nonsolvating Fluoroaromatic Cosolvent Enabled Long-Term Cycling of High-Voltage Lithium-Ion Batteries with Organosulfur Electrolytes

Here, the structure–activity relationships of nonsolvating cosolvents for organosulfur-based electrolyte systems were revealed. The performance of nonsolvating dilutant fluorobenzene (FB) was compared to various fluorinated ether dilutants in high-voltage electrolytes containing a concentration of 1.2 M LiPF 6 dissolved in fluoroethylene carbonate (FEC), ethyl methyl sulfone (EMS), and the dilutant. In a high-voltage and high-loading LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) full cell configuration, the organosulfur-based electrolyte containing FB dilutant enabled superior electrochemical performance compared to the electrolytes using other nonsolvating fluorinated ether formulations. Moreover, the FB-containing electrolyte exhibited the highest ionic conductivity and lowest viscosity among all organosulfur-based electrolytes containing nonsolvating dilutant. These improvements are attributed to the enhanced physical properties of electrolyte and lithium-ion mobility. Furthermore, by employing first-principles simulations, the observed suppression of side reactions at high voltage is linked to FB’s lower reactivity toward singlet dioxygen, which is likely produced at the NMC interface. Overall, FB is considered an excellent diluent that does not impede cell operation by mass decomposition at the cathode.

25 ENERGY STORAGE

Quantitative Modeling of High-Energy Electron Scattering in Thick Samples Using Monte Carlo Techniques

Cryo-electron microscopy (cryo-EM) is a powerful tool for imaging biological samples but is typically limited by sample thickness, which is restricted to a few hundred nanometers depending on the electron energy. However, there is a growing need for imaging techniques capable of studying biological samples up to 10 µm in thickness while maintaining nanoscale resolution. This need motivates the use of mega-electron-volt scanning transmission electron microscopy (MeV-STEM), which leverages the high penetration power of MeV electrons to generate high-resolution images of thicker samples. In this study, we employ Monte Carlo simulations to model electron–sample interactions and explore the signal decay of imaging electrons through thick specimens. By incorporating material properties, interaction cross-sections for energy loss, and experimental parameters, we investigate the relationship between the incident and transmitted beam intensities. Key factors such as detector collection angle, convergence semi-angle, and the material properties of samples were analyzed. Our results demonstrate that the relationship between incident and transmitted beam intensities follows the Beer–Lambert law over thicknesses ranging from a few microns to several tens of microns, depending on material composition, electron energy, and collection angles. The linear depth of silicon dioxide reaches 3.9 µm at 3 MeV, about 6 times higher than that at 300 keV. Meanwhile, the linear depth of amorphous ice reaches 17.9 µm at 3 MeV, approximately 11.5 times higher than that at 300 keV. These findings are crucial for advancing the study of thick biological and semiconductor samples using MeV-STEM.

36 MATERIALS SCIENCE

Q1 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q1: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $/circ$ Initially run CESOL with SOLPS and map heat flux using simple HEAT method: ▪ Generate SOLPS grid for CAT reference case, ▪ Implement and apply ‘lower’ fidelity method of using HEAT-like analytic method to map SOLPS heat flux from charged particles to first wall, and ▪ Provide initial heat flux + neutron loading for further engineering analysis. 2. Demonstrate multiphysics magnet analysis: $/circ$ Demonstrate validation of the Elmer workflow by comparing induced stresses due to EM forces generated by TF coils using multiple codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

High-Propellant Throughput Sub-kW Electric Propulsion System for Deep Space Science and Exploration

The National Aeronautics and Space Administration (NASA) is maturing high-propellant throughput sub-kilowatt electric propulsion technologies to enable small spacecraft deep space science and exploration missions with high delta-v requirements. The pathfinder model (PM) propulsion system consists of the H71M-PM Hall-effect thruster, a breadboard 1-kW power processing unit (PPU), and a propellant flow control system. The propulsion system requirements balance the needs of various high delta-v NASA and commercial industry mission concepts to achieve a design that both enables a variety of NASA small spacecraft deep space missions, while remaining viable for select commercial applications. The H71M-PM thruster has completed performance characterization and three 500-h short duration wear tests (SDWT). The propulsion system provides stable thrust generation over a wide range of operating conditions from 200 W to 1 kW, and 200 V to 400 V. The thruster has demonstrated a thrust as high as 68 mN at 300 V and 1 kW. The thruster has similarly demonstrated a specific impulse of 1850 s at 400 V and 1 kW. Key surfaces were machined between each SDWT to simulate accelerated discharge channel and pole cover erosion. Profilometry scans across masked pole cover surfaces were conducted to determine erosion rates. SDWT results support that a target thruster lifetime of 14 kh with 50% margin is achievable. Component testing has demonstrated propellant azimuthal flow uniformity better than ± 2 percent of the nominal value, azimuthal magnetic field uniformity better than ± 0.5 percent of the nominal value, and cathode heater cycle testing to greater than 30,000 cycles. Propulsion system integrated system testing is planned to use the H71M-PM and a breadboard 1-kW PPU. Pathfinder model test results are now supporting the design of the H71M-EM engineering model thruster. A second-generation breadboard PPU has been fabricated and is currently under test. NASA has made these technologies available to U.S. industry through a no cost, nonexclusive license agreement.

Hall

Control and Scaling Approach for the Emulation of Dynamic Subscale Torque Loads

Research and development of electrified aircraft propulsion powertrains are relying on the use of electromechanical systems to emulate turbomachinery/rotor loads. Replacing a physical turbomachinery/rotor with a model driving a subscale electromechanical system capable of emulating subscale torque loads and responses is a lower risk, lower cost alternative to using the full-scale turbomachinery/rotor for initial control system verification. This paper outlines a novel control and scaling approach for emulating dynamic subscale torque loads using electric machine (EM) hardware for electrified aircraft propulsion research and development purposes. The approach, known as the Sliding Mode Impedance Controller with Scaling (SMICS), drives a mechanically coupled, two-EM system to behave like a subscale hardware representation of a hybrid-electric turbomachinery shaft. One EM reflects the inertial dynamics and torque load of the subscale turbomachinery under steady-state and transient operation while the second EM represents a motor/generator connected to the shaft, which is intended to hybridize the turbomachinery. This closed loop control and scaling algorithm applies impedance and sliding mode control schemes, along with parameter scaling, to match subscale, desired dynamics in real-time and to allow this system to be driven by a full-scale hybrid-electric turbomachinery model and control. The paper elaborates on the concept of the closed loop control and scaling approach and explains the significance of using impedance and sliding mode control. It shows a derivation of the closed loop control and scaling algorithm, its implementation, and presents a comparison of theoretical and actual simulation results acquired during hardware-in-the-loop testing of a partial turboelectric propulsion concept aircraft at the NASA Electric Aircraft Testbed (NEAT).The results show that the intended dynamic responses of the hardware and the aircraft model are achieved in both the time and frequency domain. Full scale propulsion control systems can be tested using this hardware and software approach.

Emulation

SympGNNs: Symplectic Graph Neural Networks for identifying high-dimensional Hamiltonian systems and node classification

Existing neural network models to learn Hamiltonian systems, such as SympNets, although accurate in low-dimensions, struggle to learn the correct dynamics for high-dimensional many-body systems. Herein, we introduce Symplectic Graph Neural Networks (SympGNNs) that can effectively handle system identification in high-dimensional Hamiltonian systems, as well as node classification. SympGNNs combine symplectic maps with permutation equivariance, a property of graph neural networks. Specifically, we propose two variants of SympGNNs: (i) G-SympGNN and (ii) LA-SympGNN, arising from different parameterizations of the kinetic and potential energy. We demonstrate the capabilities of SympGNN on two physical examples: a 40-particle coupled Harmonic oscillator, and a 2000-particle molecular dynamics simulation in a two-dimensional Lennard-Jones potential. Furthermore, we demonstrate the performance of SympGNN in the node classification task, achieving accuracy comparable to the state-of-the-art. Finally, we also empirically show that SympGNN can overcome the oversmoothing and heterophily problems, two key challenges in the field of graph neural networks.

Deep learning

Geochemical Insights Into Volcanic and Lithospheric Evolution of Mars

Introduction: The Martian lithosphere plays a pivotal role in volcanic processes, altering surface features, climate shifts, and the planet's early habitability. Crucially, the rigid lithosphere significantly influences Mars' long-term thermal conditions. Studying various geological eons’ volcanic compositions unveils how the lithosphere has evolved, even though understanding the early Noachian volcanic chemistry remains challenging due to weathering and resurfacing. In this study, newly discovered Noachian volcanic terranes[1-2] along with Hesperian and Amazonian volcanic terranes[3] are considered(Figure 1), to infer the evolution of the Martian lithosphere. Figure 1:Mars topographic map derived from MOLA-HRSC data [4] highlights the study regions including the Noachian volcanic terranes, Arabia Terra (AT), Thaumasia Minor (TM), and the Southwest Margin of Icaria Planum (SWIP)[5].Additionally, the Hesperian volcanic terranes such as Syrtis Major (SM), Hesperia Planum (HP), and, the Amazonian volcanic terranes encompass Elysium Mons (EM), Alba Patera (AP), and Olympus Mons (OM) are marked. Methods and Datasets: We use the most recent elemental mass fraction maps derived from gamma spectroscopy data from the Mars Odyssey 2001mission [6]to ascertain the bulk composition of the Martian landscape. Utilizing the updated geochemical provinces dataset of Mars[7], we consider bulk chemistry variations within our study regions. The maps, at a resolution of 50×50, encompass elements like Al, Ca, Fe, Si, K, Th, H2O, Cl, and S, detected through Gamma-Ray Spectroscopy (GRS) with decimeter scale depth sensitivity[6].Additional major and minor elements like Mg, Na, Ti, P, and Mn, using mass balance methods [8]serve as inputs for melting and crystallization simulations. We employ the pMELTS model with FMQ-3to FMQ oxygen fugacity to estimate the Pressure(P)and melting degree(F)in various eons of volcanic terranes. We have also used conventional thermobarometric calculations to determine their formation P-T conditions, using silica activity and olivine-melt Mg-exchange thermometry[9-10]. Results: Chemical Composition of Studied Volcanic Terranes: To assess formation conditions, we have investigated if GRS-measured compositions signify primary igneous processes by evaluating weathering through the Chemical Index of Alteration (CIA), K/Th ratios, and ternary plots for aqueous alteration. Consistent K/Th ratios and low CIA (<50) across volcanic regions suggest minimal large-scale weathering, affirming well-preserved igneous material extending at decimeter depths. In addition, the calculated bulk compositions align with Martian mantle equilibrium, confirming studied volcanic terranes-Arabia Terra, Thaumasia Minor, and SWIP -as primary or less altered compositions, reinforcing their representativeness at a regional GRS scale on Mars. We have also estimated the bulk compositions of Hesperian and Amazonian volcanic terranes shown in Figure 1. All volcanic terranes represent the primary or minimally altered composition. Temporal Evolution of Martian Volcanic Terranes: Our results show P-T conditions ranging from 1.3-1.6 GPa and 1350-1390°C for Noachian terranes, while Hesperian volcanic terranes vary from 1.6-1.7 GPa and 1370-1395°C, and Amazonian volcanic terranes range between 1.9-2.8 GPa and 1380-1415°C. Our analyses correspond to lithospheric thicknesses/ depth of melting for Noachian which ranges from110-135 km, 120-235 km for Hesperian terranes, and 160-235 km for Amazonian terranes. Partial melting percentages (F) range from 8-12% for Noachian, 10-11% for Hesperian, and 10-12% for Amazonian. We have also calculated mantle potential temperatures (Tp) to show variations across epochs(Figure 2). Heat flux estimation shows variations from 51-65 mW/m2 for Noachian, 38-45 mW/m2for Hesperian, and 27-39 mW/m2for Amazonian terranes, implying temporal changes in lithospheric thickness and heat flow. Comparison of Noachian volcanic terranes with Hesperian and Amazonian reveals variations in lithospheric thickness and heat flux, indicating potential implications for Martian surface conditions, including volcanic activity, climate evolution, and early habitability. Figure 2: P-T diagram displays formation conditions of studied Martian volcanic regions using GRS data, alongside Noachian-age surface basalt samples from Gusev, Gale, and Jezero. Dashed lines represent Martian mantle potential temperatures, calculated by adding the latent heat of fusion to the equilibrium temperature at which melt and solid mantle coexist and subtracting the Martian adiabatic gradient corresponding to the depth of melt. The solid line shows the Martian mantle solidus [11].Discussion& Implication: The study explores Martian volcanic activity across eons. Uniform lithospheric thickness from the Noachian to Hesperian eons suggests persistent weak plumes until Hesperian, influencing heat flow variations. Noachian's high heat flux, likely due to concentrated radioactive elements, shows spatial heterogeneity, indicating early Mars' differentiated mantle. Volcanic outgassing plays avital role in the formation of the Martian climate, likely sustained water on the surface, influencing warm-wet or cold-dry conditions. Hesperian's favorable lithospheric heat flow (45 mW/m2) supports extended water existence, resembling Earth's conditions. Whereas, Amazonian's diminished heat flux (25 mW/m2) corresponds to water loss and arid climate. Hesperian's similar lithospheric thickness to mid-Noachian implies prolonged water retention, influencing climate and potential for astrobiological research, expanding perspectives for future Martian missions. Future work: As volcanism was active throughout its history. The emergence of Mars' extensive volcanic regions, known as large igneous provinces, which are linked to intermittent volcanic activity driven by mantle plumes, poses challenges to understanding the scale and persistence of these phenomena in Martian convection models. Therefore, understanding the regional-scale changes in eruptive processes within Martian volcanic provinces over time remains unclear and crucial for deciphering the evolution of the Martian interior without Earth-like plate tectonics. We will continue this work with petrological and thermoelastic analyses to study the compositional and spatiotemporal evolution of the Elysium Volcanic Province.

A Rani