Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mars 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Role and Lifetime of Dissociative Heterogeneous Processes in Improving Simulated Ozone on Mars

Abstract Ozone simulated in Mars Global Climate Models (MGCMs) is used to assess the underlying chemistry occurring in the atmosphere. Currently, ozone total column abundance (TCA) is under‐predicted in MGCMs by up to 120%, implying missing or inaccurate chemistry in models. Heterogeneous reactions of hydroxyl radicals (HO X ) have been offered as an explanation for some of this bias, because they cause ozone to increase at locations where it's currently under‐predicted. We use four simulations to compare modeled ozone TCA with observations from the UVIS spectrometer aboard the ExoMars Trace Gas Orbiter to improve the representation of heterogeneous processes and their impact on ozone. We use a gas‐phase only run, a dissociative scheme, an adsorbed HO X retention scheme, and a hybrid scheme that combines the dissociative mechanism with the retention of HO X on water ice. We find retention of HO X is dependent on water ice sublimation, and ozone abundance increases when water ice persists for longer periods (1–20 sols). Over time, the loss of HO X causes a depletion in H 2 O 2 concentration (HO X reservoir), and thus allows ozone concentration to increase. When adsorbed HO X are desorbed and dissociate into other by‐products, HO X are not immediately available to destroy ozone. This results in larger ozone concentrations than if desorbed HO X are released directly back into their gaseous states. When using the hybrid scheme, ozone TCA is increased up to 50% where the ozone deficit is greatest, demonstrating the best agreement with observations, and implying that HO X radicals are both retained when adsorbed and dissociate.

Brown, M. A. J.↗
DOE OSTI2024

A Neural Network-Based Power Mismatch Elimination Strategy for Integrated Solar and ESS AC/DC Systems (MARS)

The multiport autonomous reconfigurable solar power plant (MARS) is a promising concept for the integration of photovoltaic (PV) and energy storage system (ESS) to the transmission ac grid and a high-voltage direct current (HVdc) link. The presence of PV and ESS in each arm of the MARS results in uneven distribution of active power among different submodules (SMs), thereby leading to unbalanced SM capacitor voltages and potentially compromising the system stability. Moreover, in the case of partial shadings, shaded PV SMs will suffer from decreased power injections causing power mismatch in the MARS system. To address this issue, a neural-network-based power mismatch elimination (NNPME) strategy is proposed in this article. The proposed NNPME strategy optimizes ESS usage and leverages both dc and ac circulating currents to facilitate power transfer among the SMs, arms, and phases of the MARS system. Simulation and control hardware-in-the-loop (cHIL) experiments demonstrate the effectiveness of the proposed NNPME strategy. Compared with the traditional approaches, the proposed NNPME strategy can significantly enhance system efficiency and ensure stable and continuous operation, even in the presence of uneven power distribution within the MARS system.

14 SOLAR ENERGY↗
DOE OSTI2024

Overview of recent experimental results on the EAST Tokamak

Since the last IAEA-FEC in 2021, significant progress on the development of long pulse steady state scenario and its related key physics and technologies have been achieved, including the reproducible 403 s long-pulse steady-state H-mode plasma with pure radio frequency (RF) power heating. A thousand-second time scale (~1056 s) fully non-inductive plasma with high injected energy up to 1.73 GJ has also been achieved. The EAST operational regime of high β P has been significantly extended (H 98y2 > 1.3, β P ~ 4.0, β N ~ 2.4 and n e /n GW ~ 1.0) using RF and neutral beam injection (NBI). The full edge localized mode suppression using the n = 4 resonant magnetic perturbations has been achieved in ITER-like standard type-I ELMy H-mode plasmas with q 95 ≈ 3.1 on EAST, extrapolating favorably to the ITER baseline scenario. The sustained large ELM control and stable partial detachment have been achieved with Ne seeding. The underlying physics of plasma-beta effect for error field penetration, where toroidal effect dominates, is disclosed by comparing the results in cylindrical theory and MARS-Q simulation in EAST. Breakdown and plasma initiation at low toroidal electric fields (<0.3 V m -1 ) with EC pre-ionization is developed. A beneficial role on the lower hybrid wave injection to control the tungsten concentration in the NBI discharge is observed for the first time in EAST suggesting a potential way toward steady-state H-mode NBI operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Investigations of atomic and molecular processes of NBI-heated discharges in the MAST Upgrade Super-X divertor with implications for reactors

This experimental study presents an in-depth investigation of the performance of the MAST-U Super-X divertor during NBI-heated operation (up to 2.5 MW) focussing on volumetric ion sources and sinks as well as power losses during detachment. The particle balance and power loss analysis revealed the crucial role of Molecular Activated Recombination and Dissociation (MAR and MAD) ion sinks in divertor particle and power balance, which remain pronounced in the change from ohmic to higher power (NBI heated) L-mode conditions. The importance of MAR and MAD remains with double the absorbed NBI heating. MAD results in significant power dissipation (up to ${\sim} 20\%$ of $P_\textrm{SOL}$), mostly in the cold ($T_e \lt 5$ eV) detached region. Theoretical and experimental evidence is found for the potential contribution of $D^-$ to MAR and MAD, which warrants further study. These results suggest that MAR and MAD can be relevant in higher power conditions than the ohmic conditions studied previously. Post-processing reactor-scale simulations suggests that MAR and MAD can play a significant role in divertor physics and synthetic diagnostic signals of reactor-scale devices, which are currently underestimated in exhaust simulations. This raises implications for the accuracy of reactor-scale divertor simulations of particularly tightly baffled (alternative) divertor configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Detection of visible-wavelength aurora on Mars

Mars hosts various auroral processes despite the planet’s tenuous atmosphere and lack of a global magnetic field. To date, all aurora observations have been at ultraviolet wavelengths from orbit. We describe the discovery of green visible-wavelength aurora, originating from the atomic oxygen line at 557.7 nanometers, detected with the SuperCam and Mastcam-Z instruments on the Mars 2020 Perseverance rover. Near–real-time simulations of a Mars-directed coronal mass ejection (CME) provided sufficient lead-time to schedule an observation with the rover. The emission was observed 3 days after the CME eruption, suggesting that the aurora was induced by particles accelerated by the moving shock front. To our knowledge, detection of aurora from a planetary surface other than Earth has never been reported, nor has visible aurora been observed at Mars. This detection demonstrates that auroral forecasting at Mars is possible, and that during events with higher particle precipitation, or under less dusty atmospheric conditions, aurorae will be visible to future astronauts.

Science & Technology - Other Topics↗
DOE OSTI2025

Simulating evaporative wet and dry cycles in Gale crater, Mars using thermochemical modeling techniques

The aim of this work is to provide a model‐backed hypothesis for the formation of evaporites—sulfates, borates—in Gale crater using thermochemical modeling to determine constraints on their formation. We test the hypothesis that primary evaporites required multiple wet–dry cycles to form, akin to how evaporite assemblages form on Earth. Starting with a basalt‐equilibrated Mars fluid, Mars‐relevant concentrations of B and Li were added, and then equilibrated with Gale lacustrine bedrock. We simulated the cycles of evaporation followed by groundwater recharge/dilution to establish an approximate minimum number of wet–dry cycles required to form primary evaporites. We determine that a minimum of 250 wet–dry cycles may be required to start forming primary evaporites that consist of borates and Ca‐sulfates. We estimate that ~14,250 annual cycles (~25.6 k Earth years) of wet and dry periods may form primary borates and Ca‐sulfates in Gale crater. These primary evaporites could have been remobilized during secondary diagenesis to form the veins that the Curiosity rover observes in Gale crater. No Li salts form after 14,250 cycles modeled for the Gale‐relevant scenario (approximately 10 6 cycles would be needed) which implies Li may be leftover in a groundwater brine after the time of the lake. No major deposits of borates are observed to date in Gale crater which also implies that B may be leftover in the subsequent groundwater brine that formed after evaporites were remobilized into Ca‐sulfate veins.

58 GEOSCIENCES↗
DOE OSTI2026

Toroidal phase of post-disruption runaway electron loss to wall in presence of applied 3-D fields in DIII-D

The toroidal distribution of runaway electrons (REs) striking the centre post in DIII-D post-disruption RE final loss events is simulated utilizing the MARS-F code combined with the REORBIT module, by tracing the guiding-centre drift orbits of test REs in the presence of external (applied) plus internal (intrinsic instability) 3-D fields. To better recover the experimental results, three different equilibria with different safety factor profiles are adopted. Dominant resistive magnetohydrodynamic instabilities are found to be m/n=2/1, 3/2 and 1/1 modes (where m and n are the poloidal and toroidal mode numbers). Externally applied n=1 resonant magnetic perturbations (RMPs), with both even and odd parity configuration, are considered in each equilibrium. Locking of MHD instability to the applied RMPs is simulated by minimizing the total perturbed magnetic energy in the plasma. Here, the simulated toroidal locations, where most REs strike in the presence of an edge m/n= 2/1 instability and the n=1 RMPs in even parity, are found to agree well with experiments. Modeling also captures the measured n = 2 RE loss pattern on the wall for experiments where an edge m/n=3/2 instability is dominant instead of the m/n=2/1 mode. The key overall findings include (a) the measured global RE toroidal impact pattern in experiments can be explained by a locked resistive kink instability, (b) the toroidal phase of the locked mode is well predicted by minimizing the total perturbed magnetic energy, and (c) this edge locked mode determines the toroidal phase of the peak RE wall impact.

Bai, Xue [University of California San Diego, La J↗
DOE OSTI2024

Validation of SOLPS-ITER and EDGE2D-EIRENE simulations for H, D, and T JET ITER-like wall low-confinement mode plasmas

Both experiments and simulations with SOLPS-ITER and EDGE2D-EIRENE show that the onset of detachment for the low-field side (LFS) divertor – defined here as the line-averaged upstream density (〈n e 〉 edge ) at which the plasma flux to the LFS target (I LFS-plate ) starts to decrease with increasing 〈n e 〉 edge – is independent of the isotope mass. However, there are three major simulation-experiment discrepancies: (i) the absolute values of I LFS-plate and the electron density (n e ) in the LFS divertor at the onset of detachment are significantly lower in simulations, i.e., approximately a factor of 2 for I LFS-plate and a factor of 3-4 for n e ; (ii) the degree of detachment – defined here as the difference between I LFS-plate at the onset of detachment and at an 〈n e 〉 edge value close to the density limit – is smaller in simulations compared to experiments; and (iii) the experimentally observed larger degree of detachment for D and T plasmas compared to H plasmas cannot be clearly distinguished from the simulation results. There are strong indications that discrepancy (i) is to a large extent caused by neglecting Lyman-opacity effects in our simulations. The simulations predict a similar net volumetric recombination source for all isotopes due to the fact that molecule-activated recombination (MAR) compensates for the reduced electron–ion recombination (EIR) for H, whereas MAR is negligible for D and T. This similar net volumetric recombination source for all isotopes leads to an isotope-independent degree of detachment in simulations. An analysis of the Balmer-α and Balmer-$γ$ emission confirms the underestimate of MAR in simulations (especially for D and T) for the JET metallic wall, which was previously observed for devices with a carbon wall. The underestimate of MAR is an important cause for discrepancy (ii) and the fact that there is a stronger underestimate of MAR for D and T than for H explains discrepancy (iii). Extending the plasma grid to the vessel wall increases I LFS-plate and n e at the onset of detachment by 25%, and the EIR source increases by 80% in detached conditions. Hence, while the extended grid results are closer to the experimental observations, the previously described qualitative discrepancies still persist.

Balmer emission↗
DOE OSTI2025

NB4 MARS Preliminary Optics Design

This report describes the calculations and simulations to determine the optimal guide layout, dimensions and coatings for the NB4 (MARS) beam line post beryllium reflector changeout. The reported guide performance meets the goals described in the MARS scientific requirements document (1).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Optimization Studies of Radiation Shielding for PIP-II Project at Fermilab

The Proton Improvement Plan-II (PIP-II) at Fermilab represents a significant advancement in the quest to answer some of the most profound questions about our universe using the world's most intense high-energy neutrino beam. The project requires the construction of a new addition to the Fermilab accelerator complex – an 800-MeV high-intensity superconducting linear accelerator. Ensuring the safety and regulatory compliance of this ambitious project is paramount, necessitating thorough dose rate assessments under both normal operational and accidental scenarios to align with the Fermilab Radiological Control Manual (FRCM) standards. Our approach included a shielding optimization used for the simulations with the Monte Carlo code MARS [1,2,3] to incorporate new magnet and collimator designs, essential for reflecting the current state of PIP-II infrastructure. The implementation of high-resolution detector planes, despite their computational demands, enabled us to gather detailed radiation field data crucial for optimizing shielding configurations. To overcome the significant computational demands, we developed a branching code that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through geometry splitting and the application of Russian Roulette techniques, tailored to prioritize regions of interest based on predefined importances and weight limits.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Optimization Studies of Radiation Shielding for PIP-II Project at Fermilab

The Proton Improvement Plan-II (PIP-II) at Fermilab represents a significant advancement in the quest to answer some of the most profound questions about our universe using the world's most intense high-energy neutrino beam. The project requires the construction of a new addition to the Fermilab accelerator complex – an 800-MeV high-intensity superconducting linear accelerator. Ensuring the safety and regulatory compliance of this ambitious project is paramount, necessitating thorough dose rate assessments under both normal operational and accidental scenarios to align with the Fermilab Radiological Control Manual (FRCM) standards. Our approach included a shielding optimization used for the simulations with the Monte Carlo code MARS [1,2,3] to incorporate new magnet and collimator designs, essential for reflecting the current state of PIP-II infrastructure. The implementation of high-resolution detector planes, despite their computational demands, enabled us to gather detailed radiation field data crucial for optimizing shielding configurations. To overcome the significant computational demands, we developed a branching code that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through geometry splitting and the application of Russian Roulette techniques, tailored to prioritize regions of interest based on predefined importances and weight limits.

Makovec, Alajos↗
DOE OSTI2025

Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D

Post-disruption runaway electron (RE) kinetic energy K and pitch angle sin$\vartheta$ are critical parameters for determining resulting first wall material damage during wall strikes, but are very challenging to measure experimentally. During the final loss instability, confined RE K and sin$\vartheta$ are reconstructed during center-post wall strikes for both high impurity (high-Z) and low impurity (low-Z) plasmas by combining soft x-ray, hard x-ray, synchrotron emission, and total radiated power measurements. Deconfined (wall impacting) RE sin$\vartheta$ is then reconstructed for these shots by using time-decay analysis of infra-red imaging. Additionally, deconfined RE K and sin$\vartheta$ are reconstructed for a low-Z downward loss shot by analyzing resulting damage to a sacrificial graphite dome limiter. The damage analysis uses multi-step modeling simulating plasma instability, RE loss orbits, energy deposition, and finally material expansion (MARS-F, KORC, GEANT-4, and finally COMSOL). Overall, mean kinetic energies are found to be in the range $\langle$$K$$\rangle$ $≈ 3-4$ MeV for confined REs. KORC simulations indicate that the final loss instability process does not change individual RE kinetic energy K. In conclusion, confined RE pitch angles are found to be fairly low initially pre-instability, $\langle$sin$\vartheta$$\rangle$ $≈ 0.1 – 0.2$, but appear to increase roughly 2x, to $\langle$sin$\vartheta$$\rangle$ $≈ 0.3 – 0.4$ for both confined and deconfined REs during instability onset in the low-Z case; this increase is not observed in the high-Z case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Novel Two Stage Collimation Unit for Fermilab Booster

A new two-stage collimation unit (2SC) for Fermilab Booster will be installed during 2024 summer shutdown. It is a supplementary collimator for existing single stage Booster collimators. Design details of this 2SC adapted to Booster conditions are described. Results of beam dynamics simulations on collimation efficiency and evaluation of collimator shielding with MARS code on this new system are presented. The analysis on beam dynamics shows that the system will help to collimate up to about 1% of the beam in a confined way and, prompt and residual activation meet Fermilab Radiological Control limits in current as well as during PIP-II era high intensity operation of the Booster.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Proton Cyclotron Waves and Pickup Ion Ring Distribution Instabilities Upstream of Mars

Proton cyclotron waves (PCWs) upstream of Mars are thought to be associated with the instabilities of pickup ions. The instabilities of pickup ions of ring distributions have larger maximum growth rates compared with beam distributions. However, observations revealed a notably‐reduced occurrence rate of PCWs for pickup ions of ring distributions. Linear instability analysis and a corresponding two‐dimensional particle‐in‐cell (PIC) simulation are performed to investigate the instabilities of pickup ion ring distributions upstream of Mars. Linear instability analysis indicates that a pickup ion ring distribution is unstable to the ion cyclotron, ion Bernstein, and mirror instabilities. The corresponding PIC simulation confirms the linear analysis results and further demonstrates that the pickup ions are scattered toward an isotropic shell distribution by the waves excited. Interestingly, the saturation energy of the waves is much lower than that driven by a corresponding pickup ion beam distribution, and mirror waves eventually dominate the system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Monte Carlo Simulation with CAD Interface for Calculation of 3D Maps of Residual Dose (CRADA)

Objective: To develop an easy-to-use software application to predict and mitigate radiation effects in research environment, space instruments, nuclear plants and medical facilities and help nonproliferation and national security efforts. Tech-X will develop standalone software libraries and command-line tools for ( 1) translating CAD into tessellated surfaces and tetrahedral meshes in GDML (for Geant4 and MARS 15), ROOT (for MARS 15) and HDF5 (for compact representation and for the visualization) formats, (2) healing CAD geometries to make them suitable for Monte Carlo simulations; (3) creating uniform and variable Cartesian and cylindrical meshes for detailed scoring; and ( 4) efficient Monte Carlo navigation in CAD geometries. JLAB will finish automation of simulations of residual dose in CAD geometries and integrate Tech-X software into Geant4 and MARS15. Finally, Tech-X will develop a Graphical User Interface to set up and heal CAD geometries, create input files, run and visualize simulations for residual dose. This application will run on local desktops, local and remote clusters and supercomputers and will be made available through public clouds, such as Amazon Web Services.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Optimization Studies of Radiation Shielding for the PIP-II Project at Fermilab

The PIP-II project at Fermilab, which includes an 800-MeV superconducting LINAC, demands rigorous radiation shielding optimization to meet safety requirements. We updated the MARS geometry model to reflect new magnet and collimator designs and introduced high-resolution detector planes to better capture radiation field distributions. To overcome the significant computational demands, we implemented a well-known branching technique that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through particle splitting and the application of Russian Roulette techniques. Additionally, new graphical tools were created to streamline data visualization and MARS code usability.

Makovec, Alajos [Fermilab] (ORCID:0000000286157492↗
DOE OSTI2024

Characterization of SimulCam, a standoff Raman system for scientific support of SuperCam operations on Mars

During the development activities of SuperCam Calibration Target, target intended for one of the two first Raman instruments to be deployed on another planetary body, our group developed a laboratory instrument that could simulate to some extent the Raman capabilities of one of such instruments and could provide data with similar quality. The use of this kind of laboratory instruments has demonstrated its utility in the evaluation of potential calibration targets or anticipating the science outcome that an instrument could provide. The present work describes our laboratory setup to support SuperCam, evaluating similarities between both instruments, despite of differences in the hardware. Evaluation of data gathered by SuperCam on Mars and the availability of one replica of SuperCam’s Calibration Target allowed the comparison on the same set of targets, demonstrating how similar Signal-to-Noise Ratio (SNR) could be achieved from both instruments. The higher energy per pulse on SimulCam is compensated by a greater analytical footprint and the use of smaller collection optics. The results show how spectra obtained at representative distances of SuperCam are comparable. Operational principles are also comparable in terms of time resolution, and close in terms of spectral resolution. This similarity has allowed different science support works using SimulCam data, as well as the support to Mars detections using our setup. We provide examples of this support that will be shared with the community in different papers, as well as examples of possible operations activities that could benefit from experiments performed with SimulCam. We show how this setup can complement the two laboratory replicas in Los Alamos and Toulouse in providing support data to different experiments.

47 OTHER INSTRUMENTATION↗