In this study, accurate prediction of fusion performance in present and future tokamaks requires taking into account the strong interplay between core transport, pedestal structure, current profile and plasma equilibrium. An integrated modeling workflow capable of calculating the steady-state self- consistent solution to this strongly-coupled problem has been developed. The workflow leverages state-of-the-art components for collisional and turbulent core transport, equilibrium and pedestal stability. Validation against DIII-D discharges shows that the workflow is capable of robustly pre- dicting the kinetic profiles (electron and ion temperature and electron density) from the axis to the separatrix in good agreement with the experiments. An example application is presented, showing self-consistent optimization for the fusion performance of the 15 MA D-T ITER baseline scenario as functions of the pedestal density and ion effective charge Z eff.

Accurate prediction of fusion performance in present and future tokamaks requires taking into account the strong interplay between core transport, pedestal structure, current profile and plasma equilibrium. An integrated modeling workflow capable of calculating the steady-state selfconsistent solution to this strongly-coupled problem has been developed. The workflow leverages state-of-the-art components for collisional and turbulent core transport, equilibrium and pedestal stability. Validation against DIII-D discharges shows that the workflow is capable of robustly predicting the kinetic profiles (electron and ion temperature and electron density) from the axis to the separatrix in good agreement with the experiments. An example application is presented,more » showing self-consistent optimization for the fusion performance of the 15 MA D-T ITER baseline scenario as functions of the pedestal density and ion effective charge Z eff.« less

Fusion whole device modeling simulations require comprehensive models that are simultaneously physically accurate, fast, robust, and predictive. In this paper we describe the development of two neural-network (NN) based models as a means to perform a snon-linear multivariate regression of theory-based models for the core turbulent transport fluxes, and the pedestal structure. Specifically, we find that a NN-based approach can be used to consistently reproduce the results of the TGLF and EPED1 theory-based models over a broad range of plasma regimes, and with a computational speedup of several orders of magnitudes. These models are then integrated into a predictive workflowmore » that allows prediction with self-consistent core-pedestal coupling of the kinetic profiles within the last closed flux surface of the plasma. Finally, the NN paradigm is capable of breaking the speed-accuracy trade-off that is expected of traditional numerical physics models, and can provide the missing link towards self-consistent coupled core-pedestal whole device modeling simulations that are physically accurate and yet take only seconds to run.« less

In this study, a self-consistent model of the interstellar pickup protons, the slab component of the Alfvénic turbulence, and core solar wind (SW) protons is presented for r ≥ 1 along with the initial results of and comparison with the Voyager 2 (V2) observations. Two kinetic equations are used for the pickup proton distribution and Alfvénic power spectral density, and a third equation governs SW temperature including source due to the Alfvén wave energy dissipation. A fraction of the pickup proton free energy, fD , which is actually released in the waveform during isotropization, is taken from the quasi-linear considerationmore » without preexisting turbulence, whereas we use observations to specify the strength of the large-scale driving, C sh, for turbulence. The main conclusions of our study can be summarized as follows. (1) For C sh ≈ 1-1.5 and f D ≈ 0.7-1, the model slab component agrees well with the V2 observations of the total transverse magnetic fluctuations starting from ~8 AU. This indicates that the slab component at low-latitudes makes up a majority of the transverse magnetic fluctuations beyond 8-10 AU. (2) The model core SW temperature agrees well with the V2 observations for r ≳ 20 AU if f D ≈ 0.7-1. (3) A combined effect of the Wentzel-Kramers-Brillouin attenuation, large-scale driving, and pickup proton generated waves results in the energy sink in the region r ≲ 10 AU, while wave energy is pumped in the turbulence beyond 10 AU. Without energy pumping, the nonlinear energy cascade is suppressed for r ≲ 10 AU, supplying only a small energy fraction into the k-region of dissipation by the core SW protons. A similar situation takes place for the two-dimensional turbulence. (4) The energy source due to the resonant Alfvén wave damping by the core SW protons is small at heliocentric distances r ≲ 10 AU for both the slab and the two-dimensional turbulent components. As a result, adiabatic cooling mostly controls the model SW temperature in this region, and the model temperature disagrees with the V2 observations in the region r ≲ 20 AU.« less

We demonstrate a self-consistent numerical scheme for simulating an electronic device which contains active defects. As a specific case, we consider copper defects in cadmium telluride solar cells. The presence of copper has been shown experimentally to play a crucial role in predicting device performance. The primary source of this copper is migration away from the back contact during annealing, which likely occurs predominantly along grain boundaries. We introduce a mathematical scheme for simulating this effect in 2D and explain the numerical implementation of the system. Lastly, we will give numerical results comparing our results to known 1D simulations tomore » demonstrate the accuracy of the solver and then show results unique to the 2D case.« less