Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed simulation”

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 667 records · Page 37

Multidimensional electron beam-plasma instabilities in the relativistic regime

The interest in relativistic beam-plasma instabilities has been greatly rejuvenated over the past two decades by novel concepts in laboratory and space plasmas. Recent advances in this long-standing field are here reviewed from both theoretical and numerical points of view. The primary focus is on the two-dimensional spectrum of unstable electromagnetic waves growing within relativistic, unmagnetized, and uniform electron beam-plasma systems. Although the goal is to provide a unified picture of all instability classes at play, emphasis is put on the potentially dominant waves propagating obliquely to the beam direction, which have received little attention over the years. First, the basic derivation of the general dielectric function of a kinetic relativistic plasma is recalled. Next, an overview of two-dimensional unstable spectra associated with various beam-plasma distribution functions is given. Both cold-fluid and kinetic linear theory results are reported, the latter being based on waterbag and Maxwell-Juettner model distributions. The main properties of the competing modes (developing parallel, transverse, and oblique to the beam) are given, and their respective region of dominance in the system parameter space is explained. Later sections address particle-in-cell numerical simulations and the nonlinear evolution of multidimensional beam-plasma systems. The elementary structures generated by the various instability classes are first discussed in the case of reduced-geometry systems. Validation of linear theory is then illustrated in detail for large-scale systems, as is the multistaged character of the nonlinear phase. Finally, a collection of closely related beam-plasma problems involving additional physical effects is presented, and worthwhile directions of future research are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrated analysis of X-ray diffraction patterns and pair distribution functions for machine-learned phase identification

Abstract To bolster the accuracy of existing methods for automated phase identification from X-ray diffraction (XRD) patterns, we introduce a machine learning approach that uses a dual representation whereby XRD patterns are augmented with simulated pair distribution functions (PDFs). A convolutional neural network is trained directly on XRD patterns calculated using physics-informed data augmentation, which accounts for experimental artifacts such as lattice strain and crystallographic texture. A second network is trained on PDFs generated via Fourier transform of the augmented XRD patterns. At inference, these networks classify unknown samples by aggregating their predictions in a confidence-weighted sum. We show that such an integrated approach to phase identification provides enhanced accuracy by leveraging the benefits of each model’s input representation. Whereas networks trained on XRD patterns provide a reciprocal space representation and can effectively distinguish large diffraction peaks in multi-phase samples, networks trained on PDFs provide a real space representation and perform better when peaks with low intensity become important. These findings underscore the importance of using diverse input representations for machine learning models in materials science and point to new avenues for automating multi-modal characterization.

36 MATERIALS SCIENCE↗

The Fate of Saharan Dust Across the Atlantic: An Integrated Modeling and Observational Study of the TC4 Field Campaign

During the NASA TC-4 field campaign in July 2007, several Saharan dust events were observed over the Caribbean basin. A-Train observations suggest that these Saharan dust events are confined the Caribbean and rarely transported across Central America to the Pacific Ocean. We investigate the nature of this barrier to dust transport using the NASA GEOS-5 atmospheric general circulation model. Our simulations with GEOS-5 are driven by the Modern Era Retrospective-Analysis for Research and Applications (MERRA) meteorological analyses, and include online simulation of aerosol distributions using a version of the Goddard Chemistry, Aerosol, Radiation, and Transport (GOCART) model. Simulated dust distributions are evaluated using A-Train observations from MODIS and CALIOP, as well as MISR and ground-based AERONET sun photometers, and show good agreement with the observations in terms of the timing and magnitude of dust events. A component analysis of the dust transport and removal pathways is used to evaluate the relative roles of these processes in establishing the observed dust transport barrier. From this analysis, we show that while both atmospheric dynamics and wet removal contribute towards the Caribbean dust barrier, northward dust transport is the more dominant term. Additional simulations are performed to ascertain the sensitivity of our results to uncertain loss processes (i.e., wet removal) in our model.

Nowottnick, E.↗

HTS Transformer Windings Design Using Distributive Ratios for Minimization of Short Circuit Forces

High-temperature superconducting (HTS) transformers have a promising feature in reduction of total weight, total size, and the losses of large-scale distribution transformers. However, the lower leakage reactance of HTS transformers results in a higher short-circuit fault currents and electromagnetic forces. Therefore, optimization of short-circuit electromagnetic forces is one of the crucial aspects in the design of HTS transformers. In this paper, a novel analytical method is proposed for determination of optimum distributive ratios resulting in minimization of these forces for asymmetrical multi-segment windings of an HTS transformer. Employing these distributive ratios, radial and axial components of short-circuit electromagnetic forces in an HTS transformer are significantly reduced. Two- and three-dimensional (2D and 3D) finite element method (FEM) simulations are employed to verify the analytical method results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Studies for the Loss of Atomic and Molecular Species from Io

The general objective of this project has been to advance our theoretical understanding of Io's atmosphere and how various atomic and molecular species are lost from this atmosphere and are distributed in the circumplanetary environment of Jupiter. The scientific objectives of the larger collaborative program between AER, Inc., and the University of Michigan have been to undertake theoretical modeling studies to simulate the distributions of the exospheric gases in Io's corona and extended clouds, to investigate the importance of the various physical processes that shape their relative abundances, and with these tools to analyze observations of O, S and Na obtained by four observers: M.A. McGrath of the Space Telescope Science Institute and G.E. Ballester of the University of Michigan who each have obtained Hubble Space Telescope observations of O and S near Io, F. Scherb who continues an effort to obtain 6300 A OI observations as part of the University of Wisconsin Fabry-Perot program, and N.M. Schneider of the University of Colorado who obtained an extensive set of spectral and spatial observations of the Na emission near Io in the D-lines.

Combi, Michael R.↗

Transition region and chromospheric signatures of impulsive heating events. II. Modeling

Results from the Solar Maximum Mission showed a close connection between the hard X-ray (HXR) and transition region (TR) emission in solar flares. Analogously, the modern combination of RHESSI and IRIS data can inform the details of heating processes in ways that were never before possible. We study a small event that was observed with RHESSI, IRIS, SDO, and Hinode, allowing us to strongly constrain the heating and hydrodynamical properties of the flare, with detailed observations presented in a previous paper. Long duration redshifts of TR lines observed in this event, as well as many other events, are fundamentally incompatible with chromospheric condensation on a single loop. We combine RHESSI and IRIS data to measure the energy partition among the many magnetic strands that comprise the flare. Using that observationally determined energy partition, we show that a proper multithreaded model can reproduce these redshifts in magnitude, duration, and line intensity, while simultaneously being well constrained by the observed density, temperature, and emission measure. We comment on the implications for both RHESSI and IRIS observations of flares in general, namely that: (1) a single loop model is inconsistent with long duration redshifts, among other observables; (2) the average time between energization of strands is less than 10 s, which implies that for a HXR burst lasting 10 minutes, there were at least 60 strands within a single IRIS pixel located on the flare ribbon; (3) the majority of these strands were explosively heated with an energy distribution well described by a power law of slope ≈−1.6; (4) the multi-stranded model reproduces the observed line profiles, peak temperatures, differential emission measure distributions, and densities.

79 ASTRONOMY AND ASTROPHYSICS↗

Extension of moment projection method to the fragmentation process

The method of moments is a simple but efficient method of solving the population balance equation which describes particle dynamics. Recently, the moment projection method (MPM) was proposed and validated for particle inception, coagulation, growth and, more importantly, shrinkage; here the method is extended to include the fragmentation process. The performance of MPM is tested for 13 different test cases for different fragmentation kernels, fragment distribution functions and initial conditions. Comparisons are made with the quadrature method of moments (QMOM), hybrid method of moments (HMOM) and a high-precision stochastic solution calculated using the established direct simulation algorithm (DSA) and advantages of MPM are drawn.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Emittance growth in the DARHT Axis-II Downstream Transport

Using a particle-in-cell (PIC) code, we investigated the possibilities for emittance growth through the quadrupole magnets of the system used to transport the high-current electron beam from an induction accelerator to the bremsstrahlung converter target used for flash radiography. We found that even highly mismatched beams exhibited little emittance growth (< 6%), which we attribute to softening of their initial hard edge current distributions. We also used this PIC code to evaluate the accuracy of emittance measurements using a solenoid focal scan following the quadrupole magnets. If the beam is round after the solenoids, the simulations indicate that the measurement is highly accurate, but it is substantially inaccurate for elliptical beams

43 PARTICLE ACCELERATORS↗

Cost and Efficiency Analysis for Hybrid AC/DC Distribution System Planning with PV and Battery

Distribution system planning is one of the main tasks for the electric utilities. Cost and energy efficiency are two main concerns in the distribution planning process. In this work, cost and efficiency analyses are conducted for the hybrid AC and DC distribution grids for the planning purposes. The IEEE 34-bus feeder has been applied for the case studies. The investment cost, operation and maintenance cost, replacement cost, etc., have been considered in the life cycle cost analysis. Power losses of both the AC and DC section are calculated with an open source simulation tool, the Open Distribution Simulator Software (OpenDSS). The proposed framework can be applied to quantitatively evaluate the cost and efficiency of hybrid AC and DC distribution grids and help comparing different planning options.

AC/DC hybrid system planning↗

Lens Model Accuracy in the Expected LSST Lensed AGN Sample

Strong gravitational lensing of active galactic nuclei (AGN) enables measurements of cosmological parameters through time-delay cosmography (TDC). With data from the upcoming LSST survey, we anticipate using a sample of O(1000) lensed AGN for TDC. To prepare for this dataset and enable this measurement, we construct and analyze a realistic mock sample of 1300 systems drawn from the OM10 (Oguri & Marshall 2010) catalog of simulated lenses with AGN sources at $z<3.1$ in order to test a key aspect of the analysis pipeline, that of the lens modeling. We realize the lenses as power law elliptical mass distributions and simulate 5-year LSST i-band coadd images. From every image, we infer the lens mass model parameters using neural posterior estimation (NPE). Focusing on the key model parameters, $θ_E$ (the Einstein Radius) and $γ_{lens}$ (the projected mass density profile slope), with consistent mass-light ellipticity correlations in test and training data, we recover $θ_E$ with less than 1% bias per lens, 6.5% precision per lens and $γ_{lens}$ with less than 3% bias per lens, 8% precision per lens. We find that lens light subtraction prior to modeling is only useful when applied to data sampled from the training prior. If emulated deconvolution is applied to the data prior to modeling, precision improves across all parameters by a factor of 2. Finally, we combine the inferred lens mass models using Bayesian Hierarchical Inference to recover the global properties of the lens sample with less than 1% bias.

Venkatraman, Padmavathi [Illinois U., Urbana; KIPA↗

A 3D simulation of the early winter distribution of reactive chlorine in the north polar vortex

Early in December 1991, high values of ClO are seen by the Microwave Limb Sounder (MLS) on the Upper Atmosphere Research Satellite at latitudes south of areas of temperatures cold enough to form polar stratospheric clouds (PSCs). A 3D simulation shows that the heterogeneous conversion of chlorine reservoirs to reactive chlorine on the surfaces of PSCs (processing) takes place at high latitudes. Often the processed air must be transported to lower latitudes, where the reactive chlorine is photochemically converted to ClO, to be observed by MLS. In this simulation, one incidence of cold temperatures is associated with an anticyclone, and a second with a cyclone. The transport of processed air associated with the anticyclone is marked by shearing; a decrease in the maximum of the processed air is accompanied by growth of the area influenced by the processing. In contrast, the air processed in the cyclonic event spreads more slowly. This shows that transport and shearing is a crucial element to the evolution of reactive chlorine associated with a processing event. In particular, transport and shearing, as well as photochemical processes, can cause variations in observed ClO.

Douglass, A.↗

Coordinated Optimization of Multiservice Dispatch for Energy Storage Systems with Degradation Model for Utility Applications

The rapid development of energy storage technologies is motivating researchers to understand the possible value of energy storage. This paper evaluates the multiservice dispatch of distribution-sited energy storage systems with realistic degradation models by integrating a multiservice optimization model with a distribution system simulator. This work attempts to define the theoretical value energy storage can bring when multiple value streams are considered in the modeling environment. Use cases presented in this paper deliver multiple services to both distribution and transmission stakeholders that can provide positive net value to the electricity system under prevailing energy storage cost structures. Furthermore, this work was performed in collaboration with a distribution utility to help utility personnel understand the value distributed energy storage can provide by examining different utility use cases without considering implementation costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

EVIDENCE FOR NONLINEAR DEVELOPMENT OF MAGNETOHYDRODYNAMIC SCALE INTERMITTENCY IN THE INNER HELIOSPHERE

The formation of coherent structures in turbulence is a signature of a developing cascade and therefore might be observable by analyzing inner heliospheric solar wind turbulence. To test this idea, data from the Helios 2 mission, for six streams of solar wind at different heliocentric distances and of different velocities, were subjected to statistical analysis using the partial variance of increments (PVI) approach. We see a clear increase of the PVI distribution function versus solar wind age for higher PVI cutoff, indicating development of non-Gaussian coherent structures. The plausibility of this interpretation is confirmed by a similar behavior observed in two-dimensional magnetohydrodynamics simulation data at corresponding dimensionless nonlinear times.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic Field Mapping of a 2.5 T Fixed-Field HTS Gantry Magnet for Proton Therapy

We present results from testing a high-temperature superconducting (HTS) magnet prototype for proton therapy. This magnet is specifically designed for a novel rotating gantry capable of delivering the entire proton beam energy range (70225 MeV) while maintaining a fixed magnetic field in the superconducting magnets. The gantry layout simplifies the magnet design, enabling the use of straight, flat racetrack Bi-2223 (DI-BSCCO) coil technology and operation at higher temperatures (1015 K). The magnet has a non-linear field distribution for bending and focusing the proton beams. To validate this feature, we developed a system for measuring the magnetic field distribution in the magnet aperture. We present the design of this hall probe array and experimental results from two different magnet tests at 4.2 K in a liquid helium bath. These results are compared with the simulated field distribution and discussed in the context of the required field quality for the application.

Mosat, M↗

Nanoflare activity in the solar chromosphere

We use ground-based images of high spatial and temporal resolution to search for evidence of nanoflare activity in the solar chromosphere. Through close examination of more than 1 × 10{sup 9} pixels in the immediate vicinity of an active region, we show that the distributions of observed intensity fluctuations have subtle asymmetries. A negative excess in the intensity fluctuations indicates that more pixels have fainter-than-average intensities compared with those that appear brighter than average. By employing Monte Carlo simulations, we reveal how the negative excess can be explained by a series of impulsive events, coupled with exponential decays, that are fractionally below the current resolving limits of low-noise equipment on high-resolution ground-based observatories. Importantly, our Monte Carlo simulations provide clear evidence that the intensity asymmetries cannot be explained by photon-counting statistics alone. A comparison to the coronal work of Terzo et al. suggests that nanoflare activity in the chromosphere is more readily occurring, with an impulsive event occurring every ∼360 s in a 10,000 km{sup 2} area of the chromosphere, some 50 times more events than a comparably sized region of the corona. As a result, nanoflare activity in the chromosphere is likely to play an important role in providing heat energy to this layer of the solar atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical studies of the high power microwave breakdown in gas using the fluid model with a modified electron energy distribution function

A modified electron energy distribution function (EEDF) is introduced into the fluid model and its effects on the high power microwave (HPM) breakdown in air and argon are investigated. A proper numerical scheme for the finite-difference time-domain method is employed to solve the fluid model. Numerical simulations show that the HPM breakdown time in argon predicted by the fluid model with the modified EEDF agrees well with the results of Particle-in-cell-Monte Carlo collision simulations, while the Maxwellian EEDF results in faster HPM breakdown when the mean electron energy is less than 20 eV. We also confirm that the Maxwellian EEDF can be used in the fluid model for simulating the air breakdown at the low frequencies based on the reported experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development and Validation of a Two-Phase Thermal-Hydraulic CFD Code NEK-2P

A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable, high-performance Computational Fluid Dynamics (CFD) code NEK5000. The goal of this work is to verify and validate the two-phase version of the NEK5000 code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework (EBF) previously developed at Argonne but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed. The development of the NEK-2P two-phase solver and the implementation of the Extended Boiling Framework two-phase models were initially supported by Argonne National Laboratory (Argonne) through a Laboratory Directed Research and Development (LDRD) project during FY2014-2016. The development and validation of the two-phase models through analyses of selected two-phase boiling flow experiments was supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY2017-2020. This report focuses on verification and validation of the water-steam boiling model NEK-2P Two-Phase, CFD code. The NEK-2P was validated with Nuclear Power Engineering Corporation (NUPEC) Pressurized Water Reactor (PWR) Sub-channel and Bundle Test (PSBT) void distribution benchmark. Three different simulations were performed and analyzed for various operating conditions such as wall-heat flux and sub-cooled inlet temperatures. Reasonably good agreement with measured data was obtained in predicting the measured void distributions. Simulations were performed for Virginia Tech. (VT) 3x3 rod bundle geometry with and without spacers. The preliminary results were presented for Simplified Spacer Grid (SSG). In addition, the implementation of interface reconstruction model was tested with one of the Becker benchmark Critical Heat Flux (CHF) experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Simulations of emission from microcavity tandem organic light-emitting diodes

Microcavity tandem organic light-emitting diodes (OLEDs) are simulated and compared to experimental results. The simulations are based on two complementary techniques: rigorous finite element solutions of Maxwell's equations and Fourier space scattering matrix solutions. A narrowing and blue shift of the emission spectrum relative to the noncavity single unit OLED is obtained both theoretically and experimentally. In the simulations, a distribution of emitting sources is placed near the interface of the electron transport layer tris(8-hydroxyquinoline) Al (Alq{sub 3}) and the hole transport layer (N,N'-bis(naphthalen-1-yl)-N,N'-bis(phenyl)benzidine) ({alpha}-NPB). Far-field electric field intensities are simulated. The simulated widths of the emission peaks also agree with the experimental results. The simulations of the 2-unit tandem OLEDs shifted the emission to shorter wavelength, in agreement with experimental measurements. The emission spectra's dependence on individual layer thicknesses also agreed well with measurements. Approaches to simulate and improve the light emission intensity from these OLEDs, in particular for white OLEDs, are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗