Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability regime”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Long-pulse high-performance H-mode plasmas achieved on EAST

A record duration of a 310 s H-mode plasma (H 98y2 ~ 1.3, n e /n GW ~ 0.7, f BS > 50%) has been recently achieved on experimental advanced superconducting tokamak (EAST) with metal walls, exploiting the device's improved long-pulse capabilities. The experiment demonstrates good control of tungsten concentration, core/edge MHD stability, and particle and heat exhaust with an ITER-like tungsten divertor and zero injected torque, establishing a milestone on the path to steady-state long-pulse high-performance scenarios in support of ITER and CFETR. Important synergistic effects are leveraged toward this result, which relies purely on radio frequency (RF) powers for heating and current drive (H&CD). On-axis electron cyclotron heating enhances the H&CD efficiency from lower hybrid wave injection, increasing confinement quality and enabling fully non-inductive operation at high density (n e /n GW ~ 70%) and high poloidal beta (β P ~ 2.5). A small-amplitude grassy edge localized mode regime facilitates the RF power coupling to the H-mode edge and reduces divertor sputtering/erosion. The high energy confinement quality (H 98y2 ~ 1.3) is achieved with the experimental and simulated results pointing to the strong effect of Shafranov shift on turbulence. Transport analysis suggests that trapped electron modes dominate in the core region during the record discharge. The detailed physics processes (RF synergy, core-edge integration, confinement properties, etc.) of the steady-state operation will be illustrated in the content. In the future, EAST will aim at accessing more relevant dimensionless parameters to develop long-pulse high-performance plasma toward ITER and CFETR steady-state advanced operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The surface of halide perovskites from nano to bulk

The surface of a semiconductor often has a key role in determining its properties. Here, for metal halide perovskites, understanding the surface features and their impact on the materials and devices is becoming increasingly important. At length scales down to the nanoscale regime, surface features become dominant in regulating the properties of perovskite materials, owing to the high surface-to-volume ratio. For perovskite bulk films in the micrometre range, defects and structural disorder readily form at the surface and affect device performance. Through concerted efforts to optimize processing techniques, high-quality perovskite thin films can now be fabricated with monolayer-like polycrystalline grains or even single crystals. Surface defects therefore remain the major obstacle to progress, pushing surface studies to the forefront of perovskite research. In this Review, we summarize and assess recent advances in the understanding of perovskite surfaces and surface strategies towards improving perovskite materials and the efficiency and stability of perovskite devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamics and resilience of the unconventional charge density wave in ScV 6 Sn 6 bilayer kagome metal

Long-range electronic ordering descending from a metallic parent state constitutes a rich playground to study the interplay of structural and electronic degrees of freedom. In this framework, kagome metals are in the most interesting regime where both phonon and electronically mediated couplings are significant. Several of these systems undergo a charge density wave transition. However, to date, the origin and the main driving force behind this charge order is elusive. Here, we use the kagome metal ScV 6 Sn 6 as a platform to investigate this problem, since it features both a kagome-derived nested Fermi surface and van-Hove singularities near the Fermi level, and a charge-ordered phase that strongly affects its physical properties. By combining time-resolved reflectivity, first principles calculations and photo-emission experiments, we identify the structural degrees of freedom to play a fundamental role in the stabilization of charge order, indicating that ScV 6 Sn 6 features an instance of charge order predominantly originating from phonons.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Nonlinear time-series-based adaptive control applications

A control design methodology based on a nonlinear time-series reference model is presented. It is indicated by highly nonlinear simulations that such designs successfully stabilize troublesome aircraft maneuvers undergoing large changes in angle of attack as well as large electric power transients due to line faults. In both applications, the nonlinear controller was significantly better than the corresponding linear adaptive controller. For the electric power network, a flexible AC transmission system with series capacitor power feedback control is studied. A bilinear autoregressive moving average reference model is identified from system data, and the feedback control is manipulated according to a desired reference state. The control is optimized according to a predictive one-step quadratic performance index. A similar algorithm is derived for control of rapid changes in aircraft angle of attack over a normally unstable flight regime. In the latter case, however, a generalization of a bilinear time-series model reference includes quadratic and cubic terms in angle of attack.

Mohler, R. R.↗
NASA NTRS2009

Reflectance Spectra Comparison of Orbital Debris, Intact Spacecraft, and Intact Rocket Bodies in the GEO Regime

A key objective of NASA s Orbital Debris program office at Johnson Space Center (JSC) is to characterize the debris environment by way of assessing the physical properties (type, mass, density, and size) of objects in orbit. Knowledge of the geosynchronous orbit (GEO) debris environment in particular can be used to determine the hazard probability at specific GEO altitudes and aid predictions of the future environment. To calculate an optical size from an intensity measurement of an object in the GEO regime, a 0.175 albedo is assumed currently. However, identification of specific material type or types could improve albedo accuracy and yield a more accurate size estimate for the debris piece. Using spectroscopy, it is possible to determine the surface materials of space objects. The study described herein used the NASA Infrared Telescope Facility (IRTF) to record spectral data in the 0.6 to 2.5 micron regime on eight catalogued space objects. For comparison, all of the objects observed were in GEO or near-GEO. The eight objects consisted of two intact spacecraft, three rocket bodies, and three catalogued debris pieces. Two of the debris pieces stemmed from Titan 3C transtage breakup and the third is from COSMOS 2054. The reflectance spectra of the Titan 3C pieces share similar slopes (increasing with wavelength) and lack any strong absorption features. The COSMOS debris spectra is flat and has no absorption features. In contrast, the intact spacecraft show classic absorption features due to solar panels with a strong band gap feature near 1 micron. The two spacecraft are spin-stabilized objects and therefore have solar panels surrounding the outer surface. Two of the three rocket bodies are inertial upper stage (IUS) rocket bodies and have similar looking spectra. The slopes flatten out near 1.5 microns with absorption features in the near-infrared that are similar to that of white paint. The third rocket body has a similar flattening of slope but with fewer features of white paint - indicating that the surface paint on the SL-12 may be different than the IUS. This study shows that the surface materials of debris appear different spectrally than intact rocket bodies and spacecraft and therefore are not believed to be solar panel material or pristine white paint. Further investigation is necessary in order to eliminate materials as possible choices for the debris pieces.

Barker, Ed↗
DOE OSTI2023

Plasma-Coupled Flow Reactor Studies of Low-Temperature Plasma Assisted Kinetics of Methanol Blended with CO2

Ignition technologies based on low-temperature plasmas (LTP) have the potential to operate next generation engines at elevated pressures and increased dilution limits promoting higher efficiencies. From a practical standpoint, research on LTP igniters has shown to enhance combustion and ignition, improve flame stability, and extend the dilution limits of combustion. All of the aforementioned gains are complemented with higher ignition efficiencies. However, the biggest challenge with incorporating this technology into engines is the knowledge gap of how exactly plasma chemistry effects can enhance the basic combustion phenomena. This coupled with the lack of validated kinetic mechanisms for plasma-combustion chemistry is the biggest obstruction to recognize efficient ignition, especially for application relevant fuels and biofuels. In order to comprehensively evaluate the effects of LTP on an oxygenated fuel specific system, this present study examines the kinetics of methanol plasma-assisted pyrolysis and oxidation using a custom-built plasma flow reactor (PFR). Experimental regimes are also further extended to understand the effects of adding CO2 to the mixture and its consequence on reaction kinetics. The PFR is installed with a dielectric-barrier discharge (DBD) configuration to induce LTP into the fuel mixture. Non-equilibrium plasmas are generated by high-voltage pulses (nearing 20 KV) administered by a plasma pulser at high-pulse repetition rates (up to 10 kHz). In order to better understand and isolate plasma chemistry and its effect on neutral chemistry, the experiments were carried out by heavily diluting reactive mixtures in nitrogen at near isothermal conditions. This suppresses the effect of exothermic reactions on chemistry allowing stable intermediates and products to be detected and quantified using ex-situ GC/MS diagnostic methods. Experiments were carried out at 0.5 atm pressure and over a wide range of temperatures from 523 K to 1203 K. Experimental results depicted the enhancement in intermediates production as well as overall lower temperatures required for complete fuel consumption in the plasma specific cases as opposed to their pure thermal counterpart. Formation of oxygenated and nitrile compounds specific to plasma assisted pyrolysis cases illustrated the efficacy of LTP to introduce new reaction pathways accelerating fuel decomposition. Increase in reactivity at lower temperatures is sought to be an effect induced by plasma chemistry postulating the acceleration in intermediates production. The onset of thermal ignition in plasma assisted oxidation is seen 200 K earlier than thermal oxidation highlighting efficient fuel conversion to final byproducts. Enhanced collisional processes afforded by LTP is seen to perturb reaction pathways between oxygenated fuel radicals and N-atoms to alter overall chemical reactivity. New insights into kinetics are gained from this study highlighting governing plasma assisted combustion (PAC) pathways for oxygenated fuel reaction chemistry. The results from this study can be used to develop future mechanisms specific to plasma chemistry which can bridge the mechanistic knowledge gap for LTP ignition so that future engines can adopt LTPs in their design for efficient combustion.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Dependence of alpha-particle-driven Alfven eigenmode linear stability on device magnetic field strength and consequences for next-generation tokamaks

Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conven- tional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfven eigenmode (AE) stability. Energetic particles, including alphas from D-T fusion, can destabilize AEs, possibly causing loss of alpha heat and damage to the device. AEs are sensitive to device magnetic field via the field dependence of resonances, alpha particle beta, and alpha orbit width. We describe the origin and effect of these dependences analytically and by using recently- developed numerical techniques (Rodrigues et al. 2015 Nucl. Fusion 55 083003). The work suggests high-field machines where fusion-born alphas are sub-Alfvenic or nearly sub-Alfvenic may partially cut off AE resonances, reducing growth rates of AEs and the energy of alphas interacting with them. High-field burning plasma regimes have non-negligible alpha particle beta and AE drive, but faster slowing down time, provided by high electron density, and higher field strength reduces this drive relative to low-field machines with similar power densities. The toroidal mode number of the most unstable modes will tend to be higher in high magnetic field devices. The work suggests that high magnetic field devices have unique, and potentially advantageous, AE instability properties at both low and high densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Code validation for the simulation of supersonic viscous flow about the F-16XL

The viewgraphs and discussion on code validation for the simulation of supersonic viscous flow about the F-16XL are provided. Because of the large potential gains related to laminar flow on the swept wings of supersonic aircraft, interest in the applications of laminar flow control (LFC) techniques in the supersonic regime has increased. A supersonic laminar flow control (SLFC) technology program is currently underway within NASA. The objective of this program is to develop the data base and design methods that are critical to the development of laminar flow control technology for application to supersonic transport aircraft design. Towards this end, the program integrates computational investigations underway at NASA Ames-Moffett and NASA Langley with flight-test investigations being conducted on the F-16XL at the NASA Ames-Dryden Research Facility in cooperation with Rockwell International. The computational goal at NASA Ames-Moffett is to integrate a thin-layer Reynolds averaged Navier-Stokes flow solver with a stability analysis code. The flow solver would provide boundary layer profiles to the stability analysis code which in turn would predict transition on the F-16XL wing. To utilize the stability analysis codes, reliable boundary layer data is necessary at off-design cases. Previously, much of the prediction of boundary layer transition has been accomplished through the coupling of boundary layer codes with stability theory. However, boundary layer codes may have difficulties at high Reynolds numbers, of the order of 100 million, and with the current complex geometry in question. Therefore, a reliable code which solves the thin-layer Reynolds averaged Navier-Stokes equations is needed. Two objectives are discussed, the first in greater depth. The first objective is method verification, via comparisons of computations with experiment, of the reliability and robustness of the code. To successfully implement LFC techniques to the F-16XL wing, the flow about the leading edge must be maintained as laminar flow. Therefore, the second objective is to focus on a series of numerical simulations with different values of angle of attack, alpha, and Reynolds numbers. The purpose of the simulations is to study their effects on the two main factors which precipitate transition to turbulence at leading edges of highly swept wings (e.g., 'spanwise contamination' and 'crossflow instability').

Flores, Jolen↗
NASA NTRS2023

Evaluation of Transport-Equation-Based Transition Models for High-Speed Boundary Layers Using OVERFLOW

Accurate modeling of laminar-turbulent transition is crucial for the design of hypersonic flight systems. However, the current transition models used in production CFD codes are insufficient for high-speed flows. Many extensions to low-speed models have been suggested; however, a thorough verification and validation effort is needed before these models can be used in design settings. Challenges include potentially missing details of the model implementation requirements and/or a complete specification of the input parameters needed to replicate the test findings. A meaningful assessment of the generalization capability of these models is also hindered by a lack of information regarding the specific flow configurations and associated grids employed for model calibration. As a key first step toward model verification, we present an independent assessment of two recently proposed models for high-speed transition, namely, a model within the SST-$\gamma$ framework and a model based on the SST-$\gamma-\nu_L$ equations. These models are implemented in the NASA OVERFLOW 2.3e solver and their performance in predicting first mode, second mode, and crossflow transition has been evaluated for several test cases in the supersonic and hypersonic regimes. Besides the test cases employed by the model developers, which could have also been used for model calibration, the present assessment includes supplementary configurations that contribute to an unbiased assessment of the models. The outcomes presented in this study indicate the potential for the models to be applied to high-speed flight configurations. Key steps toward future improvements to these models are also outlined.

High-speed flow↗
DOE OSTI2021

Capturing in-field root system dynamics with RootTracker

Abstract Optimizing root system architecture offers a promising approach to developing stress tolerant cultivars in the face of climate change, as root systems are critical for water and nutrient uptake as well as mechanical stability. However, breeding for optimal root system architecture has been hindered by the difficulty in measuring root growth in the field. Here, we describe the RootTracker, a technology that employs impedance touch sensors to monitor in-field root growth over time. Configured in a cylindrical, window shutter-like fashion around a planted seed, 264 electrodes are individually charged multiple times over the course of an experiment. Signature changes in the measured capacitance and resistance readings indicate when a root has touched or grown close to an electrode. Using the RootTracker, we have measured root system dynamics of commercial maize (Zea mays) hybrids growing in both typical Midwest field conditions and under different irrigation regimes. We observed rapid responses of root growth to water deficits and found evidence for a “priming response” in which an early water deficit causes more and deeper roots to grow at later time periods. Genotypic variation among hybrid maize lines in their root growth in response to drought indicated a potential to breed for root systems adapted for different environments. Thus, the RootTracker is able to capture changes in root growth over time in response to environmental perturbations.

Aguilar, Jeffrey J.↗
DOE OSTI2026

Electronic Structure of Compressively Strained Bilayer Nickelate Thin Film

The discovery of superconductivity in bulk bilayer nickelates under high pressure, and its subsequent stabilization in compressively strained thin films at ambient pressure, has established a new platform for exploring high-𝑇 c superconductivity beyond the cuprates. Central to this development is the prominent role of the 3⁢𝑑$^2_𝑧$ orbital in shaping the low-energy electronic structure, imposing constraints on microscopic theories and fueling debate over the superconducting mechanism. Here we report a systematic in situ angle-resolved photoemission spectroscopy study of compressively strained bilayer nickelate thin films spanning Ca doping, oxygen stoichiometry and film thickness. Despite variations in oxygen-vacancy disorder and surface termination, the electronic structure remains robust and exhibits a systematic strain-driven evolution consistent with an intermediate-correlation regime. In particular, we demonstrate a Ca-doping-induced electronic structure evolution that is consistent with Fermi-level crossing of the 𝛾 band and is decoupled from the presence of superconductivity, suggesting that, while important for the 𝛾 band to be near the Fermi level, the debated 𝛾 crossing and the resulting Fermi pocket may not be a prerequisite for superconductivity. Together, our results establish key spectroscopic constraints on the minimal fermiology and correlation strength relevant to superconductivity in bilayer nickelates, providing an experimental foundation for microscopic theories of their pairing mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1986

On the instability of a three-dimensional attachment-line boundary layer - Weakly nonlinear theory and a numerical approach

The instability of a three-dimensional attachment-line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite-amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time-dependent Navier-Stokes equations for the attachment-line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite-amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment-line boundary layer is also investigated.

Hall, P.↗
DOE OSTI2025

Inkless, dry printing nanographene via in-Situ Coordinated laser ablation and sintering processes

Printing carbon, such as graphene and other carbon structures, for flexible and printed electronics currently relies on either ink-based printing, laser-induced forward transfer (LIFT), or laser-induced graphitization (LIG) to convert carbon-rich precursor materials, such as polymers, to graphene-like carbon structures. Liquid inks contain toxic solvents, surfactants, and stabilizing additives that degrade electrical conductivity and require high-temperature post-processing. On the other hand, LIG is limited by the substrate. Here, this study introduces an additive manufacturing method for dry-printing carbon nanomaterials, ranging from amorphous carbon to crystalline graphene-like structures, on various substrates. The system utilizes laser ablation of a solid graphite target to create pure carbon nanoparticles in situ and on demand. An inert gas carries the nanoparticles onto the substrate, where they can be deposited either as amorphous structures or laser-sintered in real time to form various graphitic structures. The study of laser processing parameters, specifically fluence and pulse repetition frequency, revealed three unique regimes of nanostructure evolution that influence the morphological and electrical properties of these printed structures. Raman spectroscopy confirmed graphitization with a resistivity slightly higher than that of the bulk graphite target. The conductivity/resistivity could be tuned as a function of sintering laser power. Scanning transmission electron microscopy (STEM) revealed that turbostratic nanographene formed with an interlayer spacing of 0.40 nm. Despite ink-based printing methods, such as screen printing, inkjet printing (IJP), and aerosol jet printing (AJP), this eco-friendly and green manufacturing technique could eliminate toxic chemicals, reduce environmental impact, and enable single-step fabrication of carbon-based devices for applications in wearable sensors, energy storage, flexible electronics, and Internet of Things (IoT) devices.

Additive nanomanufacturing↗
DOE OSTI2026

Selective Reduction of Carbon Dioxide in Water Using [M(bpy2+)(CO)3(I)]2+ (M = Mn, Re) Electrocatalysts with Pendent Cations

Manganese(I) carbonyl complexes are promising electrocatalysts for CO2 reduction, yet their application in homogeneous aqueous media remains limited by poor solubility and selectivity. Here, we report water-soluble Mn(I) and Re(I) complexes fac-[M(bpy2+)(CO)3X]2+ (X = I or Cl), featuring bipyridine ligands functionalized with -Ph-CH2-(NMe3)+ cationic ammonium groups that integrate water solubility with secondary-sphere stabilization. In bicarbonate buffer at pH 6.8, the Mn catalyst is completely selective for CO production at a low overpotential (η = 0.3 V), operating by a protonation-first mechanism with observed rates of ~10 s−1. Pulse radiolysis reveals that the one-electron reduced Mn species undergoes dimerization in the absence of CO2 but uniquely reacts competitively with CO2 through an initial pre-equilibrium followed by fast formation of a dinuclear CO2-bridged species (ΔGo = −12.4 kcal mol−1). At a higher 0.6 V over-potential, a faster reduction-first pathway (~100 s−1) is available upon reduction of the metallo-carboxylic acid intermediate, Mn-CO2H2+; however, this regime is functionally limited by the formation of a resistive, noncatalytic film on the electrode surface. Comparison to the analo-gous water-soluble Re catalyst (kobs = 440 s−1, η = 0.6 V) highlights the distinct mechanistic ad-vantages of earth-abundant Mn in low-potential catalysis. These results demonstrate how cati-onic second-sphere design enables selective, homogeneous CO2 reduction in water while reveal-ing competing radical and electrode-mediated processes that govern catalytic performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Wind tunnel test of a variable-diameter tiltrotor (VDTR) model

This report documents the results from a wind tunnel test of a 1/6th scale Variable Diameter Tiltrotor (VDTR). This test was a joint effort of NASA Ames and Sikorsky Aircraft. The objective was to evaluate the aeroelastic and performance characteristics of the VDTR in conversion, hover, and cruise. The rotor diameter and nacelle angle of the model were remotely changed to represent tiltrotor operating conditions. Data is presented showing the propulsive force required in conversion, blade loads, angle of attack stability and simulated gust response, and hover and cruise performance. This test represents the first wind tunnel test of a variable diameter rotor applied to a tiltrotor concept. The results confirm some of the potential advantages of the VDTR and establish the variable diameter rotor a viable candidate for an advanced tiltrotor. This wind tunnel test successfully demonstrated the feasibility of the Variable Diameter rotor for tilt rotor aircraft. A wide range of test points were taken in hover, conversion, and cruise modes. The concept was shown to have a number of advantages over conventional tiltrotors such as reduced hover downwash with lower disk loading and significantly reduced longitudinal gust response in cruise. In the conversion regime, a high propulsive force was demonstrated for sustained flight with acceptable blade loads. The VDTR demonstrated excellent gust response capabilities. The horizontal gust response correlated well with predictions revealing only half the response to turbulence of the conventional civil tiltrotor.

Matuska, David↗
NASA NTRS2016

Development of Semi-Empirical Damping Equation for Baffled Tank with Oblate Spheroidal Dome

Propellant slosh is a potential source of disturbance that can significantly impact the stability of space vehicles. The slosh dynamics are typically represented by a mechanical model of a spring-mass-damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. The typical parameters required by the mechanical model include natural frequency of the slosh, slosh mass, slosh mass center location, and the critical damping ratio. A fundamental study has been undertaken at NASA MSFC to understand the fluid damping physics from a ring baffle in the barrel section of a propellant tank. An asymptotic damping equation and CFD blended equation have been derived by NASA MSFC team to complement the popularly used Miles equation at different flow regimes. The new development has found success in providing a nonlinear damping model for the Space Launch System. The purpose of this study is to further extend the semi-empirical damping equations into the oblate spheroidal dome section of the propellant tanks. First, previous experimental data from the spherical baffled tank are collected and analyzed. Several methods of taking the dome curvature effect, including a generalized Miles equation, area projection method, and equalized fill height method, are assessed. CFD simulation is used to shed light on the interaction of vorticity around the baffle with the locally curved wall and liquid-gas interface. The final damping equation will be validated by a recent subscale test with an oblate spheroidal dome conducted at NASA MSFC.

Yang, H. Q.↗
DOE OSTI2022

Orbital stability analysis of hypothetical Earth-mass and Luna-mass moons in the Sagarmatha (HD 100777) star system

Abstract Six of the solar system planets have 150 confirmed moons (Earth—1, Mars—2, Jupiter—53, Saturn—53, Uranus—27, Neptune—14) and seven of them (Ganymede, Titan, Callisto, Io, the Moon, Europa, Triton) have masses >0.001 $M_\oplus$. However, no exomoons have yet been discovered despite the successful detection of ∼5000 exoplanets. We can infer, based on the solar system planets, that these exoplanets are capable of hosting one or more exomoons. In this paper, we study the possible existence of hypothetical Earth-mass and Luna-mass moons orbiting the Jupiter-mass planet, Laligurans (HD 100777b) in the Sagarmatha (HD 100777) star system by means of orbital stability. We apply long-term orbital integrations and the MEGNO (mean exponential growth of nearby orbits) chaos indicator to study the orbital stability of the moons and predict a phase-space region comprising periodic, chaotic, and unstable orbits. The phase spaces primarily constitute the moon's semimajor axis, which extends from the host planet's Roche radius to the Hill radius, and full range of eccentricity. Specific points are picked from three different regions of the MEGNO map and run as single-orbit integration for up to 10 billion periods of the innermost orbit. Furthermore, the lifetime and maximum eccentricity maps are generated from the direct integration to inspect the stable and unstable orbital configurations. The analyses of these maps, with the aid of time-series plots, show that both moons maintain stable orbits in the low-eccentricity regime and semimajor axis between the Roche limit and 28.4% of the Hill radius of the planet.

Gautam, Aayush (ORCID:0000000247257589)↗