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 595 records · Page 33

Generation of valley-polarized electron beam in bilayer graphene

We propose a method to produce valley-polarized electron beams using a bilayer graphene npn junction. By analyzing the transmission properties of electrons through the junction with zigzag interface in the presence of trigonal warping, we observe that there exist a range of incident energies and barrier heights in which transmitted electrons are well polarized and collimated. From this observation and by performing numerical simulations, it is demonstrated that valley-dependent electronic currents with nearly perfect polarization can be generated. We also show that the peak-to-peak separation angle between the polarized currents is tunable either by incident energy or by barrier height each of which is controlled by using top and back gate voltages. The results can be used for constructing an electron beam splitter to produce valley-polarized currents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermal Analysis Of Saltstone Disposal Units

SRR requested that SRNL develop a two-dimensional computer simulation model to be capable of predicting grout temperatures attained in Saltstone Disposal Units (SDUs) during the transient addition of saltstone grout and its subsequent storage period. A simulation model was developed to calculate the transient grout temperatures of the SDU facility by taking an axisymmetric and two-dimensional Computational Fluid Dynamics (CFD) approach. The final performance calculations were based on the conduction and radiation model without a convection cooling mechanism for conservative estimate of the grout temperatures and computational efficiency. While convection within the SDU is disabled in the current calculations, the simulation model has the capability of turning convection back on if desired by the user. Hence, the thermal model in its current configuration is best characterized as a finite-volume heat transfer model. However, given that the model can account for convection if desired by the user, it will be referred to as a CFD model. The grout material contains a transient heat source due to the exothermic hydration reactions of the cementitious material mixed with salt solution. The model can predict the transient temperatures at point locations within the poured grout and the overlying vapor space, and at the interior and exterior surfaces of the structure. The thermal model has been developed for two SDUs at the Savannah River Site (SRS); namely, SDU 2A and SDU 6. SDU 2A is a cylindrical unit measuring 150 feet (ft) in diameter, 22 ft tall, and holds approximately three million gallons of grout. SDU 6 is 375 ft in diameter and 43 ft high with a capacity of over thirty million gallons. Prior to performing the thermal calculations for the SDU 6 facility, the benchmarking test for the smaller SDU 2A unit was performed to establish the solution method and to verify the computational results using existing thermocouple data. Based on the Lagrangian discrete method and boundary conditions established by the SDU 2A benchmarking test, the CFD calculations were performed to predict the transient temperatures at local positions of the poured grout and air spaces within the SDU 6 structure.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of single-event upset of magnetic tunnel junction used in spintronic circuits caused by radiation-induced current

This paper describes the possibility of a switching upset of a magnetic tunnel junction (MTJ) caused by a terrestrial radiation-induced single-event-upset (SEU) current in spintronic integrated circuits. The current waveforms were simulated by using a 3-D device simulator in a basic circuit including MTJs designed using 90-nm CMOS parameters and design rules. The waveforms have a 400 -μA peak and a 200-ps elapsed time when neutron particles with a linear energy transfer value of 14 MeV cm{sup 2}/mg enter the silicon surface. The authors also found that the SEU current may cause soft errors with a probability of more than 10{sup −12} per event, which was obtained by approximate solution of the ordinary differential equation of switching probability when the intrinsic critical current (I{sub C0}) became less than 30 μA.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Radiation transport in kinetic simulations and the influence of photoemission on electron current in self-sustaining discharges

A kinetic description for electronic excitation of helium for principal quantum number n $\leqslant $ 4 has been included into a particle-in-cell (PIC) simulation utilizing direct simulation Monte Carlo (DSMC) for electron-neutral interactions. The excited electronic levels radiate state-dependent photons with wavelengths from the extreme ultraviolet (EUV) to visible regimes. Photon wavelengths are chosen according to a Voigt distribution accounting for the natural, pressure, and Doppler broadened linewidths. This method allows for reconstruction of the emission spectrum for a non-thermalized electron energy distribution function (EEDF) and investigation of high energy photon effects on surfaces, specifically photoemission. A parallel plate discharge with a fixed field (i.e. space charge neglected) is used to investigate the effects of including photoemission for a Townsend discharge. When operating at a voltage near the self-sustaining discharge threshold, it is observed that the electron current into the anode is higher when including photoemission from the cathode than without even when accounting for self-absorption from ground state atoms. As a result, the photocurrent has been observed to account for as much as 20% of the total current from the cathode under steady-state conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Estimating internal moisture generation rates in recently constructed, occupied homes in the southeastern United States

Internal moisture generation (IMG), or moisture generated by building occupants via activities such as respiration, cooking, bathing, and cleaning, is a critical input required for design, analysis, and simulation of building enclosure and heating, ventilation, and air conditioning (HVAC) systems. Based on previously published values, ASHRAE Standard 160-2021 provides guidance for estimating IMG rates for moisture control design analysis, which is based on occupied home datasets collected in the 1980s and 1990s. Residential energy use simulation software also utilizes estimates for IMG for energy rating, energy analysis, and code compliance calculations, as specified in ANSI/RESNET/ICC Standard 301. Data quantifying IMG rates in newer homes is useful in determining the continued relevance of current design and simulation guidance, and whether the guidance represents conditions found in new housing stock. ASHRAE Research Project 1844, conducted by the Florida Solar Energy Center (FSEC), a research institute of the University of Central Florida, estimated IMG rates in newer occupied homes built since 2013 in the southeastern US using a moisture balance model approach. Furthermore, the project obtained occupied home data in cooperation with a US Department of Energy Building America research project to characterize indoor air quality (IAQ) in newer US homes, along with presence, functionality, and occupant use of control measures. Full-scale laboratory homes operating with known IMG rates were used to validate a moisture balance model and quantify the accuracy of estimates obtained when the model is applied to data from occupied homes with unknown IMG rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Simulations of particle interactions in a high-current RFQ

High-intensity beams of composite ions are susceptible to particle interactions that do not affect single ions, such as stripping of electrons and charge exchange. Beam dynamics simulation techniques used in accelerator design do not often include these particle interactions, and so the transmission of the real beam can be significantly different from the design simulations. We have modified the IMPACT-T code to include particle interactions as part of the beam dynamics simulations, including space-charge effects and beam losses. We have carried out validation simulations of interactions in a molecular hydrogen ion beam through an existing high-intensity deuteron RFQ design, as beam commissioning is often carried out using molecular hydrogen beams, to avoid deuteron-deuteron interactions in the accelerator leading to neutron activation. The new code framework successfully models the interactions within the beam, tracks the different source and product ions through the accelerator, and calculates the collective effects such as space charge from all particles together. Finally, we also discuss other simulation codes that could be modified to include similar interaction effects, and possible applications to other types of accelerators.

47 OTHER INSTRUMENTATION↗

The relevance of MRI for patient modeling in head and neck hyperthermia treatment planning: A comparison of CT and CT-MRI based tissue segmentation on simulated temperature

Purpose: In current clinical practice, head and neck (H and N) hyperthermia treatment planning (HTP) is solely based on computed tomography (CT) images. Magnetic resonance imaging (MRI) provides superior soft-tissue contrast over CT. The purpose of the authors’ study is to investigate the relevance of using MRI in addition to CT for patient modeling in H and N HTP. Methods: CT and MRI scans were acquired for 11 patients in an immobilization mask. Three observers manually segmented on CT, MRI T1 weighted (MRI-T1w), and MRI T2 weighted (MRI-T2w) images the following thermo-sensitive tissues: cerebrum, cerebellum, brainstem, myelum, sclera, lens, vitreous humor, and the optical nerve. For these tissues that are used for patient modeling in H and N HTP, the interobserver variation of manual tissue segmentation in CT and MRI was quantified with the mean surface distance (MSD). Next, the authors compared the impact of CT and CT and MRI based patient models on the predicted temperatures. For each tissue, the modality was selected that led to the lowest observer variation and inserted this in the combined CT and MRI based patient model (CT and MRI), after a deformable image registration. In addition, a patient model with a detailed segmentation of brain tissues (including white matter, gray matter, and cerebrospinal fluid) was created (CT and MRI{sub db}). To quantify the relevance of MRI based segmentation for H and N HTP, the authors compared the predicted maximum temperatures in the segmented tissues (T{sub max}) and the corresponding specific absorption rate (SAR) of the patient models based on (1) CT, (2) CT and MRI, and (3) CT and MRI{sub db}. Results: In MRI, a similar or reduced interobserver variation was found compared to CT (maximum of median MSD in CT: 0.93 mm, MRI-T1w: 0.72 mm, MRI-T2w: 0.66 mm). Only for the optical nerve the interobserver variation is significantly lower in CT compared to MRI (median MSD in CT: 0.58 mm, MRI-T1w: 1.27 mm, MRI-T2w: 1.40 mm). Patient models based on CT (T{sub max}: 38.0 °C) and CT and MRI (T{sub max}: 38.1 °C) result in similar simulated temperatures, while CT and MRI{sub db} (T{sub max}: 38.5 °C) resulted in significantly higher temperatures. The SAR corresponding to these temperatures did not differ significantly. Conclusions: Although MR imaging reduces the interobserver variation in most tissues, it does not affect simulated local tissue temperatures. However, the improved soft-tissue contrast provided by MRI allows generating a detailed brain segmentation, which has a strong impact on the predicted local temperatures and hence may improve simulation guided hyperthermia.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Mechanical Analysis of W78/88-1 Life Extension Program Warhead Design Options

Life Extension Program (LEP) is a program to repair/replace components of nuclear weapons to ensure the ability to meet military requirements. The W78/88-1 LEP encompasses the modernization of two major nuclear weapon reentry systems into an interoperable warhead. Several design concepts exist to provide different options for robust safety and security themes, maximum non-nuclear commonality, and cost. Simulation is one capability used to evaluate the mechanical performance of the designs in various operational environments, plan for system and component qualification efforts, and provide insight into the survivability of the warhead in environments that are not currently testable. The simulation efforts use several Sandia-developed tools through the Advanced Simulation and Computing program, including Cubit for mesh generation, the DART Model Manager, SIERRA codes running on the HPC TLCC2 platforms, DAKOTA, and ParaView. Several programmatic objectives were met using the simulation capability including: (1) providing early environmental specification estimates that may be used by component designers to understand the severity of the loads their components will need to survive, (2) providing guidance for load levels and configurations for subassembly tests intended to represent operational environments, and (3) recommending design options including modified geometry and material properties. These objectives were accomplished through regular interactions with component, system, and test engineers while using the laboratory's computational infrastructure to effectively perform ensembles of simulations. Because NNSA has decided to defer the LEP program, simulation results are being documented and models are being archived for future reference. However, some advanced and exploratory efforts will continue to mature key technologies, using the results from these and ongoing simulations for design insights, test planning, and model validation.

42 ENGINEERING↗

Validation of Simulation Models for Wall-Connected Scrape off Layer Currents during MHD Activity in the HBT-EP Tokamak (Final Report)

This is the final report for the DOE award DE-SC0021325, titled “Validation of Simulation Models for Wall-Connected Scrape off Layer Currents during MHD Activity in the HBT-EP Tokamak” for the period September 1, 2020 – May 31, 2023. This award is part of a broader collaborative project, which supported training of a PhD student at Columbia University under the joint supervision of Dr. Hansen and Columbia project members. The award supported research into computational models for plasma and 3D conducting structures in the vicinity relevant to tokamak disruptions and other mode activity (eg. RWMs) using the NIMROD, and PSI-Tet, codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Developing a simulation tool for evaluating in-motion detector systems [Slides]

In order to combat the illicit trafficking of nuclear materials, Radiation Detection Systems (RDS) are used to scan pedestrians, vehicles, cargo, luggage, and other conveyances. Many of these systems are used in an “in-motion” configuration where the detector or target are moving past the other. Current modeling and simulation tools lack the ability to produce time-series detector response profiles for in-motion applications.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

A 3D simulation comparison of carrier transport in green and blue c-plane multi-quantum well nitride light emitting diodes

Until recently, the electrical efficiency of green nitride light-emitting diodes (LEDs) was considerably lower than that of blue LEDs. This is particularly surprising as one would expect a reduced forward voltage with increasing emission wavelength. In this paper, we theoretically investigated the impact of the number of quantum wells on the forward voltage of III-nitride LEDs with x = 0.15 (blue) and x = 0.24 (green) In x Ga1–xN QWs. The simulated dependence of current density (J) on applied diode bias (V) shows a significant increase of 1.9 V in the forward voltage between one and five quantum well (QW) c-plane green LED structures. Artificially turning off the polarization fields in the simulation does not entirely suppress this effect. Due to the large band offsets in the green LED multiple QW stack, simulations indicate a sequential band filling of the QW sequence. Furthermore, this mechanism should not be limited to c-plane LEDs and could also be present in nonpolar or semipolar devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

3D electrostatic gyrokinetic electron and fully kinetic ion simulation of lower-hybrid drift instability of Harris current sheet

The eigenmode stability properties of three-dimensional lower-hybrid-drift-instabilities (LHDI) in a Harris current sheet with a small but finite guide magnetic field have been systematically studied by employing the gyrokinetic electron and fully kinetic ion (GeFi) particle-in-cell (PIC) simulation model with a realistic ion-to-electron mass ratio m i /m e . In contrast to the fully kinetic PIC simulation scheme, the fast electron cyclotron motion and plasma oscillations are systematically removed in the GeFi model, and hence one can employ the realistic m i /m e . The GeFi simulations are benchmarked against and show excellent agreement with both the fully kinetic PIC simulation and the analytical eigenmode theory. Our studies indicate that, for small wavenumbers, ky, along the current direction, the most unstable eigenmodes are peaked at the location where $\vec{k}$• $\vec{B}$ =0, consistent with previous analytical and simulation studies. Here, $\vec{B}$ is the equilibrium magnetic field and $\vec{k}$ is the wavevector perpendicular to the nonuniformity direction. As ky increases, however, the most unstable eigenmodes are found to be peaked at $\vec{k}$ •$\vec{B}$ ≠0. Additionally, the simulation results indicate that varying m i /m e , the current sheet width, and the guide magnetic field can affect the stability of LHDI. Simulations with the varying mass ratio confirm the lower hybrid frequency and wave number scalings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Counter-Current Flow Limitation Studies in Complex Geometries Utilizing Interface Capturing Simulations Coupled with PID Flow Rate Controller

In nuclear thermal-hydraulic studies, counter-current flow limitation (CCFL) typically refers to steam rising at a fast rate such that it prevents coolant from draining down within a confined channel. CCFL is a crucial issue in nuclear reactor safety analysis. This study investigates CCFL in debris bed channels using high-resolution interface-capturing simulations. A novel proportional-integral-derivative flow rate controller is developed to efficiently achieve the CCFL conditions. Verification studies confirm that CCFL occurs under the same conditions with or without the controller, demonstrating that PID control ensures accurate prediction. Three debris bed channel geometries were examined: a cylindrical channel, a channel with small obstacles, and a channel with large obstacles. Results show that obstacles significantly impact flow behavior, interfacial shear, wall shear, and pressure gradients required for CCFL. Furthermore, the comparison with experimental data confirmed that simulations incorporating geometric complexities align more closely with experimental CCFL conditions. A pressure gradient correlation was also developed for CCFL prediction.

Counter-current flow limitation↗

Analyzing the Effects of the Interference between One- and Two-Body Currents on Neutrino-Nucleus Scattering in ACHILLES

Modeling improvements of neutrino–nucleus scattering are an essential aspect of reducing systematic uncertainties for current and future neutrino oscillation experiments. In particular, retaining quantum mechanical effects such as the interference between one-body (1p1h) and two-body (2p2h) currents in our cross-section models is an essential consideration. This work investigates the impact of this interference term on neutrino–argon scattering using the ACHILLES event generator, which directly implements the interference model developed by Lovato, Rocco, and Steinberg. Using fluxes relevant to SBND, the DUNE near detector, and MiniBooNE, charged- and neutral-current event samples are generated and predicted event rates are plotted against a variety of kinematic variables. For events simulated using the SBND flux, it is found that interference contributes approximately 10 \% to the total charged-current event rate, and serves to enhance the quasi-elastic peak. When considering events generated using the DUNE flux, resonance interactions contribute more to the predicted event rate, as the neutrino flux extends out to higher energies. The relative interference contribution subsequently decreases to roughly 7 \%. Neutral-current simulations show a higher rate of neutron knockout. Finally, a comparison of argon-to-carbon cross-section ratios is performed utilizing fluxes from SBND and MiniBooNE.

Pando, Nathan [Fermilab]↗

Analyzing the Effects of the Interference between One- and Two-Body Currents on Neutrino-Nucleus Scattering in ACHILLES

Modeling improvements of neutrino–nucleus scattering are an essential aspect of reducing systematic uncertainties for current and future neutrino oscillation experiments. In particular, retaining quantum mechanical effects such as the interference between one-body (1p1h) and two-body (2p2h) currents in our cross-section models is an essential consideration. This work investigates the impact of this interference term on neutrino–argon scattering using the ACHILLES event generator, which directly implements the interference model developed by Lovato, Rocco, and Steinberg. Using fluxes relevant to SBND, the DUNE near detector, and MiniBooNE, charged- and neutral-current event samples are generated and predicted event rates are plotted against a variety of kinematic variables. For events simulated using the SBND flux, it is found that interference contributes approximately 10 \% to the total charged-current event rate, and serves to enhance the quasi-elastic peak. When considering events generated using the DUNE flux, resonance interactions contribute more to the predicted event rate, as the neutrino flux extends out to higher energies. The relative interference contribution subsequently decreases to roughly 7 \%. Neutral-current simulations show a higher rate of neutron knockout. Finally, a comparison of argon-to-carbon cross-section ratios is performed utilizing fluxes from SBND and MiniBooNE.

Pando, Nathan [Fermilab]↗

OBSERVATIONAL SIGNATURES OF A MASSIVE DISTANT PLANET ON THE SCATTERING DISK

The orbital element distribution of trans-Neptunian objects (TNOs) with large pericenters has been suggested to be influenced by the presence of an undetected, large planet at >200 au from the Sun. To find additional observables caused by this scenario, we present here the first detailed emplacement simulation in the presence of a massive ninth planet on the distant Kuiper Belt. We perform 4 Gyr N -body simulations with the currently known solar system planetary architecture, plus a 10 M {sub ⊕} planet with similar orbital parameters to those suggested by Trujillo and Sheppard or Batygin and Brown, and 10{sup 5} test particles in an initial planetesimal disk. We find that including a distant super-Earth-mass planet produces a substantially different orbital distribution for the scattering and detached TNOs, raising the pericenters and inclinations of moderate semimajor axis (50 < a < 500 au) objects. We test whether this signature is detectable via a simulator with the observational characteristics of four precisely characterized TNO surveys. We find that the qualitatively very distinct solar system models that include a ninth planet are essentially observationally indistinguishable from an outer solar system produced solely by the four giant planets. We also find that the mass of the Kuiper Belt’s current scattering and detached populations is required to be 3–10 times larger in the presence of an additional planet. We do not find any evidence for clustering of orbital angles in our simulated TNO population. Wide-field, deep surveys targeting inclined high-pericenter objects will be required to distinguish between these different scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗

Scaling magnetized liner inertial fusion on Z and future pulsed-power accelerators

The MagLIF (Magnetized Liner Inertial Fusion) concept [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)] has demonstrated fusion–relevant plasma conditions [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)] on the Z accelerator with a peak drive current of about 18 MA. We present 2D numerical simulations of the scaling of MagLIF on Z as a function of drive current, preheat energy, and applied magnetic field. The results indicate that deuterium-tritium (DT) fusion yields greater than 100 kJ could be possible on Z when all of these parameters are at the optimum values: i.e., peak current = 25 MA, deposited preheat energy = 5 kJ, and B{sub z} = 30 T. Much higher yields have been predicted [S. A. Slutz and R. A. Vesey, Phys. Rev. Lett. 108, 025003 (2012)] for MagLIF driven with larger peak currents. Two high performance pulsed-power accelerators (Z300 and Z800) based on linear-transformer-driver technology have been designed [W. A. Stygar et al., Phys. Rev. ST Accel. Beams 18, 110401 (2015)]. The Z300 design would provide 48 MA to a MagLIF load, while Z800 would provide 65 MA. Parameterized Thevenin-equivalent circuits were used to drive a series of 1D and 2D numerical MagLIF simulations with currents ranging from what Z can deliver now to what could be achieved by these conceptual future pulsed-power accelerators. 2D simulations of simple MagLIF targets containing just gaseous DT have yields of 18 MJ for Z300 and 440 MJ for Z800. The 2D simulated yield for Z800 is increased to 7 GJ by adding a layer of frozen DT ice to the inside of the liner.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗