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,711 records · Page 95

Two-temperature effects in Hall-MHD simulations of the HIT-SI experiment

A two-temperature Hall-magnetohydrodynamic (Hall-MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet 3D MHD code and used to model plasma dynamics in the Helicity Injected Torus–Steady Inductive (HIT-SI) experiment. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5–68.5 kHz range. At all frequencies, the NIMROD two-temperature model results in increased toroidal current, lower chord-averaged density, higher average temperatures, outward radial shift of the current centroid, and axial symmetrization of the current centroid, relative to the single-temperature NIMROD simulations. The two-temperature PSI-Tet model illustrates similar trends, but at high frequency operation, it exhibits lower electron temperature, smaller toroidal current, and decreased axial symmetrization with respect to the single-temperature PSI-Tet model. With all models, average temperatures and toroidal currents increase with the injector frequency. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and illustrate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for artificial diffusivity, wall temperature, and density. Both artificial diffusivity and the density boundary condition significantly modify the plasma density profiles, leading to larger average temperatures, toroidal current, and relative density fluctuations at low densities. As a result, a low density simulation achieves sufficiently high current gain (G > 5) to generate significant volumes of closed flux lasting 1–2 injector periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Assessment of Antarctic Circumpolar Current and Southern Ocean Meridional Overturning Circulation During 1958-2007 in a Suite of Interannual CORE-2 Simulations

In the framework of the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II), we present an analysis of the representation of the Antarctic Circumpolar Current (ACC) and Southern Ocean Meridional Overturning Circulation (MOC) in a suite of seventeen global ocean-sea ice models. We focus on the mean, variability and trends of both the ACC and MOC over the 19582007 period, and discuss their relationship with the surface forcing. We aim to quantify the degree of eddy saturation and eddy compensation in the models participating in CORE-II, and compare our results with available observations, previous fine-resolution numerical studies and theoretical constraints. Most models show weak ACC transport sensitivity to changes in forcing during the past five decades, and they can be considered to be in an eddy saturated regime. Larger contrasts arise when considering MOC trends, with a majority of models exhibiting significant strengthening of the MOC during the late 20th and early 21st century. Only a few models show a relatively small sensitivity to forcing changes, responding with an intensified eddy-induced circulation that provides some degree of eddy compensation, while still showing considerable decadal trends. Both ACC and MOC interannual variability are largely controlled by the Southern Annular Mode (SAM). Based on these results, models are clustered into two groups. Models with constant or two-dimensional (horizontal) specification of the eddy-induced advection coefficient show larger ocean interior decadal trends, larger ACC transport decadal trends and no eddy compensation in the MOC. Eddy-permitting models or models with a three-dimensional time varying show smaller changes in isopycnal slopes and associated ACC trends, and partial eddy compensation. As previously argued, a constant in time or space is responsible for a poor representation of mesoscale eddy effects and cannot properly simulate the sensitivity of the ACC and MOC to changing surface forcing. Evidence is given for a larger sensitivity of the MOC as compared to the ACC transport, even when approaching eddy saturation. Future process studies designed for disentangling the role of momentum and buoyancy forcing in driving the ACC and MOC are proposed.

Southern Ocean dynamics↗

An active learning framework for the rapid assessment of galvanic corrosion

The current present in a galvanic couple can define its resistance or susceptibility to corrosion. However, as the current is dependent upon environmental, material, and geometrical parameters it is experimentally costly to measure. To reduce these costs, Finite Element (FE) simulations can be used to assess the cathodic current but also require experimental inputs to define boundary conditions. Due to these challenges, it is crucial to accelerate predictions and accurately predict the current output for different environments and geometries representative of in-service conditions. Machine learned surrogate models provides a means to accelerate corrosion predictions. However, a one-time cost is incurred in procuring the simulation and experimental dataset necessary to calibrate the surrogate model. Therefore, an active learning protocol is developed through calibration of a low-cost surrogate model for the cathodic current of an exemplar galvanic couple (AA7075-SS304) as a function of environmental and geometric parameters. The surrogate model is calibrated on a dataset of FE simulations, and calculates an acquisition function that identifies specific additional inputs with the maximum potential to improve the current predictions. This is accomplished through a staggered workflow that not only improves and refines prediction, but identifies the points at which the most information is gained, thus enabling expansion to a larger parameter space. The protocols developed and demonstrated in this work provide a powerful tool for screening various forms of corrosion under in-service conditions.

36 MATERIALS SCIENCE↗

Evaluation of Portable Acceleration Solutions for LArTPC Simulation Using Wire-Cell Toolkit

The Liquid Argon Time Projection Chamber (LArTPC) technology plays an essential role in many current and future neutrino experiments. Accurate and fast simulation is critical to developing efficient analysis algorithms and precise physics model projections. The speed of simulation becomes more important as Deep Learning algorithms are getting more widely used in LArTPC analysis and their training requires a large simulated dataset. Heterogeneous computing is an efficient way to delegate computationally intensive tasks to specialized hardware. However, as the landscape of compute accelerators quickly evolves, it becomes increasingly difficult to manually adapt the code to the latest hardware or software environments. A solution which is portable to multiple hardware architectures without substantially compromising performance would thus be very beneficial, especially for long-term projects such as the LArTPC simulations. In search of a portable, scalable and maintainable software solution for LArTPC simulations, we have started to explore high-level portable programming frameworks that support several hardware backends. In this paper, we present our experience porting the LArTPC simulation code in the Wire-Cell Toolkit to NVIDIA GPUs, first with the CUDA programming model and then with a portable library called Kokkos. Preliminary performance results on NVIDIA V100 GPUs and multi-core CPUs are presented, followed by a discussion of the factors affiecting the performance and plans for future improvements.

Yu, Haiwang↗

Knowledge-based simulation for aerospace systems

Knowledge-based techniques, which offer many features that are desirable in the simulation and development of aerospace vehicle operations, exhibit many similarities to traditional simulation packages. The eventual solution of these systems' current symbolic processing/numeric processing interface problem will lead to continuous and discrete-event simulation capabilities in a single language, such as TS-PROLOG. Qualitative, totally-symbolic simulation methods are noted to possess several intrinsic characteristics that are especially revelatory of the system being simulated, and capable of insuring that all possible behaviors are considered.

Will, Ralph W.↗

What Is Unusual About the Third Largest Geomagnetic Storm of Solar Cycle 24?

We report on the solar and interplanetary (IP) causes of the third largest geomagnetic storm (26 August 2018) in solar cycle 24. The underlying coronal mass ejection (CME) originating from a quiescent filament region becomes a 440 km/s magnetic cloud (MC) at 1 au after ∼5 days. The prolonged CME acceleration (for ∼24 hr) coincides with the time profiles of the post-eruption arcade intensity and reconnected flux. Chen et al. (2019, https://doi.org/10.3847/1538-4357/ab3f36) obtain a lower speed since they assumed that the CME does not accelerate after ∼12 hr. The presence of multiple coronal holes near the filament channel and the high-speed wind from them seem to have the combined effect of producing complex rotation in the corona and IP medium resulting in a high-inclination MC. The Dst time profile in the main phase steepens significantly (rapid increase in storm intensity) coincident with the density increase (prominence material) in the second half of the MC. Simulations using the Comprehensive Inner Magnetosphere-Ionosphere model show that a higher ring current energy results from larger dynamic pressure (density) in MCs. Furthermore, the Dst index is highly correlated with the main-phase time integral of the ring current injection that includes density, consistent with the simulations. A complex temporal structure develops in the storm main phase if the underlying MC has a complex density structure during intervals of southward IP magnetic field. We conclude that the high intensity of the storm results from the prolonged CME acceleration, complex rotation of the CME flux rope, and the high density in the 1-au MC.

N. Gopalswamy↗

ExaWind: Then and Now

The scientific goal of the ExaWind project is to advance our fundamental understanding of the flow physics governing whole wind plant performance, including wake formation, complex terrain impacts, and turbine-turbine-interaction effects. The primary application codes in the ExaWind environment are Nalu-Wind, an unstructured-grid computational fluid dynamics (CFD) code, AMR-Wind, a structured-grid CFD code, and OpenFAST, a whole-turbine simulation code. In this poster we present the current status of the ExaWind software stack in the context of the modeling and simulation capabilities when the project started in 2016.

computational fluid dynamics↗

Auroral plasma transport processes in the presence of kV potential structures

We have simulated plasma transport processes in the presence of a quasi-two-dimensional current filament, that generated kV potential structure in the auroral region. The simulation consists of a set of one-dimensional flux tube simulations with different imposed time-dependent, field-aligned currents. The model uses the 16 moment system of equations and simultaneously solves coupled continuity and momentum equations and equations describing the transport along the magnetic field lines of parallel and perpendicular thermal energy and heat flows for each species. The lower end of the simulation is at an altitude of 800 km, in the collisional topside ionosphere, while the upper end is at 10 R(sub E) in the magnetosphere. The plasma consists of hot electrons and protons of magnetospheric origin and low-energy electrons, protons, and oxygen ions of ionospheric origin. The dynamical interaction of the individual current filaments with ionospheric and magnetospheric plasma generates a potential structure in the horizontal direction and kilovolt field-aligned potential drops along the field lines. The side-by-side display exhibits the evolution of the implied potential structure in the horizontial direction. In the presence of this potential structure and parallel electric field ionospheric plasma density is depleted and velocity is reduced, while density enhancement and increased velocity is observed in magnetospheric plasma. The ionospheric and magnetospheric electron temperatures increase below 2 R(sub E) due to magnetic mirror force on converging geomagnetic field lines. The primary cross-field motion produced by the horizontal E field (E x B drift) is perpendicular to both of the significant spatial directions and is thus ignorable in this geometry. The effects of other cross-field drift processes are discussed. The simulation thus provides insight into the dynamical evolution of two-dimensional potential structures driven by an imposed finite width, field-aligned current profile.

Ganguli, Supriya B.↗

Effects of electron cyclotron current drive on magnetic islands in tokamak plasmas

The effects of the electron cyclotron current drive on magnetic islands in tokamak plasmas are studied using gyrokinetic simulations. By investigating the effects of different characteristics of the driven current, such as current density distribution and deposition location, the factors which can determine the suppression effect on the resistive tearing modes have been explored. It is found that an electron cyclotron wave (ECW) driven current with a larger peak value and more focused deposition region has a better stabilization effect. When the ECW-driven current is closer to the rational surface, it has a better stabilizing effect. These gyrokinetic toroidal code (GTC) linear simulations in the electron fluid limit of the tearing modes in the cylindrical geometry agree well with the magnetohydrodynamic codes. In conclusion, the optimal timing control of the current deposition on resistive tearing modes is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Functional Comparison of Lunar Regoliths and Their Simulants

Lunar regolith simulants are essential to the development of technology for human exploration of the Moon. Any equipment that will interact with the surface environment must be tested with simulant to mitigate risk. To reduce the greatest amount of risk, the simulant must replicate the lunar surface as well as possible. To quantify the similarities and differences between simulants, the Figures of Merit were developed. The Figures of Merit software compares the simulants and regolith by particle size, particle shape, density, and bulk chemistry and mineralogy; these four properties dictate the majority of the remaining characteristics of a geologic material. There are limitations to both the current Figures of Merit approach and simulants in general. The effect of particle textures is lacking in the Figures of Merit software, and research into this topic has only recently begun with applications to simulants. In addition, not all of the properties for lunar regolith are defined sufficiently for simulant reproduction or comparison; for example, the size distribution of particles greater than 1 centimeter and the makeup of particles less than 10 micrometers is not well known. For simulants, contamination by terrestrial weathering products or undesired trace phases in feedstock material is a major issue. Vapor deposited rims have not yet been created for simulants. Fortunately, previous limitations such as the lack of agglutinates in simulants have been addressed and commercial companies are now making agglutinate material for simulants. Despite some limitations, the Figures of Merit sufficiently quantify the comparison between simulants and regolith for useful application in lunar surface technology. Over time, the compilation and analysis of simulant user data will add an advantageous predictive capability to the Figures of Merit, accurately relating Figures of Merit characteristics to simulant user parameters.

Rickman, D.↗

Activation of the Navy's indirect effects lightning simulation laboratory

The Naval Air Test Center is currently the Navy's lead lab for electromagnetic effects testing. As part of this charter, it has been performing lightning effects testing on Navy aircraft in support of specification compliance since 1973. An overview is presented of lightning test and evaluation efforts at NATC, both past and present, as well as its plans for the future. The array of simulation capabilities presently operational are described, and a high level look is given to the test methodology now being used. The principal discussion centers on the results from the recent air launched ordnance test and testing of the Navy's A-6E all weather attack aircraft. Particular attention is paid to the NATC's test approach, including details about coaxial return construction, aircraft preparation, and the test waveforms and data acquisition systems that were used.

Whitaker, Mike↗

Automated Calibration of Parallel and Distributed Computing Simulators: A Case Study

Many parallel and distributed computing research results are obtained in simulation, using simulators that mimic real-world executions on some target system. Each such simulator is configured by picking values for parameters that define the behavior of the underlying simulation models it implements. The main concern for a simulator is accuracy: simulated behaviors should be as close as possible to those observed in the real-world target system. This requires that values for each of the simulator's parameters be carefully picked, or “calibrated,” based on ground-truth real-world executions. Examining the current state of the art shows that simulator calibration, at least in the field of parallel and distributed computing, is often undocumented (and thus perhaps often not performed) and, when documented, is described as a labor-intensive, manual process. In this work we evaluate the benefit of automating simulation calibration using simple algorithms. Specifically, we use a real-world case study from the field of High Energy Physics and compare automated calibration to calibration performed by a domain scientist. Our main finding is that automated calibration is on par with or significantly outperforms the calibration performed by the domain scientist. Furthermore, automated calibration makes it straightforward to operate desirable tradeoffs between simulation accuracy and simulation speed.

Mc donald, Jesse↗

2011.2 Revision of the Evaluated Nuclear Data Library (ENDL2011.2)

LLNL's Computational Nuclear Physics Group and Nuclear Theory and Modeling Group have col- laborated to create the 2011.2 revised release of the Evaluated Nuclear Data Library (ENDL2011.2). ENDL2011.2 is designed to support LLNL's current and future nuclear data needs and will be em- ployed in nuclear reactor, nuclear security and stockpile stewardship simulations with ASC codes. This database is currently the most complete nuclear database for Monte Carlo and deterministic transport of neutrons and charged particles. This library was assembled with strong support from the ASC PEM and Attribution programs, leveraged with support from Campaign 4 and the DOE/O ce of Science's US Nuclear Data Program. This document lists the revisions made in ENDL2011.2 compared with the data existing in the original ENDL2011.0 release and the ENDL2011.1-rc4 re- lease candidate of April 2015. These changes are made in parallel with some similar revisions for ENDL2009.2.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A three-dimensional simulation of a coronal streamer

A time-dependent MHD simulation was used to investigate the magnetic topology of an idealized streamerlike configuration in three dimensions. Starting from an initially current-free multipole field and a transonic flow, a quasi-steady, streamerlike configuration forms. The simulated streamer exhibits a closed magnetic field region bounded by a current layer and surrounded by open fields. This helmet region is surmounted by a current layer, with an enhanced density in the closed field and current sheet regions. The simulated streamer has a finite longitudinal extent, with the closed field region primarily confined to the center of the structure.

Linker, Jon A.↗

Evolution of large-scale plasma structures in comets: Kinematics and physics

Cometary and solar wind data from December 1985 through April 1986 are presented for the purpose of determining the solar wind conditions associated with comet plasma tail disconnection events (DE's). The cometary data are from The International Halley Watch Atlas of Large-Scale Phenomena (Brandt, Niedner, and Rahe, 1992). In addition, we present the kinematic analysis of 4 DE's, those of Dec. 13.5 and 31.2, 1985, and Feb. 21.7 and 28.7, 1986. The circumstances of these DE's clearly illustrate the need to analyze DE's in groups. In situ solar wind measurements from IMP-8, ICE, and PVO were used to construct the variation of solar wind speed, density, and dynamic pressure during this interval. Data from these same spacecraft plus Vega-1 were used to determine the time of 48 current sheet crossings. These data were fitted to heliospheric current sheet curves extrapolated from the corona into the heliosphere in order to determine the best-fit source surface radius for each Carrington rotation. Comparison of the solar wind conditions and 16 DE's in Halley's comet (the four DE's discussed in this paper and 12 DE's in the literature) leaves little doubt that DE's are associated primarily with crossings of the heliospheric current sheet and apparently not with any other property of the solar wind. If we assume that there is a single or primary physical mechanism and that Halley's DE's are representative, efforts at simulation should concentrate on conditions at current sheet crossings. The mechanisms consistent with this result are sunward magnetic reconnection and tailward magnetic reconnection, if tailward reconnection can be triggered by the sector boundary crossing.

Brandt, John C.↗

Building energy simulation coupled with CFD for indoor environment: A critical review and recent applications

This paper presents a comprehensive review of the open literature on motivations, methods and applications of linking stratified airflow simulation to building energy simulation (BES). First, we reviewed the motivations for coupling prediction models for building energy and indoor environment. This review classified various exchanged data in different applications as interface data and state data, and found that choosing different data sets may lead to varying performance of stability, convergence, and speed for the co-simulation. Second, our review shows that an external coupling scheme is substantially more popular in implementations of co-simulation than an internal coupling scheme. The external coupling is shown to be generally faster in computational speed, as well as easier to implement, maintain and expand than the internal coupling. Third, the external coupling can be carried out in different data synchronization schemes, including static coupling and dynamic coupling. In comparison, the static coupling that performs data exchange only once is computationally faster and more stable than the dynamic coupling. However, concerning accuracy, the dynamic coupling that requires multiple times of data exchange is more accurate than the static coupling. Furthermore, the review identified that the implementation of the external coupling can be achieved through customized interfaces, middleware, and standard interfaces. The customized interface is straightforward but may be limited to a specific coupling application. The middleware is versatile and user-friendly but usually limited in data synchronization schemes. The standard interface is versatile and promising, but may be difficult to implement. Current applications of the co-simulation are mainly energy performance evaluation and control studies. Finally, we discussed the limitations of the current research and provided an overview for future research.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Microgap breakdown with floating metal rod perturbations

Here, we report the characterization of microgap breakdown with perturbations from a metal rod floating between anode and cathode electrodes. The effects of the metal rod on the electric field distribution and the field enhancement factor are evaluated by numerical simulation and the conformal mapping method, and they indicate that the field emission regime is not reached. The breakdown voltages in the Townsend discharge regime are determined based on the voltage–current characteristics, which are obtained from two-dimensional fluid simulations. It is found that the breakdown characteristics can be significantly modulated by the floating metal rod, and the breakdown curve (breakdown voltage vs the net gap distance) is no longer U-shaped, which deviates from the conventional Paschen's law. The underlying physical mechanisms are related to the electric field enhancement, curved breakdown path, and nonuniform ion flux caused by the electric shielding effect. The results provide insights into breakdown characteristics in microscale discharges, which may promote conventional investigation of simplified clean gaps toward more complex conditions (e.g., with floating microparticles) in miniaturized plasma devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Characterization of thunderstorm induced Maxwell current densities in the middle atmosphere

Middle atmospheric transient Maxwell current densities generated by lightning induced charge perturbations are investigated via a simulation of Maxwell's equations. A time domain finite element analysis is employed for the simulations. The atmosphere is modeled as a region contained within a right circular cylinder with a height of 110 km and radius of 80 km. A composite conductivity profile based on measured data is used when charge perturbations are centered about the vertical axis at altitudes of 6 and 10 km. The simulations indicate that the temporal structure of the Maxwell current density is relatively insensitive to altitude variation within the region considered. It is also shown that the electric field and Maxwell current density are not generally aligned.

Baginski, Michael Edward↗