Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “edge inference”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Trends in Antarctic Surface Temperatures and Ice Extent from Satellite and Surface Observations

Surface air temperatures observed from stations around Antarctica have indicated predominantly positive trends which are as high as 0.4K per decade along the Antarctic Peninsula. Since the average air temperature during summer in the West Antarctic region is about -6 K, and since most of the stations are near coastal regions, it is important to know whether such trends are due to a local or large scale phenomenon. We analyzed about 19 years of infrared satellite data (1979-1997) to obtain a spatially detailed study of the variability of surface temperatures over the entire Antarctic region. We also analyzed potential impact of temperature changes over the sea ice cover using co-registered and nearly coincident microwave satellite data. The surface temperatures inferred from infrared data are shown to agree well with those observed from the limited number of Antarctic station data with a correlation coefficient of 0.98 and a standard deviation of less than 3K. The 19-year temperature data set shows positive trends of 0.12, 0.68, and 0.41 C per decade in ocean, sea ice and ice sheet surfaces, respectively. The sensitivity of the trend to record length was studied, using long term station data sets, and the results show that the trends start to stabilize at around 2 decades which is about the record length of satellite data. However, interannual fluctuations of the temperatures are large, especially in the ice sheets, and the 95% confidence level for the trends ranges from 0.12 to 1.22 for sea ice and from -0.74 to 1.55 C per decade for the ice sheet. In the sea ice regions, the position of the ice edge in winter with respect to the northernmost position of sea ice during the study period is also shown to be influenced by alternating warm and cold anomalies around the continent. This pattern is similar to that suggested by the presence of the Antarctic Circumpolar Wave (ACW) but the wavenumber appears to be predominantly mode 3 (instead of reported mode 2) for the Antarctic region. Unexpectedly, trends in the sea ice extent are found to be positive in most Antarctic sectors during the same period with the exception of the Bellingshausen and Amundsen Seas region.

Comisc, J. C.↗
DOE OSTI2025

Chemical Imaging of Atmospheric Biomass Burning Particles from North American Wildfires

The effects of biomass burning aerosols (BBA) on radiative forcing and cloud formation depend on chemical composition and the internal structures of individual particles within smoke plumes. To improve our understanding of the chemical and physical properties of BBA emitted at different times of the day and their evolution during atmospheric aging, we conducted a study as a part of the Fire Influence on Regional to Global Environments and Air Quality field campaign. Particle samples were collected onboard a research aircraft from smoke plumes from a wildfire in eastern Oregon during late afternoon and nighttime flights on August 28, 2019. A time-resolved aerosol collector was used to collect samples on substrates for offline spectromicroscopic imaging to investigate the single-particle characteristics of BBA particles. Approximately 20,400 individual particles from 10 selected samples were analyzed using computer-controlled scanning electron microscopy coupled with energy-dispersive X-ray microanalysis, revealing their elemental composition, morphology, and viscosity. Elemental microanalysis indicated that aged potassium is likely found in the form of K 2 SO 4 , KNO 3 , and possible K-organic salts. Further chemical speciation and carbon bonding mapping within individual particles were conducted using synchrotron-based scanning transmission X-ray microscopy (STXM) coupled with near edge X-ray absorption fine structure (NEXAFS) spectroscopy. Real-time, water-soluble light absorption measurements were acquired using a particle-into-liquid sampler instrument coupled to a liquid waveguide capillary cell and total organic analyzer. In the late afternoon samples, 65% of the total particle number population consisted entirely of organic components, compared to 46% in the nighttime particles. These differences were attributed to discrepancies in composition at the time of emission and to the daytime condensation and accumulation of photochemically formed secondary organic material on existing BBA particles, a process that halts at night. Microscopy images indicated that particle viscosity was lower in the nighttime particles (<10 1 Pa·s), likely due to increased relative humidity and a higher contribution from hygroscopic inorganic components. The chemical heterogeneity of individual particles was quantified using STXM-derived mixing state parameters. The nature of carbon bonding within individual particles was inferred from the extent of carbon sp 2 hybridization derived from NEXAFS spectra. Average percentages of sp 2 hybridization range between 40% and 60%, with no noticeable differences between late afternoon and nighttime flights. These findings were compared with the online optical properties of both late afternoon and nighttime smoke plumes, providing valuable insights into the complex relationship between chemical composition and optical properties of BBA particles at different times of the day.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

FPGA Architectures for Distributed ML Systems for Real-time Beam Loss De-blending

The Accelerator Real-time Edge AI for Distributed Systems (READS) project’s goal is to create a Artificial Intelligence (AI) system for real-time beam loss de-blending within the accelerator enclosure, which houses two accelerators: the Main Injector (MI) and the Recycler Ring (RR). In periods of joint operation, when both machines contain high intensity beam, radioactive beam losses from MI and RR overlap on the enclosure’s beam loss monitoring Beam Loss Monitor (BLM) system, making it difficult to attribute those losses to a single machine. Incorrect diagnoses result in unnecessary downtime that incurs both financialand experimental cost. The ML system will automatically disentangle each machine’s contributions to those measured losses, while not disrupting the existing operations-critical functions of the BLM system. This paper will focus on the evolution of the architectures, which provided the high-frequency, low-latency collection of synchronized data streams to make real-time inferences. The ML models, used for learning both local and global machine signatures and producing high quality inferences based on raw BLM loss measurements, will only be discussed at a high-level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Schmid factor crack propagation and tracking crystallographic texture markers of microstructural condition in direct energy deposition additive manufacturing of Ti-6Al-4V

Metallic additive manufacturing (AM) provides a customizable and tailorable manufacturing process for new engineering designs and technologies. The greatest challenge currently facing metallic AM is maintaining control of microstructural evolution during solidification and any solid state phase transformations during the build process. Ti-6Al-4V has been extensively surveyed in this regard, with the potential solid state and solidification microstructures explored at length. In this study, we evaluate the applicability of previously determined crystallo- graphic markers of microstructural condition observed in electron beam melting powder bed fusion (PBF-EB) builds of Ti-6Al-4V in a directed energy deposition (DED) build process. The aim of this effort is to elucidate whether or not these specific crystallographic textures are useful tools for indicating microstructural conditions in AM variants beyond PBF-EB. Parent β-Ti grain size was determined to be directly related to α-Ti textures in the DED build process, and the solid state microstructural condition could be inferred from the intensity of specific α-Ti orientations. Qualitative trends on the as-solidified β-Ti grain size were also determined to be related to the presence of a fiber texture, and proposed as a marker for as-solidified grain size in any cubic metal melted by AM. Analysis of the DED Ti-6Al-4V build also demonstrated a near complete fracture of the build volume, suspected to originate from accumulated thermal stresses in the solid state. Crack propagation was found to only appreciably occur in regions of slow cooling with large α +β colonies. Schmid factors for the basal and prismatic α-Ti systems explained the observed crack pathway, including slower bifurcation in colonies with lower Schmid factors of both slip systems. Colony morphologies and localized equiaxed β-Ti solidification were also found to originate from build pauses during production and uneven heating of the build edges during deposition. Tailoring of DED Ti-6Al-4V microstructures with the insight gained here is proposed, along with cautionary insight on preventing unplanned build pauses to maintain an informed and controlled thermal environment for microstructural control.

36 MATERIALS SCIENCE↗
NASA NTRS2012

The Swift Burst Alert Telescope Detected Seyfert 1 Galaxies: X-Ray Broadband Properties and Warm Absorbers

We present results from an analysis of the broadband, 0.3-195 keV, X-ray spectra of 48 Seyfert 1-1.5 sources detected in the very hard X-rays with the Swift Burst Alert Telescope (BAT). This sample is selected in an all-sky survey conducted in the 14-195 keV band. Therefore, our sources are largely unbiased toward both obscuration and host galaxy properties. Our detailed and uniform model fits to Suzaku/BAT and XMM-Newton/BAT spectra include the neutral absorption, direct power-law, reflected emission, soft excess, warm absorption, and narrow Fe I K[alpha] emission properties for the entire sample. We significantly detect O VII and O VIII edges in 52% of our sample. The strength of these detections is strongly correlated with the neutral column density measured in the spectrum. Among the strongest detections, X-ray grating and UV observations, where available, indicate outflowing material. The ionized column densities of sources with O VII and O VIII detections are clustered in a narrow range with Nwarm [approx] 1021 cm-2, while sources without strong detections have column densities of ionized gas an order of magnitude lower. Therefore, we note that sources without strong detections likely have warm ionized outflows present but at low column densities that are not easily probed with current X-ray observations. Sources with strong complex absorption have a strong soft excess, which may or may not be due to difficulties in modeling the complex spectra of these sources. Still, the detection of a flat [Gamma] [approx] 1 and a strong soft excess may allow us to infer the presence of strong absorption in low signal-to-noise active galactic nucleus spectra. Additionally, we include a useful correction from the Swift BAT luminosity to bolometric luminosity, based on a comparison of our spectral fitting results with published spectral energy distribution fits from 33 of our sources.

Winter, Lisa M.↗
DOE OSTI2024

Application of Modified Meshgraphnets for Subsurface Prediction during CO2 Sequestration

In the face of the increasingly dire consequences of anthropogenic climate change, capturing and storing carbon dioxide is paramount. However, several impediments exist to the safe and effective subsurface storage of CO2, such as cost of transport, identification of suitable sites for subsurface storage, and assessment of long-term risk from storage in subsurface aquifers. Accurate subsurface modeling is necessary to ensure that CO2 storage is both safe and effective. Still, such modeling has traditionally required either substantial time and computational power (numerical simulation) or a substantial amount of pre-existing data for training (machine learning models). Additionally, these models lack flexibility in dealing with both changes in discretization of the input data and generalizability beyond the data on which they are trained. In order to address these issues, this research applies graph neural networks (GNNs) to predict subsurface saturation and pressure during CO₂ injection in a model of the Illinois Basin-Decatur Project (IBDP). GNNs provide a flexible, intuitive method for representing and manipulating complex unstructured data, which is often found in many practical domain problems such as fluid flow and subsurface characterization. These unstructured grids are easily represented in GNNs by representing spatially-localized features such as permeability, porosity, saturation, and pressure as nodes in a graph and relationships between these properties as edges connecting these nodes. This research applies a specific GNN model called MeshGraphNets (MGN) to model the change in CO2 saturation and pressure over a 50-month time period (36 months of injection, 14 months post-injection). The MGN model leverages a message passing process that allows the network to learn both the spatial and temporal dynamics of this system simultaneously. Additionally, training on a limited dataset (64 realizations, 20 time points each) resulted in a high degree of accuracy in saturation prediction both within the same timeframe as the training (20 months, 0.039 average RMSE) and when projecting out to the end of injection (36 months, 0.053 average RMSE). Temporal predictions such as those generated by MGNs and other similar models are prone to accumulated error over time; in order to address this, a multi-step rollout (MSR) training process was applied to calculate training loss. This method mimics the forward prediction during inference by “rolling out” multiple time points in a single training step using the previous prediction as input to the MGN model. By calculating the loss several time steps forward from the current prediction, the model is forced to find a more stable state over time. Application of MSR to the MGN model resulted in an average 15% reduction in inference error over time during forward prediction. This study showcases the immense potential of GNNs as a game-changing methodology for predicting pressure and saturation evolution in CCS projects, ultimately paving the way for more sustainable and effective carbon storage solutions. Presentation prepared for the 2024 AiChE Annual Meeting, October 27 to November 1 2024, San Diego, CA.

Holcomb, Paul↗
NASA NTRS1998

Carbonates in fractures of Martian meteorite Allan Hills 84001: petrologic evidence for impact origin

Carbonates in Martian meteorite Allan Hills 84001 occur as grains on pyroxene grain boundaries, in crushed zones, and as disks, veins, and irregularly shaped grains in healed pyroxene fractures. Some carbonate disks have tapered Mg-rich edges and are accompanied by smaller, thinner and relatively homogeneous, magnesite microdisks. Except for the microdisks, all types of carbonate grains show the same unique chemical zoning pattern on MgCO3-FeCO3-CaCO3 plots. This chemical characteristic and the close spatial association of diverse carbonate types show that all carbonates formed by a similar process. The heterogeneous distribution of carbonates in fractures, tapered shapes of some disks, and the localized occurrence of Mg-rich microdisks appear to be incompatible with growth from an externally derived CO2-rich fluid that changed in composition over time. These features suggest instead that the fractures were closed as carbonates grew from an internally derived fluid and that the microdisks formed from a residual Mg-rich fluid that was squeezed along fractures. Carbonate in pyroxene fractures is most abundant near grains of plagioclase glass that are located on pyroxene grain boundaries and commonly contain major or minor amounts of carbonate. We infer that carbonates in fractures formed from grain boundary carbonates associated with plagiociase that were melted by impact and dispersed into the surrounding fractured pyroxene. Carbonates in fractures, which include those studied by McKay et al. (1996), could not have formed at low temperatures and preserved mineralogical evidence for Martian organisms.

Carbonates/analysis↗
DOE OSTI2022

Utah FORGE 2439: A Multi-Component Approach to Characterizing In-Situ Stress

Core-based in-situ stress estimation, Triaxial Ultrasonic Velocity (labTUV) data, and Deformation Rate Analysis (DRA) data for Utah FORGE well 16A(78)-32 using triaxial ultrasonic velocity and deformation rate analysis. Report documenting a multi-component approach to characterizing in-situ stress at the U.S. DOE FORGE EGS site: laboratory, modeling and field measurement. Core-based methods for in-situ stress estimation were applied using samples from 5 intervals within the Utah FORGE 16A(78)-32 well. At three of these locations, Triaxial Ultrasonic Velocity (labTUV) tests were performed, resulting in experimentally-determined relationships between wave velocities and stresses. Non-monotonic increase in the velocity-stress relationships are inferred provide evidence of stress history and are therefore used to estimate in-situ stress magnitudes. Additionally, Deformation Rate Analysis (DRA) tests were run on core plugs from various orientations at each of the 5 sampling locations. These, too, provide evidence of stress history based on stress-strain behavior. A novel Weight of Evidence (WoE) method was developed as a means of synthesizing in-situ stress evidence from these two types of tests. Results indicate the minimum horizontal stress gradient ranges from 0.58 psi/ft to 0.69 psi/ft, with 4 of the 5 values between 0.66 psi/ft and 0.69 psi/ft. The vertical stress gradient ranges from 1.05 psi/ft to 1.12 psi/ft, with 4 of the 5 zones given results between 1.09 psi/ft and 1.12 psi/ft. The maximum horizontal stress gradient ranges from 0.98 psi/ft to 1.34 psi/ft, with 4 of the 5 zones falling between 0.98 psi/ft and 1.24 psi/ft. The stress regime thus appears to be on the edge between normal faulting and strike-slip faulting, potentially flipping back and forth between the two regimes due to variability of rock properties, structures such as faults, and/or thermal anomalies.

15 GEOTHERMAL ENERGY↗
NASA NTRS1992

Derivation of surface properties from Magellan altimetry data

The fit of the Hagfors model to the Magellan altimetry data provides a means to characterize the surface properties of Venus. However, the derived surface properties are only meaningful if the model provides a good representation of the data. The Hagfors model provides a good representation of the data. The Hagfors model is generally a realistic fit to surface scattering properties of a nadir-directed antenna such as the Magellan altimeter; however, some regions of the surface of Venus are poorly described by the existing model, according to the goodness of fit parameter provided on the ARCDR CD-ROMs. Poorly characterized regions need to be identified and fit to new models in order to derive more accurate surface properties for use in inferring the geological processes that affect the surface in those regions. We have compared the goodness of fit of the Hagfors model to the distribution of features across the planet, and preliminary results show a correlation between steep topographic slopes and poor fits to the standard model, as has been noticed by others. In this paper, we investigate possible relations between many classes of features and the ability of the Hagfors model to fit the observed echo profiles. In the regions that are not well characterized by existing models, we calculate new models that compensate for topographic relief in order to derive improved estimates of surface properties. Areas investigated to date span from longitude 315 through 45, at all latitudes covered by Magellan. A survey of those areas yields preliminary results that suggest that topographically high regions are well suited to the current implementation of the Hagfors model. Striking examples of such large-scale good fits are Alpha Regio, the northern edges of Lada Terra, and the southern edge of Ishtar Terra. Other features that are typically well fit are the rims of coronae such as Heng-O and the peaks of volcanos such as Gula Mons. Surprisingly, topographically low regions, such as the ubiquitous plains areas, are modeled poorly in comparison. However, this generalization has has exceptions: Lakshmi Planum is an elevated region that is not well fit compared to the rest of neighboring Ishtar, while the southern parts of topographically low Guinevere Planitia are characterized quite well by the Hagfors model. Features that are candidates for improved models are impact craters, coronae, ridges of significant scale, complex ridged terrains, moderate-sized mountains, and sharp terrain boundaries. These features are chosen because the goodness of fit is likely to be most affected either by departures from normal incidence angles or by sharp changes in terrain type within a single footprint. Most large features that are elevated with respect to their surroundings will suffer from steep slope effects, and smaller coronae and impact craters will probably suffer due to rapid changes in their appearance within a single footprint (10-20 km).

Lovell, Amy J.↗
NASA NTRS1975

Eutectic superalloys by edge-defined, film-fed growth

The feasibility of producing directionally solidified eutectic alloy composites by edge-defined, film-fed growth (EFG) was carried out. The three eutectic alloys which were investigated were gamma + delta, gamma/gamma prime + delta, and a Co-base TaC alloy containing Cr and Ni. Investigations into the compatibility and wettability of these metals with various carbides, borides, nitrides, and oxides disclosed that compounds with the largest (negative) heats of formation were most stable but poorest wetting. Nitrides and carbides had suitable stability and low contact angles but capillary rise was observed only with carbides. Oxides would not give capillary rise but would probably fulfill the other wetting requirements of EFG. Tantalum carbide was selected for most of the experimental portion of the program based on its exhibiting spontaneous capillary rise and satisfactory slow rate of degradation in the liquid metals. Samples of all three alloys were grown by EFG with the major experimental effort restricted to gamma + delta and gamma/gamma prime + delta alloys. In the standard, uncooled EFG apparatus, the thermal gradient was inferred from the growth speed and was 150 to 200 C/cm. This value may be compared to typical gradients of less than 100 C/cm normally achieved in a standard Bridgman-type apparatus. When a stream of helium was directed against the side of the bar during growth, the gradient was found to improve to about 250 C/cm. In comparison, a theoretical gradient of 700 C/cm should be possible under ideal conditions, without the use of chills. Methods for optimizing the gradient in EFG are discussed, and should allow attainment of close to the theoretical for a particular configuration.

Hurley, G. F.↗
DOE OSTI2022

Ion thermal transport in the H-mode edge transport barrier on DIII-D

The power balance ion heat flux in the pedestal region on DIII-D increases and becomes increasingly anomalous (above conventional neoclassical) in experiments with higher temperature and lower density pedestals where the ion collisionality ($v^*_i$) is lowered toward values expected on ITER. Direct measurements of the main-ion temperature are shown to be essential on DIII-D when calculating the ion heat flux due to differences between the temperature of $D^+$ and the more commonly measured $C^{6+}$ impurity ions approaching the separatrix. Neoclassical transport calculations from NEO and non-linear gyrokinetic calculations using CGYRO are consistent with these observations and show that while neoclassical transport plays an important role, the turbulent ion heat flux due to ion scale electrostatic turbulence is significant and can contribute similar or larger ion heat fluxes at lower collisionality. Beam emission spectroscopy and Doppler backscattering measurements in the steep gradient region of the H-mode pedestal reveal increased broadband, long-wavelength ion scale fluctuations for the low $v^*_i$ discharges at the radius where the non-linear CGYRO simulations were run. Taken together, increased fluctuations, power balance calculations, and gyrokinetic simulations show that the above neoclassical ion heat fluxes, including the increases at lower $v^*_i$, are likely due to weakly suppressed ion scale electrostatic turbulence. These new results are based on world first inferred ion and electron heat fluxes in the pedestal region of deuterium plasmas using direct measurements of the deuterium temperature for power balance across ion collisionalities covering an order of magnitude from high $v^*_i$ values of 1.3 down to ITER relevant $v^*_i$ ~0.1.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2002

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

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

Young, Andrew J.↗
DOE OSTI2022

Diagnosing the pedestal magnetic field and magnetohydrodynamics radial structure with pedestal–scrape of layer electron cyclotron emission radiation inversion in H-mode plasma (invited)

Forward modeling is used to interpret inversion patterns of the pedestal–Scrape of Layer (SOL) Electron Cyclotron Emission (ECE) in DIII-D H-mode experiments. The modeling not only significantly improves the ECE data interpretation quality but also leads to the potential measurements of (1) the magnetic field strength |B| at the separatrix, (2) the pedestal |B| evolution during an inter-Edge Localized Mode (ELM) period, and (3) the pedestal Magnetohydrodynamics (MHD) radial structure. The ECE shine-through effect leads to three types of pedestal–SOL radiation inversions that are discussed in this paper. The first type of inversion is the non-monotonic T e,rad profile with respect to the major radius. Using the ECE frequency at the minimum T e,rad , the inversion can be applied to measure the magnetic field |B| at the separatrix and calibrate the mapping of the ECE channels with respect to the separatrix. The second type of inversion refers to the opposite phase between the radiation fluctuations δT e,rad at the pedestal and SOL. This δT e,rad phase inversion is sensitive to density and temperature fluctuations at the pedestal foot and, thus, can be used to qualitatively measure the MHD radial structure. The third type of inversion appears when the pedestal and SOL T e,rad evolve in an opposite trend, which can be used to infer the pedestal |B| field change during an inter-ELM period. Finally, the bandwidth effect on measuring δT e,rad due to pedestal MHD is also investigated in the radiation modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2016

Probabilistic Prognosis of Non-Planar Fatigue Crack Growth

Quantifying the uncertainty in model parameters for the purpose of damage prognosis can be accomplished utilizing Bayesian inference and damage diagnosis data from sources such as non-destructive evaluation or structural health monitoring. The number of samples required to solve the Bayesian inverse problem through common sampling techniques (e.g., Markov chain Monte Carlo) renders high-fidelity finite element-based damage growth models unusable due to prohibitive computation times. However, these types of models are often the only option when attempting to model complex damage growth in real-world structures. Here, a recently developed high-fidelity crack growth model is used which, when compared to finite element-based modeling, has demonstrated reductions in computation times of three orders of magnitude through the use of surrogate models and machine learning. The model is flexible in that only the expensive computation of the crack driving forces is replaced by the surrogate models, leaving the remaining parameters accessible for uncertainty quantification. A probabilistic prognosis framework incorporating this model is developed and demonstrated for non-planar crack growth in a modified, edge-notched, aluminum tensile specimen. Predictions of remaining useful life are made over time for five updates of the damage diagnosis data, and prognostic metrics are utilized to evaluate the performance of the prognostic framework. Challenges specific to the probabilistic prognosis of non-planar fatigue crack growth are highlighted and discussed in the context of the experimental results.

Leser, Patrick E.↗
DOE OSTI2025

Predicted thresholds for RMP ELM suppression access in double-null configurations

Modeling of DIII-D plasmas spanning shapes from single to double-null (DN) reveals new insights into the nature of resonant magnetic perturbation (RMP) conditions necessary for edge-localized mode (ELM) suppression. The suppression of ELMs with RMPs has proven difficult in DN configurations, where no device has thus far reported any hints of suppression. Modeling using the GPEC code finds a reduced high-field side response closer to DN shaping. The resulting synthetic diagnostic measurements are consistent with what has been observed in experiments on DIII-D, validating the plasma response model in this regime. While common metrics for suppression do not illustrate a clear distinction between single-null (SN) and DN cases, the pedestal top resonant field does show a ∼20% decrease at DN shaping in modeling. Field penetration is assessed using linear tearing theory, which demonstrates a lack of sufficient pedestal top resonant flux in the DN shape, requiring at least 1.5 × greater RMP coil currents than what was used in experiment. Analysis from drift kinetic simulations further indicate up to 2 × larger critical island widths are required at the pedestal top for tearing mode growth compared to SN cases. Effective island widths inferred from 3D ideal MHD are also analyzed, where maximum widths in lab coordinates indicate a threshold of ∼18–24× the ion gyroradius for sufficient profile flattening for ELM suppression. These results suggest that ELM suppression may be possible in DN with sufficiently large RMP coil amplitude. Future prospects of achieving RMP ELM suppression in the DN configuration may involve going to lower triangularity, which is also highlighted in this work.

RMP ELM suppression↗
DOE OSTI2024

Codiscovering graphical structure and functional relationships within data: A Gaussian Process framework for connecting the dots

Most problems within and beyond the scientific domain can be framed into one of the following three levels of complexity of function approximation. Type 1: Approximate an unknown function given input/output data. Type 2: Consider a collection of variables and functions, some of which are unknown, indexed by the nodes and hyperedges of a hypergraph (a generalized graph where edges can connect more than two vertices). Given partial observations of the variables of the hypergraph (satisfying the functional dependencies imposed by its structure), approximate all the unobserved variables and unknown functions. Type 3: Expanding on Type 2, if the hypergraph structure itself is unknown, use partial observations of the variables of the hypergraph to discover its structure and approximate its unknown functions. These hypergraphs offer a natural platform for organizing, communicating, and processing computational knowledge. While most scientific problems can be framed as the data-driven discovery of unknown functions in a computational hypergraph whose structure is known (Type 2), many require the data-driven discovery of the structure (connectivity) of the hypergraph itself (Type 3). We introduce an interpretable Gaussian Process (GP) framework for such (Type 3) problems that does not require randomization of the data, access to or control over its sampling, or sparsity of the unknown functions in a known or learned basis. Its polynomial complexity, which contrasts sharply with the super-exponential complexity of causal inference methods, is enabled by the nonlinear ANOVA capabilities of GPs used as a sensing mechanism.

Science & Technology - Other Topics↗
DOE OSTI2024

Chasing the multi-modal plasma response in MAST-U

Abstract Achieving edge localized modes (ELMs) suppression in spherical tokamaks by applying resonant magnetic perturbations (RMPs) has proven challenging. The poloidal spectrum of the applied RMP is a key parameter that has an impact on the capability to mitigate and eventually suppress ELMs. In this work the resistive magnetohydrodynamic code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is used to evaluate the possibility of directly measuring the plasma response in MAST-U, and particularly its variation as function of the applied poloidal spectrum, in order to guide the experimental validation of the predicted best RMP configuration for ELM suppression. Toroidal mode number n = 2 RMP is considered to minimize the presence of sidebands, and to avoid the deleterious core coupling of n = 1. Singular Value Decomposition is used to highlight linearly independent structures in the simulated magnetic 3D fields and how those structures can be measured at the wall where the magnetic sensors are located. Alternative ways to measure the multimodal plasma response and how they can be used to infer the best RMP configuration to achieve ELM suppression are also presented, including the plasma displacement and the 3D footprints at the divertor plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

A Remarkable Three Hour Thermonuclear Burst from 4U 1820-30

We present a detailed observational and theoretical study of an approximately three hour long X-ray burst (the "super burst") observed by the Rossi X-ray Timing Explorer (RXTE) from the low mass X-ray binary (LMXB) 4U 1820-30. This is the longest X-ray burst ever observed from this source, and perhaps one of the longest ever observed in great detail from any source. We show that the super burst is thermonuclear in origin. Its peak luminosity of approximately 3.4 x 10(exp 38) ergs s(exp -1) is consistent with the helium Eddington limit for a neutron star at approximately 7 kpc, as well as the peak luminosity of other, shorter, thermonuclear bursts from the same source. The super burst begins in the decaying tail of a more typical (approximately equal to 20 s duration) thermonuclear burst. These shorter, more frequent bursts are well known helium flashes from this source. The level of the accretion driven flux as well as the observed energy release of upwards of 1.5 x 10(exp 42) ergs indicate that helium could not be the energy source for the super burst. We outline the physics relevant to carbon production and burning on helium accreting neutron stars and present calculations of the thermal evolution and stability of a carbon layer and show that this process is the most likely explanation for the super burst. Ignition at the temperatures in the deep carbon "ocean" requires greater than 30 times the mass of carbon inferred from the observed burst energetics unless the He flash is able to trigger a deflagration from a much smaller mass of carbon. We show, however, that for large columns of accreted carbon fuel, a substantial fraction of the energy released in the carbon burning layer is radiated away as neutrinos, and the heat that is conducted from the burning layer in large part flows inward, only to be released on timescales longer than the observed burst. Thus the energy released during the event possibly exceeds that observed in X-rays by more than a factor of ten, making the scenario of burning a large mass of carbon at great depths consistent with the observed fluence without invoking any additional trigger. A strong constraint on this scenario is the recurrence time: to accrete an ignition column of 1013 g cm (exp -1) takes approximately 13/(M/3 x 10(exp 17) g s(exp -1) yr. Spectral analysis during the super burst reveals the presence of a broad emission line between 5.8 - 6.4 keV and an edge at 8 - 9 keV likely due to reflection of the burst flux from the inner accretion disk in 4U 1820-30. We believe this is the first time such a signature has been unambiguously detected in the spectrum of an X-ray burst.

Strohmayer, Tod E.↗