Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 109 records · Page 6

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗

1D model of tokamak scrape-off layer and divertor plasmas: Dimensionless parameters and the results of numerical simulations

Solutions of a 1D model of the scrape-off layer and a divertor plasma in a tokamak are analyzed on the basis of dimensionless parameters governing the physics of both plasma and neutral transport and atomic physics processes. Numerical simulations based on fluid plasma and Monte Carlo neutral modeling confirm earlier results obtained with simplified analytic considerations and numerical simulations performed with fluid neutral description and demonstrate bifurcation of plasma parameters. Further, these simulations also confirm that the major physical reasons for the rollover of the plasma flux to the target, which is the key feature of the phenomenon of divertor plasma detachment, are the impurity radiation loss and plasma recombination. It is suggested that simulations of the edge plasma based on dimensionless parameters provide a powerful tool for checking the consistency of the physics incorporated in existing complex 2D codes used for edge plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology

Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An implementation of neural simulation-based inference for parameter estimation in ATLAS

Neural simulation-based inference (NSBI) is a powerful class of machine-learning-based methods for statistical inference that naturally handles high-dimensional parameter estimation without the need to bin data into low-dimensional summary histograms. Such methods are promising for a range of measurements, including at the Large Hadron Collider, where no single observable may be optimal to scan over the entire theoretical phase space under consideration, or where binning data into histograms could result in a loss of sensitivity. This work develops a NSBI framework for statistical inference, using neural networks to estimate probability density ratios, which enables the application to a full-scale analysis. It incorporates a large number of systematic uncertainties, quantifies the uncertainty due to the finite number of events in training samples, develops a method to construct confidence intervals, and demonstrates a series of intermediate diagnostic checks that can be performed to validate the robustness of the method. As an example, the power and feasibility of the method are assessed on simulated data for a simplified version of an off-shell Higgs boson couplings measurement in the four-lepton final states. This approach represents an extension to the standard statistical methodology used by the experiments at the Large Hadron Collider, and can benefit many physics analyses.

frequentist statistics↗

Transient cycling of nitrogen, organic carbon and oxygen within the free-flowing Columbia River corridor: Linking exposure time dependent biogeochemical reactions to river stage fluctuations (Final project report)

The objective of this project was to develop next-generation techniques for representing the transformations of complex reaction networks like those impacting transient river-corridors, then apply them to develop upscaling tools. The project contributed significantly to the development of Lagrangian “mass transfer particle tracking” (MTPT) tools that are the first numerical framework that explicitly separates mixing and spreading processes; a crucial distinction for accurate representations of reaction rates across scales. The MTPT approach uses an operator splitting scheme where physical transport processes are simulated using classical random walk methods and each particle is treated as a “container” that carries with it any number of chemical species. Mixing between particles is simulated using a colocation-probability based mass transfer kernel and reactions are evaluated on each particle after mixing. The method is stable, accurate, and also allows the explicit modeling of residence time distributions of the individual containers of mass. The MTPT scheme allowed us to explore several research questions related to reactive transport in river corridors and complex reaction networks. We found that it is possible to quantify the exposure time distributions (ETDs) of reactants and that this does suggest some pathways to upscaling. However, comparable approximations can be obtained using simplified (reduced dimensionality) MTPT simulations in less time with similar uncertainty, so development of extensive ETD-based methods was not productive. The major findings are that i) the MTPT schemes are robust and highly accurate across a range of arbitrarily complex reactions, ii) the methods can be efficiently parallelized and the parallel performance characteristics of MTPT are predictable, and iii) the MTPT tools allow simultaneous tracking of residence time. The development of these research tools into efficient software packages continues. The methods are already available to the community because open-source, working examples have been included with all publications.

54 ENVIRONMENTAL SCIENCES↗

An Online Tool for Preliminary Design and Techno-Economic Analysis of District Geothermal Heating and Cooling Systems

District geothermal heating and cooling systems (DGHCS) have significant benefits for reducing energy consumption as well as building- and grid-level peak electric demand. Currently, no publicly available tools are available to effectively design and conduct techno-economic analysis of DGHCS. GeoWISE was originally developed for preliminary design and techno-economic analysis of geothermal heating and cooling systems in an individual commercial or residential building. This paper introduces recent upgrades of GeoWISE that allow users to design and conduct techno-economic analysis of DGHCS. Several new features are implemented in GeoWISE to allow selection and specification of multiple new or existing buildings. A database of information for over 125 million existing U.S. buildings was used in GeoWISE that allows users easily locate existing buildings of interest based on street addresses, and optionally edit information of the buildings (e.g., footprint, vintage, principal functions, number of floors, window-to-wall ratio). Unique energy simulation models of the selected buildings are then automatically created using the Automatic Building Energy Modeling (AutoBEM) and EnergyPlus simulations are performed to predict thermal loads of the buildings. A simplified DGHCS is then designed and simulated to predict its energy use. A central borehole heat exchanger (BHE) of the DGHCS is sized using the RowWise algorithm of GHEDesigner to meet the thermal loads within user-specified land areas for installing BHE. The upgraded GeoWISE reports the needed capacity of heating and cooling equipment in each building, design of the central BHE, energy consumption reduction, and energy cost saving resulting from using DGHCS compared with conventional HVAC systems. A case study is showcased using the upgraded GeoWISE to design and conduct techno-economic analysis of a simplified DGHCS.

Prem Anand Jayaprabha, Jyothis Anand [ORNL] (ORCID↗

On a simplified solution of climate-carbon dynamics in idealized flat10MIP simulations

Abstract. Idealized experiments with coupled climate-carbon Earth system models (ESMs) provide a basis for understanding the response of the carbon cycle to external forcing and for quantifying climate-carbon feedbacks. Here, we analyze globally-averaged results from idealized esm-flat10 experiments and show that most models exhibit a quasi-linear relationship between cumulative carbon uptake on land and in the ocean during a period of constant fossil fuel emissions of 10 Pg C yr−1. We hypothesize that this relationship does not depend on emission pathways. Further, as a simplification, we quantify the relationship between cumulative ocean carbon uptake and changes in ocean heat content using a linear approximation. In this way, changes in oceanic heat content and atmospheric CO2 concentration become interdependent variables, reducing the coupled temperature-CO2 system to just one differential equation. The equation can be solved analytically or numerically for the atmospheric CO2 concentration as a function of fossil fuel emissions. This approach leads to a simplified description of global carbon and climate dynamics, which could be used for applications beyond existing analytical frameworks.

Brovkin, Victor↗

Computer-Based Procedures Among Older Adults: Findings from the Rancor Microworld Simulator

As with many industries worldwide, nuclear energy is experiencing an aging workforce; older than other energy sources and the national average. For example, while almost one in three nuclear industry professionals are 55+ years, for oil and gas that number is one in five. The challenges brought about by an aging workforce are one of the industry’s top concerns. However, the industry lacks for empirical research that examines the effects of developmental aging in older workers, as well as the interactions between aging and new digital technologies. We present the results of an experiment that tested three different computer-based procedures in a sample of 30 older adults (55+ years). While the industry has traditionally relied on paper-based procedures, these are being modernized by digital technology. Participants were randomly assigned to one of three procedure-types that varied by the level of digitalization, based on the IEEE Standard-1786. Type 1 essentially represents a digital representation of a paper-based procedure, Type 2 adds embedded indicators, and Type 3 adds soft controls. Participants performed two different operational scenarios (startup and loss of feedwater) on a simplified nuclear power plant simulator. Results revealed a weak signal that Type 2 may produce lower workload and lower completion times in some instances. However, there were no significant effects of procedure-type across any other outcomes, including simulator log data, situation awareness, and preference ratings. We discuss our findings in terms of optimal levels of digitalization/automation for an aging nuclear workforce and suggest pathways for future directions.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Volatile Products of the Autoxidation of Poly(ethylenimine) in CO 2 Sorbents

Poly(ethyleneimine) (PEI)-based CO 2 sorbents are a promising material class for use in CO 2 capture applications, particularly direct air capture (DAC), due to their high amine content and CO 2 capacities. The sorbent lifetime is a key uncertainty in the deployment of these materials in such applications, as they oxidize under conditions relevant to candidate DAC process cycles. Here, in this study, we utilize thermogravimetric analysis/differential scanning calorimetry/FTIR spectroscopy to characterize the nature, rate, and quantity of volatile species formed from the oxidation of PEI when supported on mesoporous Al 2 O 3 . We show that NH 3 , CO 2 , and H 2 O are primary volatile species formed from PEI oxidation and that their production rate matches the overall oxidation rate as measured via heat flow from a differential scanning calorimeter. We show that the total quantity of NH 3 evolved is consistent with terminal primary amines being cleaved from the polymer chain and that approximately one hydrogen atom per PEI repeat unit reacts to form H 2 O during the reaction. Finally, metadynamics and quantum chemical simulations performed on a simplified system are used to highlight the likely importance of radicals in the elementary reactions that comprise the oxidation reaction set. The findings here represent the first characterization of the volatile products from the oxidation of PEI in CO 2 sorbents and lead to a greater understanding of the mechanisms behind the oxidative degradation of these materials.

36 MATERIALS SCIENCE↗

Bay Area Regional Energy: Network Integrated Commercial Retrofits (BRICR) Project. Final Report

The BRICR project applied large-scale building energy modeling concepts with the aim of reducing the cost of energy efficiency targeting, design, and project development, and measurement of energy savings for energy efficiency programs implemented by local governments that serve small and medium commercial buildings (SMB). The project leveraged the services and resources of existing local government energy programs serving disadvantaged and hard-to-reach SMB customers. In contrast to programs run by utilities, local government programs generally do not have direct access to energy billing records for an entire class of customers in a geographic area, which prior research demonstrated useful for large-scale building energy model baseline development and calibration. , However, local governments are rich in public records that offer important clues about physical attributes and uses that, along with behavior, determine energy use. Relying only on public records, BRICR demonstrated development of credible baseline energy models for 3,792 office, retail, and hotel buildings. Publicly disclosed annual energy use data from a local energy benchmarking program and anonymized data from the Building Performance Database, the nation’s largest dataset about energy-related characteristics of buildings, were utilized to validate and calibrate energy models via an innovative method comparing distributions of energy intensity by fuel type for portfolios of buildings of similar size, vintage, and use. Portfolio calibration does not provide certainty that an energy model fits an individual building; the method is useful when billing data is not accessible – a common situation for researchers, energy service providers and ESCOs, local governments, and any party other than a utility. A software component was developed, the BRICR gem, which automates simulation when relevant data is added or edited by the user to a file saved in the standardized BuildingSync XML schema for energy audit data. The component was demonstrated as a simplified means to generate a mass of energy models corresponding to public records containing basic attributes such as building scale, location, use, year built, and aspect ratio in combination with building energy code prototype data corresponding to use and vintage. The component was also demonstrated as a simplified means to automate energy simulation when attributes are revised; the intention was to enable iterative improvement of the baseline model and energy savings estimates for common energy conservation measures as users revise relevant attributes based on their observations. In the context of institutional change and uncertainty for the participating local government energy programs, 13 whole building retrofits were completed. Impacts were measured by applying the CalTRACK2.0 methods to standardize measurement of normalized metered energy consumption. The GRIDMeter methods of stratified sampling and individual load shape analysis were applied to adjust for impacts of the effect of COVID-19 on retrofitted buildings in the context of all local buildings of similar size and use. Excluding impacts of the pandemic, retrofitted buildings demonstrated between 1.6% and 25.1% reduction in energy use. The project contributed use cases and feedback that helped inform evolution of the software tools and data formats that were combined for the first time in the BRICR project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Super-X and conventional divertor configurations in MAST-U ohmic L-mode; a comparison facilitated by interpretative modelling

Measurements are presented, alongside corresponding interpretative SOLPS-ITER simulations, of the first MAST-U experiments comparing ohmically heated L-mode fuelling scans in Conventional divertor (CD) and Super-X divertor (SXD) configurations. In experiment, at comparable outer mid-plane separatrix electron density, $n_{e,\textrm{sep,OMP}}$, the maximum lower outer target heat load was found to be a factor 16 $\,\pm\,7$ lower in SXD compared to CD. In simulation, a factor 26.8 reduction was found (slightly higher than the experimental range), suggesting an additional reduction in SXD compared to the factor 9.3 expected from geometric considerations alone. According to the simulations, this additional reduction in the SXD is due to a net radial transport of the energy remaining downstream of the $T_e = 5$ eV location. This energy is carried out of the critical (highest heat load) flux tube by deuterium atoms, demonstrating the importance of a longer legged divertor which provides space for this to occur. Importantly, in both simulation and experiment, the SXD has minimal impact on the upstream n e and T e profiles. Spectral inferences of detachment front movement in SXD compare well between simulation and experiment. In regions of high magnetic field gradient, the parallel movement of the front towards the X-point becomes less sensitive to increasing $n_{e,\textrm{sep,OMP}}$, in qualitative agreement with simplified models and previous predictive simulations. Additional aspects, regarding the target ion flux rollover, upstream separatrix temperature and drift effects, are also presented and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Discrete element simulation of Pebble Bed Reactors on graphics processing units

Prediction of pebble positions in a Pebble Bed Reactor (PBR) is necessary for both reactor physics and thermal hydraulics simulations as the arrangement of pebbles has a significant impact on the resulting core power, coolant flow, and fuel temperature. Knowledge of pebble movement as the fuel is cycled through the core is also critical for predicting the fuel residence time and subsequently, the fuel burnup. Simulation with the Discrete Element Method (DEM) can provide knowledge of both the fuel packing and the fuel movement during cycling. Previous works that have performed 3D full-core DEM simulation of PBRs have used simplified models that neglect reflector wall features. This work employs a graphics processing unit (GPU)-enabled DEM code, Project Chrono, to analyze the differences in pebble packing and pebble velocities between a simplified smooth PBR reflector and a more realistic reflector that includes circular wall features. Additionally, a sensitivity study is performed on the depth of the wall features to ensure that crystallization is prevented. Project Chrono is also validated for PBR cycling applications using experimental data. It is found that wall features with a depth of at least 0.5 pebble diameters significantly reduce crystallization in the near-wall region, leading to discrepancies in both packing fraction and pebble velocity in this region compared to the simplified reflector models. These discrepancies are found to lead to roughly a 5–10% difference in the prediction of the near-wall porosity and a 10% difference in the prediction of the velocity of pebbles near the wall. As a result of these discrepancies, it is suggested that future DEM simulations of PBRs include wall features to reduce modeling errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Verification of diesel spray ignition phenomenon in dual-fuel diesel-piloted premixed natural gas engine

Dual-fuel (DF) engines, in which premixed natural gas and air in an open-type combustion chamber is ignited by diesel-fuel pilot sprays, have been more popular for marine use than pre-chamber spark ignition (PCSI) engines because of their superior durability. However, control of ignition and combustion in DF engines is more difficult than in PCSI engines. In this context, here we focus on the ignition stability of n-heptane pilot-fuel jets injected into a compressed premixed charge of natural gas and air at low-load conditions. To aid understanding of the experimental data, chemical-kinetics simulations were carried out in a simplified engine-environment that provided insight into the chemical effects of methane (CH4) on pilot-fuel ignition. The simulations reveal that CH4 has an effect on both stages of n-heptane autoignition: the small, first-stage, cool-flame-type, low-temperature ignition (LTI) and the larger, second-stage, high-temperature ignition (HTI). As the ratio of pilot-fuel to CH4 entrained into the spray decreases, the initial oxidization of CH4 consumes the OH radicals produced by pilot-fuel decomposition during LTI, thereby inhibiting its progression to HTI. Using imaging diagnostics, the spatial and temporal progression of LTI and HTI in DF combustion are measured in a heavy-duty optical engine, and the imaging data are analyzed to understand the cause of severe fluctuations in ignition timing and combustion completeness at low-load conditions. Images of cool-flame and hydroxyl radical (OH*) chemiluminescence serve as indicators of LTI and HTI, respectively. The cycle-to-cycle and spatial variation in ignition extracted from the imaging data are used as key metrics of comparison. The imaging data indicate that the local concentration of the pilot-fuel and the richness of the surrounding natural-gas air mixture are important for LTI and HTI, but in different ways. In particular, higher injection pressures and shorter injection durations increase the mixing rate, leading to lower concentrations of pilot-fuel more quickly, which can inhibit HTI even as LTI remains relatively robust. Decreasing the injection pressure from 80 MPa to 40 MPa and increasing the injection duration from 500 µs to 760 µs maintained constant pilot-fuel mass, while promoting robust transition from LTI to HTI by effectively slowing the mixing rate. This allows enough residence time for the OH radicals, produced by the two-stage ignition chemistry of the pilot-fuel, to accelerate the transition from LTI to HTI before being consumed by CH4 oxidation. Thus from a practical perspective, for a premixed natural gas fuel–air equivalence-ratio, it is possible to improve the “stability” of the combustion process by solely manipulating the pilot-fuel injection parameters while maintaining constant mass of injected pilot-fuel. This allows for tailoring mixing trajectories to offset changes in fuel ignition chemistry, so as to promote a robust transition from LTI to HTI by changing the balance between the local concentration of the pilot-fuel and richness of the premixed natural gas and air. This could prove to be a valuable tool for combustion design to improve fuel efficiency or reduce noise or perhaps even reduce heat-transfer losses by locating early combustion away from in-cylinder walls.

42 ENGINEERING↗

Wave-packet manipulation of He Rydberg states by a seeded free-electron laser

We report a two-dimensional pump-control-probe spectroscopic study of the dynamics of singly excited He Rydberg-state wave packets with a seeded extreme ultraviolet (XUV) free-electron laser (FEL) source. A pair of coherent XUV pulses, defined by their coarse time separation and relative phase, created and manipulated the wave packets. The He atoms were postionized by infrared (IR) pulses, and the ion yield was measured as a function of XUV phase and IR arrival time. We tagged and sorted the relative phase of the XUV pulse pair on a single-shot basis by fitting each FEL spectrum with a suitable function that accounts for nonidealities of the XUV pulse pairs, associated with the seeding process; more generally, the fit returns the time-dependent electric field of the FEL spectra. The experimental two-dimensional maps of ion yields, measured as a function of IR (probe) delay and of XUV (pump-control) phase, were compared with the solution of the first-order time-dependent Schrödinger equation for this field. Despite the fact that the experimental conditions imply strong excitation, beyond the approximations of first-order perturbation theory, the simulated map satisfactorily reproduces the experimental one for temporally well-separated pulses. We show that by selecting data at appropriate values of pump-control phase, we enhance or suppress the amplitude of chosen wave-packet components consisting of two or more Rydberg states. When the temporal overlap of the pulse pair cannot be neglected, the phase reconstruction is underdetermined, and we provide a simplified comparison between data and simulations. Published by the American Physical Society 2024

Dumergue, Mathieu (ORCID:0000000236984406)↗

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗

Towards a High Fidelity Training Environment for Autonomous Cyber Defense Agents

Cyber defenders are overwhelmed by the frequency and scale of attacks against their networks. This problem will only be exacerbated as attackers leverage AI to automate their workflows. Autonomous cyber defense capabilities could aid defenders by automating operations and adapting dynamically to novel threats. However, existing training environments fall short in areas such as generalization, explainability, scalability, and transferability, making it intractable to train agents that will be effective in real networks. In this paper we take an important step towards creating autonomous cyber defense agents — we present a high fidelity training environment called Cyberwheel that includes both simulation and emulation capabilities. Cyberwheel simplifies customization of the training network and easily allows redefining the agent’s reward function, observation space, and action space to support rapid experimentation of novel approaches to agent design. It also provides visibility into agent behaviors necessary for agent evaluation and sufficient documentation / examples to lower the barrier to entry. As an example use case of Cyberwheel, we present initial results training an autonomous agent to deploy cyber deception strategies in simulation.

Oesch, T↗