Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulations”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

ATS (Advanced Terrestrial Simulator) integrated hydrology and reactive transport model output in Copper Creek, Colorado

This dataset is generated using the ATS (Advanced Terrestrial Simulator) model at Copper Creek, Colorado, the largest catchment in the East River watershed. ATS is an integrated hydrology and reactive transport model to simulate the Concentration-Discharge (C-Q) relationship, and is used to quantify the geochemical export and understand how watersheds respond to climate disturbances. Key hydrogeochemical variables, including discharge, precipitation, evapotranspiration, and geochemical concentrations, were simulated and output at a daily frequency at the Copper Creek outlet. A group of hydrology and geochemical model variables over the spatial domain are output at a bi-monthly frequency. The model used 1-km resolution daily frequency Daymet meteorological forcing, and mineral composition from previous studies. The simulations are performed on NERSC. The datasets include the input files for running the ATS (Advanced Terrestrial Simulator) simulation, including the mesh file, input script and geochemical model files. The output files from the ATS (Advanced Terrestrial Simulator) are the daily time-series discharge and geochemical concentration simulated at the outlet of Copper Creek, and the bi-monthly spatial output over the simulation domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Path Integral Monte Carlo Simulations of Iron Plasmas (Final Technical Report)

This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Vehicle Powertrain Simulation Accuracy for Various Drive Cycle Frequencies and Upsampling Techniques

As connected and automated vehicle technologies emerge and proliferate, lower frequency vehicle trajectory data is becoming more widely available. In some cases, entire fleets are streaming position, speed, and telemetry at sample rates of less than 10 seconds. This presents opportunities to apply powertrain simulators such as the National Renewable Energy Laboratory's Future Automotive Systems Technology Simulator to model how advanced powertrain technologies would perform in the real world. However, connected vehicle data tends to be available at lower temporal frequencies than the 1-10 Hz trajectories that have typically been used for powertrain simulation. Higher frequency data, typically used for simulation, is costly to collect and store and therefore is often limited in density and geography. This paper explores the suitability of lower frequency, high availability, connected vehicle data for detailed powertrain simulation. A large data set of 1 Hz trajectories is used to quantify the accuracy loss when simulating energy consumption for conventional, hybrid, and battery electric powertrains using less than 1 Hz data. Techniques to upsample lower frequency drive cycle data in order to increase accuracy are also explored. Median energy consumption errors when simulating energy consumption for a 1/10 Hz trajectory are found to be 3-6% when compared to 1 Hz trajectories. Applying upsampling and interpolation techniques are shown to reduce the simulation errors by roughly 50%. The findings in this work can guide connected vehicle data collection specifications and processing techniques applied when using collected data for powertrain simulation.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

High Performance Computing Traffic Simulations for Real-Time Traffic Control of Mobility in Chattanooga Region

In 2019, highway congestion wasted over 3 billion gallons of fuel and caused 8.8 billion hours of lost productivity.1 Research has shown that introducing near-real time traffic controls can significantly reduce congestion. Validated and calibrated traffic simulations enable the modeling of transportation systems and the evaluation of different traffic control actions and schemes given a variety of circumstances that represent likely future scenarios. The developed scenarios can inform the deployment of controls in near real-time to improve freight and passenger vehicle congestion and energy use. In this work, we present simulations used to model the traffic in the Chattanooga, Tennessee, metropolitan area. Simulations were constructed and calibrated using a variety of local, data science enhanced, data sources utilizing open source software including the Simulation of Urban Mobility (SUMO) simulator. High-Performance Computing (HPC) provides a scalable platform with enough computing for the high-fidelity simulation of many scenarios and the application of advance data science especially for large-scale systems. Our simulations include microscopic simulations at a corridor level for traffic signal control, and mesoscopic simulations to evaluate regional operational controls and infrastructure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗
DOE OSTI2026

Coupling Microstructural Evolution Simulations to Material Property Degradation Predictions for Plasma-Facing Materials

Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers part I: Computational models and workflow

Computational simulations have become central to the seismic analysis and design of major infrastructure over the past several decades. Most major structures are now “proof tested” virtually through representative simulations of earthquake-induced response. More recently, with the advancement of high-performance computing (HPC) platforms and the associated massively parallel computational ecosystems, simulation is beginning to play a role in increased understanding and prediction of ground motions for earthquake hazard assessments. However, the computational requirements for regional-scale geophysics-based ground motion simulations are extreme, which has restricted the frequency resolution of direct simulations and limited the ability to perform the large number of simulations required to numerically explore the problem parametric space. In this article, recent developments toward an integrated, multidisciplinary earth science-engineering computational framework for the regional-scale simulation of both ground motions and resulting structural response are described with a particular emphasis on advancing simulations to frequencies relevant to engineered systems. This multidisciplinary computational development is being carried out as part of the US Department of Energy (DOE) Exascale Computing Project with the goal of achieving a computational framework poised to exploit emerging DOE exaflop computer platforms scheduled for the 2022–2023 timeframe.

58 GEOSCIENCES↗
DOE OSTI2025

An open-access simulated earthquake ground-motion database for an M7 Hayward Fault earthquake in the San Francisco Bay Region

Comprehensive understanding of earthquake ground motions, particularly in the near-fault region of large-magnitude events, is limited by gaps in strong-motion data. This challenge is prominent in areas with high seismic hazard but infrequent large earthquakes where data is sparse and difficult to interpret. These data limitations lead to uncertainties in the development of site-specific ground motions, which are crucial for engineering risk assessments. To address these challenges, physics-based regional-scale ground-motion simulations have been developed. With the emergence of exaflop-scale computing ecosystems, it is now possible to simulate regional earthquake processes at unprecedented fidelity and generate the large number of fault rupture realizations necessary to characterize both intra- and inter-event ground-motion variability. This article introduces a new database of simulated earthquake ground motions, created for applications in earthquake engineering, earthquake planning, and emergency response. The inaugural version of the database features simulated ground motions for a magnitude 7 Hayward Fault earthquake in the San Francisco Bay Region (SFBR), using the EarthQuake SIMulation (EQSIM) simulation framework and the Graves–Pitarka kinematic rupture model. The aim is to provide high-fidelity, spatially dense, three-component motions generated on the Department of Energy’s (DOE) newest generation of graphics processing unit (GPU)-accelerated supercomputers. These motions are being made openly available to the engineering, scientific, and disaster planning communities. In addition, this work develops protocols for the efficient dissemination of these large data sets and emphasizes community engagement to build confidence in their application. This article discusses the methodology behind the data, underlying software verification and validation, scalable data management, and a user interface for data access. The goal is to facilitate widespread use and elicit expert feedback to maximize the utility and exploitation of simulated motions. While the initial focus is on the San Francisco Region, simulations for additional regions will be added as the DOE program progresses.

Simulated ground-motion database↗
DOE OSTI2022

Diurnal cycle of precipitation over global monsoon systems in CMIP6 simulations

Accurate simulation of the diurnal cycle of precipitation over different climate regimes, including global monsoon systems, is an ongoing challenge for current state-of-the-art climate models. Here, in this study, the diurnal cycles of summertime precipitation from 21 CMIP6 historical simulations are evaluated with the NASA IMERG data via the diurnal harmonic amplitude and phase from the Fourier analysis in seven regional monsoon systems: Northern and Southern Africa monsoons, South Asia monsoon, East Asia monsoon, Australia monsoon, and North and South America monsoons. Particularly, we focus on four subregions including ocean, coastal ocean, land, and coastal land considering the different underlying physical processes. It is found that the CMIP6 historical simulations perform fairly well on reproducing the diurnal cycle of precipitation over ocean and coastal ocean, while the simulated diurnal harmonic phase and amplitude over land and coastal land show larger biases and model spread compared to biases over ocean and coastal ocean. The simulation of precipitation diurnal cycle is substantially improved in the multi-model mean of both CMIP6 historical and CMIP6 AMIP simulations compared with that of CMIP5 historical simulations. But the model bias of precipitation diurnal cycle is in general similar between CMIP6 AMIP and CMIP6 historical simulations, which suggests that the impact of the interactive ocean on the simulated diurnal precipitation averaged over a large domain is generally minor. The influence of model resolution on the diurnal cycle is shown in mixed directions across various subregions, and significant improvement from the low- to high-resolution models is only notable over coastal ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Myna: Connecting powder bed fusion build data to simulation tools for digital twin applications

Additive manufacturing (AM), as a digital process, can generate a detailed digital thread linking a part’s design and manufacturing to its operational performance. As AM systems advance, an increasing amount of process data is stored in manufacturing databases. In principle, this data can be utilized by simulation-based digital twin approaches, such as real-time process control and asynchronous post-processing guidance. However, few tools currently exist for systematically integrating digital thread data with computational tools. Here, in this study, we propose a software package, called Myna, for connecting data from powder bed fusion processes to simulation tools. The utility of such a platform is demonstrated using build data from the Oak Ridge National Laboratory Manufacturing Demonstration Facility “Peregrine v2023-10” public dataset to automatically configure and run 54 semi-analytical 3DThesis melt pool simulations, 78 numerical Additive FOAM melt pool simulations, and 3 ExaCA microstructure simulations. The simulated, spatially registered microstructures are then compared directly with electron backscatter diffraction characterization of the corresponding as-built part locations. The resulting simulated microstructure showed variation as a function of process parameters, particularly stripe width; however, the experimental data had little variation between the microstructure texture and grain size resulting from different processing conditions. Analysis of the discrepancies suggest that it is possible a two-phase ferritic-austenitic solidification model is needed to accurately predict grain size and texture for certain stainless steel 316L feedstock compositions under powder bed fusion conditions, providing direction for future research. As illustrated here, due to the number and complexity of the simulations involved in AM process-structure–property predictions, automated methods to connect process data and simulations will remain necessary tools for testing hypotheses and implementing digital twin applications.

Knapp, Gerald L. [Oak Ridge National Laboratory (O↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Validation of SOLPS-ITER and EDGE2D-EIRENE simulations for H, D, and T JET ITER-like wall low-confinement mode plasmas

Both experiments and simulations with SOLPS-ITER and EDGE2D-EIRENE show that the onset of detachment for the low-field side (LFS) divertor – defined here as the line-averaged upstream density (〈n e 〉 edge ) at which the plasma flux to the LFS target (I LFS-plate ) starts to decrease with increasing 〈n e 〉 edge – is independent of the isotope mass. However, there are three major simulation-experiment discrepancies: (i) the absolute values of I LFS-plate and the electron density (n e ) in the LFS divertor at the onset of detachment are significantly lower in simulations, i.e., approximately a factor of 2 for I LFS-plate and a factor of 3-4 for n e ; (ii) the degree of detachment – defined here as the difference between I LFS-plate at the onset of detachment and at an 〈n e 〉 edge value close to the density limit – is smaller in simulations compared to experiments; and (iii) the experimentally observed larger degree of detachment for D and T plasmas compared to H plasmas cannot be clearly distinguished from the simulation results. There are strong indications that discrepancy (i) is to a large extent caused by neglecting Lyman-opacity effects in our simulations. The simulations predict a similar net volumetric recombination source for all isotopes due to the fact that molecule-activated recombination (MAR) compensates for the reduced electron–ion recombination (EIR) for H, whereas MAR is negligible for D and T. This similar net volumetric recombination source for all isotopes leads to an isotope-independent degree of detachment in simulations. An analysis of the Balmer-α and Balmer-$γ$ emission confirms the underestimate of MAR in simulations (especially for D and T) for the JET metallic wall, which was previously observed for devices with a carbon wall. The underestimate of MAR is an important cause for discrepancy (ii) and the fact that there is a stronger underestimate of MAR for D and T than for H explains discrepancy (iii). Extending the plasma grid to the vessel wall increases I LFS-plate and n e at the onset of detachment by 25%, and the EIR source increases by 80% in detached conditions. Hence, while the extended grid results are closer to the experimental observations, the previously described qualitative discrepancies still persist.

Balmer emission↗
DOE OSTI2021

Simulated Dust Transport in the Convective Boundary Layer

Soil dusts are an important source of aerosol in agricultural regions and can affect the Earth's radiation budget through the modification of cloud properties, and in particular, through their ability to act as ice nucleating particles. In order to impact cloud properties, agricultural soil dusts need to be transported from the point of emission to cloud–relevant altitudes. Vertical transport within the planetary boundary layer is strongly controlled by turbulence and is challenging to represent accurately in regional and global models. Large–eddy simulations (LES) are run at resolutions capable of resolving most of the turbulent energy directly and can thus better simulate vertical transport. In this study, we leverage the LES ARM Symbiotic Simulation and Observation (LASSO) large–eddy simulations to simulate vertical transport of agricultural dust within the turbulent boundary layer using a modified version of the stochastic dispersion model FLEXPART–WRF. We find that the modified model is better capable of simulating particle transport due to turbulence, and that particle size was the greatest factor in determining particle lifetime. Individual meteorology and particle density had intermediate effects upon particle transport, while release height had little effect upon simulation results. Finally, we utilize a quasi–single column model (QSCM) approach to determine how our results compare to a parameterized treatment of turbulence. The QSCM simulations led to greater tracer transport out of the boundary layer, with ramifications for any studies utilizing a Lagrangian stochastic model to understand tracer dispersion. Furthermore, these results highlight the importance of accurately simulating turbulence for understanding particle transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Enhancing traffic safety analysis with digital twin technology: integrating vehicle dynamics and environmental factors into microscopic traffic simulation

Traffic safety is a critical concern in transportation engineering and urban planning. Traditional traffic safety analysis requires trained observers to collect data in the field, which is time-consuming, labor-intensive, and sometimes inaccurate. In recent years, microscopic traffic simulation, which simulates individual vehicles’ movements within a transportation network, have been utilized to study traffic safety. However, microscopic traffic simulation only focuses on traffic-related factors, such as traffic volume, traffic signals, and lane configurations, neglecting vehicle dynamics and environment-related factors like weather and lighting conditions, which can significantly impact traffic safety. In light of this, this paper explores the application of digital twin technology in traffic safety analysis, integrating vehicle simulators, which consider vehicle dynamics and environmental factors, and microscopic traffic simulators, which simulate the operations of traffic flow, for enhanced safety evaluations. Various scenarios, including different weather conditions and visibility levels, are simulated using a digital twin of a road segment in Tuscaloosa, Alabama. The simulations employ Surrogate Safety Measures (SSMs) like Time to Collision (TTC) and Deceleration Rate to Avoid a Crash (DRAC) to assess safety under varying conditions. The results demonstrate that traffic digital twin can identify potential safety issues that traditional microscopic simulation cannot, providing insights for improving traffic control strategies and transportation infrastructure to enhance traffic safety.

Civil engineering↗
DOE OSTI2020

Fully compressible simulations of waves and core convection in main-sequence stars

Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Evaluating Connected and Automated Vehicles in Co-Simulation Environment of Traffic Microsimulation and Vehicle Dynamics

Connected and automated vehicles (CAVs) have the potential to improve many aspects of the current transportation systems such as safety, mobility, and energy efficiency. In order to evaluate the benefits and impacts of a CAV, the CAV control algorithm is typically implemented on vehicles simulated in a traffic microsimulation environment. However, traffic microsimulation usually lacks detailed vehicle and powertrain dynamics, making it challenging to fully understand how a CAV control algorithm will perform and respond on an actual vehicle. Whether the same benefits measured in the simulation will also be observed in real-world remains an open question. One potential approach to fill in this gap is to conduct a co-simulation of traffic microsimulation with detailed vehicle and powertrain dynamics models, often developed in MATLAB Simulink. However, current microsimulation tools such as VISSIM and SUMO do not have a ready-to-use interface for co-simulation with vehicle dynamics and Simulink. Also, even if such an interface exists, it will be tool-specific, making it challenging to shift from one tool to another or test CAV controls in different tools. There are needs for tool-agnostic co-simulation as different microsimulation tools have their pros and cons, and researchers often need to use different tools based on the purposes of the simulation, project needs, and applications. In this work, Flexible Interface for X-in-the-loop Simulation (FIXS) is developed that can support the co-simulation of microsimulation, CAV control algorithm, and vehicle dynamics model in Simulink. Enabled by the FIXS, the benefit and performance of a CAV control algorithm can be better understood with the consideration of vehicle responses and dynamics. The connection to VISSIM and SUMO is handled internally by the interface, and users can easily switch tools by changing a configuration file. The co-simulation capability is demonstrated for a VISSIM eco-approach and departure CAV scenario and a SUMO cooperative merging scenario for both a passenger CAV and a class 8 heavy-duty connected and automated trucks.

Shao, Yunli↗
DOE OSTI2021

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗
DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗