Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation framework”

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 397 records · Page 22

Resolved photons in S HERPA

We present the first complete simulation framework, in the SHERPA event generator, for resolved photon interactions at next-to leading order accuracy. It includes photon spectra obtained through the equivalent-photon approximation, parton distribution functions to parametrize the hadronic structure of quasi-real photons, the matching of the parton shower to next-to leading order QCD calculations for resolved photon cross sections, and the modelling of multiple-parton interactions. We validate our framework against a wide range of photo-production data from LEP and HERPA experiments, observing good overall agreement. We identify important future steps relevant for high-quality simulations at the planned Electron–Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Situational Awareness of Grid Anomalies (SAGA)

The modern power industry becomes more vulnerable to cyber events due to the growing interconnectivity, interdependence, and complexity of the electric power grid. High-fidelity modeling and simulation tools that support the preventative risk analysis on potential cyber-relevant events is essential for ensuring the situational awareness of the system operator as it provides an inexpensive and risk-free environment to test the system responses under various cyber-relevant events and hereby can support research on cyber anomaly detection, optimal protective resource allocation, and mitigation measures. In this webinar, we will share NREL's cybersecurity research capabilities by highlighting the development of a scalable cyber-physical event test bed and demonstration with real hardware in the loop. The developed cyber-physical event test bed is backboned by an integrated transmission, distribution, and communication dynamic co-simulation framework and a plug-and-play cyber event generation module. It is designed to be modular and compatible with parallel computing, and thereby supports large-scale system simulations at an affordable computation cost. The test bed can capture millisecond-to-minutes dynamic frequency and voltage responses under cyber events from the bulk transmission system to the active distribution systems and distributed energy resources at the grid edge.

co-simulation↗

Benchmarking and Exploring Parameter Space of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation

High intensity proton pulses strike the Spallation Neutron Source (SNS)’s mercury target to provide bright neutron beams. These strikes deposit extensive energy into the mercury and its steel vessel. Prediction of the resultant loading on the target is difficult when helium gas is intentionally injected into the mercury to reduce the loading and to mitigate the pitting damage on the vessel. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. We present a study comparing the measured target strains in the SNS target station with the simulation results of the solid mechanics simulation framework. We investigate a wide range of various physical model parameters, including the number of bubble families, bubble size distribution, viscosity, surface tension, etc. to understand their impact on simulation accuracy. Our initial findings reveal that using 8-10 bubble families in the model renders a simulation strain envelope that covers the experimental ones. Further optimization studies are planned to predict the strain response more accurately.

Lin, Lianshan↗

Low-diffusion Xe-He gas mixtures for rare-event detection: electroluminescence yield

High pressure xenon Time Projection Chambers (TPC) based on secondary scintillation (electroluminescence) signal amplification are being proposed for rare event detection such as directional dark matter, double electron capture and double beta decay detection. The discrimination of the rare event through the topological signature of primary ionisation trails is a major asset for this type of TPC when compared to single liquid or double-phase TPCs, limited mainly by the high electron diffusion in pure xenon. Helium admixtures with xenon can be an attractive solution to reduce the electron diffusion significantly, improving the discrimination efficiency of these optical TPCs. We have measured the electroluminescence (EL) yield of Xe-He mixtures, in the range of 0 to 30% He and demonstrated the small impact on the EL yield of the addition of helium to pure xenon. For a typical reduced electric field of 2.5 kV/cm/bar in the scintillation region, the EL yield is lowered by ~ 2%, 3%, 6% and 10% for 10%, 15%, 20% and 30% of helium concentration, respectively. This decrease is less than what has been obtained from the most recent simulation framework in the literature. The impact of the addition of helium on EL statistical fluctuations is negligible, within the experimental uncertainties. The present results are an important benchmark for the simulation tools to be applied to future optical TPCs based on Xe-He mixtures.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

User-Oriented Improvements in the MOOSE framework in support of Multiphysics Simulation

The MOOSE Framework is a foundational capability used by NEAMS to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms including optimization, improved user experience, streamlined APIs, parallelism, and new capability. The work transcribed in the report was in direct support of the simulation tools and is already deployed or will be deployed in the coming months. The capabilities implemented were, in the same order as this report, chainable execution objects or executors, support for transfers between applications at the same level in a coupling scheme, support for boundary/subdomain restricted transfers, support for transfers between applications with different coordinate or unit systems, support for MOOSE applications in the NEAMS workbench, deployment of MOOSE application of the INL HPC OnDemand platform, addition of a triangular meshing library in libMesh and increased support of face variables.

97 MATHEMATICS AND COMPUTING↗

Modeling transient edge plasma transport with dynamic recycling

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A-LEAF

A-LEAF is an electricity market simulation framework. It includes models for short-term market operations and long-term investment planning of generation, storage, and transmission. Model formulations include least cost and two-stage market equilibrium. Simulations can be flexibly structured over one or more years, including options for periodic samples or jumping to a single future year. Years can be simulated as 365 consecutive days at hourly or 5-minute sample rates, or they can be subsampled using a representative day group selection algorithm.

MANN, WILLIAM↗

Assessing the feasibility of a spaceborne 3D lightning observing concept

The distribution of electrical charge in thunderclouds results from thermodynamic, microphysical, and kinematic processes, which also modulate thunderstorm evolution. It is no surprise that the connection between lightning and these physical processes is so strong that the increase and vertical growth of lightning activity closely follows the vertical growth of the thundercloud, but unraveling these connections is not trivial and requires observations of the three-dimensional (3D) structure of electrical activity in a cloud. Ground-based 3D lightning mapping networks give excellent 3D flash-level detail but are limited to regional coverage. Satellite-based optical lightning mappers give excellent global coverage but are largely limited to 2D summaries of flash rate and radiant intensity, albeit new flash products and stereographic techniques are chipping away this limitation. New observing strategies are needed to expand and diversify the corpus of 3D lightning datasets and motivate studies that unravel connections lightning has with these key physical processes and the surrounding environment. This study examines the feasibility of using a distributed network of orbing satellites with VHF-based lightning detectors to obtain global maps of 3D lightning activity and assess efficacy of this approach for use in a new, small satellite mission concept called CubeSpark. CubeSpark combines new VHF and high-resolution, bispectral optical instruments on a constellation of low-Earth orbiting (LEO) satellites to globally map the 3D electrical structure of thunderstorms and study how it relates to thunderstorm evolution, extreme weather, nitrogen oxide production and distribution, upper atmospheric electrical phenomena, and how 3D flash observations can complement existing satellite-based lightning mappers and improve decision support tools. To locate lightning discharges, CubeSpark seeks to use the VHF time-of-arrival technique, similar to ground-based total lightning mapping networks. The vertical location accuracy of these satellite retrievals will be of poorer quality compared to a ground-based network, which has non-trivial implications for lightning flash reconstruction and lightning-based interpretations of deep convection. We adapt a Lightning Mapping Array (LMA) simulation framework to an orbiting network and use it to address feasibility of 3D lightning detection from space with particular attention to location accuracy of VHF detections of lightning in the vertical. These simulations inform a constellation design study that defines a realistic orbital configuration and depicts the global coverage for CubeSpark. Results indicate that a 3D location accuracy of <1-2 km for each dimension can be achieved across 300-500 km wide swaths, which suggests that CubeSpark can resolve the charge structure of thunderclouds from the tropics to the mid- and high- latitudes.

Lightning↗

Simulation of Heavy-Duty Vehicles in Platooning Scenarios

Connectivity and automation provide the potential to use information about the environment and future driving to minimize energy consumption. In order to achieve this goal, the designers of control strategies need to simulate a wide range of driving situations that can interact with other vehicles and the infrastructure to account for the specific powertrain characteristics of each automated vehicle. We present here a simulation framework called RoadRunner, which aims to facilitate the design of powertrain-aware eco-driving algorithms and a better understanding of the interactions between automation and powertrain technology. RoadRunner allows users to simulate both powertrain and longitudinal dynamics within a simulated environment. The user defines the scenario to be simulated by providing a route profile, intersection control types, number of vehicles, vehicle class, powertrain technology, connectivity, and automation level. RoadRunner then builds a Simulink diagram of the scenario, including the information flows between vehicles, and between vehicles and the road. After the simulation, the user can analyze the driving, component operations, and detailed energy consumption rates for each simulated vehicle. We present a case study on heavy-duty vehicles platooning. Using RoadRunner and detailed road data from digital maps, we quantify the impact of gap setting on fuel consumption, for a real-world route.

Kim, Namdoo↗

Experimental validation for thermal transmittances of window shading systems with perimeter gaps

Virtually all residential and commercial windows in the U.S. have some form of window attachment, but few have been designed for energy savings. ISO 15099 presents a simulation framework to determine thermal performance of window attachments, but the model has not been validated for these products. This paper outlines a review and validation of the ISO 15099 centre-of-glass heat transfer correlations for perimeter gaps (top, bottom, and side) in naturally ventilated cavities through measurement and simulation. The thermal transmittance impact due to dimensional variations of these gaps is measured experimentally, simulated using computational fluid dynamics, and simulated utilizing simplified correlations from ISO 15099. Results show that the ISO 15099 correlations produce a mean error between measured and simulated heat flux of 2.5 ± 7%. These tolerances are similar to those obtained from sealed cavity comparisons and are deemed acceptable within the ISO 15099 framework.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Linking constituent phase properties to ductility and edge stretchability of two DP 980 steels

Two DP 980 steels were characterized with multiple methods to determine its microstructure characteristics, tensile properties and edge stretchability. The first (DP980-B) is a commercial grade DP980 steel produced by AK Steel, while the second (DP980-T) is a modification of the first one by a tempering process. Both contain two individual phases: ferrite and martensite with body centered cubic (BCC) and tetragonal (BCT) crystal structure respectively. ASTM standard tensile tests were performed for both materials, revealing that DP980-B has higher work hardening and larger total elongation than DP980-T. Standard hole piercing and subsequent hole expansion and extrusion tests were also performed. The critical hole expansion ratio (HER) of hole pierced sheet at different nominal clearances, on the other hand, show an opposite trend: the HER of DP980-T almost doubles that of DP980-B. This indicates the tensile ductility and edge stretchability are determined by different mechanisms due to different deformation modes. The ductility of sheet metal during uniaxial tension is determined by deformation instability, i.e. necking, due to less constraint at the region of deformation and fracture will occur quickly after necking. The contribution of post-necking deformation to the total elongation is insignificant for the high strength grades such as DP980 steels. The ductility during hole-expansion, however, is mainly determined by material fracture behavior that is dependent on phase property disparity and material intrinsic fracture parameters. With the individual phase properties for both DP steels determined by in situ tensile tests under high energy X-ray diffraction, the results show that DP980-T steel has lower property disparities between the two phases than DP980-B steel, due to the tempering process. This explains why the DP980-T steel has higher HER than DP980-B steel: lower phase disparity will lead to less local deformation during loading. The higher work hardening rate for DP980-B steel contributes to its higher uniform elongation compared with DP980-T steel based on the maximum load condition of deformation instability. An integrated finite element simulation framework for studying hole expansion is also presented here, based on the calculated individual phase properties from the combined high energy x-ray diffraction (HEXRD) and elastic plastic self-consistent modeling. The simulation results correlate well with the experimental results on the HER difference between the two materials.

36 MATERIALS SCIENCE↗

Circular Economy Modeling Efforts at NREL: Circular Economy Lifecycle Analysis and Visualization (CELAVI) and the Circular Economy Agent-Based Model (ABM)

Circular Economy (CE) aims at decoupling human activities from economic growth and resource use. While CE economic and environmental benefits are often touted by its proponents, increased circularity and improved sustainability do not necessarily go hand-in-hand. The Circular Economy Lifecycle Assessment and Visualization (CELAVI) framework simulates the transition of renewable energy technology supply chains towards circularity and quantifies spatially explicit environmental impacts and supply chain costs. Using wind blade end-of-life (EOL) management as a case study, this presentation will highlight CELAVI's features. Furthermore, envisioned extension to the framework could improve the stakeholder decision model by accounting for logistical and other factors. The CELAVI framework could be leveraged to answer crucial questions regarding renewables EOL management.

circular economy↗

Influence of Pore Water on the Fracture of Silica Sand at Particle Scale

The fracture of particles has a significant influence on the engineering behavior of granular materials. There are no reported conclusive answers in the literature to explain the influence of pore water on the fracture of sand, gravel, railroad ballast rock, aggregates, rockfill, and so on. Here, in this paper, two novel techniques were adopted to investigate the influence of pore water on the fracture of natural silica sand at the particle scale. The three-dimensional (3D) synchrotron microcomputed tomography (SMT) imaging technique was used to acquire and analyze multiple 3D images of specimens composed of wetted silica sand that were subjected to confined one-dimensional (1D) compression loading up to the fracture stage. A few representative particles were identified and digitally separated from the liquid water and gas phases. The 3D SMT images offer clear experimental evidence of the phase transition of water from liquid to gas within the opening cracks. In addition, a computational fluid dynamics numerical simulation framework for water flow into an opening crack was adopted and the results show that the fast opening of the crack induced cavitation, water phase transition, and water hammer effects. The results of the numerical simulations agree with the confined 1D compression experiments where the onset of the bubbles (gas phase) within the crack is mainly generated during the secondary cavitation stage, and the generated bubble near the crack mouth opening resulted from the primary cavitation, which was directly caused by the crack expansion. The results reported in this paper present a new cavitation phenomenon that occurs during particle fracture, which offers a physics-based explanation of why water-wetted sand fractures at smaller stresses than dry and can pave the way for more in-depth future studies on the fracture of water-wetted sand.

cavitation↗

Image-Driven Hybrid Structural Analysis Based on Continuum Point Cloud Method with Boundary Capturing Technique

Conventional approaches for the structural health monitoring of infrastructures often rely on physical sensors or targets attached to structural members, which require considerable preparation, maintenance, and operational effort, including continuous on-site adjustments. This paper presents an image-driven hybrid structural analysis technique that combines digital image processing (DIP) and regression analysis with a continuum point cloud method (CPCM) built on a particle-based strong formulation. Polynomial regressions capture the boundary shape change due to the structural loading and precisely identify the edge and corner coordinates of the deformed structure. The captured edge profiles are transformed into essential boundary conditions. This allows the construction of a strongly formulated boundary value problem (BVP), classified as the Dirichlet problem. Capturing boundary conditions from the digital image is novel, although a similar approach was applied to the point cloud data. It was shown that the CPCM is more efficient in this hybrid simulation framework than the weak-form-based numerical schemes. Unlike the finite element method (FEM), it can avoid aligning boundary nodes with regression points. A three-point bending test of a rubber beam was simulated to validate the developed technique. The simulation results were benchmarked against numerical results by ANSYS and various relevant numerical schemes. The technique can effectively solve the Dirichlet-type BVP, yielding accurate deformation, stress, and strain values across the entire problem domain when employing a linear strain model and increasing the number of CPCM nodes. In addition, comparative analysis with conventional displacement tracking techniques verifies the developed technique’s robustness. The proposed technique effectively circumvents the inherent limitations of traditional monitoring methods resulting from the reliance on physical gauges or target markers so that a robust and non-contact solution for remote structural health monitoring in real-scale infrastructures can be provided, even in unfavorable experimental environments.

Chemistry↗

Surface roughness scattering of electrons in bulk mosfets

Surface-roughness scattering of electrons at the Si-SiO 2 interface is a very important consideration when analyzing Si metal-oxide-semiconductor field-effect transistors (MOSFETs). Scattering reduces the mobility of the electrons and degrades the device performance. 250-nm and 50-nm bulk MOSFETs were simulated with varying device parameters and mesh sizes in order to compare the effects of surface-roughness scattering in multiple devices. The simulation framework includes the ensemble Monte Carlo method used to solve the Boltzmann transport equation coupled with a successive over-relaxation method used to solve the two-dimensional Poisson's equation. Four methods for simulating the surface-roughness scattering of electrons were implemented on both devices and compared: the constant specularity parameter, the momentum-dependent specularity parameter, and the real-space-roughness method with both uniform and varying electric fields. The specularity parameter is the probability of an electron scattering speculariy from a rough surface. It can be chosen as a constant, characterizing partially diffuse scattering of all electrons from the surface the same way, or it can be momentum dependent, where the size of rms roughness and the normal component of the electron wave number determine the probability of electron-momentum randomization. The real-space rough surface method uses the rms roughness height and correlation length of an actual MOSFET to simulate a rough interface. Due to their charge, electrons scatter from the electric field and not directly from the surface. If the electric field is kept uniform, the electrons do not perceive the roughness and scatter as if from a at surface. However, if the field is allowed to vary, the electrons scatter from the varying electric field as they would in a MOSFET. These methods were implemented for both the 50-nm and 250-nm MOSFETs, and using the rms roughness heights and correlation lengths for real devices. The current-voltage and mobility-electric field curves were plotted for each method on the two devices and compared. The conclusion is that the specularity-parameter methods are valuable as simple models for relatively smooth interfaces. However, they have limitations, as they cannot accurately describe the drastic reduction in the current and the electron mobility that occur in MOSFETs with very rough Si-SiO 2 interfaces.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electron energization in reconnection: Eulerian vs Lagrangian perspectives

Particle energization due to magnetic reconnection is an important unsolved problem for myriad space and astrophysical plasmas. Electron energization in magnetic reconnection has traditionally been examined from a particle, or Lagrangian, perspective using particle-in-cell (PIC) simulations. Guiding-center analyses of ensembles of PIC particles have suggested that Fermi (curvature drift) acceleration and direct acceleration via the reconnection electric field are the primary electron energization mechanisms. However, both PIC guiding-center ensemble analyses and spacecraft observations are performed in an Eulerian perspective. For this work, we employ the continuum Vlasov–Maxwell solver within the Gkeyll simulation framework to reexamine electron energization from a kinetic continuum, Eulerian, perspective. We separately examine the contribution of each drift energization component to determine the dominant electron energization mechanisms in a moderate guide-field Gkeyll reconnection simulation. In the Eulerian perspective, we find that the diamagnetic and agyrotropic drifts are the primary electron energization mechanisms away from the reconnection x-point, where direct acceleration dominates. We compare the Eulerian (Vlasov Gkeyll) results with the wisdom gained from Lagrangian (PIC) analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advancing coupled fire-atmosphere research with HIGRAD/FIRETEC (w20_firetec) (CY 2021 Institutional Computing Annual Progress Report) [Slides]

LANL’s HIGRAD and FIRETEC models were used in a multi-scale simulation framework to investigate impacts of small-scale fuel heterogeneities on large-scale Pyro-Cb formation. Forest structural diversity is found to increase canopy fire spread, and thus risk of tree mortality in FIRETEC simulations. Structural diversity is associated with an increase in ladder fuels, which facilitate fire spread into the canopy. Low-intensity fire burning in marginal conditions is highly sensitive to small perturbation in atmospheric turbulence. These results highlight predictive limitations of individual fire simulations, and emphasize the need for ensemble simulations, especially for low-intensity, prescribed fire.

36 MATERIALS SCIENCE↗

Java-based Graphical User Interface for MAVERIC-II

A computer program entitled "Marshall Aerospace Vehicle Representation in C II, (MAVERIC-II)" is a vehicle flight simulation program written primarily in the C programming language. It is written by James W. McCarter at NASA/Marshall Space Flight Center. The goal of the MAVERIC-II development effort is to provide a simulation tool that facilitates the rapid development of high-fidelity flight simulations for launch, orbital, and reentry vehicles of any user-defined configuration for all phases of flight. MAVERIC-II has been found invaluable in performing flight simulations for various Space Transportation Systems. The flexibility provided by MAVERIC-II has allowed several different launch vehicles, including the Saturn V, a Space Launch Initiative Two-Stage-to-Orbit concept and a Shuttle-derived launch vehicle, to be simulated during ascent and portions of on-orbit flight in an extremely efficient manner. It was found that MAVERIC-II provided the high fidelity vehicle and flight environment models as well as the program modularity to allow efficient integration, modification and testing of advanced guidance and control algorithms. In addition to serving as an analysis tool for techno logy development, many researchers have found MAVERIC-II to be an efficient, powerful analysis tool that evaluates guidance, navigation, and control designs, vehicle robustness, and requirements. MAVERIC-II is currently designed to execute in a UNIX environment. The input to the program is composed of three segments: 1) the vehicle models such as propulsion, aerodynamics, and guidance, navigation, and control 2) the environment models such as atmosphere and gravity, and 3) a simulation framework which is responsible for executing the vehicle and environment models and propagating the vehicle s states forward in time and handling user input/output. MAVERIC users prepare data files for the above models and run the simulation program. They can see the output on screen and/or store in files and examine the output data later. Users can also view the output stored in output files by calling a plotting program such as gnuplot. A typical scenario of the use of MAVERIC consists of three-steps; editing existing input data files, running MAVERIC, and plotting output results.

Seo, Suk Jai↗