Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape constraints”

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 919 records · Page 51

Refinements to the Graves–Pitarka Kinematic Rupture Generator, Including a Dynamically Consistent Slip-Rate Function, Applied to the 2019 Mw 7.1 Ridgecrest Earthquake

The main objective of this study is to develop physics–based constraints on the spatiotemporal variation of the slip–rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth–dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip–rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface–rupturing, planar strike–slip fault that includes a weak (negative to low–stress–drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high–stress–drop patches designed to mirror the large–slip areas seen in various rupture model inversions of the event. The locations of the high–stress–drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip–rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth–dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth–dependent parametric variation of the slip–rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth–dependent variation in the shape of the slip–rate function improves the fit to the observed near–fault ground motions in the 0.5–3 s period range.

58 GEOSCIENCES↗

Sustaining Thermal Protection Systems Needed for Uranus Probes

Summary: Achieving Uranus in-situ science objectives through deploying probes in the atmosphere requires thermal protection systems capable of withstanding extreme entry environments. NASA’s 3-D Woven based thermal protection systems (TPS), developed specifically to meet this challenge, are mature, efficient, and capable of withstanding such entry. These science-enabling capabilities need to be sustained, otherwise mission implementation risk may increase to the point that the mission is untenable. The heatshield technology has been developed in concert with industry, and manufacturing processes have been technology transferred. Current technology readiness does not guarantee future availability. Sustainment requires maintaining expertise within NASA as well as ensuring that proven industrial expertise and capabilities will be readily available in the future. A constant awareness and risk assessment followed by risk mitigation are required. The scientific community needs to be aware of these challenges. It needs to engage NASA to ensure these capabilities will be available when needed. As was the case with the atrophy of heritage carbon phenolic that led to developing the 3-D Woven capabilities, our vigilance is again needed to sustain these enabling TPS capabilities for Uranus and for other missions. Background: NASA invested in and developed Heatshield for Extreme Entry Environment Technology (HEEET) TPS that was matured to TRL 6 in 2019 as a tiled system with seams and gap-fillers. Manufacturing readiness included development of specialized looms, molding, and infusion processes using commercial partners. In support of the Mars Sample Return mission, a single-layer, seamless, single-piece variant referred to as 3MDCP (3-D woven Mid-Density Carbon-Phenolic) was developed and is limited to a 1.3m maximum diameter scale due to weaving width limitations of 80”. Molding techniques developed to transform a flat, 2”-thick 3-D woven preform into a sphere-cone shape prior to phenolic resin infusion has been demonstrated. A team at NASA Ames is currently working with industry to demonstrate infusion at full-scale. By 2026, 3MDCP will be at a high level of maturity in technical, manufacturing, and integration readiness as well as material characterization. HEEET and 3MDCP for Uranus Entry: Several system studies have evaluated and baselined 3-D Woven TPS for Uranus entry, both HEEET as well as 3MDCP. If the aeroshell carrying the descent probe is 1.3m or less, then 3MDCP is a very capable and the most mass-efficient choice. If the aeroshell is bigger than 1.3m diameter, then HEEET can meet the mission need as it allows for an aeroshell of any diameter. HEEET and 3MDCP are capable of Saturn in-situ science missions prioritized in the New Frontiers-5 draft Announcement of Opportunity (AO) and if Saturn is the mission of choice, there is a pathway for sustaining the capability. Capability Sustainment: Given NASA’s current budgetary constraints, the postponement of the next New Frontiers AO to no earlier than 2026 and the unknown delay in starting the UOP flagship effort point to potential gap years for the HEEET and 3MDCP capabilities. The first step in capability sustainment is to understand the risks as they emerge, assess those risks, and develop plans to mitigate them. Keeping the science community, who are interested in UOP as well as the larger Giant Planet missions, informed of the emerging challenges requires continuous risk assessment/mitigation. Without these steps, atrophy is more likely within the industrial partners with which NASA has partnered. This presentation will provide more details as well as steps NASA can take to minimize the impact to the UOP mission.

Ethiraj Venkatapathy↗

Supportability Technologies for Future Exploration Missions

Future long-duration human exploration missions will be challenged by resupply limitations and mass and volume constraints. Consequently, it will be essential that the logistics footprint required to support these missions be minimized and that capabilities be provided to make them highly autonomous from a logistics perspective. Strategies to achieve these objectives include broad implementation of commonality and standardization at all hardware levels and across all systems, repair of failed hardware at the lowest possible hardware level, and manufacture of structural and mechanical replacement components as needed. Repair at the lowest hardware levels will require the availability of compact, portable systems for diagnosis of failures in electronic systems and verification of system functionality following repair. Rework systems will be required that enable the removal and replacement of microelectronic components with minimal human intervention to minimize skill requirements and training demand for crews. Materials used in the assembly of electronic systems (e.g. solders, fluxes, conformal coatings) must be compatible with the available repair methods and the spacecraft environment. Manufacturing of replacement parts for structural and mechanical applications will require additive manufacturing systems that can generate near-net-shape parts from the range of engineering alloys employed in the spacecraft structure and in the parts utilized in other surface systems. These additive manufacturing processes will need to be supported by real-time non-destructive evaluation during layer-additive processing for on-the-fly quality control. This will provide capabilities for quality control and may serve as an input for closed-loop process control. Additionally, non-destructive methods should be available for material property determination. These nondestructive evaluation processes should be incorporated with the additive manufacturing process - providing an in-process capability to ensure that material deposited during layer-additive processing meets required material property criteria.

Watson, Kevin↗

Mineralogy, fabric and deformation domains in D" across the southwestern border of the African LLSVP

Recent advances in seismic anisotropy studies that jointly use reflections and shear wave splitting have proven to place tight constraints on the plausible anisotropic and deformation scenarios in the D" region. We apply this novel methodology to a large area of the D" region beneath the South Atlantic, in proximity to and within the African large low seismic velocity province (LLSVP). This area of the mantle is characterized by a transition from fast to slow seismic velocity anomalies and it is thought to be the location of deep-seated plumes responsible for hotspot volcanism. Attempting to probe mantle composition and deformation along the LLSVP borders may provide key information on mantle dynamics. By analysing seismic phases sampling this region, we detect a D" discontinuity over a large area beneath the South Atlantic, with inferred depth ranges ~170 to ~240 km above the core–mantle boundary. We find evidence for a D" reflector within the area of the LLSVP. Shear wave splitting observations suggest that anisotropy is present in this region of the mantle, in agreement with previous studies that partially sampled this region. We model the observations considering lattice- and shape-preferred orientation of materials expected in the D" region. A regional variation of mineralogy, phase transition boundaries, and deformation direction is required to explain the data. We infer two distinct domains of mineralogy and deformation: aligned post-perovskite outside the LLSVP and aligned bridgmanite within the LLSVP. Finally, the scenario depicted by this study agrees well with the current hypotheses for the composition of the LLSVP and with the prevalence of vertical deformation directions expected to occur along the LLSVPs borders.

58 GEOSCIENCES↗

Parametric Modeling and Economic Analysis of a 2MW th 3-Stream sCO 2 Heat Exchanger

Here, this paper presents the design and cost optimization of a novel 2MW th 3-stream sCO2 plate-fin heat exchanger. This heat exchanger design is unique in that it uses reduced metal oxide particle-to-sCO2 heat exchanger for cost-effective energy storage applications. The design uses low velocity, laminar air as the re-oxidizing reactant to transfer the heat of the re-oxidizing reaction to a sCO2 power loop. The design of the heat exchanger is based on a 2-D, 3-fluid plate/fin heat transfer model. The model parameterizes the size, shape, and number of passages of the heat exchanger to calculate the temperature profile, pressure drop, and fluid velocities of all three fluids. Global heat exchanger parameters such as the effectiveness and total heat transferred to the sCO2 are then calculated for overall performance. Due to the value and increased use of sCO2 heat exchangers in power cycles, a cost model of the system based on the unique high temperature/high pressure operating conditions was created using quotes from reference projects and market analysis. These quoted air-to-sCO2 heat exchangers are then processed using multiple weighting factors pertinent to heat exchanger design, including heat exchanger type, maximum temperature, differential pressures, fluids, duty, and more. These factors are then used in an exponential function in order to generate a parameterized cost curve. The design and cost of the heat exchanger are then optimized using the SMPSO genetic algorithm in Python. The optimization objectives for the system are to maximize the overall system effectiveness, including an air recuperator for preheating, and to minimize unit costs. Additional constraints are added to the system for the sCO2 and air pressure drops, air velocity to reduce particle entrainment, and the length and volume of the heat exchanger.

Cost Model↗

Fomalhaut's Debris Disk and Planet: Constraining the Mass of Formalhaut B from Disk Morphology

Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M(sub pl) < 3M(sub J), an orbital semimajor axis a(sub pl) > 101.5AU, and an orbital eccentricity e(sub pl) = 0.11 - 0.13. These conclusions are independent of Fom b's photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom b's orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a approximately equals 133AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of e approximately equals 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of approximately 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to approximately 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut b's nominal space velocity does not bear this out, but the astrometric uncertainties are difficult to quantify. Even if the apsidal misalignment proves real, our calculated upper mass limit of 3 M(sub J) still holds. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M(sub Earth) of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.

Chiang, E.↗

Using Cassini CIRS Data to Constrain Enceladus' Libration State

Given the non-spherical shape of Enceladus, the satellite may experience gravitational torques that will cause it to physically librate as it orbits Saturn. Physical libration would produce a diurnal oscillation in the longitude of Enceladus' tidal bulge, which could have a profound effect on the diurnal stresses experienced by the surface of the satellite. Although Cassini ISS has placed an observational upper limit on Enceladus' libration amplitude, small amplitude librations may have geologically significant consequences. For example, a physical libration will affect heat production along the tiger stripes as produced by tidal shear heating. We have modeled the expected power emitted along the tiger stripes for various types of physical libration and have quantified which types of physical libration best reproduce the observed power f1ux as detailed in Cassini CIRS data. We find that including a physical libration does allow better fits to the observations and we have identified regions of the libration phase space that where these fits are optimized. A physical libration has important implications for tidal dissipation within Ence1adus and if identified may provide an additional constraint on its interior mass distribution.

Hurford, T. A.↗

Cross-correlation of Dark Energy Survey Year 3 lensing data with ACT and P l a n c k thermal Sunyaev-Zel’dovich effect observations. II. Modeling and constraints on halo pressure profiles

Hot, ionized gas leaves an imprint on the cosmic microwave background via the thermal Sunyaev-Zel’dovich (tSZ) effect. The cross-correlation of gravitational lensing (which traces the projected mass) with the tSZ effect (which traces the projected gas pressure) is a powerful probe of the thermal state of ionized baryons throughout the Universe and is sensitive to effects such as baryonic feedback. In a companion paper (Gatti et al. Phys. Rev. D 105, 123525 (2022)), we present tomographic measurements and validation tests of the cross-correlation between Galaxy shear measurements from the first three years of observations of the Dark Energy Survey and tSZ measurements from a combination of Atacama Cosmology Telescope and Planck observations. In this work, we use the same measurements to constrain models for the pressure profiles of halos across a wide range of halo mass and redshift. We find evidence for reduced pressure in low-mass halos, consistent with predictions for the effects of feedback from active Galactic nuclei. We infer the hydrostatic mass bias (B ≡ M 500c /M SZ from our measurements, finding B = 1.8 ± 0.1 when adopting the Planck-preferred cosmological parameters. We additionally find that our measurements are consistent with a nonzero redshift evolution of B, with the correct sign and sufficient magnitude to explain the mass bias necessary to reconcile cluster count measurements with the Planck-preferred cosmology. Our analysis introduces a model for the impact of intrinsic alignments (IAs) of galaxy shapes on the shear-tSZ correlation. We show that IA can have a significant impact on these correlations at current noise levels.

79 ASTRONOMY AND ASTROPHYSICS↗

Using Cassini CIRS Data to Constrain Enceladus' Libration State

Given the non-spherical shape of Enceladus, the satellite may experience gravitational torques that will cause it to physically librate as it orbits Saturn. Physical Libration would produce a diurnal oscillation in the longitude of Enceladus' tidal bulge, which could have a profound effect on the diurnal stresses experienced by the surface of the satellite. Although Cassini ISS has placed an observational upper limit on Enceladus' libration amplitude, small amplitude librations may have geologically significant consequences. For example, a physical libration will affect heat production. along the tiger stripes as produced by tidal shear heating. We have modeled the expected power en-litted along the tiger stripes for various types of physical libration and have quantified which types of physical libration best reproduce the observed power flux as detailed in Cassini CIRS data. We find that including a physical libration does allow better fits to the observations and we have identified regions of the libration phase space that where these fits are optimized. A physical libration has important implications for tidal dissipation within Enceladus and if identified may provide an additional constraint on its interior mass distribution.

Hurford, T. A.↗

Implications of Wind-Assisted Aerial Navigation for Titan Mission Planning and Science Exploration

The recent Titan Saturn System Mission (TSSM) proposal incorporates a montgolfiere (hot air balloon) as part of its architecture. Standard montgolfiere balloons generate lift through heating of the atmospheric gases inside the envelope, and use a vent valve for altitude control. A Titan aerobot (robotic aerial vehicle) would have to use radioisotope thermoelectric generators (RTGs) for electric power, and the excess heat generated can be used to provide thermal lift for a montgolfiere. A hybrid montgolfiere design could have propellers mounted on the gondola to generate horizontal thrust; in spite of the unfavorable aerodynamic drag caused by the shape of the balloon, a limited amount of lateral controllability could be achieved. In planning an aerial mission at Titan, it is extremely important to assess how the moon-wide wind field can be used to extend the navigation capabilities of an aerobot and thereby enhance the scientific return of the mission. In this paper we explore what guidance, navigation and control capabilities can be achieved by a vehicle that uses the Titan wind field. The control planning approach is based on passive wind field riding. The aerobot would use vertical control to select wind layers that would lead it towards a predefined science target, adding horizontal propulsion if available. The work presented in this paper is based on aerodynamic models that characterize balloon performance at Titan, and on TitanWRF (Weather Research and Forecasting), a model that incorporates heat convection, circulation, radiation, Titan haze properties, Saturn's tidal forcing, and other planetary phenomena. Our results show that a simple unpropelled montgolfiere without horizontal actuation will be able to reach a broad array of science targets within the constraints of the wind field. The study also indicates that even a small amount of horizontal thrust allows the balloon to reach any area of interest on Titan, and to do so in a fraction of the time needed by the unpropelled balloon. The results show that using the Titan wind field allows an aerobot to significantly extend its scientific reach, and that a montgolfiere (unpropelled or propelled) is a highly desirable architecture that can very significantly enhance the scientific return of a future Titan mission.

aerobot↗

Supra Arcade Downflows in the Earth's Magnetotail

Pinpointing the location of a single reconnection event in the corona is difficult due to observational constraints, although features directly resulting from this rapid reconfiguration of the field lines can be observed beyond the reconnection site. One set of such features are outflows in the form of post-reconnection loops, which have been linked to observations of supra-arcade downflows (SADs). SADs appear as sunward-traveling, density-depleted regions above flare arcades that develop during long duration eruptions. The limitations of remote sensing methods inherently results in ambiguities regarding the interpretation of SAD formation. Of particular interest is how these features are related to post-reconnection retracting magnetic field lines. In planetary magnetospheres, similar events to solar flares occur in the form of substorms, where reconnection in the anti-sunward tail of the magnetosphere causes field lines to retract toward the planet. Using data from the Time History of Events and Macroscopic Interactions during Substorms (THEMIS), we compare one particular aspect of substorms, dipolarization fronts, to SADs. Dipolarization fronts are observed as rapid but temporary changes in the magnetic field of the magnetotail plasma sheet into a more potential-like dipolar shape. These dipolarization fronts are believed to be retracting post-reconnection field lines. We combine data sets to show that the while the densities and magnetic fields involved vary greatly between the regimes, the plasma betas and Alfvén speeds are similar. These similarities allow direct comparison between the retracting field lines and their accompanying wakes of rarified plasma observed with THEMIS around the Earth to the observed morphological density depletions visible with XRT and AIA on the Sun. These results are an important source of feedback for models of coronal current sheets.

Kobelski, Adam↗

Intrinsic and environmental drivers of pairwise cohesion in wild Canis social groups

Animals within social groups respond to costs and benefits of sociality by adjusting the proportion of time they spend in close proximity to other individuals in the group (cohesion). Variation in cohesion between individuals, in turn, shapes important group-level processes such as subgroup formation and fission–fusion dynamics. Although critical to animal sociality, a comprehensive understanding of the factors influencing cohesion remains a gap in our knowledge of cooperative behavior in animals. We tracked 574 individuals from six species within the genus Canis in 15 countries on four continents with GPS telemetry to estimate the time that pairs of individuals within social groups spent in close proximity and test hypotheses regarding drivers of cohesion. Pairs of social canids (Canis spp.) varied widely in the proportion of time they spent together (5%–100%) during seasonal monitoring periods relative to both intrinsic characteristics and environmental conditions. The majority of our data came from three species of wolves (gray wolves, eastern wolves, and red wolves) and coyotes. For these species, cohesion within social groups was greatest between breeding pairs and varied seasonally as the nature of cooperative activities changed relative to annual life history patterns. Across species, wolves were more cohesive than coyotes. For wolves, pairs were less cohesive in larger groups, and when suitable, small prey was present reflecting the constraints of food resources and intragroup competition on social associations. Pair cohesion in wolves declined with increased anthropogenic modification of the landscape and greater climatic variability, underscoring challenges for conserving social top predators in a changing world. We show that pairwise cohesion in social groups varies strongly both within and across Canis species, as individuals respond to changing ecological context defined by resources, competition, and anthropogenic disturbance. Our work highlights that cohesion is a highly plastic component of animal sociality that holds significant promise for elucidating ecological and evolutionary mechanisms underlying cooperative behavior.

59 BASIC BIOLOGICAL SCIENCES↗

Particle transport constraints via Bayesian spectral fitting of multiple atomic lines

Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.

47 OTHER INSTRUMENTATION↗

PWR loading pattern optimization with reinforcement learning

The core loading pattern optimization problem belongs to the class of combinatorial optimization problem and has been studied since the dawn of commercial nuclear energy industry. It is characterized by multiple objectives and constraints, with a very high number of candidate patterns, which makes it impossible to solve explicitly. Stochastic optimization methodologies including Genetic Algorithms and Simulated Annealing are used by different nuclear utilities and vendors to perform fuel cycle reload design. Nevertheless, hand-designed solutions continue to be the prevalent method in the industry. To improve the state-of-the-art core reload patterns, we aim to create a method as scalable as possible, that agrees with the designer's goal of performance and safety. To help in this task Deep Reinforcement Learning (DRL), in particular Proximal Policy Optimization is leveraged. DRL has recently experienced a strong impetus from its successes applied to games, sometimes even reaching 'super-human' performances. This paper lays out the foundation of this method and proposes to study the behavior of several hyper-parameters that influence the DRL algorithm. The algorithm is highly dependent on multiple factors such as an exploration/exploitation trade-off that manifests through different parameters such as the number of loading patterns seen and the number of samples collected before a policy update, but also the shape of the objective function derived for the core design. Experimental results also demonstrate the effectiveness of the method in finding high-quality solutions from scratch within a reasonable amount of time. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Constraints on the formation history of the elliptical galaxy NGC 3923 from the colors of its globular clusters

We present a study of the colors of globular clusters associated with the elliptical galaxy NGC 3923. Our final sample consists of Wasington system C and T(sub 1) photometry for 143 globular cluster candidates with an expected contamination of no more than 10%. We find that the color distribution of the NGC 3923 globular cluster system (GCS) is broad and appears to have at least two peaks. A mixture modeling analysis of the color distribution indicates that a two-component model is favored over a single-component one at a high level of confidence (greater than 99%). This evidence for more than one population in the GCS of NGC 3923 is similar to that previously noted for the four other elliptical galaxies for which similar data have been published. Furthermore, we find that the NGC 3923 GCS is redder than the GCSs of previously studed elliptical galaxies of similar luminosity. The median metallicity inferred from our (C-(T(sub 1)))(sub 0) colors is (Fe/H)(sub med) = -0.56, with an uncertainty of 0.14 dex arising from all sources of uncertainty in the mean color. This is more metal rich than the median metallicity found for the GCS of M87 using the same method, (Fe/H)(sub med) = -0.94. Since M87 is more luminous than NGC 3923, this result points to significant scatter about any trend of higher GCS metallicity with increasing galaxy luminosity. We also show that there is a color gradient in the NGC 3923 GCS corresponding to about -0.5 dex in Delta(Fe/H)/Delta(log r). We conclude that the shape of the color distribution of individual GCSs and the variation in mean color among the GCSs of ellipticals are difficult to understand if elliptical galaxies are formed in a single protogalactic collapse. Models in which ellipticals and their globular clusters are formed in more than one event, such as a merger scenario, are more successful in accounting for these observations.

Zepf, Stephen E.↗

Stability and transport of gyrokinetic critical pedestals

Abstract A gyrokinetic threshold model for pedestal width–height scaling prediction is applied to multiple devices. A shaping and aspect ratio scan is performed on National Spherical Torus Experiment (NSTX) equilibria, finding Δ ped = 0.92 A 1.04 κ − 1.24 0.38 δ β θ , ped 1.05 for the wide-pedestal branch with pedestal width Δ ped , aspect ratio A , elongation κ , triangularity δ , and normalized pedestal height β θ , ped . The width–transport scaling is found to vary significantly if the pedestal height is varied either with a fixed density or fixed temperature, showing how fueling and heating sources affect the pedestal density and temperature profiles for the kinetic-ballooning-mode (KBM) limited profiles. For an NSTX equilibrium, at fixed density, the wide branch is Δ ped = 0.028 ( q e / Γ e − 1.7 ) 1.5 ∼ η e 1.5 and at fixed temperature Δ ped = 0.31 ( q e / Γ e − 4.7 ) 0.85 ∼ η e 0.85 , where q e and Γ e are turbulent electron heat and particle fluxes and η e = ∇ ln ⁡ T e / ∇ ln ⁡ n e for an electron temperature T e and density n e . Pedestals close to the KBM limit are shown to have modified turbulent transport coefficients compared to the strongly driven KBMs. The role of flow shear is studied as a width–height scaling constraint and pedestal saturation mechanism for a standard and lithiated wide pedestal discharge. Finally, the stability, transport, and flow shear constraints are combined and examined for an NSTX experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensitivity of magnetic islands in permanent magnet stellarators using the gradient and Hessian methods

Stellarator plasmas are known to be very sensitive to perturbations in the magnetic field. The permanent magnet stellarator was in part developed as a solution to high machining tolerances placed on the shape properties of electromagnetic coils in traditional stellarators. However, as a consequence of this high sensitivity to the field structure, sensitivities of permanent magnet stellarator plasmas to perturbations of permanent magnet properties must necessarily be well-understood. The gradient and Hessian matrix methods have been previously demonstrated to be useful sensitivity analysis methods for modular coils. We apply these two methods to the study of island width sensitivities in both the MUSE and PM4STELL permanent magnet stellarator projects. These sensitivity methods were used to determine the relative impacts of permanent magnet parameter perturbations on island widths in the vacuum field approximation of both stellarator equilibria. The square of resonant magnetic field perturbation is used here as a proxy for island width. In particular, gradients of magnetizations of individual magnets were examined in MUSE, as well as gradients of magnet group displacements informed by device design. Three different forms of permanent magnet magnetization perturbations are investigated for MUSE, and the flux surface response to perturbations is demonstrated. The Hessian matrix method is applied to PM4STELL, illustrating the sensitivity of dominant island widths to displacements of toroidal wedge structures. These methods allow for selective direction of experimental resources toward regions of heightened sensitivity, while constraints on less impactful permanent magnet parameters can be relaxed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations

Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.

42 ENGINEERING↗