Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computing”

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

Surface-Driven Evolution of the Anomalous Hall Effect in Magnetic Topological Insulator MnBi 2 Te 4 Thin Films

We report understanding the effects of the interfacial modification to the functional properties of magnetic topological insulator thin films is crucial for developing novel technological applications from spintronics to quantum computing. Here, a large electronic and magnetic response is reported to be induced in the intrinsic magnetic topological insulator MnBi 2 Te 4 by controlling the propagation of surface oxidation. It is shown that the formation of the surface oxide layer is confined to the top 1–2 unit cells but drives large changes in the overall magnetic response. Specifically, a dramatic reversal of the sign of the anomalous Hall effect is observed to be driven by finite thickness magnetism, which indicates that the film splits into distinct magnetic layers each with a unique electronic signature. These data reveal a delicate dependence of the overall magnetic and electronic response of MnBi 2 Te 4 on the stoichiometry of the top layers. This study suggests that perturbations resulting from surface oxidation may play a non-trivial role in the stabilization of the quantum anomalous Hall effect in this system and that understanding targeted modifications to the surface may open new routes for engineering novel topological and magnetic responses in this fascinating material.

36 MATERIALS SCIENCE↗

Suppression of one-dimensional weak localization by band asymmetry

Herein we investigate disorder-induced localization in metals that break time-reversal and inversion symmetries through their energy dispersion, ε k ≠ε– k , but lack Berry phases. In the perturbative regime of disorder, we show that weak localization is suppressed due to a mismatch of the Fermi velocities of left and right movers. To substantiate this analytical result, we perform quench numerics on chains shorter than the Anderson localization length ξ—the latter computed and verified to be finite using the recursive Green's function method—and find a sharp rise in the saturation value of the participation ratio due to band asymmetry, indicating a tendency to delocalize. Interestingly, for weak disorder strength η, we see a better fit to the scaling behavior ξ∝1/η 2 for asymmetric bands than conventional symmetric ones.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the prediction of pressure effects for the combination kinetics of two alkyl radicals with the geometric mean rule

Here, alkyl combination reactions, along with their reverse alkane unimolecular decomposition reactions, are widely present in pyrolysis and combustion systems. These reactions occur mainly at temperatures at which their rate coefficients depend on pressure, and only a few studies have explored this dependence for large size alkyls. The accurate calculation of these kinetic data is limited by the size of the alkyl fragments involved, since the theoretical methods used are computationally expensive. Therefore, little reliable data is available for the modeling of real fuels alkane decomposition kinetics. In this paper, we propose an alternative approach that allows an accurate estimation of pressure-dependent kinetic rates for the combination of long-chain alkyl radicals, based on tabulated reaction models. The rates of combination of linear alkyls radicals forming ethane to pentane were calculated using the CASPT2-F12/cc-pVDZ-F12 level within VRC-TST theory combined with master-equation simulations. The computed high-pressure limit rate coefficients showed that n-propyl could be a minimum-size representative structure for larger n-alkyl chain, but this approach fails to capture the pressure dependence of the reaction that remains size-dependent for larger alkyls. The Geometric Mean Rule (GMR), a simple rule that links the self and cross-reactions of three alkyl radicals, is used to address this issue. Our theoretical calculations show for the first time that the GMR applies not only at the high-pressure limit, but also for every pressure explored. The GMR is shown to be accurate, with a maximum 22 % deviation in the falloff region for the estimation of rates for small alkyl combinations. Using this method, we propose a tabulated reaction rate rule for the combination of alkyl radicals up to C5+C5. The use of computed rate coefficients with the GMR approximation opens ways to develop pressure-dependent reaction rate rules for large n-alkyl combinations, based on accurate and computationally expensive theoretical calculations on smaller alkyl fragments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A review of the shear rheology of carbon black suspensions

The microstructural link to the rheology of carbon black suspensions has recently become clear as a result of advances in computational and experimental methods. This understanding reveals the important role of the restructuring, build-up, and break-up of carbon black agglomerates in simple shear, rationalized by a dimensionless balance of the hydrodynamic forces acting to break the agglomerates apart against the cohesive forces holding them together (i.e., the Mason number). The Mason number not only can predict the origin of reversible thixotropy seen in carbon black suspensions observed at higher shear intensities, but can also be used to rationalize the evolution of microstructure at lower shear intensities. This review focuses on carbon black suspension behavior, but the insights derived from carbon black suspensions are broadly applicable to a diverse class of soft matter including colloidal gels relevant to a variety of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advancing Thermal-Hydraulic Modeling Capabilities for Molten Salt Reactors in Pronghorn

This is an M2 report for the NEAMS program. This work focuses on advancing the modeling and simulation capabilities for Molten Salt Reactors (MSRs) within Pronghorn, addressing key thermal-hydraulic challenges that must be addressed for the safe and efficient operation of these advanced reactors. This report summarizes the development and validation of the following new models in Pronghorn: (i) two-equation RANS turbulence models, (ii) near-wall turbulence corrections due to curvature and roughness, and (iii) multidimensional two-phase drift-flux models. For each of these cases, along with the experimental validation, an application study of interest to MSRs is presented. In addition, this report also presents recent improvements in the domain-overlapping coupling between Pronghorn and SAM to address the cases of weakly compressible flow and flow reversal. Finally, an MSR action that greatly simplifies the setup of MSR models using NEAMS tools is described.

42 ENGINEERING↗

Strain Fluctuations Unlock Ferroelectricity in Wurtzites

Ferroelectrics are of practical interest for non-volatile data storage due to their reorientable, crystallographically defined polarization. Yet efforts to integrate conventional ferroelectrics into ultrathin memories have been frustrated by film-thickness limitations, which impede polarization reversal under low applied voltage. Wurtzite materials, including magnesium-substituted zinc oxide (Zn,Mg)O, have been shown to exhibit scalable ferroelectricity as thin films. In this work, the origins of ferroelectricity in (Zn,Mg)O are explained, showing that large strain fluctuations emerge locally in (Zn,Mg)O and can reduce local barriers to ferroelectric switching by more than 40%. Concurrent experimental and computational evidence of these effects are provided by demonstrating polarization switching in ZnO/(Zn,Mg)O/ZnO heterostructures featuring built-in interfacial strain gradients. These results open up an avenue to develop scalable ferroelectrics by controlling strain fluctuations atomistically.

36 MATERIALS SCIENCE↗

CO 2 Hydrogenation to Methanol over Inverse ZrO 2 /Cu(111) Catalysts: The Fate of Methoxy under Dry and Wet Conditions

Understanding the surface chemistry of CH 3 O species is essential for the production of methanol by CO 2 hydrogenation over Cu-based heterogeneous catalysts, as it facilitates the rational design of more efficient conversion processes. Recent research has identified inverse ZrO 2 /Cu catalysts as highly active and selective systems for the transformation of CO 2 to methanol with a performance that can be better than that of commercial Cu/ZnO catalysts. Here, we employed synchrotron-based ambient pressure X-ray photoelectron spectroscopy (AP-XPS) and calculations based on density functional theory (DFT) to understand the fate of CH 3 O groups under dry and wet environments. AP-XPS spectra revealed that under CO 2 hydrogenation conditions, formate and methoxy are two key intermediates to produce methanol. Furthermore, there are three different types of reactive sites on the surface: One is active for methoxy adsorption, which is stable and responsible for the methanol synthesis; Another one transforms CO 2 into CO; and a third one is active for CO 2 and methoxy dissociation, leading to C and methane formation. The theoretical calculations indicate that CH 3 OH readily dissociates to CH 3 O species following a highly exothermic (ΔE = -20.99 kcal/mol) and barrierless process. The water produced by the reverse water-gas shift reaction (CO 2 + H 2 → H 2 O + CO) can prevent the decomposition of CH 3 O species. We discovered that by introducing a tiny amount of water vapor (2 × 10 -6 Torr) into the reaction chamber, the energy barrier for the reaction CH 3 O(ads) + H(ads) → CH 3 OH(gas) is dramatically reduced. AP-XPS and computational modelling showed that water is quite capable of extracting adsorbed methoxy to form gaseous methanol. With this in mind, one could boost the methanol selectivity by adding appropriate amounts of water or steam, which is an inexpensive and feasible solution for industrial operations.

36 MATERIALS SCIENCE↗

Global Characteristics of Observable Foreshocks for Large Earthquakes

Abstract Foreshocks are the only currently widely identified precursory seismic behavior, yet their utility and even identifiability are problematic, in part because of extreme variation in behavior. Here, we establish some global trends that help identify the expected frequency of foreshocks as well the type of earthquake most prone to foreshocks. We establish these tendencies using the global earthquake catalog of the U.S. Geological Survey National Earthquake Information Center with a completeness level of magnitude 5 and mainshocks with Mw≥7.0. Foreshocks are identified using three clustering algorithms to address the challenge of distinguishing foreshocks from background activity. The methods give a range of 15%–43% of large mainshocks having at least one foreshock but a narrower range of 13%–26% having at least one foreshock with magnitude within two units of the mainshock magnitude. These observed global foreshock rates are similar to regional values for a completeness level of magnitude 3 using the same detection conditions. The foreshock sequences have distinctive characteristics with the global composite population b-values being lower for foreshocks than for aftershocks, an attribute that is also manifested in synthetic catalogs computed by epidemic-type aftershock sequences, which intrinsically involves only cascading processes. Focal mechanism similarity of foreshocks relative to mainshocks is more pronounced than for aftershocks. Despite these distinguishing characteristics of foreshock sequences, the conditions that promote high foreshock productivity are similar to those that promote high aftershock productivity. For instance, a modestly higher percentage of interplate mainshocks have foreshocks than intraplate mainshocks, and reverse faulting events slightly more commonly have foreshocks than normal or strike-slip-faulting mainshocks. The western circum-Pacific is prone to having slightly more foreshock activity than the eastern circum-Pacific.

Geochemistry & Geophysics↗

Miniaturized Magnetoelastic Sensor System

This article describes the design, assembly, and implementation of a hand-held, magnetic-field-based sensor system that can be adapted for a variety of sensing applications. The miniaturized system is based on Chemical Identification by Magneto-Elastic Sensing (ChIMES) technology, which uses three concentric solenoid coils to wirelessly interrogate a sensor body comprised of a response material coupled to a magnetoelastic wire. The response material expands when it encounters a target, imposing mechanical stress on the wire and altering its magnetic permeability. The sensor bodies are passive, requiring no external power source, and they are small, measuring about 15 mm in length and 3.0 mm in diameter. Up to four sensor bodies can be configured as an evenly-spaced linear array. The sensor system operates by applying a low-frequency, current-stabilized, filtered triangle wave to a uniform-density excitation coil to switch the magnetic domains within the wire. Further, the responses from the sensors are picked up by a detection coil as stress-induced changes in the Faraday voltage, and the strong magnetic field induced by the excitation coil in the detection coil is nullified by a cancellation coil reverse-wound in series with the detection coil. The responses of the sensors in an array are separated in time by a linear gradient dc biasing coil. The sensors can be interrogated through metallic and nonmetallic barriers. The signals from the detection coil and the excitation coil are digitized by a pair of bipolar analog-to-digital converters (ADCs). A Raspberry Pi single-board computer (SBC) and associated software perform data acquisition and control all aspects of the sensor system hardware. The program allows the user to select the number of sensors in the array, the type of signal that is being collected, and the number of samples to take. The program also allows for signal processing of the sensor data, such as baseline correction. The program can differentiate sensor peaks from each other and calculate the magnitude of each sensor response with less than 1% error. The data are then displayed along with a graph of the signal.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Programmable Dynamic Self-Assembly of DNA Nanostructures (Final Technical Report)

This technical report summarizes our research toward the synthesis of active biomolecular materials that take advantage of the biological properties of DNA and RNA, but can operate in non-biological contexts. Our approach takes inspiration from the dynamic assembly and disassembly of cytoskeletal filaments inside cells, and focuses on the generation of artificial filaments with comparable adaptation and responsiveness. We report the construction of DNA nanotube systems that assemble and disassemble reversibly in response to a variety of molecular stimuli including other nucleic acids (DNA and RNA), enzymes, and pH, as well as our progress on developing a reaction network framework to regulate self-assembly. Our DNA polymers are active in that they rely on the programmable energetics of hybridization and on the kinetics of enzymatic reactions to perform assembly and disassembly. Results include experiments and computational models. The research described here was performed between June 15, 2016 and June 15, 2019 at the University of California at Riverside with the support of the collaborative award DE-SC0010595 to PI Elisa Franco at UC Riverside and PI Rebecca Schulman at Johns Hopkins. The report focuses on the research performed by the team of PI Franco. The report also details work performed during a no-cost extension (June 15, 2019 - December 15, 2019) awarded to PI Guillermo Aguilar and PI Franco.

36 MATERIALS SCIENCE↗

Revolutionizing Turbine Cooling with Micro-Architectures Enabled by Direct Metal Laser Sintering

The objective of this research effort was to explore innovative cooling architectures enabled by additive manufacturing techniques for improved turbine cooling performance. The ability to create complex internal geometries was leveraged to better distribute coolant as well as to integrate inherently unstable flow devices to enhance internal and external heat transfer. This was accomplished with a multi-faceted approach including analytical, experimental, and computational components. This final report documents progress during the total 4.25 year effort that was extended (at no additional cost) from the original 3-year cooperative agreement. Prior to starting the project, significant effort was invested in investigating innovative cooling designs from the open literature. When the project started, we quickly down-selected to 4 configurations: impingement jet fluidic oscillators (for internal leading edge cooling), reverse film cooling (for pressure surface), sweeping film cooling jets (for suction surface), and trailing edge slot cooling with microchannels. During the 2nd year, the reverse film cooling design was shown to be sub-optimal and work on that topic was halted. The three remaining technologies were integrated into a large scale nozzle guide vane and installed into a low-speed linear cascade facility for further interrogation. Initial results showed the sweeping film cooling jets (on the suction surface) to be the most promising technology compared to a baseline diffusion shaped film hole (777 design). The benefit was particularly evident at high blowing ratios (>1.5) when the 777 coolant film separated from the downstream surface. Benefits were also evident at elevated freestream turbulence levels. The leading edge cooling with unsteady jets was less effective in terms of peak or average cooling – however it was superior to round hole impingement cooling in terms of spatial uniformity. Finally, the trailing edge cooling design with micro-channels was scrapped in favor of a pinned arrangement with centerbody. During the 3rd year, the same three technologies were integrated into a transonic linear cascade for an assessment of compressibility effects. Again, the sweeping film cooling jets proved superior at high blowing ratios while the unsteady leading edge impingement and trailing edge pin-fins with centerbody designs yielded mixed results. Finally, during the 4th and final year, these 3 technologies were integrated into a direct metal laser sintered (DMLS) nozzle guide vane for testing in a high temperature, transonic annular vane cascade. The facility matches the flow temperature, Mach number, and coolant pressure ratios of an actual gas turbine. Before doing this, an additional series of tests were completed to validate a redesigned trailing edge cooling architecture. This final design included a centerbody with triangular pins between the centerbody and the vane external skin. The majority of these triangular pins are fabricated with a 30% gap to the centerbody – reducing their pressure drop considerably while still providing excellent heat transfer augmentation. This additional test campaign required an additional 3-month extension request to complete testing in the high temperature NGV test facility (Turbine Reacting Flow Rig – TuRFR). Testing of the DMLS vane in TuRFR pitted the 3 innovative cooling technologies against more traditional technologies (shaped 777 film hole, round impingement jet, and full pin arrangement in the trailing edge. The sweeping film cooling jets (on the suction surface) showed approximately 15% improvement in overall effectiveness compared to the 777 film hole while the trailing edge showed up to 20% improvement. The leading edge impingement was disappointing with a substantially lower effectiveness than traditional direct impingement. In summary, the study was successful in demonstrating that DMLS-enabled cooling technologies can yield significant gains in cooling performance.

03 NATURAL GAS↗

Coupling Anionic Oxygen Redox with Selenium for Stable High‐Voltage Sodium Layered Oxide Cathodes

Utilizing anion redox reaction is crucial for developing the next generation of high-energy density, low-cost sodium-ion batteries. However, the irreversible oxygen redox reaction in Na-ion layered cathodes, which leads to voltage fading and reduced overall lifespan, has hindered their practical application. In this study, selenium is incorporated as a synergistic redox active center of oxygen to improve the stability of Na-ion cathodes. The redesigned cathode maintains stable voltage by demonstrating reversible oxygen redox while significantly suppressing the redox activity of manganese. The anionic redox contribution capacity of the selenium-doped Na 0.6 Li 0.2 Mn 0.8 O 2 cathode remains as high as 84% after 50 cycles, while the pristine Na 0.6 Li 0.2 Mn 0.8 O 2 cathode experiences a reduction to 39% of its initial capacity. The X-ray photoelectron spectroscopy data and computational analysis further revealed that selenium doping participates in redox as Se +4/5 which stabilizes the charged state and increases the energy step for O─O dimerization, thus improving the stability and lifespan of Na 0.6 Li 0.2 Mn 0.8 O 2 cathodes. In conclusion, the findings highlight the potential of redox coupling design to address the issue of voltage fade caused by irreversible anionic redox.

25 ENERGY STORAGE↗

Photoresponsive Rotaxanes Switch Lipid Bilayer Neuromorphic Behavior with Light

A rotaxane consisting of a macrocycle ring with two azobenzene units mechanically interlocked onto a bolaamphiphilic axle was incorporated into droplet interface bilayers (DIBs). The azobenzene groups on the ring underwent quasi-reversible, photoisomerization-induced cycling between 1-E and 1-Z configurations when irradiated with 370 and 467 nm light, respectively, enabling programmable access to different history-dependent electrical behaviors from the same membrane. In the 1-E configuration, bilayers exhibited type-IIactive memristance that coincided with increasingly elevated ionic conduction, associated with progressively enhanced bilayer permeability during voltage cycling. In the 1-Z configuration, bilayers displayed type-I, passive memcapacitive behavior, reflecting tighter lipid packing and reduced ionic permeability. Photoswitching also yielded a nonvolatile, photoresponsive memcapacitor that could be modulated repetitively with negligible loss, likely via reversible changes in membrane thickness. Concurrent ohmic leakage currents across the membrane were less than 0.3%. These results agree with previous studies with increasing membrane permeability using photoswitchable rotaxanes and provide new insights into the coupling between volatile and nonvolatile memcapacitance during photoisomerization. More broadly, they demonstrate a new strategy for the manipulation of neuromorphic behaviors in soft materials using light, with implications for brain-inspired computation and sensing.

droplet interface bilayers↗

Spy in the GPU-box: Covert and Side Channel Attacks on Multi-GPU System

The deep learning revolution has been enabled in large part by GPUs, and more recently accelerators, which make it possible to carry out computationally demanding training and inference in acceptable times. As the size of machine learning networks and workloads continues to increase, multi-GPU machines have emerged as an important platform offered on High Performance Computing and cloud data centers. Since these machines are shared among multiple users, it becomes increasingly important to protect applications against potential attacks. In this paper, we explore the vulnerability of Nvidia's DGX multi-GPU machines to covert and side channel attacks. These machines consist of a number of discrete GPUs that are interconnected through a combination of custom interconnect (NVLink) and PCIe connections. We reverse engineer the interconnected cache hierarchy and show that it is possible for an attacker on one GPU to cause contention on the L2 cache of another GPU. We use this observation to first develop a covert channel attack across two GPUs, achieving the best bandwidth of around 4 MB/s. We also develop a prime and probe attack on a remote GPU allowing an attacker to recover the cache access pattern of another workload. This access pattern can be used in any number of side channel attacks: we demonstrate a proof of concept attack that fingerprints the application running on the remote GPU, with high accuracy. We also develop a proof of concept attack to extract hyperparameters of a machine learning workload. Our work establishes for the first time the vulnerability of these machines to microarchitectural attacks and can guide future research to improve their security.

Dutta, Sankha↗

TBASS: A Robust Adaptation of Bayesian Adaptive Spline Surfaces

The R package TBASS is an extension of the BASS package created by Francom and Sansó (2019). The package is used to fit a Bayesian multivariate adaptive spline to a dataset that either follows a Student’s t-distribution or has outliers. Much of the framework for TBASS is adapted from the concepts of Bayesian Multivariate Adaptive Regression Splines (BMARS), specifically the work done by Denison, Mallick, and Smith (1998). The spline function is fit using a Reversible-Jump Markov Chain Monte Carlo algorithm,. By including this more robust generalization, a dataset with outliers can be accurately fit using the BMARS model, without the possibility of overfitting or variance inflation.

97 MATHEMATICS AND COMPUTING↗

Dehydroxylation kinetics of kaolinite and montmorillonite examined using isoconversional methods

The use of calcined clays as supplementary cementitious materials (SCMs) in concrete is a promising strategy towards decarbonizing the cement and concrete industry. This is especially relevant considering the ever-increasing demand for concrete. Comprehensive understanding of the kinetics of calcination is essential towards maximizing the potential reactivity of clay minerals while ensuring energy efficiency. In this study, the kinetics of the dehydroxylation of kaolinite and montmorillonite are investigated under non-isothermal conditions at constant heating rate. Activation energies ( E a ) are determined via Friedman differential and advanced Vyazovkin incremental methods over the isoconversional range; these are devoid of computational approximations, thus allowing kinetic analysis without assuming a specific reaction model. Kinetic equations—in the differential form as well as a combination of differential and integral forms are compared against the experimentally determined reaction models to identify the most probable dehydroxylation mechanism for kaolinite and montmorillonite. A reaction order mechanism is established for dehydroxylation of kaolinite, while montmorillonite is noted to undergo dehydroxylation via a single-step reversible diffusion-controlled process. Kinetic triplet—comprising activation energy, reaction model and pre-exponential factor—is used to predict isothermal calcination conditions, which is further verified using analytical techniques. Heat release rates of clay-portlandite blends from isothermal calorimetry are used within a thermodynamic framework to quantify reactivity of the calcined clays. Here, the study demonstrates a general approach based on isoconversional methods to predict calcination conditions for different clays that can be used in efficient and optimized production of blended cements or SCMs.

36 MATERIALS SCIENCE↗

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗

Energetic particle marginal stability profile for HL-2M integrated simulation based on neural network module

Abstract A critical gradient model is employed to develop a module of energetic particle (EP) marginal stability profiles in OMFIT integrated simulations for studying EP transport. Currently, each iteration of transport evolution is approximately 10 min in the integrated simulation, whereas, the EP marginal stability profile, which serves as an input in the integrated simulation could take much longer; the reason being a combination of the TGLFEP and EPtran codes is employed in our previous investigation. To reduce the simulation time, the critical gradient is predicted by a neural network instead of the TGLFEP code, and the EPtran code is revised with parallel computing, so that the running time of this module can be controlled to within 5 min. The predictions are in good agreement with previous approaches. The integrated simulation of HL-2M with Alfven eigenmodes transported by neutral beam EP profiles indicates that EP transport reduces the total pressure and current as expected, but could also under some conditions raise the safety factor in the core, which is favorable for reversed magnetic shear and high-performance plasmas.

Physics↗