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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Impact of divertor plasma parameters on characteristics of current-convective turbulence under DIII-D-like detached conditions

Recent observations of turbulence in plasma of the DIII-D inner divertor revealed formation of poloidally localized fluctuations with frequencies in the low-kHz range, when the machine was operating in the fluctuating state of detachment. This phenomenon can be related to the onset of the current-convective instability (CCI). In this contribution, we use the basic physical model of the instability to simulate the saturated current-convective turbulence under the DIII-D-like detached conditions. Parameter scans of turbulence characteristics are performed. The spatial and temporal spectra of fluctuations are obtained and compared with the available experimental data. It is demonstrated that simulations are able to qualitatively reproduce most of the distinct features of turbulence found in the experiments, indicating that the CCI can be a plausible mechanism responsible for plasma activity recently observed in DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

An assessment of using variable blade pitch for moored ocean current turbine flight control

This study investigates the possibility of using active individual blade pitch control for positioning moored ocean current turbines within an array. Using a numerical simulation of a representative ocean current turbine it is shown that harmonic blade oscillations with amplitudes of 1.5° can be used to displace an ocean current turbine laterally by a distance equal to one rotor diameter, when the mooring cable length is approximately 30 times the diameter of the rotor blade. For current flow directions that would result in a downstream turbine operating in the wake of an upstream system at a distance of 10 diameters, this technique could be used to reduce the power loss of the downstream system from about 50% (caused by the reduced flow speed available in the turbine wake) to the 1.2% power loss associated with the suggested control approach.

16 TIDAL AND WAVE POWER↗

DOE OSTI2016

Mapping the current–current correlation function near a quantum critical point

The current–current correlation function is a useful concept in the theory of electron transport in homogeneous solids. The finite-temperature conductivity tensor as well as Anderson’s localization length can be computed entirely from this correlation function. Based on the critical behavior of these two physical quantities near the plateau–insulator or plateau–plateau transitions in the integer quantum Hall effect, we derive an asymptotic formula for the current–current correlation function, which enables us to make several theoretical predictions about its generic behavior. For the disordered Hofstadter model, we employ numerical simulations to map the current–current correlation function, obtain its asymptotic form near a critical point and confirm the theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

On the applicability of surrogate-based Markov chain Monte Carlo-Bayesian inversion to the Community Land Model: Case studies at flux tower sites: SURROGATE-BASED MCMC FOR CLM

The Community Land Model (CLM) has been widely used in climate and Earth system modeling. Accurate estimation of model parameters is needed for reliable model simulations and predictions under current and future conditions, respectively. In our previous work, a subset of hydrological parameters has been identified to have significant impact on surface energy fluxes at selected flux tower sites based on parameter screening and sensitivity analysis, which indicate that the parameters could potentially be estimated from surface flux observations at the towers. To date, such estimates do not exist. In this paper, we assess the feasibility of applying a Bayesian model calibration technique to estimate CLM parameters at selected flux tower sites under various site conditions. The parameters are estimated as a joint probability density function (PDF) that provides estimates of uncertainty of the parameters being inverted, conditional on climatologically-average latent heat fluxes derived from observations. We find that the simulated mean latent heat fluxes from CLM using the calibrated parameters are generally improved at all sites when compared to those obtained with CLM simulations using default parameter sets. Further, our calibration method also results in credibility bounds around the simulated mean fluxes which bracket the measured data. The modes (or maximum a posteriori values) and 95% credibility intervals of the site-specific posterior PDFs are tabulated as suggested parameter values for each site. Analysis of relationships between the posterior PDFs and site conditions suggests that the parameter values are likely correlated with the plant functional type, which needs to be confirmed in future studies by extending the approach to more sites.

Huang, Maoyi↗

DOE OSTI2019

Experimental Determination of PV Inverter Response to Grid Phase Shift Events: Preprint

With the continued growth of renewable energy resources which interface to the electric grid via inverters, the understanding of such devices becomes ever more important to the safe and reliable operation of the bulk power system. This work investigates the specific response of a utility-scale PV inverter to grid voltage phase shift-type disturbances which sometimes occur during grid fault events. The role of the PV inverter's phase-locked-loop (PLL) is identified as important to modeling the response. Switching-level simulations of a utility-scale PV inverter with a modeled PLL show a characteristic response when phase shift disturbances require the PLL to track what appear as fast frequency changes. Additionally, in this work a full-scale laboratory testing was carried out using the Opal real time (RT) OP5142 real time simulator and a large grid simulator to create phase shift disturbances with a high degree of repeatability. A photovoltaic (PV) inverter was connected to a grid simulator, and phase shifts were instantaneously implemented on the simulated grid, the results of the currents were then obtained. The experimental results obtained were validated with simulation results obtained from MATLAB/Simulink.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Runaway electron plateau current profile reconstruction from synchrotron imaging and Ar-II line polarization angle measurements in DIII-D

Abstract Current profile reconstructions are obtained for high current ( I p ≃ 550 kA) post-disruption runaway electron (RE) plateau plasmas in DIII-D. Two novel methods of measuring the RE current profile in high-current RE plateaus are introduced and compared: localization of the q = 2 rational surface using visible synchrotron emission (SE) imaging and the measurement of the polarization angle of line-integrated Ar-II line emission. The two methods are found to be consistent with each other within the data uncertainties. Different simulations of the RE current profile are compared with the measurements: the toroidal fluid RE model is found to best fit the data, within the measurement uncertainties. In addition to introducing two novel methods to measure the RE current profile and validating present simulation capabilities, this work demonstrates that instabilities can grow at q = 2 and q = 1 surfaces without necessarily causing a RE final loss instability. Numerical simulations are also presented to elucidate the role of these instabilities on synchrotron emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Integration of RAM-SCB into the Space Weather Modeling Framework

We present that numerical simulations of the ring current are a challenging endeavor. They require a large set of inputs, including electric and magnetic fields and plasma sheet fluxes. Because the ring current broadly affects the magnetosphere-ionosphere system, the input set is dependent on the ring current region itself. This makes obtaining a set of inputs that are self-consistent with the ring current difficult. To overcome this challenge, researchers have begun coupling ring current models to global models of the magnetosphere-ionosphere system. This paper describes the coupling between the Ring current Atmosphere interaction Model with Self-Consistent Magnetic field (RAM-SCB) to the models within the Space Weather Modeling Framework. Full details on both previously introduced and new coupling mechanisms are defined. Finally, the impact of self-consistently including the ring current on the magnetosphere-ionosphere system is illustrated via a set of example simulations.

58 GEOSCIENCES↗

DOE OSTI2018

Calculations of Beam-Beam Effect and Luminosity for Crab Dynamics Simulations in JLEIC

Crab crossing is an integral part of the Jefferson Lab Electron-Ion Collider (JLEIC) design to achieve high luminosity while meeting the detection and physics pro-gram requirements. The crab crossing scheme provides a head-on beam-beam collision for beams with a nonzero crossing angle. Simulations of crabbing dynamics currently do not include beam-beam effects. We describe a framework for accurate simulation of beam-beam effects on crabbing dynamics by applying a numerical calculation of the Bassetti-Erskine analytic solution to symplectic particle tracking codes. The numerical calculation is benchmarked against the analytic solution by calculating the luminosity reduction for several colliding beam scenarios. Benchmarking results show good agreement be-tween the numerical calculation and analytic solution, paving the way for implementation of the beam-beam kick to Elegant tracking simulations.

Huang, He↗

DOE OSTI2015

Integrated Waste Treatment Unit (IWTU) Input Coal Analyses and Off-Gass Filter (OGF) Content Analyses

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2015

Integrated Waste Treatment Unit (IWTU) Input Coal Analyses and Off-Gass Filter (OGF) Content Analyses

A full engineering scale Fluidized Bed Steam Reformer (FBSR) system is being used at the Idaho Nuclear Technology and Engineering Center (INTEC) to stabilize acidic Low Activity Waste (LAW) known as Sodium Bearing Waste (SBW). The INTEC facility, known as the Integrated Waste Treatment Unit (IWTU), underwent an Operational Readiness Review (ORR) and a Technology Readiness Assessment (TRA) in March 2014. The IWTU began non-radioactive simulant processing in late 2014 and by January, 2015 ; the IWTU had processed 62,000 gallons of simulant. The facility is currently in a planned outage for inspection of the equipment and will resume processing simulated waste feed before commencing to process 900,000 gallons of radioactive SBW. The SBW acidic waste will be made into a granular FBSR product (carbonate based) for disposal in the Waste Isolation Pilot Plant (WIPP). In the FBSR process calcined coal is used to create a CO2 fugacity to force the waste species to convert to carbonate species. The quality of the coal, which is a feed input, is important because the reactivity, moisture, and volatiles (C,H,N,O, and S) in the coal impact the reactions and control of the mineralizing process in the primary steam reforming vessel, the Denitration and Mineralizing Reformer (DMR). Too much moisture in the coal can require that additional coal be used. However since moisture in the coal is only a small fraction of the moisture from the fluidizing steam this can be self-correcting. If the coal reactivity or heating value is too low then the coal feed-rate needs to be adjusted to achieve the desired heat generation. Too little coal and autothermal heat generation in the DMR cannot be sustained and/or the carbon dioxide fugacity will be too low to create the desired carbonate mineral species. Too much coal and excess S and hydroxide species can form. Excess sulfur from coal that (1) is too rich in sulfur or (2) from overfeeding coal can promote wall scale and contribute to corrosion in process piping and materials, in excessive off-gas absorbent loading, and in undesired process emissions. The ash content of the coal is important as the ash adds to the DMR and other vessel products which affect the final waste product mass and composition. The amount and composition of the ash also affects the reaction kinetics. Thus ash content and composition contributes to the mass balance. In addition, sodium, potassium, calcium, sulfur, and maybe silica and alumina in the ash may contribute to wall-scale formation. Sodium, potassium, and alumina in the ash will be overwhelmed by the sodium, potassium, and alumina from the feed but the impact from the other ash components needs to be quantified. A maximum coal particle size is specified so the feed system does not plug and a minimum particle size is specified to prevent excess elutriation from the DMR to the Process Gas Filter (PGF). A vendor specification was used to procure the calcined coal for IWTU processing. While the vendor supplied a composite analysis for the 22 tons of coal (Appendix A), this study compares independent analyses of the coal performed at the Savannah River National Laboratory (SRNL) and at the National Energy Technology Laboratory (NETL). Three supersacks a were sampled at three different heights within the sack in order to determine within bag variability and between bag variability of the coal. These analyses were also compared to the vendor’s composite analyses and to the coal specification. These analyses were also compared to historic data on Bestac coal analyses that had been performed at Hazen Research Inc. (HRI) between 2004-2011.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2017

Adaptive Individual Blade Pitch Control of an Ocean Current Turbine for Load Reduction

A Direct Adaptive Disturbance Rejection (DADR) control algorithm is evaluated using a numer-ical simulation of an Ocean Current Turbine (OCT) for its ability to minimize the cyclic loadings on the rotor blades caused by current shear. This control algorithm reduces these cyclic loadings by individually altering the pitch angles of the rotor blades as a function of rotor rotation angle. The OCT numerical model represents a system with a single 3 bladed 20 meter diameter rotor operated in the Gulf Stream off the coast of Ft. Lauderdale, Florida for mean and maximum measured current shear values of 0.004 and 0.03 (m/s)/m. The results indicate that the adaptive controller reduced the amplitude of the cyclic axial blade loads on a rigidly mounted and moored turbine by 91.18% and 92.3% in the mean current shear and by 94.12% and 81.88% in maximum current shear demonstrating the value of this controls approach.

adaptive control↗

DOE OSTI2025

Accelerating Neutrino Event Generation in MARLEY Using CUDA-Based RNG and GPU Parallelization

MARLEY is a simulation tool that helps scientists study how low-energy neutrinos interact with matter. To work properly, MARLEY uses random numbers thousands of times in each simulation. These random numbers are important for modeling things like how neutrinos collide with atoms and what particles they produce. Right now, MARLEY runs on a regular computer processor (CPU) and uses a built-in random number generator called the Mersenne Twister. This setup works, but it can be slow, especially when trying to simulate many events. This research focuses on making MARLEY run faster by moving the random number generation and some of the repetitive calculations from the CPU to a graphics processing unit (GPU), which can handle many tasks at the same time. We use CUDA (a tool for programming NVIDIA GPUs) and cuRAND (a GPU-based random number library) to test faster alternatives to the current random number system. We compare different GPU-based generators, like curand_mtgp32, xorwow, and philox, to see which ones are the quickest and still give reliable results. Early tests show that using the GPU can make MARLEY simulations much faster. This project not only helps improve current simulation performance but also moves closer to a full simulation chain where all stages can run on modern GPU hardware.

Dunkley, Kimieka [Florida A-M]↗

DOE OSTI2016

Dynamic phase-control of a rising sun magnetron using modulated and continuous current

Phase-control of a magnetron is studied via simulation using a combination of a continuous current source and a modulated current source. The addressable, modulated current source is turned ON and OFF at the magnetron operating frequency in order to control the electron injection and the spoke phase. Prior simulation work using a 2D model of a Rising Sun magnetron showed that the use of 100% modulated current controlled the magnetron phase and allowed for dynamic phase control. In this work, the minimum fraction of modulated current source needed to achieve a phase control is studied. The current fractions (modulated versus continuous) were varied from 10% modulated current to 100% modulated current to study the effects on phase control. Dynamic phase-control, stability, and start up time of the device were studied for all these cases showing that with 10% modulated current and 90% continuous current, a phase shift of 180° can be achieved demonstrating dynamic phase control.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

3D modeling of the Staged Z-pinch with the FLASH code

The Staged Z-pinch (SZP) fusion concept is a magneto-inertial compression scheme developed by Magneto-Inertial Fusion Technologies, Inc. (MIFTI), in which small amounts of fusion fuel are brought to fusion-relevant conditions by passing multi-million amperes strong current through a cylindrical shell of high atomic number material. One- and two-dimensional modeling performed by MIFTI with the MACH2 code, suggests that net fusion energy gain can be achieved when currents in the 10 million amperes range compress a 50%-50% mixture of deuterium and tritium gas. In this project we enhance and use the FLASH code to execute high-fidelity, simulations of various SZP configurations, a collaboration between the Flash Center for Computational Science at the University of Rochester and the Laboratory for Laser Energetics, and MIFTI. FLASH is a high-performance computing, multi-physics, radiation magnetohydrodynamic (MHD) code with extended physics capabilities, which is developed at the Flash Center. The goal of this project is to assess the shell/fuel stability of the pinch to two- and three-dimensional MHD instabilities and to utilize FLASH’s extended physics capabilities to understand how extended-MHD effects impact implosion dynamics and plasma conditions at stagnation. The project is a natural extension of an ongoing collaborative effort between MIFTI and the Flash Center through the U.S. DOE Advanced Research Projects Agency-Energy (ARPA-E) BETHE program and can provide MIFTI with simulation capabilities that are currently beyond their reach with MACH2. Ultimately, MIFTI wants FLASH to become one of their simulation workhorses for reliable, high-fidelity SZP platform design. During the project, the Flash Center team developed a suite of in one- and two-dimensional FLASH simulations to model different variants of the SZP platform, performed code-to-code comparisons with MACH2, and used the FLASH code to design and model SZP experiments like the Double Eagle experiments that MIFTI recently performed. Also, the Flash Center team developed the capabilities of FLASH to be able to do three-dimensional simulations of pulsed-power experiments with the code for the first time. Two publications from this effort are currently under review and two more are in preparation. MIFTI has also committed to the use of FLASH for future SZP simulation efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Symmetry Parameter Constraints from a Lower Bound on Neutron-matter Energy

We propose the existence of a lower bound on the energy of pure neutron matter (PNM) on the basis of unitary-gas considerations. We discuss its justification from experimental studies of cold atoms as well as from theoretical studies of neutron matter. We demonstrate that this bound results in limits to the density-dependent symmetry energy, which is the difference between the energies of symmetric nuclear matter and PNM. In particular, this bound leads to a lower limit to the volume symmetry energy parameter S {sub 0}. In addition, for assumed values of S {sub 0} above this minimum, this bound implies both upper and lower limits to the symmetry energy slope parameter L , which describes the lowest-order density dependence of the symmetry energy. A lower bound on neutron-matter incompressibility is also obtained. These bounds are found to be consistent with both recent calculations of the energies of PNM and constraints from nuclear experiments. Our results are significant because several equations of state that are currently used in astrophysical simulations of supernovae and neutron star mergers, as well as in nuclear physics simulations of heavy-ion collisions, have symmetry energy parameters that violate these bounds. Furthermore, below the nuclear saturation density, the bound on neutron-matter energies leads to a lower limit to the density-dependent symmetry energy, which leads to upper limits to the nuclear surface symmetry parameter and the neutron-star crust–core boundary. We also obtain a lower limit to the neutron-skin thicknesses of neutron-rich nuclei. Above the nuclear saturation density, the bound on neutron-matter energies also leads to an upper limit to the symmetry energy, with implications for neutron-star cooling via the direct Urca process.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

SU-E-T-08: A Beam Source Model for Monte Carlo Simulations of a Double-Scattering Proton Beam Delivery System Using Beam Current Modulation

Purpose: To reconstruct phase-space information upstream of patient specific collimators for Monte Carlo simulations using only radiotherapy planning system data. Methods: The proton energies are calculated based on residual ranges, e.g., sum of prescribed ranges in a patient and SSD. The Kapchinskij and Vladimirskij (KV) distribution was applied to sample proton’s x-y positions and momentum direction and the beam shape was assumed to be a circle. Free parameters, e.g., the initial energy spread and the emittance of KV distribution were estimated from the benchmarking with commissioning data in a commercial treatment planning system for an operational proton therapy center. The number of histories, which defines the height of individual pristine Bragg peaks (BP) of Spread-out Bragg peak (SOBP), are weighted based on beam current modulation and a correction factor is applied to take into account the fluence reduction as the residual range decreases due to the rotation of the range modulator wheel. The timedependent behaviors, e.g., the changes of the residual range and histories per a pristine BP, are realized by utilizing TOPAS (Tool for Particle Simulation). Results: Benchmarking simulations for selected SOBPs ranging 7.5 cm to 15.5 cm matched within 2 mm in range and up to 5 mm in SOBP width against measurement data in water phantom. We found this model tends to underestimate entrance dose by about 5 % in comparison to measurement. This was attributed to the situation that the energy distribution used in the model was limited in its granularity at the limit of single energy spectrum for the narrow angle modulator steps used in the proximal pull back region of the SOBPs. Conclusion: Within these limitations the source modeling method proved itself an acceptable alternative of a full treatment head simulation when the machine geometry and materials information are not available.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Nonlinear dynamics and dissipation of a curvilinear vortex driven by a strong time-dependent Meissner current

Here, we report numerical simulations of large-amplitude oscillations of a trapped vortex line under a strong ac magnetic field $\textit{H(t) = H}$ sin$\textit{ωt}$ parallel to the surface. The power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was calculated by taking into account the nonlinear vortex line tension, vortex mass, and a nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient $\textit{η(v)}$. We show that the LO decrease of $\textit{η(v)}$ with the vortex velocity $\textit{v}$ can radically change the field dependence of the surface resistance R i ( H ) caused by trapped vortices. At low frequencies, $R_i(H)$ exhibits a conventional increase with $\textit{H}$, but as $\textit{ω}$ increases, the surface resistance becomes a nonmonotonic function of $\textit{H}$ which decreases with $\textit{H}$ at higher fields. The effects of frequency, pin spacing, and the mean-free path $l_i$ on the field dependence of $R_i(H)$ were calculated. It is shown that, as the surface gets dirtier and $l_i$ decreases, the anomalous drop of $R_i(H)$ with $\textit{H}$ shifts to lower fields which can be much smaller than the lower critical magnetic field. Our numerical simulations also show that the LO decrease of $\textit{η(v)}$ with $\textit{v}$ can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along the vortex. Measurements of $R_i(H)$ caused by sparse vortices trapped perpendicular to the surface can offer opportunities to investigate an extreme nonlinear dynamics of vortices driven by strong current densities up to the depairing limit at low temperatures. The behavior of $R_i(H)$, which can be tuned by varying the rf frequency or concentration of nonmagnetic impurities, is not masked by strong heating effects characteristic of dc or pulse transport measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2013

Classical trajectory Monte Carlo code for simulating ion beam focusing or defocusing with magnetic elements modeled as current loops or current lines

A computational tool is described that can be used for designing magnetic focusing or defocusing systems. A fully three-dimensional classical trajectory Monte Carlo simulation has been developed. Ion trajectories are simulated in the presence of magnetic elements that can be modeled as any combination of current loops and current lines. Each current loop or line may be located anywhere in the system and oriented along any of the three coordinate axes. The configuration need not be axisymmetric. The solutions are obtained using normalized parameters, which can be used for easily scaling the results. Examples are provided of the utility of the code.

43 PARTICLE ACCELERATORS↗