Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated current”

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 955 records · Page 53

Coupling characteristics of the ITER relevant lower hybrid antenna in Tore Supra: experiments and modelling

A new concept of lower hybrid antenna for current drive has been proposed for ITER [Bibet et al, Nuclear Fusion 1995]: the Passive Active Multijunction (PAM) antenna that relies on a periodic combination of active and passive waveguides. An actively cooled PAM antenna at 3.7 GHz has been recently installed on the tokamak Tore Supra. The paper summarizes the comprehensive experimental characterization of the linear coupling properties of the PAM antenna to the Tore Supra plasmas. These experimental results are systematically compared with the linear wave coupling theory via the linear ALOHA code. Good agreement between experimental results and ALOHA have been obtained. The detailed validation of the coupling modelling is an important step toward the validation of the PAM concept in view of further optimizing the electromagnetic properties of the future ITER antenna.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlinear evolution of three-dimensional instabilities of thin and thick electron scale current sheets: Plasmoid formation and current filamentation

Nonlinear evolution of three dimensional electron shear flow instabilities of an electron current sheet (ECS) is studied using electron-magnetohydrodynamic simulations. The dependence of the evolution on current sheet thickness is examined. For thin current sheets (half thickness =d{sub e}=c/ω{sub pe}), tearing mode instability dominates. In its nonlinear evolution, it leads to the formation of oblique current channels. Magnetic field lines form 3-D magnetic spirals. Even in the absence of initial guide field, the out-of-reconnection-plane magnetic field generated by the tearing instability itself may play the role of guide field in the growth of secondary finite-guide-field instabilities. For thicker current sheets (half thickness ∼5 d{sub e}), both tearing and non-tearing modes grow. Due to the non-tearing mode, current sheet becomes corrugated in the beginning of the evolution. In this case, tearing mode lets the magnetic field reconnect in the corrugated ECS. Later thick ECS develops filamentary structures and turbulence in which reconnection occurs. This evolution of thick ECS provides an example of reconnection in self-generated turbulence. The power spectra for both the thin and thick current sheets are anisotropic with respect to the electron flow direction. The cascade towards shorter scales occurs preferentially in the direction perpendicular to the electron flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamics of electrical double layer formation in room-temperature ionic liquids under constant-current charging conditions

We report detailed simulation results on the formation dynamics of an electrical double layer (EDL) inside an electrochemical cell featuring room-temperature ionic liquids (RTILs) enclosed between two planar electrodes. Under relatively small charging currents, the evolution of cell potential during charging can be suitably predicted by the Landau-Ginzburg-type continuum model proposed recently (M. Z. Bazant, B. D. Storey, and A. A. Kornyshev, Phys. Rev. Lett., 106, 046102, 2011). Under very large charging currents, the cell potential shows pronounced oscillation during the initial stage of charging, a feature not captured by the continuum model. Such oscillation originates from the sequential growth of the ionic space charge layers near the electrode surface, allowing the evolution of EDLs in RTILs with time, an atomistic process difficult to visualize experimentally, to be studied by analyzing the cell potential under constant current charging conditions. While the continuum model cannot predict the potential oscillation under such far-from-equilibrium charging conditions, it can nevertheless qualitatively capture the growth of cell potential during the later stage of charging. Furthermore, improving the continuum model by introducing frequency-dependent dielectric constant and density-dependent ion diffusion coefficients may help to further extend the applicability of the model. Keywords: ionic

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reactor Pressure Vessel Fracture Mechanics Development and Concrete Application Testing for Grizzly

The Grizzly code is being developed to address degradation issues in nuclear reactor structures and components. For light-water reactors, Grizzly currently has capabilities to simulate degradation processes and their effects on structural integrity in two key areas: reactor pressure vessels (RPVs) and reinforced concrete structures. This report documents improvements made to Grizzly’s ability to address both of these structural systems. For RPVs, the reduced-order models (ROMs) used in fracture mechanics calculations have been expanded to allow their application over a broader range of the parameter space than was permitted by the previous models. The ROMs currently used in Grizzly for the evaluation of flaws that are fully embedded within the RPV (as opposed to surface-breaking flaws) are based on a model that is known to be conservative, indicating higher stress intensity factors than would be obtained from direct simulations. A more accurate model that eliminates these excess conservatisms has been recently included in the American Society of Mechanical Engineers Boiler and Pressure Vessel Code but was not applicable for flaws near the RPV surface, which is where the most critical flaws are usually located. That model has recently been extended for increased applicability in this near-surface region. The ROMs for embedded flaws in the Grizzly code have been expanded to include these recent extensions, which permit their use in a much broader set of cases than previously possible. Direct 3D simulations have been used to check these ROMs and have shown good agreement in most cases, although there are still some cases that need further investigation. There are considerable benefits to using these these more accurate and less conservative ROMs for embedded flaws. On a benchmark probabilistic fracture mechanics problem tested here, the conditional probability of fracture initiation computed for a population of flaws in a single plate in an RPV decreased by over a factor of 3. To address aging in reinforced concrete structures, a capability to simulate multiple degradation mechanisms, including alkali-silica reaction and radiation-induced volumetric expansion has been developed in Grizzly over the past several years. This had previously been demonstrated on laboratory-scale specimens but not on full-scale nuclear concrete structures with reinforcement. To demonstrate the applicability of Grizzly to the analysis of large-scale structures of interest, a full 3D model of a representative reinforced concrete structure, including a complex arrangement of reinforcing bars, was developed and demonstrated in Grizzly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Digital Full-Scope Simulation of a Conventional Nuclear Power Plant Control Room, Phase 2: Installation of a Reconfigurable Simulator to Support Nuclear Plant Sustainability

The U.S. Department of Energy’s Light Water Reactor Sustainability program has developed a control room simulator in support of control room modernization at nuclear power plants in the U.S. This report highlights the recent completion of this reconfigurable, full-scale, full-scope control room simulator buildout at the Idaho National Laboratory. The simulator is fully reconfigurable, meaning it supports multiple plant models developed by different simulator vendors. The simulator is full-scale, using glasstop virtual panels to display the analog control boards found at current plants. The present installation features 15 glasstop panels, uniquely achieving a complete control room representation. The simulator is also full-scope, meaning it uses the same plant models used for training simulators at actual plants. Unlike in the plant training simulators, the deployment on glasstop panels allows a high degree of customization of the panels, allowing the simulator to be used for research on the design of new digital control systems for control room modernization. This report includes separate sections discussing the glasstop panels, their layout to mimic control rooms at actual plants, technical details on creating a multi-plant and multi-vendor reconfigurable simulator, and current efforts to support control room modernization at U.S. utilities. The glasstop simulator provides an ideal testbed for prototyping and validating new control room concepts. Equally importantly, it is helping create a standardized and vetted human factors engineering process that can be used across the nuclear industry to ensure control room upgrades maintain and even improve current reliability and safety.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Three-dimensional spontaneous magnetic reconnection

Magnetic reconnection is best known from observations of the Sun where it causes solar flares. Observations estimate the reconnection rate as a small, but non-negligible fraction of the Alfvén speed, so-called fast reconnection. Until recently, the prevailing pictures of reconnection were either of resistivity or plasma microscopic effects, which was contradictory to the observed rates. Alternative pictures were either of reconnection due to the stochasticity of magnetic field lines in turbulence or the tearing instability of the thin current sheet. In this paper we simulate long-term three-dimensional nonlinear evolution of a thin, planar current sheet subject to a fast oblique tearing instability using direct numerical simulations of resistive-viscous magnetohydrodynamics. The late-time evolution resembles generic turbulence with a −5/3 power spectrum and scale-dependent anisotropy, so we conclude that the tearing-driven reconnection becomes turbulent reconnection. The turbulence is local in scale, so microscopic diffusivity should not affect large-scale quantities. This is confirmed by convergence of the reconnection rate toward ∼0.015v{sub A} with increasing Lundquist number. In this spontaneous reconnection, with mean field and without driving, the dissipation rate per unit area also converges to ∼0.006ρv{sub A}{sup 3}, and the dimensionless constants 0.015 and 0.006 are governed only by self-driven nonlinear dynamics of the sheared magnetic field. Remarkably, this also means that a thin current sheet has a universal fluid resistance depending only on its length to width ratio and to v{sub A}/c.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum Algorithm for High Energy Physics Simulations

Simulating quantum field theories is a flagship application of quantum computing. However, calculating experimentally relevant high energy scattering amplitudes entirely on a quantum computer is prohibitively difficult. It is well known that such high energy scattering processes can be factored into pieces that can be computed using well established perturbative techniques, and pieces which currently have to be simulated using classical Markov chain algorithms. These classical Markov chain simulation approaches work well to capture many of the salient features, but cannot capture all quantum effects. To exploit quantum resources in the most efficient way, we introduce a new paradigm for quantum algorithms in field theories. This approach uses quantum computers only for those parts of the problem which are not computable using existing techniques. In particular, we develop a polynomial time quantum final state shower that accurately models the effects of intermediate spin states similar to those present in high energy electroweak showers with a global evolution variable. The algorithm is explicitly demonstrated for a simplified quantum field theory on a quantum computer.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Latitudinal dependence of nonlinear interaction between electromagnetic ion cyclotron wave and terrestrial ring current ions

Electromagnetic ion cyclotron (EMIC) waves can lead to the rapid decay (on a timescale of hours) of the terrestrial ring current. Such decay process is usually investigated in the framework of quasi-linear theory. Here, both theoretical analysis and test-particle simulation are performed to understand the nonlinear interaction between ring current ions and EMIC waves. In particular, the dependence of the nonlinear wave-particle interaction processes on the ion initial latitude is investigated in detail. These nonlinear processes are classified into the phase trapping and phase bunching, and the phase bunching is further divided into the channel and cluster effects. Compared to the prediction of the quasi-linear theory, the ring current decay rate can be reduced by the phase trapping, increased by the channel effect phase bunching, but non-deterministically influenced by the cluster effect phase bunching. The ion initial latitude changes the occurrence of the phase trapping, modulates the transport direction and strength of the cluster effect phase bunching, and only slightly affects the channel effect phase bunching. The current results suggest that the latitudinal dependence of these nonlinear processes should be considered in the evaluation of the ring current decay induced by EMIC waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pressure-induced emergence of unusually high-frequency transverse excitations in a liquid alkali metal: Evidence of two types of collective excitations contributing to the transverse dynamics at high pressures

Unlike phonons in crystals, the collective excitations in liquids cannot be treated as propagation of harmonic displacements of atoms around stable local energy minima. The viscoelasticity of liquids, reflected in transition from the adiabatic to elastic high-frequency speed of sound and in absence of the long-wavelength transverse excitations, results in dispersions of longitudinal (L) and transverse (T) collective excitations essentially different from the typical phonon ones. Practically, nothing is known about the effect of high pressure on the dispersion of collective excitations in liquids, which causes strong changes in liquid structure. Here dispersions of L and T collective excitations in liquid Li in the range of pressures up to 186 GPa were studied by ab initio simulations. Two methodologies for dispersion calculations were used: direct estimation from the peak positions of the L/T current spectral functions and simulation-based calculations of wavenumber-dependent collective eigenmodes. It is found that at ambient pressure, the longitudinal and transverse dynamics are well separated, while at high pressures, the transverse current spectral functions, density of vibrational states, and dispersions of collective excitations yield evidence of two types of propagating modes that contribute strongly to transverse dynamics. Emergence of the unusually high-frequency transverse modes gives evidence of the breakdown of a regular viscoelastic theory of transverse dynamics, which is based on coupling of a single transverse propagating mode with shear relaxation. The explanation of the observed high-frequency shift above the viscoelastic value is given by the presence of another branch of collective excitations. With the pressure increasing, coupling between the two types of collective excitations is rationalized within a proposed extended viscoelastic model of transverse dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Study of the effects of guide field on Hall reconnection

The results from guide field studies on the Magnetic Reconnection Experiment (MRX) are compared with results from Hall magnetohydrodynamic (HMHD) reconnection simulation with guide field. The quadrupole field, a signature of two-fluid reconnection at zero guide field, is modified by the presence of a finite guide field in a manner consistent with HMHD simulation. The modified Hall current profile contains reduced electron flows in the reconnection plane, which quantitatively explains the observed reduction of the reconnection rate. The present results are consistent with the hypothesis that the local reconnection dynamics is dominated by Hall effects in the collisionless regime of the MRX plasmas. While very good agreement is seen between experiment and simulations, we note that an important global feature of the experiments, a compression of the guide field by the reconnecting plasma, is not represented in the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improving Fission Products at CARIBU: Near Field Detection (Q4/FY21 Quarterly Progress Report)

Since the last quarterly, we have been working on improving the GEANT4 simulation model, which we currently use for the 4-pi gas counter efficiency calculations, to streamline the work when calculating multi-level decay properties. We have also been preparing for the next branching-ratio measurement, which will be for 161 Tb. Because travel is still limited and restricted, we have worked with a summer student on performing simulations to estimate our detector’s efficiency for different thresholds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of bootstrap current and Landau-fluid closure on ELM crashes and transport

Results presented here are from 6-field Landau-Fluid simulations using shifted circular cross-section tokamak equilibria on BOUT++ framework. Linear benchmark results imply that the collisional and collisionless Landau resonance closures make a little difference on linear growth rate spectra which are quite close to the results with the flux limited Spitzer-Härm parallel flux. Both linear and nonlinear simulations show that the plasma current profile plays dual roles on the peeling-ballooning modes that it can drive the low-n peeling modes and stabilize the high-n ballooning modes. For fixed total pressure and current, as the pedestal current decreases due to the bootstrap current which becomes smaller when the density (collisionality) increases, the operational point is shifted downwards vertically in the J ped – α diagram, resulting in threshold changes of different modes. In conclusion, the bootstrap current can slightly increase radial turbulence spreading range and enhance the energy and particle transports by increasing the perturbed amplitude and broadening cross-phase frequency distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stray RF Power Estimates From EC Exploitation During ITER Plasma Operations

The EC H and CD system of ITER tokamak is an essential tool for all the phases of ITER operation. Different levels of EC power are required through all the plasma discharge: up to 6.7 MW for assisting the breakdown and burn through, up to 20 MW for current drive and saw-teeth control from the equatorial launcher and up to 20 MW for NTM stabilization from the upper launchers. The assistance to breakdown and burn through is characterized by a very low (if not negligible) RF power absorption by the plasma. A significant level of stray radiation may also arise from partial absorption due to non-optimal plasma parameters and/or wrong injected polarization. The stray power radiated in the vacuum chamber is estimated as a first step toward mitigating potential harmful consequence to in-vessel structures and diagnostics. Power loading of the chamber walls (peek power and average power density for straight beam propagation in the empty chamber) and diffuse stray radiation effects are simulated to infer suitable strategies to avoid damage to first wall and to microwave sensitive components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electrical de-poling and re-poling of relaxor-PbTiO 3 piezoelectric single crystals without heat treatment

Re-poling of unexpected partially depoled piezoelectric materials conventionally needs to be first fully depoled through annealing above their Curie temperature to revive piezoelectric performances. Here, we investigated de-poling and re-poling of Pb(In 1/2 Nb 1/2 )O 3 -Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 single crystals under electric fields at room temperature. We found that alternating current electric fields with amplitudes near the coercive field at low frequencies (<10 Hz) can be employed to successfully depolarize poled crystals at room temperature. We also demonstrated a reversible polarization switching process with a relaxor-PbTiO 3 single crystal ultrasound transducer without device performance degradations. This experimental observation is supported by phase-field simulation, showing that alternating current electric fields can readily induce de-poling at room temperature, while direct current electric fields induce a transient depoled state only within an uncontrollable short period of time. The findings suggest new strategies for unprecedented in-device tailoring of the polarization states of ferroelectric materials.

36 MATERIALS SCIENCE↗

Incorporating bubble growth volume feedback to improve simulation of the response of a structure containing liquid and gas to sudden energy input

The SNS target module is a stainless-steel structure that contains and directs mercury. The mercury is struck with short, intense proton pulses to create neutrons. The pulses deposit energy and cause loads on the target structure and mercury cavitation. Helium bubbles are introduced into the mercury to reduce loads and cavitation. This gas complicates the prediction of the target module’s physical response. The current method used to simulate the structural response to the proton pulse incorporates a simple but effective approach to approximate mercury cavitation effects. This method cannot account for the changes from injected non-condensable gas. The ability to model the target vessel response is a prerequisite for estimating its life. Known approaches for simulating a bubbly fluid mixture’s response are computationally expensive and impractical for detailed engineering models. A constitutive model is proposed for target mercury with injected gas bubbles, with the intent of providing more accurate mercury–vessel response behavior without the computing costs required to track individual bubbles. The model’s assumptions and a computational model embodiment are described. The bubbly mercury computational model is subjected to test problems to assess its potential utility. Here, the model was able to better predict the effect of gas on wave propagation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Study of electron transport in a Hall thruster by axial–radial fully kinetic particle simulation

Electron transport across a magnetic field in a magnetic-layer-type Hall thruster was numerically investigated for the future predictive modeling of Hall thrusters. The discharge of a 1-kW-class magnetic-layer-type Hall thruster designed for high-specific-impulse operation was modeled using an r-z two-dimensional fully kinetic particle code with and without artificial electron-diffusion models. The thruster performance results showed that both electron transport models captured the experimental result within discrepancies less than 20% in thrust and discharge current for all the simulated operation conditions. The electron cross-field transport mechanism of the so-called anomalous diffusion was self-consistently observed in the simulation without artificial diffusion models; the effective electron mobility was two orders of magnitude higher than the value obtained using the classical diffusion theory. To account for the self-consistently observed anomalous transport, the oscillation of plasma properties was speculated. It was suggested that the enhanced random-walk diffusion due to the velocity oscillation of low-frequency electron flow could explain the observed anomalous diffusion within an order of magnitude. The dominant oscillation mode of the electron flow velocity was found to be 20 kHz, which was coupled to electrostatic oscillation excited by global ionization instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the possibility of the multiple inductively coupled plasma and helicon plasma sources for large-area processes

In this study, we attempted to determine the possibility of multiple inductively coupled plasma (ICP) and helicon plasma sources for large-area processes. Experiments were performed with the one and two coils to measure plasma and electrical parameters, and a circuit simulation was performed to measure the current at each coil in the 2-coil experiment. Based on the result, we could determine the possibility of multiple ICP sources due to a direct change of impedance due to current and saturation of impedance due to the skin-depth effect. However, a helicon plasma source is difficult to adapt to the multiple sources due to the consistent change of real impedance due to mode transition and the low uniformity of the B-field confinement. As a result, it is expected that ICP can be adapted to multiple sources for large-area processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accounting for herbaceous communities in process‐based models will advance our understanding of “grassy” ecosystems

Abstract Grassland and other herbaceous communities cover significant portions of Earth's terrestrial surface and provide many critical services, such as carbon sequestration, wildlife habitat, and food production. Forecasts of global change impacts on these services will require predictive tools, such as process‐based dynamic vegetation models. Yet, model representation of herbaceous communities and ecosystems lags substantially behind that of tree communities and forests. The limited representation of herbaceous communities within models arises from two important knowledge gaps: first, our empirical understanding of the principles governing herbaceous vegetation dynamics is either incomplete or does not provide mechanistic information necessary to drive herbaceous community processes with models; second, current model structure and parameterization of grass and other herbaceous plant functional types limits the ability of models to predict outcomes of competition and growth for herbaceous vegetation. In this review, we provide direction for addressing these gaps by: (1) presenting a brief history of how vegetation dynamics have been developed and incorporated into earth system models, (2) reporting on a model simulation activity to evaluate current model capability to represent herbaceous vegetation dynamics and ecosystem function, and (3) detailing several ecological properties and phenomena that should be a focus for both empiricists and modelers to improve representation of herbaceous vegetation in models. Together, empiricists and modelers can improve representation of herbaceous ecosystem processes within models. In so doing, we will greatly enhance our ability to forecast future states of the earth system, which is of high importance given the rapid rate of environmental change on our planet.

59 BASIC BIOLOGICAL SCIENCES↗