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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Fast QSTS for Distributed PV Impact Studies Using Vector Quantization and Variable Time-Steps

As the penetration level of PV systems increases, utilities are in great need of methods to quickly and effectively analyze PV's impact on electric distribution systems. Compared to other conventional methods, quasi-static time-series (QSTS) based analysis has demonstrated its advantages of accurately capturing the impact of interest, e.g. control device operations, and min and max voltages. However, yearlong QSTS simulation with 1-second granularity is currently too computationally expensive for regular use. Therefore, vector quantization (VQ) and variable time-step (VTS) solution methods have been developed to accelerate QSTS simulation. This paper proposes the combination of VQ with VTS. We investigate the effectiveness of implementing VQ in combination with 2 VTS algorithms, predetermined time-step (PT) and VTS with backtrack downsampling (VTS-BD) to demonstrate a runtime reduction of 99.55% with error metrics below acceptable thresholds.

14 SOLAR ENERGY↗
DOE OSTI2023

Enhance BISON Metal Fuel Transient Models and Enable Assessment Based on Out-of-Pile Transient Experiments

In this report, the activities and achievements made by Argonne National Laboratory for the Nuclear Energy Advanced Modeling and Simulation (NEAMS) BISON code metallic fuel validation and verification (V&V) project in FY2023, with a focus on the capability development and assessment case development for the out-of-pile transient experiments, are summarized. The liquid cladding penetration model was improved based on a comprehensive set of fuel behavior test apparatus (FBTA) test results. The new model was implemented into BISON along with a simple model to account for fuel melting effects. The metallic fuel meshing capability was also upgraded to allow local refinement using unstructured triangular mesh. Based on the previous success in developing the BISON-FIPD integration powered BISON assessment case, this database enabled BISON evaluation and V&V capability was expanded to the out-of-pile transient experiments by introducing BISON-OPTD integration. A series of initial out-of-pile transient assessment cases were implemented into BISON to provide powerful platforms to evaluate current BISON’s capabilities in simulating metallic fuel performance during out-of-pile transients. The established capabilities will also benefit the calibration and V&V of future advanced mechanistic models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Applying the Gyrokinetic Formulation of Magnetized Plasma Dynamics to Z-pinch plasmas

This project adapts and applies the gyrokinetic Eulerian code COGENT to an emerging area of highly-magnetized Z-pinch plasma systems. Simulation efforts for these plasmas currently use either magneto-hydrodynamic (MHD) fluid codes, which suffer from fidelity being restricted to narrow plasma conditions, or fully kinetic particle-in-cell (PIC) methods, which suffer from the run-time needed to resolve the very small spatial and temporal scales required for the plasma conditions. The progress with achieving detailed theoretical understanding of Z-pinch plasma systems can be greatly expedited using some combination of gyrokinetic and drift-ordered MHD models, developed in the tokamak community and implemented in the COGENT code. These formulations are applicable to a much larger parameter space than standard MHD models, while providing a substantial speed-up in run-time over fully kinetic codes. In particular, the gyrokinetic formalism is successfully applied to a flow-stabilized Z pinch (FSZP) system – a promising fusion concept where a sheared axial flow is used to provide stability against large-scale ideal-MHD modes. Results of gyrokinetic simulations with COGENT agree well with those from fully kinetic simulations, while requiring ~500x less run-time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

Hybrid MD-Nernst Planck Model of Alpha-hemolysin Conductance Properties

Motivated by experiments in which an applied electric field translocates polynucleotides through an alpha-hemolysin protein channel causing ionic current transient blockade, a hybrid simulation model is proposed to predict the conductance properties of the open channel. Time scales corresponding to ion permeation processes are reached using the Poisson-Nemst-Planck (PNP) electro-diffusion model in which both solvent and local ion concentrations are represented as a continuum. The diffusion coefficients of the ions (K(+) and Cl(-)) input in the PNP model are, however, calculated from all-atom molecular dynamics (MD). In the MD simulations, a reduced representation of the channel is used. The channel is solvated in a 1 M KCI solution, and an external electric field is applied. The pore specific diffusion coefficients for both ionic species are reduced 5-7 times in comparison to bulk values. Significant statistical variations (17-45%) of the pore-ions diffusivities are observed. Within the statistics, the ionic diffusivities remain invariable for a range of external applied voltages between 30 and 240mV. In the 2D-PNP calculations, the pore stem is approximated by a smooth cylinder of radius approx. 9A with two constriction blocks where the radius is reduced to approx. 6A. The electrostatic potential includes the contribution from the atomistic charges. The MD-PNP model shows that the atomic charges are responsible for the rectifying behaviour and for the slight anion selectivity of the a-hemolysin pore. Independent of the hierarchy between the anion and cation diffusivities, the anionic contribution to the total ionic current will dominate. The predictions of the MD-PNP model are in good agreement with experimental data and give confidence in the present approach of bridging time scales by combining a microscopic and macroscopic model.

Cozmuta, Ioana↗
DOE OSTI2019

Fabrication and Tests of a Resistive-Type Superconducting Fault Current Limiter Module Based on Coated Conductors

The resistivetype superconducting fault current limiter (R-SFCL) is simple in structure, near zero impedance at normal state and can be quickly triggered by the fault current. Thus, it is a possible solution to limit the fault current in the power grid. In this paper, we report the fabrication and test for an R-SFCL module The R-SFCL module consists of seven bifilar coils connected in parallel, which were wound with two parallel 2G tapes. A 5mm spacer was co-wound between adjacent turns to act as the highvoltage insulation. In the tests, we test the current distribution among the bifilar coils, and then simulated short-circuit tests were conducted on the module to obtain its performance under fault current surge.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Battery state of health monitoring by estimation of the number of cyclable Li-ions

Herein this paper introduces a method to monitor battery state of health (SOH) by estimating the number of cyclable Li-ions, a health-relevant electrochemical variable. SOH monitoring is critical to battery management in balancing the trade-off between maximizing system performance and minimizing battery degradation. The decrease of cyclable Li-ions indicates the effect on the SOH of degradation mechanisms that consume cyclable Li-ions. The unavailability of the number of cyclable Li-ions through non-invasive measurements makes its estimation necessary for in-situ SOH monitoring. In this paper, the extended Kalman filter (EKF) is used to estimate the number of cyclable Li-ions as an unknown battery parameter. The single particle model (SPM), a simplified battery electrochemical model, is used as the model in the EKF to achieve a computational complexity suitable for on-line estimation. Simulations are performed under typical electric vehicle current trajectories using an example parameter set for a hybrid-electric-vehicle battery. In the simulations, the battery is represented by the Doyle–Fuller–Newman (DFN) model, an electrochemical model with higher fidelity than the SPM. To comply with the practice, instead of using the same parameters as the DFN model in the SPM, parameterization of the SPM is performed before estimation of the number of cyclable Li-ions. The simulations show high estimation accuracy of the number of cyclable Li-ions using the EKF, even with the structural and parametric differences between the DFN model and the SPM, under both the ideal conditions and various non-ideal conditions (i.e., SOC estimation error, additional modeling error, and measurement noise).

25 ENERGY STORAGE↗
DOE OSTI2017

Comparison of the Three NIF Ablators

Indirect drive implosion experiments on NIF have now been performed using three different ablator materials: glow discharge polymer (GDP) or CH, high density carbon (HDC, which we also refer to as diamond), and sputtered beryllium (Be). It has been appreciated for some time that each of these materials has specific advantages and disadvantages as an ICF ablator.[1-4] In light of experiments conducted on NIF in the last few years, how do these ablators compare? Given current understanding, is any ablator more or less likely to reach ignition on NIF? Has the understanding of their respective strengths and weaknesses changed since NIF experiments began? How are those strengths and weaknesses highlighted by implosion designs currently being tested or planned for testing soon? This document aims to address these questions by combining modern simulation results with a survey of the current experimental data base. More particularly, this document is meant to fulfill an L2 Milestone for FY17 to “Document our understanding of the relative advantages and disadvantages of CH, HDC, and Be designs.” Note that this document does not aim to recommend a down-selection of the current three ablator choices. It is intended only to gather and document the current understanding of the differences between these ablators and thereby inform the choices made in planning future implosion experiments. This document has two themes: (i) We report on a reanalysis project in which post-shot simulations were done on a common basis for layered shots using each ablator. This included data from keyholes, 2D ConA, and so forth, from each campaign, leading up to the layered shots. (“Keyholes” are shots dedicated to measuring the shock timing in a NIF target, as described in Ref. 5. “2DConAs” are backlit implosions in which the symmetry of the implosion is measured between about half and full convergence, as described in Ref. 6.) This set of common-basis postshot simulations is compared to the respective shots. Each was then scaled to a “full NIF” experiment that could be done using the respective ablators at full NIF power and/or energy, and these scaled-up designs were simulated in detail. (ii) The report also contains a general survey of experimental and simulated results as pertinent to comparing and evaluating the three ablators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Temperature coefficients of multijunction solar cells

Temperature coefficients measured in solar simulators with those measured under AM0 solar illumination are compared to illustrate the challenges in making these measurements. It is shown that simulator measurements of the short-circuit current (delta Jsc/delta T) are inaccurate due to the mismatch between the solar spectrum and the simulators at the bandgaps of the solar cells. Especially susceptible to error is the delta Jsc/delta T of cells which are components in monolithic multijunction solar cells, such as GaAs filtered by 1.93-eV AlGaAs, which has an AM0 coefficient of 6.82 micro-A/sq cm/deg C, compared to a Xenon simulator coefficient of 22.2 micro-A/sq cm/deg C.

Virshup, G. F.↗
NASA NTRS2009

Large Eddy Simulation of Transonic Flow Field in NASA Rotor 37

The current paper reports on numerical investigations on the flow characteristics in a transonic axial compressor, NASA Rotor 37. The flow field was used previously as a CFD blind test case conducted by American Society of Mechanical Engineers in 1994. Since the CFD blind-test exercise, many numerical studies on the flow field in the NASA Rotor 37 have been reported. Although steady improvements have been reported in both numerical procedure and turbulence closure, it is believed that all the important aspects of the flow field have not been fully explained with numerical studies based on the Reynolds Averaged Navier-Stokes (RANS) solution. Experimental data show large dip in total pressure distribution near the hub at downstream of the rotor at 100% rotor speed. Most original numerical solutions from the blind test exercise did not predict this total pressure deficit correctly. This total pressure deficit at the rotor exit was attributed to a hub corner flow separation by the author. Several subsequent numerical studies with different turbulence closure model also calculated this dip in total pressure rise. Also, several studies attributed this total pressure deficit to a small leakage flow coming from the hub in the test article. As the experimental study cannot be repeated, either explanation cannot be validated. The primary purpose of the current investigation is to investigate the transonic flow field with both RANS and a Large Eddy Simulation (LES). The RANS approach gives similar results presented at the original blind test exercise. Although the RANS calculates higher overall total pressure rise, the total pressure deficit near the hub is calculated correctly. The numerical solution shows that the total pressure deficit is due to a hub corner flow separation. The calculated pressure rise from the LES agrees better with the measured total pressure rise especially near the casing area where the passage shock interacts with the tip clearance vortex and flow becomes unsteady due to this interaction. The LES simulation also calculates the total pressure rise deficit near the hub and it agrees well with the measured data.

Hah, Chunill↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

Radar Autofocus Algorithm Incorporating a priori Terrain Knowledge for Correction of Mars’ Ionospheric Distortion in MARSIS Observations

Low-frequency subsurface radar observations of Mars’ polar ice deposits by MARSIS (Mars Advanced Radar for Subsurface and Ionosphere Sounding) are heavily impacted by the electron content of Mars’ ionosphere. The resulting ionospheric distortion can be represented as attenuation and bulk delay, in addition to higher-order frequency dispersion effects. Baseline, uncorrected images are often unusable when the solar zenith angle is less than 90◦ (day side). In this work, a radar autofocus algorithm is developed that estimates and inverts ionospheric distortion, producing a focused radargram of the ice deposit subsurface. Previously published autofocus algorithms have sought to maximize peak-to-noise contrast, which may yield sub-optimal results for complex terrain. Instead, a maximum likelihood approach is developed that incorporates simulated surface clutter returns for the current spacecraft position, based on the Mars Orbiter Laser Altimeter (MOLA) elevation model of the Martian surface. An ancillary product is a surface-only clutter simulation for each orbit, which is necessary to identify true subsurface features.

McMichael, Joseph G.↗
DOE OSTI2018

Analysis of Dynamics of Plasma Current Quench in the Globus-M Spherical Tokamak

Data on the dynamics of the plasma current quench in the Globus-M tokamak are presented. The main current quench characteristics at different toroidal magnetic fields are compared. The distribution of the toroidal current induced in the vessel wall is determined from magnetic measurements, and the electromagnetic loads acting on the vessel wall during the current quench are calculated. By extrapolating the experimental data, the additional pressure on the vessel wall during the current quench in the upgraded Globus-M2 tokamak is estimated. It is shown that the current quench results in the appearance of bending stresses in the vessel domes. Using numerical simulations, it is shown that the best agreement between the measured and calculated plasma current dynamics during the current quench corresponds to the linear (in time) influx of the carbon impurity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Safety Analysis for the Versatile Test Reactor Conceptual Design

Safety analysis simulations have been carried out using the SAS4A/SASSYS-1 safety analysis code system to update the baseline set of results for the Versatile Test Reactor conceptual design. Results from a previous revision of this work provided input for the Conceptual Safety Design Report summarizing the overall safety basis for the VTR. This work will also support future analyses that will be required for Chapter 13, “Accident Analysis,” of a preliminary safety analysis report. A number of transients have been evaluated based on categories of postulated accident sequences suggested for consideration in NUREG-1537. Ultimately, the full spectrum of event initiators and subsequent accident sequences will be defined through a probabilistic risk assessment of the plant design following the implementation of the Licensing Modernization Project guidelines as described in the VTR Safety Design Strategy. Since that process is ongoing, the present analysis addresses several transients that are expected to be bounding accident scenarios that represent the three key ways to perturb a reactor: through changes to the core inlet temperature, mass flow rate, or reactivity. These correspond to a loss of heat sink (LOHS), loss of flow (LOF) or station blackout (SBO), and transient overpower (TOP), respectively. Both Protected and Unprotected versions of these transients have been evaluated, along with additional transients that include transient overpower with one stuck rod, pump coast-down failures, seismic events, and overcooling events. At the current stage of design, transient simulation results for the Versatile Test Reactor indicate that large safety margins exist for many event initiators. However, several “enabling” assumptions have been made, in terms of both design features and design limits, in order to perform the analyses; and these assumptions will need to be revised as the design matures. In addition, a larger spectrum of events needs to be evaluated as part of the on-going probabilistic risk assessment process.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Modeling and analysis of heat emissions from buildings to ambient air

Heat emissions from buildings is a significant source of anthropogenic heat influencing the urban microclimate; however, they are usually oversimplified in urban climate and microclimate modeling. This study developed a bottom-up physics-based approach to calculate heat emissions from buildings to the ambient air and implemented the approach in EnergyPlus. In this work, a simple result verification was conducted by comparing the EnergyPlus simulated results against the spreadsheet calculations. Simulations covering 16 commercial building types, four climates, and two energy efficiency levels were conducted to understand and evaluate the building heat emissions and their temporal patterns as well as three major components: (1) building envelope (convective heat transfer to ambient air), (2) zones (air exfiltration and exhaust air), and (3) HVAC systems (relief air and heat rejection from condensers or cooling towers). The main findings are: (1) heat emissions are usually higher than the site energy use (about 2.5 times), and their dynamics should be considered; (2) building characteristics and their energy systems lead to differences in heat emission contributions from the three components, and their dynamics, for example, in the warehouse models, the envelope component accounts for 90.4%, while it is 12.7% for the large office models; (3) for most building typologies, the climate has a strong impact on heat emissions, for example, buildings with dominant heat emissions from the zone exhaust air and/or the HVAC reject heat, a general decrease in heat emissions in hotter climates is observed, while envelope-dominated buildings show the opposite; and (4) building technologies that reduce energy use in buildings may perform differently in reducing heat emissions. The developed heat emissions calculation method can be adopted in EnergyPlus and most other building energy modeling programs. It can provide dynamic building heat emissions as an input to urban climate computational fluid dynamics (CFD) models at a higher spatial and temporal resolution than is currently available, to improve the simulation accuracy of the urban microclimate and capture the urban heat island effect and urban overheating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2014

Aeromechanical Evaluation of Smart-Twisting Active Rotor

An investigation of Smart-Twisting Active Rotor (STAR) was made to assess potential benefits of the current active twist rotor concept for performance improvement, vibration reduction, and noise alleviation. The STAR rotor is a 40% Mach-scaled, Bo105 rotor with an articulated flap-lag hinge at 3.5%R and no pre-cone. The 0-5 per rev active twist harmonic inputs were applied for various flight conditions including hover, descent, moderate to high speed level flights, and slowed rotor high advance ratio. For the analysis, the STAR partners used multiple codes including CAMRAD II, S4, HOST, rFlow3D, elsA, and their associated software. At the high thrust level in hover, the 0 per rev active twist with 80% amplitude increased figure of merit (FM) by 0.01-0.02 relative to the baseline. In descent, the largest BVI noise reduction was on the order of 2 to 5 dB at the 3 per rev active twist. In the high speed case (mu = 0.35), the 2 per rev actuation was found to be the most effective in achieving a power reduction as well as a vibration reduction. At the 2 per rev active twist, total power was reduced by 0.65% at the 60 deg active twist phase, and vibration was reduced by 47.6% at the 45 deg active twist phase. The use of the 2 per rev active twist appears effective for vibration reduction. In the high advance ratio case (mu = 0.70), the 0 per rev actuation appeared to have negligible impact on performance improvement. In summary, computational simulations successfully demonstrated that the current active twist concept provided a significant reduction of the maximum BVI noise in descent, a significant reduction of the vibration in the high speed case, a small improvement on rotor performance in hover, and a negligible impact on rotor performance in forward flight.

Lim, Joon W.↗
DOE OSTI2021

Gate-free state preparation for fast variational quantum eigensolver simulations

Abstract The variational quantum eigensolver is currently the flagship algorithm for solving electronic structure problems on near-term quantum computers. The algorithm involves implementing a sequence of parameterized gates on quantum hardware to generate a target quantum state, and then measuring the molecular energy. Due to finite coherence times and gate errors, the number of gates that can be implemented remains limited. In this work, we propose an alternative algorithm where device-level pulse shapes are variationally optimized for the state preparation rather than using an abstract-level quantum circuit. In doing so, the coherence time required for the state preparation is drastically reduced. We numerically demonstrate this by directly optimizing pulse shapes which accurately model the dissociation of H 2 and HeH + , and we compute the ground state energy for LiH with four transmons where we see reductions in state preparation times of roughly three orders of magnitude compared to gate-based strategies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

Investigating Galactic Structure with COBE/DIRBE and Simulation

In this work I applied the current version of the SKY model of the point source sky to the interpretation of the diffuse all-sky emission observed by COBE/DIRBE (Cosmic Background Explorer Satellite/Diffuse Infrared Background Experiment). The goal was to refine the SKY model using the all-sky DIRBE maps of the Galaxy, in order that a search could be made for an isotropic cosmic background."Faint Source Model" [FSM] was constructed to remove Galactic fore ground stars from the ZSMA products. The FSM mimics SKY version 1 but it was inadequate to seek cosmic background emission because of the sizeable residual emission in the ZSMA products after this starlight subtraction. At this point I can only support that such models are currently inadequate to reveal a cosmic background. Even SKY5 yields the same disappointing result.

Cohen, Martin↗