Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “G CODES”

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 703 records · Page 39

The Anatomy of an Unusual Edge-on Protoplanetary Disk. I. Dust Settling in a Cold Disk

As the earliest stage of planet formation, massive, optically thick, and gas-rich protoplanetary disks provide key insights into the physics of star and planet formation. When viewed edge-on, high-resolution images offer a unique opportunity to study both the radial and vertical structures of these disks and relate this to vertical settling, radial drift, grain growth, and changes in the midplane temperatures. In this work, we present multi-epoch Hubble Space Telescope and Keck scattered light images, and an Atacama Large Millimeter/submillimeter Array 1.3 mm continuum map for the remarkably flat edge-on protoplanetary disk SSTC2DJ163131.2–242627, a young solar-type star in ρ Ophiuchus. We model the 0.8 μm and 1.3 mm images in separate Markov Chain Monte Carlo (MCMC) runs to investigate the geometry and dust properties of the disk using the MCFOST radiative transfer code. In scattered light, we are sensitive to the smaller dust grains in the surface layers of the disk, while the submillimeter dust continuum observations probe larger grains closer to the disk midplane. An MCMC run combining both data sets using a covariance-based log-likelihood estimation was marginally successful, implying insufficient complexity in our disk model. The disk is well characterized by a flared disk model with an exponentially tapered outer edge viewed nearly edge-on, though some degree of dust settling is required to reproduce the vertically thin profile and lack of apparent flaring. A colder than expected disk midplane, evidence for dust settling, and residual radial substructures all point to a more complex radial density profile to be probed with future, higher-resolution observations.

47 OTHER INSTRUMENTATION↗

An extended moments model of quantum efficiency for metals and semiconductors

The complexity of photocathode designs and detector materials, and the need to model their performance for short pulse durations, response to high frequency photons, presence of coatings and/or thinness of the absorptive layer, necessitates modifications to three-step and moments models of photoemission that are used in simulation codes. Here, methods to include input from computationally intensive approaches, such as density functional theory to model optical properties and transfer matrix approaches to treat emission from the surface or transport past coatings, by means of parametric models are demonstrated. First, a technique to accurately represent optical behavior so as to model reflectivity and penetration depth is given. Second, modifications to bulk models arising from the usage of thin film architectures, and a means to rapidly calculate them, is provided. Third, a parameterization to model the impact of wells associated with coatings and surface layers on the transmission probably is given. In all cases, the methods are computationally efficient and designed to allow for including input from numerically intensive approaches that would otherwise be unavailable for simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spectroscopic quadrupole moments in 124 $\mathrm{Xe}$

The Xe isotopic chain with four valence protons above the Z=50 shell closure is an ideal laboratory for the study of the evolution of nuclear deformation. At the N=82 shell closure, 136 Xe presents all characteristics of a doubly closed shell nucleus with a spherical shape. In the very neutron-deficient isotopes close to N=50, the α-decay chain of Xe was investigated to probe the radioactive decay properties near the drip-line and the magicity of 100 Sn. Additionally, the Xe isotopes present higher order symmetries in the nuclear deformation such as the octupole degree of freedom near N=60 and N=90 or O(6) symmetry in stable isotopes. The relevance of the O(6) symmetry has been investigated by measuring the spectroscopic quadrupole moment of the first excited states in 124 Xe. In the O(6) symmetry limit, the spectroscopic quadrupole moment of collective states is expected to be null. A stable 124 Xe beam with energies of 4.03A MeV and 4.11 A MeV was used to bombard a nat W target at the GANIL facility. Excited states were populated via the safe Coulomb excitation reaction. The collision of the heavy ions with a large Z at low energy make this reaction sensitive to the diagonal E2 matrix element of the excited states. The recoils were detected in the VAMOS++ magnetic spectrometer and the γ rays in the AGATA tracking array. The least squares fitting code gosia was used for the analysis to extract both E2 and M1 transitional and E2 diagonal matrix elements. The rotational ground state band was populated up to the 8$^{+}_{1}$ state as well as the 2$^{+}_{2}$ and 4$^{+}_{2}$ states. Using high precision spectroscopic data to constrain the gosia fit, the spectroscopic quadrupole moments of the 2$^{+}_{1}$, 4$^{+}_{1}$, and 6$^{+}_{1}$ states were determined for the first time. The spectroscopic quadrupole moments were found to be negative, large, and constant in the ground state band underlining the prolate axially deformed ground state band of 124 Xe. Here, the present experimental data confirm that the O(6) symmetry is substantially broken in 124 Xe .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Energy transfer and scale dynamics in 2D and 3D laser-driven jets

We demonstrate a methodology for diagnosing the multiscale dynamics and energy transfer in complex HED flows with realistic driving and boundary conditions. The approach separates incompressible, compressible, and baropycnal contributions to energy scale-transfer and quantifies the direction of these transfers in (generalized) wavenumber space. We use this to compare the kinetic energy (KE) transfer across scales in simulations of 2D axisymmetric vs fully 3D laser-driven plasma jets. Using the FLASH code, we model a turbulent jet ablated from an aluminum cone target in the configuration outlined. We show that, in addition to its well known bias for underestimating hydrodynamic instability growth, 2D modeling suffers from significant spurious energization of the bulk flow by a turbulent upscale cascade. In 2D, this arises as vorticity and strain from instabilities near the jet's leading edge transfer KE upscale, sustaining a coherent circulation that helps propel the axisymmetric jet farther (⁠ ≈25% by 3.5 ns) and helps keep it collimated. In 3D, the coherent circulation and upscale KE transfer are absent. Here, the methodology presented here may also help with inter-model comparison and validation, including future modeling efforts to alleviate some of the 2D hydrodynamic artifacts highlighted in this study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Harmonic generation in a terawatt x-ray free-electron laser

Abstract Terawatt x-ray free-electron lasers (XFELs) represent the frontier in further development of x-ray sources and require high current densities with strong transverse focusing. In this paper, we investigate the implications/potentialities of TW XFELs on the generation of harmonics at still shorter wavelengths and higher photon energies. The simulations indicate that significant power levels are possible at high harmonics of the XFEL resonance and that these XFELs can be an important coherent source of hard x-rays through the gamma ray spectrum. For this purpose, we use the MINERVA simulation code which self-consistently includes harmonic generation. Both helical and planar undulators are discussed in which the fundamental is at 1.5 Å and study the associated harmonic generation. While tapered undulators are needed to reach TW powers at the fundamental, the taper does not enhance the harmonics because the taper must start before saturation of the fundamental, while the harmonics saturate before this point is reached. Nevertheless, the harmonics reach substantial powers. Simulations indicate that, for the parameters under consideration, peak powers of the order of 180 MW are possible at the fifth harmonic with a photon energy of about 41 keV and still high harmonics may also be generated at substantial powers. Such high harmonic powers are certain to enable a host of enhanced applications.

Physics↗

Investigating high-energy proton-induced reactions on spherical nuclei: Implications for the preequilibrium exciton model

A number of accelerator-based isotope production facilities utilize $100-200$ MeV proton beams due to the high production rates enabled by high-intensity beam capabilities and the greater diversity of isotope production brought on by the long range of high-energy protons. However, nuclear reaction modeling at these energies can be challenging because of the interplay between different reaction modes and a lack of existing guiding cross section data. A Tri-lab collaboration has been formed between the Lawrence Berkeley, Los Alamos, and Brookhaven National Laboratories to address these complexities by characterizing charged-particle nuclear reactions relevant to the production of established and novel radioisotopes. In the inaugural collaboration experiments, stacked-targets of niobium foils were irradiated at the Brookhaven Linac Isotope Producer ($E_p=200$ MeV) and the Los Alamos Isotope Production Facility ($E_p=100$ MeV) to measure $^{93}$Nb(p,x) cross sections between $50-200$ MeV. The results were compared with literature data as well as the default calculations of the nuclear model codes TALYS, CoH, EMPIRE, and ALICE. The default code predictions largely failed to reproduce the measurements. Therefore, we developed a standardized procedure, which determines the reaction model parameters that best reproduce the most prominent reaction channels in a physically justifiable manner. Overall, the primary focus of the procedure was to determine the best parameterization for the pre-equilibrium two-component exciton model. This modeling study revealed a trend towards a relative decrease for internal transition rates at intermediate proton energies ($E_p=20-60$ MeV) in the current exciton model as compared to the default values. The results of this work are instrumental for the planning, execution, and analysis essential to isotope production.

43 PARTICLE ACCELERATORS↗

Three dimensions, two microscopes, one code: Automatic differentiation for x-ray nanotomography beyond the depth of focus limit

Conventional tomographic reconstruction algorithms assume that one has obtained pure projection images, involving no within-specimen diffraction effects nor multiple scattering. Advances in x-ray nanotomography are leading toward the violation of these assumptions, by combining the high penetration power of x-rays, which enables thick specimens to be imaged, with improved spatial resolution that decreases the depth of focus of the imaging system. We describe a reconstruction method where multiple scattering and diffraction effects in thick samples are modeled by multislice propagation and the 3D object function is retrieved through iterative optimization. We show that the same proposed method works for both full-field microscopy and for coherent scanning techniques like ptychography. Our implementation uses the optimization toolbox and the automatic differentiation capability of the open-source deep learning package TensorFlow, demonstrating a straightforward way to solve optimization problems in computational imaging with flexibility and portability.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Experimental studies of plasma-antenna coupling with the JET Alfven Eigenmode Active Diagnostic

This paper presents a dedicated study of plasma-antenna (PA) coupling with the Alfven Eigenmode Active Diagnostic (AEAD) in JET. Stable AEs and their resonant frequencies f, damping rates gamma < 0, and toroidal mode numbers n are measured for various PA separations and limiter versus X-point magnetic configurations. Two stable AEs are observed to be resonantly excited at distinct low and high frequencies in limiter plasmas. The values of f and n do not vary with PA separation. However, |gamma| increases with PA separation for the low-f, but not high-f, mode, yet this may be due to slightly different edge conditions. The high-f AE is detected throughout the transition from limiter to X-point configuration, though its damping rate increases; the low-f mode, on the other hand, becomes unidentifiable. The linear resistive MHD code CASTOR is used to simulate the frequency scan of an AEAD-like external antenna. For the limiter pulses, the high-f mode is determined to be an n = 0 GAE, while the low-f mode is likely an n = 2 TAE. During the transition from limiter to X-point configuration, CASTOR indicates that n = 1 and 2 EAEs are excited in the edge gap. These results extend previous experimental studies in JET and Alcator C-Mod; validate the computational work performed by Dvornova et al 2020 Phys. Plasmas 27 012507; and provide guidance for the optimization of PA coupling in upcoming JET energetic particle experiments, for which the AEAD will aim to identify the contribution of alpha particles to AE drive during the DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational fluid dynamics investigations of flow, heat transfer, and oxidation in heat recovery steam generator

Modern heat recovery steam generators (HRSGs) operate at elevated temperatures, leading to the formation of oxides inside the tubes of heat exchangers (HXs). This oxide growth reduces the heat recovery efficiency. Moreover, after reaching a certain critical thickness, some oxide scales detach from the tube surface (exfoliation), causing erosion damage to the components downstream. Predicting the metal temperature distribution and associated oxide thickness in the HXs of an HRSG can aid in mitigating these problems. Here, a computational fluid dynamics (CFD) model was developed within the commercial code STAR-CCM+ for the prediction of fluid flow, conjugate heat transfer, and associated oxidation in HRSGs. Moreover, a new Porous Media Model (PMM) method was developed to model the fin effect on the heat transfer in HX, which can substantially reduce the prohibitive computational costs of fin meshes. The developed CFD model was used to conduct a high-fidelity simulation of a real-scale HRSG to investigate flow, heat transfer, and oxide growth. The calculated oxide thickness on different tubes can be used to identify HX regions that require oxide-resistant coatings to prevent exfoliation and ensuing damages. Furthermore, this CFD framework can serve as a reference for future studies that intend to model and investigate high-temperature oxidation in HXs used for any applications.

42 ENGINEERING↗

Propagation of ion cyclotron emission in the DIII-D tokamak

Experiments in L- and H-mode plasmas performed on the DIII-D tokamak explored ion cyclotron emission (ICE) propagation via the recently upgraded ICE diagnostic. The distance between the plasma and the outer wall was scanned to alter the evanescent region, which was seen to impact the detection of runaway electron-driven whistler modes in a comparable frequency range to typical ICE harmonics (f≳10fci) [Heidbrink et al., Plasma Phys. Controlled Fusion 61, 014007 (2018)]. In L-mode plasmas, core ICE did not exhibit clear trends as the plasma-wall distance was altered. Instead, inadvertent changes to the fast ion density correlated with different ICE fine structures; the most spectral bands (spaced ∼150–200 kHz apart) were observed at the highest fast ion density, and a just single band when nfast decreased by ∼ 40%. The full-wave Petra-M code simulated core ICE propagation through the plasma and to ICE probes, agreeing with experimental findings that these are likely compressional fast waves, which propagate to probe locations with no dependence on plasma-wall gap. Edge ICE in H-mode plasmas was similarly unaffected by plasma-wall distance. Unlike in core ICE cases, there were negligible changes to the edge fast ion distribution and edge bulk plasma profiles. However, changes in ELM frequency seemed most closely tied to ICE behavior. Ultimately, core and edge ICE did not exhibit clear dependencies on the plasma-wall distance, and ICE harmonics were observed with even the largest plasma-wall gaps. This independence bodes well for similar diagnostics in future tokamaks, which might be placed further from high radiation areas without adversely impacting measurement capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measuring the mass of the black widow PSR J1555-2908

ABSTRACT Accurate measurements of the masses of neutron stars are necessary to test binary evolution models, and to constrain the neutron star equation of state. In pulsar binaries with no measurable post-Keplerian parameters, this requires an accurate estimate of the binary system’s inclination and the radial velocity of the companion star by other means than pulsar timing. In this paper, we present the results of a new method for measuring this radial velocity using the binary synthesis code Icarus. This method relies on constructing a model spectrum of a tidally distorted, irradiated star as viewed for a given binary configuration. This method is applied to optical spectra of the newly discovered black widow PSR J1555–2908. By modeling the optical spectroscopy alongside optical photometry, we find that the radial velocity of the companion star is 397 ± 4 km s−1 (errors quoted at 95 per cent confidence interval), as well as a binary inclination of >75°. Combined with γ-ray pulsation timing information, this gives a neutron star mass of 1.67$^{+0.15}_{-0.09}$ M⊙ and a companion mass of 0.060$^{+0.005}_{-0.003}$ M⊙, placing PSR J1555–2908 at the observed upper limit of what is considered a black widow system.

Kennedy, M. R.↗

The Ristra Project: FY20/21 Milestone Report

The ASC Advanced Technology Development and Mitigation (ATDM) sub-program was established in 2014 to develop new simulation tools operating on exascale-class computers to serve NNSA (see Appendix B). Over the course of ATDM, LANL management have set a strategy for exascale-class application codes that follows two supportive and mutually risk-mitigating paths: evolution for established production integrated design codes (IDCs) – with a strong pedigree within the user community – based upon existing programming paradigms(MPI+X); and a new start ATDM project, Ristra, a high-risk/high-reward push for a next-generation multi-physics, multi-scale simulation toolkit based on emerging advanced programming systems(with an initial focus on data-flow task-based models exemplified by Legion). The role of Ristra as the high-risk/high-reward path for LANL’s codes was fully consistent with the goals of ATDM as described in Appendix B, in particular its emphasis on evolving ASC capabilities through novel computing programming models and computing technologies.

97 MATHEMATICS AND COMPUTING↗

Lifetimes of two- and four-neutron-hole states in 204 Pb and 206 Pb

The centroid-shift method was used to determine half-lives of excited levels in 204 Pb and 206 Pb. For the states with 𝐼 𝜋 = 12 + , 17 − , and 19 − in 204 Pb, values of 𝑇 1/2 = 4.9⁢(3), 3.9(3), and 0.9(3) ns, respectively, were determined. The half-life of the 𝐼 𝜋 = 9 − level in 206 Pb was measured to be 1.1(3) ns. In 204 Pb, the 12 + state arises from two neutron holes in the 𝑖 13/2 subshell, while the 17 − and 19 − levels are built on 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑝$^{−1}_{1/2}$) and 𝜈⁡(𝑖$^{−3}_{13/2}$, 𝑓$^{−1}_{5/2}$) four-neutron-hole configurations, respectively. The 9 − state in 206 Pb results from a 𝜈⁡(𝑖$^{−1}_{13/2}$, 𝑓${−1}_{5/2}$) excitation. As a result, the transition probabilities for the decay of these isomeric levels are reasonably consistent with the expected single-particle values, and these states are described well by shell-model calculations using the kshell code with the KHHE interaction.

190 ≤ A ≤ 219↗

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a 5-dimensional (5D) continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular 2D velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. We describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma by using 1st- and 2nd-order elements with the pseudo-inverse method and comparing to the bilinear mapping.

Fusion Plasma↗

WHONDRS Summer 2019 Sampling Campaign: Global River Corridor Sediment FTICR-MS, Dissolved Organic Carbon, Aerobic Respiration, Elemental Composition, Grain Size, Total Nitrogen and Organic Carbon Content, Bacterial Abundance, and Stable Isotopes (v10)

The WHONDRS Summer 2019 Sampling (S19S) Campaign collected samples in 97 globally distributed river corridor systems between July and September 2019. Surficial streambed sediments were collected at three locations within each site (upstream, midstream, and downstream). Surface water was collected at the downstream site. This dataset includes a portion of the data types produced from the sediment samples and does not include any results from the surface water. Surface water data can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/1603775. The S19S campaign was designed with the science community to ask questions associated with links among core/transient metabolomes, microbial metabolism, biogeochemical function, and physical properties of watershed and river corridor systems. This data package was originally published in December 2020. It was updated in April 2021 (v2; modified files), May 2022 (v3; new and modified files), November 2022 (v4; new and modified files), January 2023 (v5 and v6; new and modified files), February 2023 (v7; new and modified files), August 2023 (v8; new and modified files), July 2025 (v9; new and modified files), and June 2026 (v10; new and modified files). Version 10 includes major restructuring of the files and their formats. See the change history in the readme for details. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset contains (1) readme, (2) data dictionary, (3) file level metadata, (4) sampling protocol, (5) a data checks report, (6) field metadata, (7) international generic sample number (IGSN) mapping file, (8) joint genome institute (JGI) mapping file, (9) a folder of site photos, and (10) a folder of sample data. The sample data folder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) and FTICR methods; (2) respiration rates and normalized respiration rates calculated from laboratory sediment incubations, (3) X-ray fluorescence (XRF) elemental composition, (4) grain size distribution from the < 2 millimeter fraction of sediment, (5) total organic carbon and total nitrogen content, (6) bacterial abundance flow cytometry data, (7) sediment stable isotopes (C and N), (8) laboratory methods, (9) water and sediment mass and volume used in incubations, (10) a folder with code, inputs, and plots for the respiration rates, (11) a folder containing XRF data, (12) a folder containing unaligned, peak-picked high resolution characterization of dissolved organic matter from sediment via 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory), (12) and a folder containing processed FTICR-MS data. All files are .msa, .txt, .spx, .xml, .csv, .pdf, .txt, .jpg, or .jpeg.

13C↗

Variational Quantum Chemistry Programs in JaqalPaq

We present example quantum chemistry programs written with JaqalPaq, a python meta-programming language used to code in Jaqal (Just Another Quantum Assembly Language). These JaqalPaq algorithms are intended to be run on the Quantum Scientific Computing Open User Testbed (QSCOUT) platform at Sandia National Laboratories. Our exemplars use the variational quantum eigensolver (VQE) quantum algorithm to compute the ground state energies of the H2, HeH+, and LiH molecules. Since the exemplars focus on how to program in JaqalPaq, the calculations of the second-quantized Hamiltonians are performed with the PySCF python package, and the mappings of the fermions to qubits are obtained from the OpenFermion python package. Using the emulator functionality of JaqalPaq, we emulate how these exemplars would be executed on an error-free QSCOUT platform and compare the emulated computation of the bond-dissociation curves for these molecules with their exact forms within the relevant basis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discovery of widespread transcription initiation at microsatellites predictable by sequence-based deep neural network

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

59 BASIC BIOLOGICAL SCIENCES↗

Investigating the role of edge neutrals in exciting tearing mode activity and achieving flat temperature profiles in LTX-β

Abstract We present observations, numerical simulations, and analysis from experiments in the Lithium Tokamak Experiment-Beta (LTX- β ) in which the electron temperature profile ( T e ( r )) shifts from flat to peaked and a tearing mode is also destabilized when the average density ( n e ave ) exceeds ∼10 19 m −3 . Flat T e ( r ) is obtained routinely in LTX- β , with a lithium coated, low-recycling first wall, once the external fueling is stopped and density decays [Boyle et al 2023 Nucl. Fusion 63 056020]. In the present experiment, flat T e profiles can be sustained while maintaining constant n e ave below a line averaged density threshold ( n e ave th ) of ∼10 19 m −3 . Above n e ave th , T e ( r ) shifts from flat to peaked and a tearing mode is destabilized. Due to low recycling, the achieved n e ave can be controlled precisely by external fueling and hence, a certain threshold of the edge neutral inventory from the external fueling is experimentally manifested through n e ave th . The goal of the present work is to investigate the role of edge neutrals in determining T e ( r ) and MHD stability in the unique low-recycling regime of LTX- β . Our hypothesis is that the peaking of T e ( r ) beyond n e ave th is due ultimately to the edge cooling by the cold neutrals beyond a critical fueling flux. At lower fueling flux, flat T e ( r ) results in broader pressure profile and lower resistivity, which in turn stabilizes the tearing mode. This hypothesis is supported by edge neutral density estimation by DEGAS 2 code. Mode analysis by singular value decomposition confirms the tearing mode structure to be m / n = 2/1 ( m and n being the poloidal and toroidal mode numbers). Linear tearing stability analysis with M3D-C1 predicts that plasmas with n e ave > 10 19 are highly susceptible to a n = 1 tearing mode. ORBIT simulations, however, confirmed that the tearing modes do not contribute to the loss of fast ions from neutral beam injection. This study shows for the first time that the neutral inventory at the edge could be one of the deciding factors for the achievability of the unique operation regime of flat T e ( r ) and the excitation of tearing activity that could be disruptive for the plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗