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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Application of chord length distributions and principal component analysis for quantification and representation of diverse polycrystalline microstructures

Highlights: • Microstructure quantification based on directional chord length distribution and principal component analysis is presented • A boundary-based algorithm of calculations of chord length distributions is introduced • Utility of the quantification framework is shown in three case studies on synthetic, experimental, and simulation datasets • Angularly resolved chord length distribution captures essential details of grain size and shape in diverse microstructures • Low-dimensional representations of directional chord length distributions are obtained using principal component analysis - Abstract: Quantification of mesoscale microstructures of polycrystalline materials is important for a range of practical tasks of materials design and development. The current protocols of quantifying grain size and morphology often rely on microstructure metrics (e.g., mean grain diameter) that overlook important details of the mesostructure. In this work, we present a quantification framework based on directionally resolved chord length distribution and principal component analysis as a means of extracting additional information from 2-D microstructural maps. Towards this end, we first present in detail a method for calculating chord length distribution based on boundary segments available in modern digital datasets (e.g., from microscopy post-processing) and their low-rank representations by principal component analysis. The utility of the proposed framework for capturing grain size, morphology, and their anisotropy for efficient visualization, representation, and specification of polycrystalline microstructures is then demonstrated in case studies on datasets from synthetic generation, experiments (on Ni-base superalloys), and simulations (on steel during recrystallization).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Data-driven predictive probability density function control of fiber length stochastic distribution shaping in refining process

Pulp is the most important raw material in paper in paper industries, whose Fiber Length Stochastic Distribution (FLSD) shaping directly determines the energy consumption and paper quality of the subsequent papermaking processes. However, the mean and variance are insufficient to describe the output FLSD shaping, which display non-Gaussian distribution properties. Therefore, the traditional control, method based on the mean and variance of fiber length is difficult to control the FLSD shaping effectively. In this paper, a novel data-driven predictive probability density function (PDF) control method is proposed for output FLSD shaping in refining process. Primarily, in order to improve the approximation accuracy of the PDF employing RBF neural network, the parameters (center value and width) of each RBF basis functions are tuned via utilizing iterative learning control (ILC) law, and the corresponding estimations of weights law can be obtained. Secondly, considering that the conventional linear model of weights vector has drawbacks of low accuracy and weak generalization ability, random vector functional link networks (RVFLNs) based data driven nonlinear modelling method is employed to characterize the prediction model between the input variables and weights vector. Finally, in order to reduce the randomness of the output FLSD, a minimum entropy control method under mean constraint is employed to design the predictive PDF controller for the output FLSD shaping with the help of the established nonlinear stochastic distribution model. Both simulations and experiments are given to illustrate the effectiveness and practicability of the proposed method.

Probability debsity function, Fiber length distrib↗
DOE OSTI2016

Three-dimensional stochastic model of actin–myosin binding in the sarcomere lattice

The effect of molecule tethering in three-dimensional (3-D) space on bimolecular binding kinetics is rarely addressed and only occasionally incorporated into models of cell motility. The simplest system that can quantitatively determine this effect is the 3-D sarcomere lattice of the striated muscle, where tethered myosin in thick filaments can only bind to a relatively small number of available sites on the actin filament, positioned within a limited range of thermal movement of the myosin head. Here we implement spatially explicit actomyosin interactions into the multiscale Monte Carlo platform MUSICO, specifically defining how geometrical constraints on tethered myosins can modulate state transition rates in the actomyosin cycle. The simulations provide the distribution of myosin bound to sites on actin, ensure conservation of the number of interacting myosins and actin monomers, and most importantly, the departure in behavior of tethered myosin molecules from unconstrained myosin interactions with actin. In addition, MUSICO determines the number of cross-bridges in each actomyosin cycle state, the force and number of attached cross-bridges per myosin filament, the range of cross-bridge forces and accounts for energy consumption. At the macroscopic scale, MUSICO simulations show large differences in predicted force-velocity curves and in the response during early force recovery phase after a step change in length comparing to the two simplest mass action kinetic models. The origin of these differences is rooted in the different fluxes of myosin binding and corresponding instantaneous cross-bridge distributions and quantitatively reflects a major flaw of the mathematical description in all mass action kinetic models. Consequently, this new approach shows that accurate recapitulation of experimental data requires significantly different binding rates, number of actomyosin states, and cross-bridge elasticity than typically used in mass action kinetic models to correctly describe the biochemical reactions of tethered molecules and their interaction energetics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Local and average structure of Mn- and La-substituted BiFeO{sub 3}

The local and average structure of solid solutions of the multiferroic perovskite BiFeO{sub 3} is investigated by synchrotron X-ray diffraction (XRD) and electron density functional theory (DFT) calculations. The average experimental structure is determined by Rietveld refinement and the local structure by total scattering data analyzed in real space with the pair distribution function (PDF) method. With equal concentrations of La on the Bi site or Mn on the Fe site, La causes larger structural distortions than Mn. Structural models based on DFT relaxed geometry give an improved fit to experimental PDFs compared to models constrained by the space group symmetry. Berry phase calculations predict a higher ferroelectric polarization than the experimental literature values, reflecting that structural disorder is not captured in either average structure space group models or DFT calculations with artificial long range order imposed by periodic boundary conditions. Only by including point defects in a supercell, here Bi vacancies, can DFT calculations reproduce the literature results on the structure and ferroelectric polarization of Mn-substituted BiFeO{sub 3}. The combination of local and average structure sensitive experimental methods with DFT calculations is useful for illuminating the structure-property-composition relationships in complex functional oxides with local structural distortions. - Graphical abstract: The experimental and simulated partial pair distribution functions (PDF) for BiFeO{sub 3}, BiFe{sub 0.875}Mn{sub 0.125}O{sub 3}, BiFe{sub 0.75}Mn{sub 0.25}O{sub 3} and Bi{sub 0.9}La{sub 0.1}FeO{sub 3}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Toward Verifying Fossil Fuel CO2 Emissions with the CMAQ Model: Motivation, Model Description and Initial Simulation

Motivated by the urgent need for emission verification of CO2 and other greenhouse gases, we have developed regional CO2 simulation with CMAQ over the contiguous U.S. Model sensitivity experiments have been performed using three different sets of inputs for net ecosystem exchange (NEE) and two fossil fuel emission inventories, to understand the roles of fossil fuel emissions, atmosphere-biosphere exchange and transport in regulating the spatial and diurnal variability of CO2 near the surface, and to characterize the well-known ‘signal-to-noise’ problem, i.e. the interference from the biosphere on the interpretation of atmospheric CO2 observations. It is found that differences in the meteorological conditions for different urban areas strongly contribute to the contrast in concentrations. The uncertainty of NEE, as measured by the difference among the three different NEE inputs, has notable impact on regional distribution of CO2 simulated by CMAQ. Larger NEE uncertainty and impact are found over eastern U.S. urban areas than along the western coast. A comparison with tower CO2 measurements at Boulder Atmospheric Observatory (BAO) shows that the CMAQ model using hourly varied and high-resolution CO2 emission from the Vulcan inventory and CarbonTracker optimized NEE reasonably reproduce the observed diurnal profile, whereas switching to different NEE inputs significantly degrades the model performance. Spatial distribution of CO2 is found to correlate with NOx, SO2 and CO, due to their similarity in emission sources and transport processes. These initial results from CMAQ demonstrate the power of a state-of-the art CTM in helping interpret CO2 observations and verify fossil fuel emissions. The ability to simulate CO2 in CMAQ will also facilitate investigations of the utility of traditionally regulated pollutants and other species as tracers to CO2 source attribution.

CMAQ↗
DOE OSTI2020

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Cross-correlation measurements with the EJ-299-33 plastic scintillator

New organic–plastic scintillation compositions have demonstrated pulse-shape discrimination (PSD) of neutrons and gamma rays. We present cross-correlation measurements of 252 Cf and mixed uranium–plutonium oxide (MOX) with the EJ-299-33 plastic scintillator. For comparison, equivalent measurements were performed with an EJ-309 liquid scintillator. Offline, digital PSD was applied to each detector. These measurements show that EJ-299-33 sacrifices a factor of 5 in neutron–neutron efficiency relative to EJ-309, but could still utilize the difference in neutron–neutron efficiency and neutron single-to-double ratio to distinguish 252 Cf from MOX. These measurements were modeled with MCNPX-PoliMi, and MPPost was used to convert the detailed collision history into simulated cross-correlation distributions. MCNPX-PoliMi predicted the measured 252 Cf cross-correlation distribution for EJ-309 to within 10%. Greater photon uncertainty in the MOX sample led to larger discrepancy in the simulated MOX cross-correlation distribution. In conclusion, the modeled EJ-299-33 plastic also gives reasonable agreement with measured cross-correlation distributions, although the MCNPX-PoliMi model appears to under-predict the neutron detection efficiency.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-T-196: Heat Diffusion Modeling for Digital Holographic Interferometry Dosimetry

Purpose: We have previously demonstrated that with Digital Holographic Interferometry (DHI) 2D spatial calorimetric measurements of high dose rate radiation sources can be obtained. The impact of heat transfer must be considered when undertaking any form of calorimetric measurement, as the radiation induced temperature distributions are subject to degradation due to heat diffusion. Unaccounted for, this limits the accuracy of the approach especially for long delivery times. Methods: 3D modelling of the heat diffusion in water was undertaken, and two different approaches developed to account for this effect. The mathematical framework to describe heat diffusion in 3D was applied, with the differential equations solved numerically using an implicit method. The first approach involved the comparison of the DHI measurements to an independent dose model of the source. The model was forward modeled to account for the heat diffusion during irradiation, allowing a direct comparison to validate the measured results. The second approach involved the correction of the measured data directly, by comparing the temperature distribution of two instances and subtracting the effects of heat diffusion of the first distribution from the second instance. This required the use of the Abel transform to approximate the 3D dose distribution from the 2D DHI results, thus limiting the approach to radiation applications possessing cylindrical symmetry. Results: The first approach resulted in higher accuracy and was more straightforward, but has a major limitation in that the measured results are only able to be utilized in comparison with an independent dose model. The applicability of the second approach is affected by noise in the measurement data and introduces higher uncertainties, but results in higher usability of the final data. Conclusion: Both approaches were implemented, and if used in conjunction would provide the most utility for the interpretation and use of DHI measurements.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2018

Single crystal to polycrystal neutron transmission simulation

A collection of routines for calculation of the total cross section that determines the attenuation of neutrons by crystalline solids is presented. The total cross section is calculated semi-empirically as a function of crystal structure, neutron energy, temperature, and crystal orientation. The semi-empirical formula includes the contribution of parasitic Bragg scattering to the total cross section using both the crystal’s mosaic spread value and its orientation with respect to the neutron beam direction as parameters. These routines allow users to enter a distribution of crystal orientations for calculation of total cross sections of user defined powder or pseudo powder distributions, which enables simulation of non-uniformities such as texture and strain. In conclusion, the spectra for neutron transmission simulations in the neutron thermal energy range (2 meV–100 meV) are presented for single crystal and polycrystal samples and compared to measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ASPIRES: Airport Shuttle Planning and Improved Routing Event-driven Simulation

Most of the existing traffic simulation packages require significant calibration work to be able to reflect reality. To evaluate special operations including emerging technologies, a microscopic simulation that tracks detailed interactions of all the elements of the traffic systems is usually needed. This type of simulation is usually computationally demanding. This work developed an Airport Shuttle Planning and Improved Routing Event-driven Simulation (ASPIRES) package to simulate and evaluate current, potential, and future airport shuttle operations. Here, the simulation was driven by data and thus did not require much calibration effort. The discrete-event simulation nature of ASPIRES makes the simulation computationally efficient. Simulating 1 day of shuttle operations takes less than 2s. The study site of this work is the Dallas/Fort Worth International Airport in the U.S. The shuttle service that connects the five terminals of the airport and the rental car center was studied. Travel times, dwell times, and passenger arrivals were simulated using empirical distributions derived mainly from real data to capture the stochastic nature of the rental car center shuttle bus operations. Data on bus miles traveled, bus energy consumption, passenger wait times, and passengers left behind at stops were collected to study the trade-off between energy use and passenger experience. Electric bus and on-demand bus operations were also included. The simulation outputs can show passengers statistics at terminals, shuttles statistics, and charging station statistics. ASPIRES cannot be used to model a generic traffic system but is well-suited for fleet systems.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Impact of Open Communication Networks on Load Frequency Control with Plug-In Electric Vehicles by Cyber-Physical Dynamic Co-Simulation

With the increasing electrification of the transportation sector to achieve the carbon neutrality objective, despite the challenges of charging electric vehicles (EV), there are also opportunities through smart charging EVs to improve system frequency stability; however, EV control technologies might require nontraditional communication support. This paper investigates the impacts of communication variations of EV on power system load frequency control through a cyber-physical dynamic system (CPDS) co-simulation. Here, the CPDS is built upon our previously developed transmission-and-distribution dynamic co-simulation model with the added communication variation functions (i.e., delay and packet loss). The case studies consider multiple communication variation scenarios when the system experiences an N-1 generation trip contingency. The scenarios include communication delays and packet loss using both homogeneous and heterogeneous assumptions. The outcomes of this work can help improve EV frequency regulation services and provide robust and effective tests for different load frequency control algorithms of the future power systems.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗
DOE OSTI2011

PROBABILITY DISTRIBUTION FUNCTIONS OF COSMOLOGICAL LENSING: CONVERGENCE, SHEAR, AND MAGNIFICATION

We perform high-resolution ray-tracing simulations to investigate probability distribution functions (PDFs) of lensing convergence, shear, and magnification on distant sources up to the redshift of z{sub s} = 20. We pay particular attention to the shot noise effect in N-body simulations by explicitly showing how it affects the variance of the convergence. We show that the convergence and magnification PDFs are closely related to each other via the approximate relation {mu} = (1 - {kappa}){sup -2}, which can reproduce the behavior of PDFs surprisingly well up to the high magnification tail. The mean convergence measured in the source plane is found to be systematically negative, rather than zero as often assumed, and is correlated with the convergence variance. We provide simple analytical formulae for the PDFs, which reproduce simulated PDFs reasonably well for a wide range of redshifts and smoothing sizes. As explicit applications of our ray-tracing simulations, we examine the strong-lensing probability and the magnification effects on the luminosity functions of distant galaxies and quasars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Impact of Open Communication Networks on Load Frequency Control with Plug-in Electric Vehicles by Cyber-Physical Dynamic Co-Simulation: Preprint

With the increasing electrification of the transportation sector, vehicle to grid (V2G) control technologies are promising to improving system frequency stability. However, these control technologies might require non-traditional communication supports. This paper investigates the impacts of electric vehicles (EVs) V2G on power system frequency regulation considering communication variations through transmission-distribution-communication (TDC) dynamic simulation. The simulation model is built based on our previously developed open-source transmission-and-distribution (T&D) dynamic co-simulation framework. Here, we add the communication variation functions (i.e. delay and packet loss) to complete the TDC simulation framework. In the case studies, practical communication variation scenarios are considered when the system experiencing an N-1 generation trip contingency. The considered scenarios include communication delay and packet loss using both homogeneous and heterogeneous assumptions. The delay and packet loss rate are the same in all communication channels in the homogeneous case but varies in the heterogeneous case. The results show the communication delay has an obvious impact on the frequency recovery compared to the communication packet drop. The impacts of the standard deviations of communication delay time and packet drop rate are not significant. The outcomes of this work can help improve the EV frequency regulation services and provide robust and effect supports to the grid.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗
DOE OSTI2020

Identifying factors controlling the selective ethane dehydrogenation on Pt-based catalysts from DFT based micro-kinetic modeling

Herein, a detailed EDH reaction mechanism on Pt(111), Pt(211), Pt 3 Sn(111) and Pt 3 Sn(211) surfaces has been systematically computed with DFT. A mean-field micro-kinetic model has been established to simulate the product distribution and selectivity on each surface at experimentally relevant conditions. We have uncovered a number of different factors controlling the catalytic performance of Pt-based EDH catalysts, thus leading to a better understanding of experimentally detected phenomenon. In this study, the (111) and (211) model surfaces of Pt and Pt 3 Sn catalysts were chosen to represent the terrace and step-edge sites when simulating the EDH reaction mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Distributed Generation Market Demand Model (dGen): Documentation

The Distributed Generation Market Demand model (dGen) is a geospatially rich, bottom-up, market-penetration model that simulates the potential adoption of distributed energy resources (DERs) for residential, commercial, and industrial entities in the continental United States through 2050. The National Renewable Energy Laboratory (NREL) developed dGen to analyze the key factors that will affect future market demand for distributed solar, wind, storage, and other DER technologies in the United States. The new model builds off, extends, and replaces NREL's SolarDS model (Denholm et al. 2009a), which simulates the market penetration of distributed PV only. Unlike the SolarDS model, dGen can model various DER technologies under one platform--it currently can simulate the adoption of distributed solar (the dSolar module) and distributed wind (the dWind module) and link with the ReEDS capacity expansion model (Appendix C). The underlying algorithms and datasets in dGen, which improve the representation of customer decision making as well as the spatial resolution of analyses (Figure ES-1), also are improvements over SolarDS.

17 WIND ENERGY↗
DOE OSTI2015

STRUCTURE IN THE 3D GALAXY DISTRIBUTION. II. VOIDS AND WATERSHEDS OF LOCAL MAXIMA AND MINIMA

The major uncertainties in studies of the multi-scale structure of the universe arise not from observational errors but from the variety of legitimate definitions and detection methods for individual structures. To facilitate the study of these methodological dependencies, we have carried out 12 different analyses defining structures in various ways. This has been done in a purely geometrical way by utilizing the HOP algorithm as a unique parameter-free method of assigning groups of galaxies to local density maxima or minima. From three density estimation techniques (smoothing kernels, Bayesian blocks, and self-organizing maps) applied to three data sets (the Sloan Digital Sky Survey Data Release 7, the Millennium simulation, and randomly distributed points) we tabulate information that can be used to construct catalogs of structures connected to local density maxima and minima. We also introduce a void finder that utilizes a method to assemble Delaunay tetrahedra into connected structures and characterizes regions empty of galaxies in the source catalog.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Function reconsideration of indium bump in InSb IRFPAs

Indium bumps are generally accepted to possess functions of electrical connection, mechanical support and heat transfer in flip-chip devices. After comparing the distribution of stress components along different paths in InSb infrared focal plane arrays (IRFPAs), we ascertain that local enhancement effects of the stress components are remarkable in regions where the indium bumps are. More specifically, the local enhanced tensile stress in InSb chip connected with the indium bumps can lead to the local fracture of the InSb chip, and the locally enhanced shear and peeling stresses may give rise to the local interfacial delamination between the InSb chip and the indium bumps. These inferences are confirmed by the observed local failure characteristics, such as the distribution of the local delamination, the distribution of the crack widths, and the distribution of the cracks. In addition, the simulated Z-component strain distribution in InSb IRFPAs is also consistent with the backside surface profile of the InAs/GaSb IRFPAs fabricated in America with the identical structure, that is, the InSb chip glued with the indium bumps is concave downward, and the InSb chip glued with the underfill is convex upward. Judging from all these confirmed simulation results, we are confident that the indium bumps play a pivotal role in inducing the local failure of InSb IRFPAs. So the role of the indium bumps in causing the local failure of InSb IRFPAs needs to be supplemented and emphasized to comprehensively evaluate the structural reliability of InSb IRFPAs.

42 ENGINEERING↗
DOE OSTI2018

Wave Resource Characterization Using an Unstructured Grid Modeling Approach

This paper presents a modeling study conducted on the central Oregon coast for wave resource characterization using the unstructured-grid SWAN model coupled with a nested-grid WWIII model. The flexibility of models of various spatial resolutions and the effects of open- boundary conditions simulated by a nested-grid WWIII model with different physics packages were evaluated. The model results demonstrate the advantage of the unstructured-grid modeling approach for flexible model resolution and good model skills in simulating the six wave resource parameters recommended by the International Electrotechnical Commission in comparison to the observed data in Year 2009 at National Data Buoy Center Buoy 46050. Notably, spectral analysis indicates that the ST4 physics package improves upon the model skill of the ST2 physics package for predicting wave power density for large waves, which is important for wave resource assessment, device load calculation, and risk management. In addition, bivariate distributions show the simulated sea state of maximum occurrence with the ST4 physics package matched the observed data better than that with the ST2 physics package. This study demonstrated that the unstructured-grid wave modeling approach, driven by the nested-grid regional WWIII outputs with the ST4 physics package, can efficiently provide accurate wave hindcasts to support wave resource characterization. Our study also suggests that wind effects need to be considered if the dimension of the model domain is greater than approximately 100 km, or O (10^2 km).

Wu, Wei-Cheng↗