Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Stable compression of a spherical tokamak plasma

In this study, an analysis of toroidal plasma compression by a collapsing flux conserver reveals linearly stable scenarios of operation to high compression ratios. The resistive and ideal MHD stability is calculated in full toroidal geometry using the asymptotic matching method in realistic conditions. A time dependent MHD simulation of the compression is conducted to confirm the conservation principles used to calculate the equilibrium states for the analysis. The near edge current profile, controlled by toroidal field ramping during compression, is shown to be critical to stability due to coupling between poloidal components of the least stable mode. Two cases with slight differences in edge current are examined, one with a stable corridor to high compression and one without.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling of ExB effects on tungsten re-deposition and transport in the DIII-D divertor

Mixed-material DIVIMP-WallDYN modelling, now incorporating ExB drifts, is presented that simultaneously reproduces tungsten (W) erosion and deposition patterns observed during the DIII-D Metal Rings Campaign, in which a toroidally symmetric set of W-coated tiles were installed in the carbon (C) DIII-D divertor. Since most reactor plasma facing component (PFC) designs call for mixed-material environments, including ITER’s W/Be enviroment, the divertor targets will quickly evolve into reconstitued surfaces of multiple elements. This work identifies controlling physics that affects material migration patterns in the divertor, which impact PFC lifetimes and impurity leakage from the divertor to the core. These simulations indicate that radial and poloidal ExB transport dominates over parallel force balance for high-Z impurities such as W in the divertor region of DIII-D. It is demonstrated that ExB drifts are required to reproduce the experimental observation of non-local W and C co-accumulation in a band ~7-9 cm outboard of the outer-strike-point (OSP) W source, for attached Lmode conditions in the unfavorable ion grad-B drift direction. In addition, W gross erosion is localized to the region outboard of the OSP, as the formation of C co-deposits suppresses W erosion at the strike point. Time-dependent simulations with scaled ExB impurity drifts (60% of the OEDGE-calculated drift velocity) and W re-erosion quantitatively reproduce these features, including depth-resolved W/C ratios, within a factor of 2 over ~115 seconds of accumulated plasma exposure. The location of co-deposition regions is shown to be well represented by an analytic leakage model, driven largely by poloidal ExB drifts. Qualitative agreement is also found between campaign-integrated W deposition measurements and simulations for the favorable ion grad-B drift direction, the standard mode of operation for most tokamaks. Furthermore, these results imply that a longterm inward radial migration of material from the outer divertor through the private flux region may occur in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Extension of the energetic particle transport kick model in TRANSP to multiple fast ion species

Abstract Alfvénic instabilities (AEs) are well known to cause enhanced transport of energetic particles (EPs) in fusion devices. Most studies until now have focused on characterizing and understanding AE stability in single-species plasmas heated by neutral beams (NB), where deuterium is typically used as both main plasma species and NB fuel. As the fusion community moves toward fusion reactors that target burning plasma conditions, such as ITER, the single-species picture breaks down. Burning plasmas, which will use a mix of deuterium and tritium (DT) as main fuel, also feature the presence of several supra-thermal fusion products such as alpha particles, protons, helium isotopes and high-energy tritium ions. This work presents the extension of the EP transport kick model implemented in the TRANSP time-dependent tokamak transport code to study the combined effect of multiple EP species on AE stability and, in turn, the response of different EP species to plasma instabilities in terms of their redistribution and losses. Further validation of the enhanced model is planned based on experimental results expected from the JET DT campaign scheduled for 2021, in preparation for ITER plasmas and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Boundary condition effects on runaway electron mitigation coil modeling for the SPARC and DIII-D tokamaks

Abstract Extended-MHD modeling of planned Runaway Electron Mitigation Coils (REMC) for SPARC and DIII-D is performed with the NIMROD code. A coil has been designed for each machine, with the two differing in shape and location, but both having n = 1 symmetry (with n the toroidal mode number). Compared to previous modeling efforts, three improvements are made to the simulations boundary conditions. First a resistive wall model is used in place of an ideal wall. Second, the ThinCurr code is used to compute the time-dependent 3D fields used as magnetic boundary conditions for the simulations. Third, the simulation boundary is moved from the first-wall location to the Vacuum Vessel (VV), which extends the boundary past the location of the internal REMC. To remove the 3D coil from the simulation domain, an equivalent set of 3D fields is calculated at the VV boundary that produce approximately the same field distribution at the last closed flux surface assuming vacuum between the two. Each of these three boundary condition improvements leads to an improvement in the predicted performance of the REMC for both machines. The resistive wall alone primarily effects the resonance of the coil with the plasma after the TQ, affecting the q-profile evolution in the SPARC modeling, and allowing the applied spectrum to be modified in response to the plasma in the DIII-D modeling. The movement of the simulation boundary has the most significant effect on the RE confinement overall, including in the early stages, particularly for a DIII-D inner wall limited equilibrium, where the RE loss fraction increases from 90% to > 99%, with SPARC RE losses also occurring much earlier when the boundary is placed at the VV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Efficient ECCD non-inductive plasma current start-up, ramp-up, and sustainment for an ST fusion reactor

The elimination of the need for an Ohmic heating solenoid may be the most impactful design driver for the realization of economical compact fusion tokamak reactor systems. However, this would require fully non-inductive start-up and current ramp-up from zero plasma current and low electron temperature of sub-keV to the full plasma current of ~10–15 MA at 20–30 keV electron temperature. To address this challenge, an efficient solenoid-free start-up and ramp-up scenario utilizing a low-field-side-launched extraordinary mode at the fundamental electron cyclotron harmonic frequency (X–I) is proposed, which has more than two orders of magnitude higher electron cyclotron current drive (ECCD) efficiency than the conventional ECCD for the sub-keV start-up regime. A time dependent model was developed to simulate the start-up scenarios. For the Spherical Tokamak Advanced Reactor (STAR) (Menard et al 2023 Next-Step Low-Aspect-Ratio Tokamak Design Studies (IAEA)), it was found that to fully non-inductively ramp-up to 15 MA, it would take about 25 MW of EC power at 170 GHz. Because of the relatively large plasma volume of STAR, radiation losses must be considered. It is important to make sure that high Z impurities are kept sufficiently low during the early current start-up phase where the temperature is sub-keV range. Since the initial current ramp up takes place at a factor of ten lower density compared to the sustained regimes, it is important to transition into a higher bootstrap fraction discharge at lower density to minimize the ECCD power requirement during the densification. For the sustainment phase an array of eight gyrotron launchers with a total of about 60 MW of fundamental O-mode was found to be sufficient to provide the required axis-peaked external current drive. High efficiencies between 19–57 kA MW –1 were found with optimal aiming, and these were resilient to small changes in aiming angles and density and temperature profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evaluation of SPARC divertor conditions in H-mode operation using SOLPS-ITER

The predicted divertor conditions for the SPARC tokamak are calculated using SOLPS-ITER for a range of scrape-off-layer (SOL) heat flux widths λ q , input powers, and particle fueling locations. Under H-mode scenario conditions with an upstream separatrix density of 1 x 10 20 m -3 , the most conservative range of λ q extrapolations ( 0.15 mm) results in extremely high unmitigated particle and energy fluxes to the divertor, both under full field (12.2 T) and power (P SOL = 29 MW) conditions, and 2/3 field with P SOL = 10 MW. Increasing the cross-field SOL diffusivities by 2–10× reduces the magnitude of the mitigation challenge, however strategies such as impurity seeding or strike-point-sweeping will likely still be required. A combination of steady-state and time-dependent SOLPS-ITER simulations are used to map out phase space diagrams of upstream and divertor conditions. The simulations include parallel currents but neglect cross-field drifts. At low upstream density the inner and outer divertor conditions are highly asymmetric, with a large temperature difference and significant heat fluxes driven by parallel currents. The solution has sharp bifurcations with a region of hysteresis, depending on whether the initial state is at a low or high density. This behavior is observed even when the fueling location, cross-field diffusivity, and impurity level is changed, although the density window with asymmetry is reduced with increasing diffusivity. The addition of neon impurity seeding reduces the divertor heat fluxes, but also causes a drop in the upstream electron density with fixed particle throughput. This drop can be counteracted by increased main ion throughput, however too much neon results in a back transition into the asymmetric divertor regimes suggesting a need for control of both main ion and impurity seeding levels to achieve a desired divertor state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Combining physics-based and data-driven models for quantitatively accurate plasma profile prediction that extrapolates well; with application to DIII-D, AUG, and ITER tokamaks

For design, scenario planning, and control, ITER and all other envisioned tokamaks rely on a variety of statistical and physics-based models to extrapolate to unseen regimes; most notably from low plasma current to high. A 'meta-learning' methodology for combining the accuracy of data-driven models with the generalizability of physics-based models is described and tested, yielding a 5–10 percent improvement in performance beyond either alone for the task of extrapolating time-dependent plasma profile prediction from low- to high- plasma current DIII-D tokamak discharges. Meanwhile, it is shown that both machine learning models extrapolated far-distribution and state-of-the-art 'physics-based' profile predictors fare worse than merely assuming plasma profiles do not change from their initial values. Finally, a variety of other mechanisms for helping data-driven models generalize—transfer learning, adding contextual information from physics simulators, and adding data from the ASDEX Upgrade tokamak—are attempted for similar extrapolation tasks but, in the methodology used in this paper, yield no significant improvement beyond simple data-driven models. Results are summarized in figures 15 and 16.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Polarization Tensor for Neutron Optics

A rank 2 dimension 3 tensor formulation for polarized neutron transport in a magnetic field is derived. Using this formulation, the rotating magnetic field problem is revisited to produce a single 3 × 3 vector-transformation matrix that describes the polarization transport. That result is then used to fashion a new strategy for a numerical approximation of polarization transport in an arbitrary magnetic field using time-dependent perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Shell effects in fission and quasi-fission reactions

Quantum shell effects are responsible for asymmetric fission. They are also expected to affect the formation of fission fragments in quasi-fission reactions occurring in heavy-ion collisions. Systematic time-dependent Hartree-Fock simulations of 40–56 Ca+ 176 Yb collisions show that quasi-fission fragment properties share strong similarities with fragments formed in fission of the compound nuclei. This is an indication that similar shell effects are responsible for the final asymmtery in both mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A hybrid 3D/2D field response calculation for liquid argon detectors with PCB based anode plane

Liquid Argon Time Projection Chamber (LArTPC) technology is commonly utilized in neutrino detector designs. It enables detailed reconstruction of neutrino events with high spatial precision and low energy threshold. Its field response (FR) model describes the time-dependent electric currents induced in the anode-plane electrodes when ionization electrons drift nearby. An accurate and precise FR is a crucial input to LArTPC detector simulations and charge reconstruction. Established LArTPC designs have been based on parallel wire planes. It allows accurate and computationally economic two-dimensional (2D) FR models utilizing the translational symmetry along the direction of the wires. Recently, novel LArTPC designs utilize electrodes formed on printed circuit board (PCB) in the shape of strips with through holes. The translational symmetry is no longer a good approximation near the electrodes and a new FR calculation that employs regions with three dimensions (3D) has been developed. Extending the 2D models to 3D would be computationally expensive. Fortuitously, the nature of strips with through holes allows for a computationally economic approach based on the finite-difference method (FDM). In this paper, we present a new software package pochoir that calculates LArTPC field response for these new strip-based anode designs. This package combines 3D calculations in the volume near the electrodes with 2D far-field solutions to achieve fast and precise field response computation. We apply the resulting FR to simulate and reconstruct samples of cosmic-ray muons and 39 Ar decays from a Vertical Drift (VD) detector prototype operated at CERN. We find the difference between real and simulated data within 5%. Current state-of-the-art LArTPC software requires a 2D FR which we provide by averaging over one dimension and estimate that variations lost in this average are smaller than 7%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Proton reconstruction with the TOTEM Roman pot detectors for high- β * LHC data

The TOTEM Roman pot detectors are used to reconstruct the transverse momentum of scattered protons and to estimate the transverse location of the primary interaction. This paper presents new methods of track reconstruction, measurements of strip-level detection efficiencies, cross-checks of the LHC beam optics, and detector alignment techniques, along with their application in the selection of signal collision events. The track reconstruction is performed by exploiting hit cluster information through a novel method using a common polygonal area in the intercept-slope plane. The technique is applied in the relative alignment of detector layers with μm precision. A tag-and-probe method is used to extract strip-level detection efficiencies. The alignment of the Roman pot system is performed through time-dependent adjustments, resulting in a position accuracy of 3 μm in the horizontal and 60 μm in the vertical directions. The goal is to provide an optimal reconstruction tool for central exclusive physics analyses based on the high-β* data-taking period at $\sqrt({s})$ = 13 TeV in 2018.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Parametrically driven nonlinear Dirac equation with arbitrary nonlinearity

Here, the damped and parametrically driven nonlinear Dirac equation with arbitrary nonlinearity parameter $\kappa$ is analyzed, when the external force is periodic in space and given by $\mathcal{f(x) = r}$ cos(Kx), both numerically and in a variational approximation using five collective coordinates (time dependent shape parameters of the wave function). Our variational approximation satisfies exactly the low-order moment equations. Because of competition between the spatial period of the external force λ = 2π/K , and the soliton width $l_s$, which is a function of the nonlinearity $\kappa$ as well as the initial frequency ω 0 of the solitary wave, there is a transition (at fixed ω 0 ) from trapped to unbound behavior of the soliton, which depends on the parameters $\mathcal{r}$ and K of the external force and the nonlinearity parameter $\kappa$. We previously studied this phenomena when $\kappa$ = 1 where we showed that for λ $\gg l_s$ the soliton oscillates in an effective potential, while forλ $\ll l_s$ it moves uniformly as a free particle. In this paper we focus on the $\kappa$ dependence of the transition from oscillatory to particle behavior and explicitly compare the curves of the transition regime found in the collective coordinate approximation as a function of $\mathcal{r}$ and K when $\kappa$ = 1/2,1,2 at fixed value of the frequency ω 0 . Since the solitary wave gets narrower for fixed ω 0 as a function of $\kappa$, we expect and indeed find that the regime where the solitary wave is trapped is extended as we increase $\kappa$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Integrable multistate Landau–Zener models with parallel energy levels

In this paper, we discuss solvable multistate Landau–Zener (MLZ) models whose Hamiltonians have commuting partner operators with ~1/τ-time-dependent parameters. Many already known solvable MLZ models belong precisely to this class. We derive the integrability conditions on the parameters of such commuting operators, and demonstrate how to use such conditions in order to derive new solvable cases. We show that MLZ models from this class must contain bands of parallel diabatic energy levels. The structure of the scattering matrix and other properties are found to be the same as in the previously discussed completely solvable MLZ Hamiltonians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Stability and response of trapped solitary wave solutions of coupled nonlinear Schrödinger equations in an external, $\mathcal{PT}$- and supersymmetric potential

In this work, we present trapped solitary wave solutions of a coupled nonlinear Schrödinger (NLS) system in 1 + 1 dimensions in the presence of an external, supersymmetric and complex $\mathcal{PT}$-symmetric potential. The Schrödinger system this work focuses on possesses exact solutions whose existence, stability, and spatio-temporal dynamics are investigated by means of analytical and numerical methods. Two different variational approximations are considered where the stability and dynamics of the solitary waves are explored in terms of eight and twelve time-dependent collective coordinates (CCs). We find regions of stability for specific potential choices as well as analytic expressions for the small oscillation frequencies in the CC approximation. Our findings are further supported by performing systematic numerical simulations of the NLS system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Revealing room temperature ferromagnetism in exfoliated Fe 5 GeTe 2 flakes with quantum magnetic imaging

Van der Waals (vdW) material Fe 5 GeTe 2 , with its long-range ferromagnetic ordering near room temperature, has significant potential to become an enabling platform for implementing novel spintronic and quantum devices. To pave the way for applications, it is crucial to determine the magnetic properties when the thickness of Fe 5 GeTe 2 reaches the few-layers regime. However, this is highly challenging due to the need for a characterization technique that is local, highly sensitive, artifact-free, and operational with minimal fabrication. Prior studies have indicated that Curie temperature T C can reach up to close to room temperature for exfoliated Fe 5 GeTe 2 flakes, as measured via electrical transport; there is a need to validate these results with a measurement that reveals magnetism more directly. Here, we investigate the magnetic properties of exfoliated thin flakes of vdW magnet Fe 5 GeTe 2 via quantum magnetic imaging technique based on nitrogen vacancy centers in diamond. Through imaging the stray fields, we confirm room-temperature magnetic order in Fe 5 GeTe 2 thin flakes with thickness down to 7 units cell. The stray field patterns and their response to magnetizing fields with different polarities is consistent with previously reported perpendicular easy-axis anisotropy. Furthermore, we perform imaging at different temperatures and determine the Curie temperature of the flakes at ≈300 K. These results provide the basis for realizing a room-temperature monolayer ferromagnet with Fe 5 GeTe 2 . This work also demonstrates that the imaging technique enables rapid screening of multiple flakes simultaneously as well as time-resolved imaging for monitoring time-dependent magnetic behaviors, thereby paving the way towards high throughput characterization of potential two-dimensional (2D) magnets near room temperature and providing critical insights into the evolution of domain behaviors in 2D magnets due to degradation.

2D magnet↗
DOE OSTI2022

Optical atomic clock aboard an Earth-orbiting space station (OACESS): enhancing searches for physics beyond the standard model in space

We present a concept for a high-precision optical atomic clock (OAC) operating on an Earth-orbiting space station. This pathfinder science mission will compare the space-based OAC with one or more ultra-stable terrestrial OACs to search for space-time-dependent signatures of dark scalar fields that manifest as anomalies in the relative frequencies of station-based and ground-based clocks. This opens the possibility of probing models of new physics that are inaccessible to purely ground-based OAC experiments where a dark scalar field may potentially be strongly screened near Earth's surface. This unique enhancement of sensitivity to potential dark matter candidates harnesses the potential of space-based OACs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Numerical methodology for treating static and dynamic dislocation problems near a free surface

Abstract Simulation of dislocation dynamics enables researchers and scientists to explore the plastic behavior of crystalline materials under loading. Analytic solutions for the stress field due to a linear dislocation segment near a free surface are case-specific, e.g. dealing with either a horizontal segment or a vertical segment, and therefore hard to implement in time-dependent dislocation dynamics simulations as different dislocation segments could be differently oriented. This article presents a generalized numerical framework to find the stress field beneath a free surface due to the presence of a dislocation segment. The framework can be expanded to non-flat surfaces. Also, three-dimensional discrete dislocation dynamics simulations are performed here, which clearly show the effect of free surfaces on the flow stress of a material.

Siddique, Abu Bakar↗