Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled problems”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Low-cost, Highly Efficient and Fast Thermally Pressed Scalable Carbon-based Planar Perovskite Solar Cells

This Final Technical Report summarizes work completed under Award Number DEEE0009833, titled "Low-cost, Highly Efficient and Fast Thermally Pressed Scalable Carbon-based Planar Perovskite Solar Cells." The project addressed the central technical problem that carbon-electrode perovskite solar cells (C-PSCs) can reduce electrode costs and improve environmental durability, but they have historically lagged noble-metalcontact devices due to poor carbon/perovskite contact, high carbon-electrode resistance, imperfect hole extraction, and insufficient operational stability. The Statement of Project Objectives (SOPO), therefore, required systematic engineering of thermally pressed carbon-based planar perovskite solar cells by improving the carbon/perovskite interface quality, improving carbon-electrode conductivity, and further increasing device stability. The work advanced understanding of carbon-based perovskite photovoltaics by connecting device performance to three coupled mechanisms: interfacial recombination at the perovskite/hole-transport-layer/carbon stack; lateral and vertical transport losses in carbon contacts; and instability driven by moisture, heat, illumination, ion migration, and slow/variable manufacturing steps.

14 SOLAR ENERGY↗
DOE OSTI2026

Accelerating kinetic plasma simulations with machine-learning-generated initial conditions

Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.

Artificial neural networks↗
DOE OSTI2025

Time-dependent phenomena in correlated materials

Understanding time-dependent processes and light-matter interaction in strongly correlated materials, and the interplay between electronic, orbital, vibrational, and spin degrees of freedom, is a cornerstone of condensed matter. These mechanisms can be proven by measuring the response of the systems to time-dependent perturbations. The corresponding time scales are dictated by the way light couples to the different excitations, and how these excitations exchange energy and momentum. Our research advances our understanding of these processes, and the interpretation of different equilibrium and time-resolved spectroscopies. Our project encompasses two main themes: (i) developing and refining computational techniques to study non-equilibrium spectroscopies including non-perturbative effects and (ii) applications to non-equilibrium phenomena. We have developed a new computational approach that works directly in the time domain: by including all the degrees of freedom involved in the scattering process (e.g. electrons, photons, neutrons), we solve the time dependent problem: a faithful numerical simulation of the experiment. By measuring the energy and momentum of the outgoing particles, we can extract information about the energy and momentum absorbed by the system. Prior to our work, people attempting to model and calculate non-equilibrium spectral functions relied on a description of the scattering cross section based on a formulation in the frequency domain, a treatment that is extremely cumbersome and complex. Our technique works in and out of equilibrium and can reproduce spectra by several spectroscopic techniques, such as time-resolved photoemission, neutron scattering, Raman, X-ray spectroscopies (RIXS, Auger, XAS, XMCD), and, by not relying on analytical approximations, yields results that reveal novel overlooked transient mechanisms. These tools provide sorely needed intuition for understanding the phenomenology of strongly correlated materials and will help experimentalists in identifying signatures of relevant excitations in pump-probe experiments, such as those conducted in DOE supported facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗
DOE OSTI2024

Identifying a novel pathway for extracellular electron uptake in Methanosarcina barkeri using shotgun

The overarching goal of this project was to identify the genetic and/or enzymatic basis for extracellular electron uptake in Methanosarcina barkeri—an organism we previously characterized to be capable of this process using poised potential electrodes and a range of electrochemical techniques. In our first aim, we will utilized proteomic techniques that labels extracellular proteins with free amine groups with biotin and perform a selective purification of these proteins from the rest of the proteome. These were identified using mass spectrometry, providing us insight into the identity of the cell surface proteins, including redox active proteins present during extracellular electron uptake. We will also used an isobaric tagging technique (iTRAQ labeling) of proteomes from different experiment conditions to look at differential proteins expression and better understand the total pathway involved in extracellular electron uptake coupled to methanogenesis. In our final aim, we will investigated the electrochemical nature of electron uptake in co-cultures, compared to our pure culture data sets. We predict the process of electron uptake from electrodes replicates what goes on in nature when these organisms for co-cultures with electrogenic microbes such as Geobacter, and to test this hypothesis we have begun to probe co-cultures electrochemically during direct interspecies electron transfer. Our final aim took much longer than expected due to challenges replicating co-cultures between Geobacter metallireducens and Methanosarcina barkeri. We found the major problem with this work was the inability of the M. barkeri strain we used (and we tried strains from three different labs) to successfully grow on acetate in the absence of hydrogen. As such, we began to use the strain Methanosarcina acetivorans—another close relative that specializes on acetate growth. We are currently in the process of comparing electrochemical observations of this microorganism between pure and co-culture experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗
DOE OSTI2024

A mechanistic interpretation of Nelson curves for PVP failures under high temperature hydrogen attack

As an empirically established design criterion, Nelson curves that relate the service temperature and the allowable hydrogen partial pressure have been developed and utilized for more than sixty years in pressure vessels and piping (PVP) safety design. Despite a relatively clear thermodynamic understanding of the high-temperature-hydrogen-attack (HTHA) problem, the detailed fracture process on the microstructural length scales, however, remains elusive, and a quantitative assessment of the PVP lifetime under HTHA from the available creep fracture dataset is still not possible. This work develops a microstructure-informed and micromechanics-based model by incorporating a synergy between hydrogen transport and intergranular-cavity-based fracture process. Based on the available creep lifetime data of C-0.5Mo steels, we are able to calibrate material constitutive parameters, and then conduct nonlinear finite element simulations that reveal a real-time stress-induced hydrogen diffusional transport along grain boundaries, coupled with a microstructure-explicit failure process, from which Nelson curves can be computed. Such failure analyses allow us to delineate two distinct regimes on the Nelson curves, i.e., dislocation-creep-controlled or grain boundary diffusion-assisted cavity growth. More importantly, we found that a small change of the pipe thickness and applied stresses can significantly shift these lifetime curves. However, these two parameters are usually not provided in Nelson curves, thus limiting their usage in material selection and safety design. In conclusion, this discrepancy can clearly be mitigated by extensive parametric studies from our micromechanical modeling/simulation framework.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Erbium-Implanted WS 2 Flakes with Room-Temperature Photon Emission at Telecom Wavelengths

Optically addressable spin impurities in crystals along with device engineering provide an attractive route to realizing quantum technologies in the solid state, but reconciling disparate emitter and host material constraints for a given target application is often challenging. Rare-earth ions in two-dimensional (2D) materials could mitigate this problem given the atomic-like transitions of the emitters and the versatile nature of van der Waals systems. Here we combine ion implantation, confocal microscopy, and ab initio calculations to examine the photon emission of Er-doped WS 2 flakes. Optical spectroscopy reveals narrow, long-lived photoluminescence lines in the telecom band, which we activate after low-temperature thermal annealing. Spectroscopic and polarization-selective measurements show a uniform response across the ensemble, while the fluorescence brightness remains mostly unchanged with temperature, suggesting nonradiative relaxation channels are inefficient. Our results create opportunities for novel solid state devices coupling 2D-hosted, telecom-band emitters to photonic heterostructures separately optimized for photon manipulation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

RANGE: A robust adaptive nature-inspired global explorer of potential energy surfaces

With the growing demand for realistic representations of chemical structures and the advent of exascale computing, the intelligent sampling of potential energy surfaces and efficient identification of global minima have become more essential but also more feasible. Building on prior studies demonstrating the efficiency of the Artificial Bee Colony (ABC) swarm intelligence algorithm, we report a hybrid metaheuristic framework that integrates the adaptive exploration capabilities of ABC coupled with the exploitation strengths of genetic algorithms (GA) in a scalable, Python-based implementation. The resulting tool, RANGE (Robust Adaptive Nature-inspired Global Explorer), provides seamless interfaces to multiple potential energy evaluators, either directly or via widely used Python libraries, and is designed for high-performance computing environments. We describe the implementation details of RANGE and evaluate its performance, relative to ABC- or GA-alone based algorithms, on a variety of chemical systems, including molecular clusters and heterogeneous surfaces. In conclusion, our results demonstrate RANGE’s efficiency, robustness, and broad applicability in addressing challenging global optimization problems in computational chemistry and materials science.

Algorithms and data structure↗
DOE OSTI2025

Proton internal pressure from deeply virtual Compton scattering on collider kinematics

Abstract The unique experimental connection to the QCD energy–momentum tensor offered by generalised parton distributions has been strongly highlighted in the past few years with attempts to extract the pressure and shear forces distributions within the nucleon. If, in principle, this can be performed in a model independent way from experimental data, in practice, the current limited precision and kinematic coverage make such an extraction very challenging. Moreover, the limitation to a leading-order description in the strong coupling of the data has provided only an indirect and weakly sensitive access to gluon degrees of freedom, solely through their mixing to quarks via evolution. In this paper we address this issue by providing a next-to-leading order formalism allowing a reanalysis of global fits with genuine gluonic degrees of freedom. In addition, we provide an estimate of the reduction in uncertainty that could stem from the extended kinematic range relevant for the future Electron Ion Collider. Finally, we stress the connection between the analysis of the dispersion relation in terms of generalised parton distributions and the deconvolution problem.

Dutrieux, H. (ORCID:0000000183344885)↗
DOE OSTI2025

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Spontaneous symmetry breaking, gauge hierarchy, and electroweak vacuum metastability

The so-called metastability bound asserts that an unnaturally small Higgs mass is a necessary condition for electroweak vacuum metastability, offering a new approach toward solving the hierarchy problem. So far, this result relies on the assumption of a negative Higgs mass parameter, or equivalently, on electroweak spontaneous symmetry breaking. We derive a new, corresponding bound for the case of a positive mass parameter. When the new bound is significantly more restrictive than or comparable to its established counterpart, it may offer an explanation for the sign of the Higgs mass parameter, and thus, spontaneous symmetry breaking itself. New physics at scales O ( 1 – 10 ) TeV can lower these bounds as far as the TeV scale. As an illustration, we consider vacuum stability in the presence of additional TeV-scale fermions with Yukawa couplings to the Higgs, as well as a dimension-six term parametrizing new physics in the UV. This scenario requires new physics that couples strongly to the Higgs, and can potentially be probed at future colliders. Finally, to allow for comparison with concrete mechanisms predicting metastability, we provide the mass-dependent lifetime of the electroweak vacuum for this model. Published by the American Physical Society 2024

Benevedes, Sean (ORCID:000000016390198X)↗
DOE OSTI2026

Simultaneous Overexpression of FERULOYL‐CoA 6′‐HYDROXYLASE 1 and COUMARIN SYNTHASE Leads to Coumarin‐Enriched Lignin and Improved Saccharification in Greenhouse‐ and Field‐Grown Poplar

ABSTRACT The urgent need for renewable resources has increased the interest in woody biomass to manufacture bio‐based products. However, lignin recalcitrance limits the enzymatic conversion of wood into fermentable sugars, posing a major challenge for biomass deconstruction. To address this problem, we aimed at incorporating the coumarin scopoletin into the lignin polymer of poplar ( Populus tremula × P . alba ) by expressing FERULOYL‐CoA 6′‐HYDROXYLASE 1 ( F6′H1 ) and COUMARIN SYNTHASE ( COSY ) in lignifying cells. Three constructs were evaluated: two bicistronic constructs, SCOP1 ( COSY followed by F6′H1 ) and SCOP2 ( F6′H1 followed by COSY ), and one monocistronic, SCOP3 (only F6′H1 ). SCOP1 poplars produced most free scopoletin without altering overall lignin, cellulose or hemicellulose content. SCOP2 poplars were overall less efficient in scopoletin production and most of these lines showed a severe biomass yield penalty, whereas SCOP3 caused plant lethality. NMR and metabolic analyses confirmed that scopoletin cross‐coupled with G and S monomers during lignification in SCOP1 lines. In addition to scopoletin, the detection of benzodioxane structures revealed the incorporation of dihydroxycoumarins. Overall coumarin incorporation in lignin amounted up to 2.3%. After alkaline pretreatment, wood from greenhouse‐grown SCOP1 poplars released up to 29% more glucose compared to the wild type upon limited saccharification. Field‐testing of three SCOP1 lines showed a 6 to 11% increase in saccharification efficiency, with the line containing the lowest scopoletin levels maintaining normal growth. These results demonstrate that engineering lignin composition in poplar can improve saccharification, and emphasize the importance of construct design, translational research and field validation.

alternative lignin monomers↗
DOE OSTI2026

Mechanistic Insights into Dinitrogen Reduction to Ammonia in Light-Controlled Nanocrystal:Nitrogenase Complexes

Developing systems that can efficiently capture photon energy and convert this energy into fuels and chemicals requires understanding how to assemble molecular components with diverse functions into complete systems possessing selectivity and efficiency in directing charge carriers to catalytic reactions. There are many challenges to achieving this goal. One promising approach is the development of hybrid systems that combine semiconductor nanocrystals (NCs) for light capture and enzymes as efficient catalysts. Such biohybrid systems capitalize on the tunable electronic and optical properties of NCs while leveraging the unmatched specificity and efficiency of enzymes in catalyzing chemical reactions, thereby offering opportunities to surpass the limitations of each component alone. Here, we focus on recent progress in developing a biohybrid system that combines CdS NCs for photon capture with the enzyme nitrogenase to accomplish light-driven dinitrogen (N 2 ) reduction to ammonia (NH 3 ). Integrating light-harvesting materials with biological catalysts requires a deep understanding of NC properties, protein stability, and electron transfer (ET), making it an inherently multidisciplinary problem. The reduction of N 2 to NH 3 is a challenging reaction, with a high demand in both agriculture and industrial chemical production. This reaction is intrinsically energy intensive, due to the need to activate the N≡N triple bond. The current standard industrial approach to N 2 reduction, the Haber−Bosch reaction, obtains the necessary energy input from fossil fuels, whereas biological systems capable of N 2 reduction utilize the hydrolysis of ATP as their energy source. Replacing these costly, energy-intensive inputs with renewable light energy represents a critical step toward sustainable NH 3 production. Recent progress has demonstrated that semiconductor CdS NCs can be coupled to the catalytic component of nitrogenase, the MoFe protein, to form a biohybrid CdS NC:MoFe protein complex, enabling light-driven N 2 reduction rather than energy input from fossil fuels or ATP. This illustrates how inorganic NCs can functionally replace the natural Fe protein partner, yielding a biohybrid catalyst that enables controlled electron delivery and provides not only light-driven NH 3 production but also new approaches for probing enzyme catalytic function. The CdS NC:MoFe protein biohybrid system enables light-initiated electron delivery at ambient temperature, as well as temperatures below freezing, allowing for stabilization and spectroscopic characterization of key reaction intermediates. These findings highlight how photochemical biohybrids can serve as both functional catalysts and mechanistic probes. Beyond studies of the nitrogenase mechanism, studies of the CdS NC:MoFe system reveal how variables such as NC size, electrostatic binding interactions, and sacrificial electron donors (SEDs) govern complex stability, charge transfer efficiency, and catalytic performance. In addition, studies of nitrogenase and the high activation barrier for N 2 reduction are enabling investigations of new and interesting questions regarding the properties and limitations of NC biocatalysis. In this Account, we describe the key features of CdS NC:MoFe protein biohybrids and the parameters for optimal light-driven N 2 reduction, and how controlling ET with light illuminates the path to new insights into the nitrogenase mechanism.

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

Fluctuations in Hill’s equation parameters and application to cosmic reheating

Cosmic inflation provides a compelling framework for explaining several observed features of our Universe, but its viability depends on an efficient reheating phase that converts the inflaton’s energy into Standard Model particles. This conversion often proceeds through nonperturbative mechanisms such as parametric resonance, which is described by Hill’s equation. In this work, we investigate how stochastic fluctuations in the parameters of Hill’s equation can influence particle production during reheating. We show that such fluctuations can arise from couplings to light scalar fields and can significantly alter the stability bands in the resonance structure, thereby enhancing the growth of fluctuations and broadening the region of efficient energy transfer. Using random matrix theory and stochastic differential equations, we decompose the particle growth rate into deterministic and noise-induced components and demonstrate analytically and numerically that even modest noise leads to substantial particle production in otherwise stable regimes. Furthermore, these results suggest that stochastic effects can robustly enhance the efficacy of reheating across a wide swath of parameter space, with implications for early Universe cosmology, UV completions involving multiple scalar fields, and the resolution of the cosmological moduli problem.

Cosmology↗
DOE OSTI2026

Monte Carlo Event Generation with Continuous Normalizing Flows

We apply continuous normalizing flows trained with the flow matching method to the problem of phase-space sampling in Monte Carlo event generation for high-energy collider physics. Focusing on lepton-pair and top-quark pair production with multiple jets, the two computationally most expensive processes at the Large Hadron Collider, we train helicity-conditioned continuous normalizing flows to remap the random numbers used in matrix element evaluation. Compared to standard methods, we achieve unweighting efficiency improvements by factors of up to 184 and 25 for the two processes at their respective highest jet number, at the cost of an increased evaluation time. When combining the advantages of continuous normalizing flows with the fast evaluation times of coupling-layer-based flows, using the RegFlow approach, we find parton-level unweighted event generation walltime gains of about a factor of 10 at the highest jet numbers. These substantial gains highlight the promise of samplers based on machine learning for next-generation collider experiments.

Bothmann, Enrico [CERN; Gottingen U.] (ORCID:00000↗
DOE OSTI2024

Direct microstability optimization of stellarator devices

Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks and stellarators. Here, in this work, we show that a significant decrease in a microstability-based proxy, as opposed to a geometric one, for the turbulent heat flux, namely the quasilinear heat flux, can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasilinear heat flux at each step of the optimization process, as well as the deviation from quasisymmetry, and minimizing their sum, leading to a balance between neoclassical and the turbulent transport proxy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗