Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability properties”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simultaneous access to high normalized density, current, pressure, and confinement in strongly-shaped diverted negative triangularity plasmas

Strongly-shaped diverted negative triangularity (NT) plasmas in the DIII-D tokamak demonstrate simultaneous access to high normalized density, current, pressure, and confinement. NT plasmas are shown to exist across an expansive parameter space compatible with high fusion power production, revealing surprisingly good core stability properties that compare favorably to conventional positive triangularity plasmas in DIII-D. Non-dimensionalizing the key parameters, expanded operating spaces featuring edge safety factors below 3, normalized betas above 3, Greenwald density fractions above 1, and high-confinement mode (H-mode) confinement qualities above 1 are observed, even simultaneously, and all with a robustly stable edge free from deleterious edge-localized mode instabilities. Scaling of the confinement time with engineering parameters reveals at least a linear dependence on plasma current although with significant power degradation, both in excess of expected H-mode scalings. These results increase confidence that NT plasmas are a viable approach to realize fusion power and open directions for future detailed study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Stability of trapped solutions of a nonlinear Schrödinger equation with a nonlocal nonlinear self-interaction potential

This work focuses on the study of the stability of trapped soliton-like solutions of a (1 + 1)-dimensional nonlinear Schrödinger equation (NLSE) in a nonlocal, nonlinear, self-interaction potential of the form [|Ψ(x,t)| 2 +|Ψ(-x,t)| 2 ] κ where κ is an arbitrary nonlinearity parameter. Although the system with κ = 1 (i.e. fully integrable case) was first reported by Yang (2018 Phys. Rev. E 98 042202), here in the present work, we extend this model to the one in which κ is arbitrary. This allows us to compare the stability properties of the now trapped solutions to previously found solutions of the more usual NLSE with κ ≠ 1 which are moving soliton solutions. We show that there is a simple, one-component, nonlocal Lagrangian and corresponding action governing the dynamics of the system. Using a collective coordinate method derived from the action as well as assuming the validity of Derrick's theorem, we find that these trapped solutions are stable for 0 < κ < 2 and unstable when κ > 2. At the critical value of κ, i.e. κ = 2, the solution can either collapse or blowup linearly in time when q 0 = 0, where q 0 is the center of the initial density ρ(x, t = 0) = ψ*ψ of the solution. For q 0 ≠ 0 the displaced solution collapses. When κ > 2 initial small displacements from the origin also lead to collapse of the wave function. This phenomenon is not seen in the usual NLSE.

collective coordinates↗
NASA NTRS1990

On the Goertler instability in hypersonic flows: Sutherland law fluids and real gas effects

The Goertler vortex instability mechanism in a hypersonic boundary layer on a curved wall is investigated. The precise roles of the effects of boundary layer growth, wall cooling, and gas dissociation is clarified in the determination of stability properties. It is first assumed that the fluid is an ideal gas with viscosity given by Sutherland's law. It is shown that when the free stream Mach number M is large, the boundary layer divides into two sublayers: a wall layer of O(M sup 3/2) thickness over which the basic state temperature is O(M squared) and a temperature adjustment layer of O(1) thickness over which the basic state temperature decreases monotonically to its free stream value. Goertler vortices which have wavelengths comparable with the boundary layer thickness are referred to as wall modes. It is shown that their downstream evolution is governed by a set of parabolic partial differential equations and that they have the usual features of Goertler vortices in incompressible boundary layers. As the local wavenumber increases, the neutral Goertler number decreases and the center of vortex activity moves towards the temperature adjustment layer. Goertler vortices with wavenumbers of order one or larger must necessarily be trapped in the temperature adjustment layer and it is this mode which is most dangerous. For this mode, it was found that the leading order term in the Goertler number expansion is independent of the wavenumber and is due to the curvature of the basic state. This term is also the asymptotic limit of the neutral Goertler numbers of the wall mode. To determine the higher order corrections terms in the Goertler number expansion, two wall curvature cases are distinguished. Real gas effects were investigated by assuming that the fluid is an ideal dissociating gas. It was found that both gas dissociation and wall cooling are destabilizing for the mode trapped in the temperature adjustment layer, but for the wall mode trapped near the wall the effect of gas dissociation can be either destabilizing or stabilizing.

Fu, Yibin B.↗
NASA NTRS2012

Linear Estimation of Particle Bulk Parameters from Multi-Wavelength Lidar Measurements

An algorithm for linear estimation of aerosol bulk properties such as particle volume, effective radius and complex refractive index from multiwavelength lidar measurements is presented. The approach uses the fact that the total aerosol concentration can well be approximated as a linear combination of aerosol characteristics measured by multiwavelength lidar. Therefore, the aerosol concentration can be estimated from lidar measurements without the need to derive the size distribution, which entails more sophisticated procedures. The definition of the coefficients required for the linear estimates is based on an expansion of the particle size distribution in terms of the measurement kernels. Once the coefficients are established, the approach permits fast retrieval of aerosol bulk properties when compared with the full regularization technique. In addition, the straightforward estimation of bulk properties stabilizes the inversion making it more resistant to noise in the optical data. Numerical tests demonstrate that for data sets containing three aerosol backscattering and two extinction coefficients (so called 3 + 2 ) the uncertainties in the retrieval of particle volume and surface area are below 45% when input data random uncertainties are below 20 %. Moreover, using linear estimates allows reliable retrievals even when the number of input data is reduced. To evaluate the approach, the results obtained using this technique are compared with those based on the previously developed full inversion scheme that relies on the regularization procedure. Both techniques were applied to the data measured by multiwavelength lidar at NASA/GSFC. The results obtained with both methods using the same observations are in good agreement. At the same time, the high speed of the retrieval using linear estimates makes the method preferable for generating aerosol information from extended lidar observations. To demonstrate the efficiency of the method, an extended time series of observations acquired in Turkey in May 2010 was processed using the linear estimates technique permitting, for what we believe to be the first time, temporal-height distributions of particle parameters.

linear estimation↗
DOE OSTI2021

Characteristics of Ultralow Aspect Ratio Toroidal Plasmas and Non-Solenoidal Startup and Edge Stability Studies at Near-Unity Aspect Ratio (Final Report)

This grant supported a long-term research and development program to explore the physics and technology of spherical tokamak (ST) plasma behavior at near-unity aspect ratio ( A ). The PEGASUS Toroidal Experiment is a university-scale ( R ) ~ 0.40 m, α ~ .35 m, I p ≤ 0.3 MA, B T ≤ 0.15 T) experiment that accesses A → 1, to provide large values of the normalized plasma current, I N = I p /α B T , or equivalently the toroidal field utilization factor I p / I TF . The research activities focused on: 1) developing the technique of Local Helicity Injection (LHI) as a viable means of initiating and growing an ST plasma without the use of a central solenoid; 2) exploring the equilibrium and stability properties of spherical tokamak plasmas at extremely low aspect ratio; and 3) obtaining new insights in the properties of the enhanced confinement H-mode regime at A ~ 1. The initial (Phase I) facility employed a unique high-stress central solenoid Ohmic for current drive and heating with simple magnet current waveforms. A subsequent major upgrade (Phase II) offered more flexible control of the plasma evolution along with a powerful LHI system for plasma initiation, heating and sustainment. The LHI technique consists of injection of high-density current streams into vacuum toroidal and vertical magnetic fields. These streams become unstable and support magnetic reconnection to allow the system to relax into a tokamak-like magnetic configuration. This program focused on: 1) development of electron current sources that provide sufficiently high injection currents in the challenging edge plasma region; 2) validation of theoretical limits to achievable I p arising from magnetic helicity and energy conservation plus absolute limits arising from relaxation to a minimal energy states; 3) tests of a power balance model to describe the evolution of the plasma current I p (t); and 4) study of magnetic fluctuations and their role in the I p evolution. This culminated in the production of toroidal currents of I p ≥ 0.2 MA with only ~8 kA of injected current. Operation with Ohmic induction alone in Phase I created plasmas with soft limits on the achievable I p and field utilization ( I p / I TF ~ 1) due to the onset of large low-order (2/1 or 3/2) internal resonant tearing modes. With the expanded capabilities in Phase II operation, methods employed to mitigate those modes and achieve I p / I TF > 1 were: 1) I TF ramp-downs; 2) non-solenoidal startup via LHI; and 3) extremely low I TF Ohmic operations assisted by LHI current injector pre-ionization. Extensive use of LHI-driven startup and drive for non-solenoidal operation produced favorable current profiles to stably operate tokamak plasmas with unprecedentedly high I N = I p / aB T ~ 15 and I p / I TF > 2. With strong reconnection-driven anomalous ion heating, this in turn provided access to tokamak plasmas with B t up to 100%, a world record 2 times higher than that previously achieved in advanced or spherical tokamaks. This high B t regime had a large fraction (~50%) of the plasma volume in a region of an absolute minimum- B configuration. The MHD stability limit in this regime was the ideal external kink mode without a stabilizing wall. These results provided the first demonstration of access to tokamak plasmas with very high B t , high i N , high elongation κ, and low internal inductance ℓ i , and as such served as a proof-ofprinciple for access to this long-sought ST regime. The naturally low B T needed for stability at A ≤ 1.2 reduces the L-H mode power threshold, P LH , so that H-mode plasmas were achieved with Ohmic heating alone. The power threshold for the transition to the H-mode regime was found to be 10× the accepted values from the ITER international scaling, confirming initial trends indicated by experiments on MAST and NSTX, and suggesting that the physics behind this transition is not yet fully understood. The density dependence of the transition and the observed insensitivity to whether the plasma is bounded by a limiter or magnetic divertor are consistent with the recent FM 3 model. Multiple unstable toroidal magnetic modes were measured during ELM (Edge Localized Mode) crashes, with generally lower toroidal mode numbers than seen at higher A , which was attributed to the larger peeling mode instability drive at low A . First-ever edge current profile measurements through an ELM crash showed current-carrying filaments were generated, as seen in ELM simulations and in electromagnetic blob transport theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Implications of carbon content for the processing, stability, and mechanical properties of cast and wrought Ni-based superalloy Nimonic 105

As the temperature and pressure requirements in land-based turbines for power generation increase, the shift towards more advanced alloys, e.g., Ni-based superalloys, requires improved understanding of their long-term stability and mechanical properties. Additionally, the larger size of the components also presents fabrication and cost-reduction problems. In this study, we investigated Nimonic 105 as a potential rotor material at two carbon content levels – the commonly used maximum of the alloy specification and a lower carbon content variant. Reducing the carbon content presented fabrication challenges, resulting in surface crack formation during hot working. Additionally, a lower M 6 C carbide fraction was present after fabrication/solutioning, leading to a significantly larger grain size. After the standard aging treatment, similar γ' populations formed in the two variants, however the M 6 C carbides in the high-carbon variant already started transforming to intergranular M 23 C 6 . Although we did not observe a noticeable effect on the tensile properties of the two variants, the lower carbon content alloy possessed a superior creep property – mainly due to the larger grain size as the microstructural stability was inferior. Further, upon prolonged thermal exposure, the M 6 C carbides in the low-carbon variant further decomposed to intergranular and intragranular σ and μ precipitates, as opposed to mainly M 23 C 6 as in the high-carbon variant.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Production of High Yield, Ultra‐Shallow Nitrogen‐Vacancy Centre Ensembles Through High Purity Forming gas Annealing

Dense, ultra‐shallow (<10 nm from the surface) ensembles of nitrogen‐vacancy (NV) centers in diamond are sought‐after for realising practical quantum sensing at the nanoscale. However, these ensembles typically suffer from reduced NV creation efficiency and spin coherence compared to their deeper and bulk counterparts. Consequently, most demonstrations in the literature to date have been restricted to using single NV centers with atypical coherence and charge stability properties not achievable across entire ensembles. Here, in this work, we demonstrate that controlled annealing in a high purity forming gas atmosphere improves near‐surface NV formation efficiency to over 20%. This environment preserves a high quality hydrogen termination throughout the anneal, avoiding material etching over many hours at temperatures up to 1100°C and possibly helping to suppress vacancy out‐diffusion. We show that high temperature annealing is especially vital for improving spin properties in high density, ultra‐shallow implanted samples due to the localised damage profiles produced. By allowing the formation of high yield NV ensembles with improved spin properties, the optimised forming gas anneal can facilitate the reliable production of more sensitive diamond sensors than have previously been accessible.

Healey, Alexander J. [Royal Melbourne Institute of↗
DOE OSTI2023

Epitaxial Design of Complex Nickelates as Electrocatalysts for the Oxygen Evolution Reaction

Abstract The oxygen evolution reaction (OER) is a crucial process in electrochemical water splitting, a promising technology to renewably yield hydrogen gas from water. Designing and developing earth‐abundant, efficient, and stable OER electrocatalysts to replace the most widely used but scarce RuO 2 and IrO 2 are thus of critical interest. Recently, ABO 3 ‐structured perovskite oxides, especially rare‐earth nickelates, are extensively studied for their potential use as OER electrocatalysts. In particular, the epitaxial synthesis of complex oxide thin films allows flexible and precise control over the materials so that their structure–stability–property relationships can be established. Using nickelate thin films as model systems, this review illustrates how epitaxial design allows researchers to test different hypotheses and proposed descriptors, as well as formulate new design principles. Following a brief introduction to the background of OER mechanisms, proposed activity descriptors, and synthesis methods, various epitaxial design strategies are surveyed including strain tuning, composition control, surface termination/orientation selection, defect engineering, and interface design. These have led to precise control over the atomic structures and electronic properties of nickelates which in turn determine their electrochemical performance. Finally, the remaining challenges and perspectives toward a deeper understanding and use of complex oxides as OER catalysts are discussed.

Choi, Min‐Ju↗
DOE OSTI2021

Numerical analyses of exponential time-differencing schemes for the solution of atmospheric models

In high-resolution numerical weather prediction models, fast-moving acoustic waves must be treated in a stable manner. Among the implicit–explicit (IMEX) schemes used for the solution of these models, horizontally explicit, vertically implicit (HEVI) methods show good stability and scalability on massively parallel machines. In this work, we present two classes of exponential time-differencing (ETD) methods for atmospheric models that use a HEVI splitting strategy, one being a three-stage method with the addition of artificial diffusion, the second based on a Strang splitting approach. Overall, the stability properties of the methods are analyzed and numerical examples are provided, which compare time-step restrictions and cost-to-accuracy ratios of the new methods with those for existing approaches.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Entropy–Preserving and Entropy–Stable Relaxation IMEX and Multirate Time–Stepping Methods

In this work, we propose entropy-preserving and entropy-stable partitioned Runge–Kutta (RK) methods. In particular, we extend the explicit relaxation Runge–Kutta methods to IMEX–RK methods and a class of explicit second-order multirate methods for stiff problems arising from scale-separable or grid-induced stiffness in a system. The proposed approaches not only mitigate system stiffness but also fully support entropy-preserving and entropy-stability properties at a discrete level. The key idea of the relaxation approach is to adjust the step completion with a relaxation parameter so that the time-adjusted solution satisfies the entropy condition at a discrete level. The relaxation parameter is computed by solving a scalar nonlinear equation at each timestep in general; however, as for a quadratic entropy function, we theoretically derive the explicit form of the relaxation parameter and numerically confirm that the relaxation parameter works the Burgers equation. Several numerical results for ordinary differential equations and the Burgers equation are presented to demonstrate the entropy-conserving/stable behavior of these methods. We also compare the relaxation approach and the incremental direction technique for the Burgers equation with and without a limiter in the presence of shocks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Implicit extensions of an explicit multirate Runge–Kutta scheme

In this work, we propose a new method that extends conservative explicit multirate methods to implicit explicit-multirate methods. We develop extensions of order one and two with different stability properties on the implicit side. The method is suitable for time-stepping adaptive mesh refinement PDE discretizations with different degrees of stiffness. A numerical example with an advection-diffusion problem illustrates the new method's properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Hamiltonian kinetic-Hall magnetohydrodynamics with fluid and kinetic ions in the current and pressure coupling schemes

We present two generalized hybrid kinetic-Hall magnetohydrodynamics (MHD) models describing the interaction of a two-fluid bulk plasma, which consists of thermal ions and electrons, with energetic, suprathermal ion populations described by Vlasov dynamics. The dynamics of the thermal components are governed by standard fluid equations in the Hall MHD limit with the electron momentum equation providing an Ohm's law with Hall and electron pressure terms involving a gyrotropic electron pressure tensor. The coupling of the bulk, low-energy plasma with the energetic particle dynamics is accomplished through the current density (current coupling scheme; CCS) and the ion pressure tensor appearing in the momentum equation (pressure coupling scheme; PCS) in the first and the second model, respectively. The CCS is a generalization of two well-known models, because in the limit of vanishing energetic and thermal ion densities, we recover the standard Hall MHD and the hybrid kinetic-ions/fluid-electron model, respectively. This provides us with the capability to study in a continuous manner, the global impact of the energetic particles in a regime extending from vanishing to dominant energetic particle densities. The noncanonical Hamiltonian structures of the CCS and PCS, which can be exploited to study equilibrium and stability properties through the energy-Casimir variational principle, are identified. As a first application here, we derive a generalized Hall MHD Grad–Shafranov–Bernoulli system for translationally symmetric equilibria with anisotropic electron pressure and kinetic effects owing to the presence of energetic particles using the PCS.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Multifunctional Thiol-Containing Additives for Improved Photoluminescence and Photovoltaic Performance of Cs 0.15 FA 0.85 PbI 3 Perovskites

Thiol containing molecules as both interfacial surface ligands and additives are promising modulators for enhancing photoluminescence (PL) properties, stability, and photovoltaic (PV) performance of metal halide perovskites. However, alkanethiols are much more effective for improving photoluminescence (PL) intensity and stability in metal halide perovskite nanocrystals than in their thin film analogues. Herein, we investigate how additional functional groups on a pyrimidine core can alter thiol reactivity and influence the PL, stability, and PV performance of organic metal halide perovskites. Through an investigation of five different pyrimidine derivatives, it is shown that all derivatives containing thiol groups form thiolates in the presence of the perovskite precursors and increase the photoluminescence intensity of the perovskite film. The largest all-around improvement to the PL intensity, stability, and photovoltaic performance of Cs 0.15 FA 0.85 PbI 3 perovskites is realized through the addition of a hydroxyl group combined with a trifluoromethyl group to form 4-hydroxy-6-(trifluoromethyl)pyrimidine-2-thiol. Furthermore, this investigation helps illuminate how combinations of functional groups can be employed to further increase the beneficial effects over monofunctional additives in organic metal halide perovskites.

14 SOLAR ENERGY↗
DOE OSTI2021

Single Iridium Atom Doped Ni 2 P Catalyst for Optimal Oxygen Evolution

Single-atom catalysts (SACs) with 100% active sites have excellent prospects for application in the oxygen evolution reaction (OER). However, further enhancement of the catalytic activity for OER is quite challenging, particularly for the development of stable SACs with overpotentials <180 mV. Here, we report an iridium single atom on Ni 2 P catalyst (Ir SA -Ni 2 P) with a record low overpotential of 149 mV at a current density of 10 mA·cm –2 in 1.0 M KOH. The Ir SA -Ni 2 P catalyst delivers a current density up to ∼28-fold higher than that of the widely used IrO 2 at 1.53 V vs RHE. Both the experimental results and computational simulations indicate that Ir single atoms preferentially occupy Ni sites on the top surface. The reconstructed Ir–O–P/Ni–O–P bonding environment plays a vital role for optimal adsorption and desorption of the OER intermediate species, which leads to marked enhancement of the OER activity. Additionally, the dynamic “top-down” evolution of the specific structure of the Ni@Ir particles is responsible for the robust single-atom structure and, thus, the stability property. In conclusion, this Ir SA -Ni 2 P catalyst offers novel prospects for simplifying decoration strategies and further enhancing OER performance.

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

A high-temperature mixed potential CO gas sensor for in situ combustion control

An in situ, accurate and robust sensor that can sustain a high temperature of above 1000 °C is needed for on-site combustion monitoring because it can give real-time and local data to the control system to adjust the overall combustion efficiency. In this work, we find that nickel oxide (NiO) is a promising sensing material for CO which is a direct indicator of the status of the combustion process of a power plant. Under the conditions of 0.5–3% O 2 and 1000 °C, the fabricated yttrium-stabilized zirconia (YSZ)-based mixed potential sensor using porous NiO demonstrates good sensitivity to CO, showing a signal as high as 36 mV to 1000 ppm CO. The effects of gas transport, structure and geometry of the NiO electrode on the sensing performance are studied. Results show that fast gas transport is much beneficial to improved sensitivity. NiO having a porous structure is much more sensitive to CO than that having a dense structure. But the t90 (time to achieve 90% final signal magnitude) of the former is much longer due to slower gas diffusion inside the pores. It's worth noting that the NiO sensor exhibits a positive relationship with the CO content, opposite to other reported results of mixed potential sensors to detect CO. We find that this might be due to the electrochemical reduction, instead of oxidation of CO during the interaction with NiO at 1000 °C. Selectivity tests on how CO 2 , CH 4 and steam affect CO sensing are also demonstrated. NiO is insensitive to even 10% CO 2 . CH 4 does not shift the average value of the CO sensing response. However, it makes the sensing signal fluctuate more intensively. 2% steam exerts a great influence on NiO sensitivity to CO: it magnifies the sensitivity of the porous NiO electrode to a low CO range of 0–100 ppm, but inhibits its sensitivity to a CO range of 100 ppm to 1000 ppm. Lastly, an 11 day stability test demonstrates the promising stability property of the proposed NiO-based YSZ sensor.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

On the convergence of an IEQ-based first-order semi-discrete scheme for the Beris-Edwards system

We present a convergence analysis of an unconditionally energy-stable first-order semi-discrete numerical scheme designed for a hydrodynamic Q-tensor model, the so-called Beris-Edwards system, based on the Invariant Energy Quadratization Method (IEQ). The model consists of the Navier–Stokes equations for the fluid flow, coupled to the Q-tensor gradient flow describing the liquid crystal molecule alignment. By using the Invariant Energy Quadratization Method, we obtain a linearly implicit scheme, accelerating the computational speed. However, this introduces an auxiliary variable to replace the bulk potential energy and it is a priori unclear whether the reformulated system is equivalent to the Beris-Edward system. In this work, we prove stability properties of the scheme and show its convergence to a weak solution of the coupled liquid crystal system. We also demonstrate the equivalence of the reformulated and original systems in the weak sense.

Weber, Franziska↗
DOE OSTI2020

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗