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 415 records · Page 23

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↗

Functionalization of SBA-15 mesoporous silica by Cu-phosphonate units: Probing of synthesis route

Mesoporous silica SBA-15 containing propyl-copper phosphonate units was investigated. The structure of mesoporous samples was tested by N{sub 2} isothermal sorption (BET and BHJ analysis), TEM microscopy and X-Ray scattering. Quantitative analysis EDX has given information about proportions between component atoms in the sample. Quantitative elemental analysis has been carried out to support EDX. To examine bounding between copper atoms and phosphonic units the Raman spectroscopy was carried out. As a support of Raman scattering, the theoretical calculations were made based on density functional theory, with the B3LYP method. By comparison of the calculated vibrational spectra of the molecule with experimental results, distribution of the active units inside silica matrix has been determined. - Graphical abstract: The present study is devoted to mesoporous silica SBA-15 containing propyl-copper phosphonate units. The species were investigated to confirm of synthesis procedure correctness by the micro-Raman technique combined with DFT numerical simulations. Complementary research was carried out to test the structure of mesoporous samples. - Highlights: • SBA-15 silica functionalized with propyl-copper phosphonate units was synthesized. • Synthesis efficiency probed by Raman study supported with DFT simulations. • Homogenous distribution of active units was proved. • Synthesis route enables precise control of distance between copper ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Short-range ordering of ion-implanted nitrogen atoms in SiC-graphene

We perform a structural analysis of nitrogen-doped graphene on SiC(0001) prepared by ultra low-energy ion bombardment. Using scanning tunneling microscopy, we show that nitrogen atoms are incorporated almost exclusively as graphitic substitution in the graphene honeycomb lattice. With an irradiation energy of 25 eV and a fluence of approximately 5 × 10{sup 14 }cm{sup −2}, we achieve a nitrogen content of around 1%. By quantitatively comparing the position of the N-atoms in the topography measurements with simulated random distributions, we find statistically significant short-range correlations. Consequently, we are able to show that the dopants arrange preferably at lattice sites given by the 6 × 6-reconstruction of the underlying substrate. This selective incorporation is most likely triggered by adsorbate layers present during the ion bombardment. This study identifies low-energy ion irradiation as a promising method for controlled doping in epitaxial graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

GridLAB-D Technical Support Document: Residential Equivalent Thermal Parameter Model

GridLAB-D is a power distribution systems simulation and analysis tool developed by Pacific Northwest National Laboratory. Although GridLAB-D and its underlying residential load model has been used in a variety of power systems analyses, comprehensive documentation of the derivation of the residential house model has not been formally published. The purpose of this technical support document is to serve as a comprehensive resource for the thermal dynamics modeling and house definition implemented in GridLAB-D. It covers model derivation, parameter definitions, and implementation steps in the codebase.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SURVEY DESIGN FOR SPECTRAL ENERGY DISTRIBUTION FITTING: A FISHER MATRIX APPROACH

The spectral energy distribution (SED) of a galaxy contains information on the galaxy's physical properties, and multi-wavelength observations are needed in order to measure these properties via SED fitting. In planning these surveys, optimization of the resources is essential. The Fisher Matrix (FM) formalism can be used to quickly determine the best possible experimental setup to achieve the desired constraints on the SED-fitting parameters. However, because it relies on the assumption of a Gaussian likelihood function, it is in general less accurate than other slower techniques that reconstruct the probability distribution function (PDF) from the direct comparison between models and data. We compare the uncertainties on SED-fitting parameters predicted by the FM to the ones obtained using the more thorough PDF-fitting techniques. We use both simulated spectra and real data, and consider a large variety of target galaxies differing in redshift, mass, age, star formation history, dust content, and wavelength coverage. We find that the uncertainties reported by the two methods agree within a factor of two in the vast majority ({approx}90%) of cases. If the age determination is uncertain, the top-hat prior in age used in PDF fitting to prevent each galaxy from being older than the universe needs to be incorporated in the FM, at least approximately, before the two methods can be properly compared. We conclude that the FM is a useful tool for astronomical survey design.

79 ASTRONOMY AND ASTROPHYSICS↗

THE DEPENDENCE OF PRESTELLAR CORE MASS DISTRIBUTIONS ON THE STRUCTURE OF THE PARENTAL CLOUD

The mass distribution of prestellar cores is obtained for clouds with arbitrary internal mass distributions using a selection criterion based on the thermal and turbulent Jeans mass and applied hierarchically from small to large scales. We have checked this methodology by comparing our results for a log-normal density probability distribution function with the theoretical core mass function (CMF) derived by Hennebelle and Chabrier, namely a power law at large scales and a log-normal cutoff at low scales, but our method can be applied to any mass distributions representing a star-forming cloud. This methodology enables us to connect the parental cloud structure with the mass distribution of the cores and their spatial distribution, providing an efficient tool for investigating the physical properties of the molecular clouds that give rise to the prestellar core distributions observed. Simulated fractional Brownian motion (fBm) clouds with the Hurst exponent close to the value H = 1/3 give the best agreement with the theoretical CMF derived by Hennebelle and Chabrier and Chabrier's system initial mass function. Likewise, the spatial distribution of the cores derived from our methodology shows a surface density of companions compatible with those observed in Trapezium and Ophiucus star-forming regions. This method also allows us to analyze the properties of the mass distribution of cores for different realizations. We found that the variations in the number of cores formed in different realizations of fBm clouds (with the same Hurst exponent) are much larger than the expected root N statistical fluctuations, increasing with H.

79 ASTRONOMY AND ASTROPHYSICS↗

Traveling Wave Relays for Distribution Feeder Protection with High Penetrations of Distributed Energy Resources

Increased penetration of power electronics interfaced Distributed Energy Resources (DER) like PV, electric vehicle and battery storage in the distribution system (both in numbers and sizes) may cause bi-directional power flows under normal operation. In addition, they contribute low fault current under short circuit conditions. These two changes can impact the proper operation of legacy protection systems. This can add additional challenges during weak grid and islanded operation. In this research, traveling wave (high frequency) signature based protection scheme is proposed for future distribution systems. Simulations are performed under different transient scenarios to provide insight into the high frequency signatures to show the traveling wave behavior. Finally, this paper discusses the importance of sensing and digital processing requirements for traveling-wave protection in distribution system.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Action-based Dynamical Modeling for the Milky Way Disk: The Influence of Spiral Arms

RoadMapping is a dynamical modeling machinery developed to constrain the Milky Way’s (MW) gravitational potential by simultaneously fitting an axisymmetric parametrized potential and an action-based orbit distribution function (DF) to discrete 6D phase-space measurements of stars in the Galactic disk. In this work, we demonstrate RoadMapping's robustness in the presence of spiral arms by modeling data drawn from an N-body simulation snapshot of a disk-dominated galaxy of MW mass with strong spiral arms (but no bar), exploring survey volumes with radii 500 pc⩽r{sub max}⩽5 kpc. The potential constraints are very robust, even though we use a simple action-based DF, the quasi-isothermal DF. The best-fit RoadMapping model always recovers the correct gravitational forces where most of the stars that entered the analysis are located, even for small volumes. For data from large survey volumes, RoadMapping finds axisymmetric models that average well over the spiral arms. Unsurprisingly, the models are slightly biased by the excess of stars in the spiral arms. Gravitational potential models derived from survey volumes with at least r{sub max}=3 kpc can be reliably extrapolated to larger volumes. However, a large radial survey extent, r{sub max}∼5 kpc, is needed to correctly recover the halo scale length. In general, the recovery and extrapolability of potentials inferred from data sets that were drawn from inter-arm regions appear to be better than those of data sets drawn from spiral arms. Our analysis implies that building axisymmetric models for the Galaxy with upcoming Gaia data will lead to sensible and robust approximations of the MW’s potential.

79 ASTRONOMY AND ASTROPHYSICS↗

Traveling Wave Relays for Distribution Feeder Protection with High Penetrations of Distributed Energy Resources: Preprint

Increased penetration of power electronics interfaced Distributed Energy Resources (DER) like PV, electric vehicle and battery storage in the distribution system (both in numbers and sizes) may cause bi-directional power flows under normal operation. In addition, they contribute low fault current under short circuit conditions. These two changes can impact the proper operation of legacy protection systems. This can add additional challenges during weak grid and islanded operation. In this research, traveling wave (high frequency) signature based protection scheme is proposed for future distribution systems. Simulations are performed under different transient scenarios to provide insight into the high frequency signatures to show the traveling wave behavior. Finally, this paper discusses the importance of sensing and digital processing requirements for traveling-wave protection in distribution system.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Traveling Wave Relays for Distribution Feeder Protection with High Penetrations of Distributed Energy Resources

Increased penetration of power electronics interfaced Distributed Energy Resources (DER) like PV, electric vehicle and battery storage in the distribution system (both in numbers and sizes) may cause bi-directional power flows under normal operation. In addition, they contribute low fault current under short circuit conditions. These two changes can impact the proper operation of legacy protection systems. This can add additional challenges during weak grid and islanded operation. In this research, traveling wave (high frequency) signature based protection scheme is proposed for future distribution systems. Simulations are performed under different transient scenarios to provide insight into the high frequency signatures to show the traveling wave behavior. Finally, this paper discusses the importance of sensing and digital processing requirements for traveling-wave protection in distribution system.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Effect of Particle Size Distribution on Slurry Rheology: Nuclear Waste Simulant Slurries

Controlling the rheological properties of slurries has been of great interest in various industries such as cosmetics, ceramic processing, and nuclear waste treatment. Many physicochemical parameters, such as particle size, pH, ionic strength, and mass/volume fraction of particles, can influence the rheological properties of slurry. Among such parameters, the particle size distribution of slurry would be especially important for nuclear waste treatment because most nuclear waste slurries show a broad particle size distribution. We studied the rheological properties of several different low activity waste nuclear simulant slurries having different particle size distributions under high salt and high pH conditions. Using rheological and particle size analysis, it was found that the percentage of colloid-sized particles in slurry appears to be a key factor for rheological characteristics and the efficiency of rheological modifiers. This behavior was shown to be coupled with an existing electrostatic interaction between particles under a low salt concentration. Our study suggests that one may need to implement the particle size distribution as a critical factor to understand and control rheological properties in nuclear waste treatment plants, such as the U.S. Department of Energy’s Hanford and Savannah River sites, because the particle size distributions significantly vary over different types of nuclear waste slurries.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

ANALYZING METEOROID FLIGHTS USING PARTICLE FILTERS

Fireball observations from camera networks provide position and time information along the trajectory of a meteoroid that is transiting our atmosphere. The complete dynamical state of the meteoroid at each measured time can be estimated using Bayesian filtering techniques. A particle filter is a novel approach to modeling the uncertainty in meteoroid trajectories and incorporates errors in initial parameters, the dynamical model used, and observed position measurements. Unlike other stochastic approaches, a particle filter does not require predefined values for initial conditions or unobservable trajectory parameters. The Bunburra Rockhole fireball, observed by the Australian Desert Fireball Network (DFN) in 2007, is used to determine the effectiveness of a particle filter for use in fireball trajectory modeling. The final mass is determined to be 2.16±1.33 kg with a final velocity of 6030±216 m s{sup −1}, similar to previously calculated values. The full automatability of this approach will allow an unbiased evaluation of all events observed by the DFN and lead to a better understanding of the dynamical state and size frequency distribution of asteroid and cometary debris in the inner solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-I-28: Evaluating the Organ Dose From Computed Tomography Using Monte Carlo Calculations

Purpose: To evaluate organ doses from computed tomography (CT) using Monte Carlo (MC) calculations. Methods: A Philips Brilliance CT scanner (64 slice) was simulated using the GMctdospp (IMPS, Germany) based on the EGSnrc user code. The X-ray spectra and a bowtie filter for MC simulations were determined to coincide with measurements of half-value layer (HVL) and off-center ratio (OCR) profile in air. The MC dose was calibrated from absorbed dose measurements using a Farmer chamber and a cylindrical water phantom. The dose distribution from CT was calculated using patient CT images and organ doses were evaluated from dose volume histograms. Results: The HVLs of Al at 80, 100, and 120 kV were 6.3, 7.7, and 8.7 mm, respectively. The calculated HVLs agreed with measurements within 0.3%. The calculated and measured OCR profiles agreed within 3%. For adult head scans (CTDIvol) =51.4 mGy), mean doses for brain stem, eye, and eye lens were 23.2, 34.2, and 37.6 mGy, respectively. For pediatric head scans (CTDIvol =35.6 mGy), mean doses for brain stem, eye, and eye lens were 19.3, 24.5, and 26.8 mGy, respectively. For adult chest scans (CTDIvol=19.0 mGy), mean doses for lung, heart, and spinal cord were 21.1, 22.0, and 15.5 mGy, respectively. For adult abdominal scans (CTDIvol=14.4 mGy), the mean doses for kidney, liver, pancreas, spleen, and spinal cord were 17.4, 16.5, 16.8, 16.8, and 13.1 mGy, respectively. For pediatric abdominal scans (CTDIvol=6.76 mGy), mean doses for kidney, liver, pancreas, spleen, and spinal cord were 8.24, 8.90, 8.17, 8.31, and 6.73 mGy, respectively. In head scan, organ doses were considerably different from CTDIvol values. Conclusion: MC dose distributions calculated by using patient CT images are useful to evaluate organ doses absorbed to individual patients.

60 APPLIED LIFE SCIENCES↗

X-ray Scattering and Coarse-Grained Simulations for Clustering and Interactions of Monoclonal Antibodies at High Concentrations

Attractive protein-protein interactions (PPI) in concentrated monoclonal antibody (mAb) solutions may lead to reversible oligomers (clusters) that impact colloidal stability and viscosity. Herein, the PPI are tuned for two mAbs via the addition of arginine (Arg), NaCl, or ZnSO 4 as characterized by the structure factor ( S eff (q)) with small angle x-ray scattering (SAXS). The SAXS data are fit with molecular dynamics simulations by placing a physically relevant short-ranged attractive interaction on selected beads in coarse-grained 12 bead models of the mAb shape. The optimized 12-bead models are then used to differentiate key microstructural properties, including center of mass radial distribution functions ( g COM (r)), coordination numbers and cluster size distributions (CSD). The addition of co-solutes results in more attractive Seff(q) relative to the no co-solute control for all systems tested, with the most attractive systems showing an upturn at low q . Only the All1 model with an attractive site in each Fab and Fc region (possessing Fab-Fab, Fab-Fc and Fc-Fc interactions) can reproduce this upturn, and the corresponding CSDs show the presence of larger clusters compared to the control. In general, for models with similar net attractions, i.e., second osmotic virial coefficients, the size of the clusters increases as the attraction is concentrated on a smaller number of evenly distributed beads. Here, the cluster size distributions from simulations are used to improve the understanding and prediction of experimental viscosities. The ability to discriminate between models with bead interactions at particular Fab and Fc bead sites from SAXS simulations, and to provide real space properties (CSD and g COM (r)), will be of interest in engineering protein sequence and formulating protein solutions for weak PPI to minimize aggregation and viscosities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Decomposition of multi-particle azimuthal correlations in Q-cumulant analysis*

Abstract The method of Q -cumulants is a powerful tool for studying the fine details of azimuthal anisotropies in high energy nuclear collisions. This paper presents a new method, based on mathematical induction, to evaluate the analytical form of high-order Q -cumulants. The capability of this method is demonstrated via a toy model that uses the elliptic power distribution to simulate the anisotropic emission of particles, quantified in terms of Fourier flow harmonics . The method can help in studying the large amount of event statistics that can be collected in the future and allow measurements of the very high central moments of the distribution. This can, in turn, facilitate progress in understanding the initial geometry, the input to the hydrodynamic calculations of medium expansion in high energy nuclear collisions, and the constraints on it.

Physics↗