Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic Timing Analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Ultraviolet relaxation dynamics of aniline, N, N-dimethylaniline and 3,5-dimethylaniline at 250 nm

Time-resolved photoelectron imaging was used to investigate the electronic relaxation dynamics of gas-phase aniline, N, N-dimethylaniline, and 3,5-dimethylaniline following ultraviolet excitation at 250 nm. Our analysis was supported by ab initio coupled-cluster calculations evaluating excited state energies and (in aniline) the evolution of a range of excited state physical properties as a function of N–H bond extension. Due to a lack of consistency between several earlier studies undertaken in aniline, the specific aim of this present work was to gain new insight into the previously proposed non-adiabatic coupling interaction between the two lowest lying singlet excited states S{sub 1}(ππ{sup ∗}) and S{sub 2}(3s/πσ{sup ∗}). The methyl-substituted systems N, N-dimethylaniline and 3,5-dimethylaniline were included in order to obtain more detailed dynamical information about the key internal molecular coordinates that drive the S{sub 1}(ππ{sup ∗})/S{sub 2}(3s/πσ{sup ∗}) coupling mechanism. Our findings suggest that in all three systems, both electronic states are directly populated during the initial excitation, with the S{sub 2}(3s/πσ{sup ∗}) state then potentially decaying via either direct dissociation along the N–X stretching coordinate (X = H or CH{sub 3}) or internal conversion to the S{sub 1}(ππ{sup ∗}) state. In aniline and N, N-dimethylaniline, both pathways most likely compete in the depletion of S{sub 2}(3s/πσ{sup ∗}) state population. However, in 3,5-dimethylaniline, only the direct dissociation mechanism appears to be active. This is rationalized in terms of changes in the relative rates of the two decay pathways upon methylation of the aromatic ring system.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Recovering missing CFD data for high-order discretizations using deep neural networks and dynamics learning

Data I/O poses a significant bottleneck in large-scale CFD simulations; thus, practitioners would like to significantly reduce the number of times the solution is saved to disk, yet retain the ability to recover any field quantity (at any time instance) a posteriori . The objective of this work is therefore to accurately recover missing CFD data a posteriori at any time instance, given that the solution has been written to disk at only a relatively small number of time instances. We consider in particular high-order discretizations (e.g., discontinuous Galerkin), as such techniques are becoming increasingly popular for the simulation of highly separated flows. To satisfy this objective, this work proposes a methodology consisting of two stages: 1) dimensionality reduction and 2) dynamics learning. For dimensionality reduction, we propose a novel hierarchical approach. First, the method reduces the number of degrees of freedom within each element of the high-order discretization by applying autoencoders from deep learning. Second, the methodology applies principal component analysis to compress the global vector of encodings. This leads to a low-dimensional state, which associates with a nonlinear embedding of the original CFD data. For dynamics learning, we propose to apply regression techniques (e.g., kernel methods) to learn the discrete-time velocity characterizing the time evolution of this low-dimensional state. In conclusion, a numerical example on a large-scale CFD example characterized by nearly 13 million degrees of freedom illustrates the suitability of the proposed method in an industrial setting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

OPTICAL AND DYNAMICAL CHARACTERIZATION OF COMET-LIKE MAIN-BELT ASTEROID (596) SCHEILA

We present observations and a dynamical analysis of the comet-like main-belt object, (596) Scheila. V-band photometry obtained on UT 2010 December 12 indicates that Scheila's dust cloud has a scattering cross-section {approx}1.4 times larger than that of the nucleus, corresponding to a dust mass of M{sub d} {approx} 3 Multiplication-Sign 10{sup 7} kg. V-R color measurements indicate that both the nucleus and dust are redder than the Sun, with no significant color differences between the dust cloud's northern and southern plumes. We also undertake an ultimately unsuccessful search for CN emission, where we find CN and H{sub 2}O production rates of Q{sub CN} < 9 Multiplication-Sign 10{sup 23} s{sup -1} and Q{sub H{sub 2O}}<10{sup 27} s{sup -1}. Numerical simulations indicate that Scheila is dynamically stable for >100 Myr, suggesting that it is likely native to its current location. We also find that it does not belong to a dynamical asteroid family of any significance. We consider sublimation-driven scenarios that could produce the appearance of multiple plumes of dust emission, but reject them as being physically implausible. Instead, we concur with previous studies that the unusual morphology of Scheila's dust cloud is most simply explained by a single oblique impact, meaning that this object is likely not a main-belt comet but is instead the second disrupted asteroid after P/2010 A2 (LINEAR) to be discovered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Synthesizing Virtual Oscillators to Control Islanded Inverters

Virtual oscillator control (VOC) is a decentralized control strategy for islanded microgrids where inverters are regulated to emulate the dynamics of weakly nonlinear oscillators. Compared to droop control, which is only well defined in sinusoidal steady state, VOC is a time-domain controller that enables interconnected inverters to stabilize arbitrary initial conditions to a synchronized sinusoidal limit cycle. However, the nonlinear oscillators that are elemental to VOC cannot be designed with conventional linear-control design methods. We address this challenge by applying averaging- and perturbation-based nonlinear analysis methods to extract the sinusoidal steady-state and harmonic behavior of such oscillators. The averaged models reveal conclusive links between real- and reactive-power outputs and the terminal-voltage dynamics. Similarly, the perturbation methods aid in quantifying higher order harmonics. The resultant models are then leveraged to formulate a design procedure for VOC such that the inverter satisfies standard ac performance specifications related to voltage regulation, frequency regulation, dynamic response, and harmonic content. Experimental results for a single-phase 750 VA, 120 V laboratory prototype demonstrate the validity of the design approach. They also demonstrate that droop laws are, in fact, embedded within the equilibria of the nonlinear-oscillator dynamics. This establishes the backward compatibility of VOC in that, while acting on time-domain waveforms, it subsumes droop control in sinusoidal steady state.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

A Dynamic Trajectory Fit to Multisensor Fireball Observations

Meteorites with known orbital origins are key to our understanding of solar system formation and the source of life on Earth. Fireball networks have been developed globally in a unified effort to record and ultimately retrieve these cosmic samples. However, the accuracy of the determined orbit and the likelihood of meteorite recovery depend directly on the accuracy of the chosen meteoroid triangulation method. There are three leading techniques for meteoroid triangulation discussed in the literature: the method of planes, the straight-line least-squares method, and the multiparameter fit method. Here we describe an alternative method to meteoroid triangulation, called the dynamic trajectory fit. This approach uses the meteoroid’s 3D dynamic equations of motion to fit a realistic trajectory directly to multisensor line-of-sight observations. This method has the ability to resolve fragmentation events, fit systematic observatory timing offsets, and determine mass estimates of the meteoroid along its observable trajectory. Through a comprehensive Monte Carlo analysis of over 100,000 trajectory simulations, we find this new method to more accurately estimate meteoroid trajectories of slow entry events (<25 km s{sup −1}) and events observed from low convergence angles (<10°) compared to existing meteoroid triangulation techniques. Additionally, we triangulate an observed fireball event with visible fragmentation using the various triangulation methods to show that the proposed dynamic trajectory fit implementing fragmentation to best match the captured multisensor line-of-sight data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Examining the role of fluctuations in the early stages of homogenous polymer crystallization with simulation and statistical learning

Here, we propose a relationship between the dynamics in the amorphous and crystalline domains during polymer crystallization: the fluctuations of ordering-rate about a material-specific value in the amorphous phase drive those fluctuations associated with the increase in percent crystallinity. This suggests a differential equation that satisfies the three experimentally observed time regimes for the rate of crystal growth. To test this postulated expression, we applied a suite of statistical learning tools to molecular dynamics simulations to extract the relevant phenomenology. This study shows that the proposed relationship holds in the early time regime. It illustrates the effectiveness of soft computing tools in the analysis of coarse-grained simulations in which patterns exist, but may not easily yield to strict quantitative evaluation. This ability assists us in characterizing the critical early time molecular arrangement during the primary nucleation phase of polymer melt crystallization. In addition to supporting the validity of the proposed kinetics expression, the simulations show that (i) the classical nucleation and growth mechanism is active in the early stages of ordering; (ii) the number of nuclei and their masses grow linearly during this early time regime; and (iii) a fixed inter-nuclei distance is established.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Material coherence from trajectories via Burau eigenanalysis of braids

In this paper, we provide a numerical tool to study a material’s coherence from a set of 2D Lagrangian trajectories sampling a dynamical system, i.e., from the motion of passive tracers. We show that eigenvectors of the Burau representation of a topological braid derived from the trajectories have levelsets corresponding to components of the Nielsen–Thurston decomposition of the dynamical system. One can thus detect and identify clusters of space–time trajectories corresponding to coherent regions of the dynamical system by solving an eigenvalue problem. Unlike previous methods, the scalable computational complexity of our braid-based approach allows the analysis of large amounts of trajectories.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Instantaneous difference frequency locking observed during toroidicity-induced Alfvén eigenmode coupling in the DIII-D tokamak

In magnetic confinement fusion, toroidicity-induced Alfvén eigenmodes (TAEs) are well-studied, weakly stable solutions of the linearized ideal magnetohydrodynamic equations. Driven unstable by suprathermal populations of energetic particles, TAE pose a key vulnerability to the confinement of high-energy alpha particles generated by fusion reactions. Hence, it is paramount to understand TAE dynamics if a working reactor is to be realized. In this work, we detect and characterize signatures of nonstationary nonlinear coupling between TAE using a novel, time-resolved bispectral analysis; results are supported by analytic signal of band-passed data. Crucially, a stationary phase relationship between two TAE and a nascent low frequency fluctuation is observed precisely when the triple product of magnetic fluctuation amplitudes is enhanced. Local mode number and frequency spectrum, gleaned from beam-emission spectroscopy, corroborates simultaneous satisfaction of nonlinear matching conditions, and provides a tool to identify theorized pathways of energy transfer, e.g. TAE parametric instability.

bispectral analysis↗
DOE OSTI2021

Protist diversity and community complexity in the rhizosphere of switchgrass are dynamic as plants develop

Abstract Background Despite their widespread distribution and ecological importance, protists remain one of the least understood components of the soil and rhizosphere microbiome. Knowledge of the roles that protists play in stimulating organic matter decomposition and shaping microbiome dynamics continues to grow, but there remains a need to understand the extent to which biological and environmental factors mediate protist community assembly and dynamics. We hypothesize that protists communities are filtered by the influence of plants on their rhizosphere biological and physicochemical environment, resulting in patterns of protist diversity and composition that mirror previously observed diversity and successional dynamics in rhizosphere bacterial communities. Results We analyzed protist communities associated with the rhizosphere and bulk soil of switchgrass (SG) plants ( Panicum virgatum ) at different phenological stages, grown in two marginal soils as part of a large-scale field experiment. Our results reveal that the diversity of protists is lower in rhizosphere than bulk soils, and that temporal variations depend on soil properties but are less pronounced in rhizosphere soil. Patterns of significantly prevalent protists groups in the rhizosphere suggest that most protists play varied ecological roles across plant growth stages and that some plant pathogenic protists and protists with omnivorous diets reoccur over time in the rhizosphere. We found that protist co-occurrence network dynamics are more complex in the rhizosphere compared to bulk soil. A phylogenetic bin-based null model analysis showed that protists’ community assembly in our study sites is mainly controlled by homogenous selection and dispersal limitation, with stronger selection in rhizosphere than bulk soil as SG grew and senesced. Conclusions We demonstrate that environmental filtering is a dominant determinant of overall protist community properties and that at the rhizosphere level, plant control on the physical and biological environment is a critical driver of protist community composition and dynamics. Since protists are key contributors to plant nutrient availability and bacterial community composition and abundance, mapping and understanding their patterns in rhizosphere soil is foundational to understanding the ecology of the root-microbe-soil system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Dynamic Analysis Methods for Detecting Anomalies in Asynchronously Interacting Systems

Detecting modifications to digital system designs, whether malicious or benign, is problematic due to the complexity of the systems being analyzed. Moreover, static analysis techniques and tools can only be used during the initial design and implementation phases to verify safety and liveness properties. It is computationally intractable to guarantee that any previously verified properties still hold after a system, or even a single component, has been produced by a third-party manufacturer. In this paper we explore new approaches for creating a robust system design by investigating highly-structured computational models that simplify verification and analysis. Our approach avoids the need to fully reconstruct the implemented system by incorporating a small verification component that dynamically detects for deviations from the design specification at run-time. The first approach encodes information extracted from the original system design algebraically into a verification component. During run-time this component randomly queries the implementation for trace information and verifies that no design-level properties have been violated. If any deviation is detected then a pre-specified fail-safe or notification behavior is triggered. Our second approach utilizes a partitioning methodology to view liveness and safety properties as a distributed decision task and the implementation as a proposed protocol that solves this task. Thus the problem of verifying safety and liveness properties is translated to that of verifying that the implementation solves the associated decision task. We develop upon results from distributed systems and algebraic topology to construct a learning mechanism for verifying safety and liveness properties from samples of run-time executions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

stochprop Documentation (Release 1.0)

Simulations of infrasonic propagation in the atmosphere typically utilize a single atmospheric specification describing the acoustic sound speed, ambient winds, and density as a function of altitude. Due to the dynamic and sparsely sampled nature of the atmosphere, there is a notable amount of uncertainty in the atmospheric state at a given location and time so that a more robust analysis of infrasonic propagation requires inclusion of this uncertainty. This Python library, stochprop, has been implemented using methods developed jointly by infrasound scientists at Los Alamos National Laboratory (LANL) and the University of Mississippi’s National Center for Physical Acoustics (NCPA). This software library includes methods to quantify variability in the atmospheric state, identify typical seasonal variability in the atmospheric state and generate suites of representative atmospheric states during a given season, as well as perform uncertainty analysis on a specified atmospheric state given some level of uncertainty. These methods have been designed to interface between propagation modeling capabilities such as InfraGA/GeoAc and NCPAprop and signal analysis methods in the LANL InfraPy tool.

58 GEOSCIENCES↗
DOE OSTI2026

Reconstructing the Stripping History of the Sagittarius Stream with Neural Networks

The Sagittarius (Sgr) Stream is produced by the ongoing disruption of the Sgr dwarf spheroidal (dSph) galaxy and is thought to contain multiple wraps that were stripped during different pericentric passages. In this study, we introduce a neural-network–based method trained on N-body simulations to infer the stripping time of Sgr Stream stars directly from their phase-space coordinates. We combine spectroscopic data from SEGUE, APOGEE DR17, and LAMOST DR7 low-resolution spectroscopic (LRS) survey with Gaia EDR3 astrometry and distance estimates from the latest StarHorse catalog to identify high-quality Sgr Stream members. Applying our method to these stars, we measure a clear metallicity gradient with stripping time, well described by a linear relation with slope ∼0.3 dex Gyr −1 . We further predict the stripping times of globular clusters previously suggested to originate from the Sgr dSph. M 54, Terzan 7, Terzan 8, and Arp 2 exhibit stripping times consistent with being currently bound to the Sgr remnant. Pal 12, Whiting 1, and NGC 2419 are inferred to have been stripped 0.9 ± 0.1, 1.1 ± 0.2, and 2.1 ± 0.2 Gyr ago, respectively. For NGC 4147 and NGC 5634, whose membership in the Sgr system remains uncertain, our analysis suggests stripping times of 1.1 ± 0.4 and 1.1 ± 0.1 Gyr, respectively, if they are ultimately confirmed as genuine Sgr members. These results demonstrate that data-driven models of dynamical stripping histories offer a promising approach for reconstructing the formation and chemical evolution of the Sgr Stream.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Direct observation of coherent femtosecond solvent reorganization coupled to intramolecular electron transfer

It is well known that the solvent plays a critical role in ultrafast electron-transfer reactions. However, solvent reorganization occurs on multiple length scales, and selectively measuring short-range solute–solvent interactions at the atomic level with femtosecond time resolution remains a challenge. Here we report femtosecond X-ray scattering and emission measurements following photoinduced charge-transfer excitation in a mixed-valence bimetallic (Fe ii Ru iii ) complex in water, and their interpretation using non-equilibrium molecular dynamics simulations. Combined experimental and computational analysis reveals that the charge-transfer excited state has a lifetime of 62 fs and that coherent translational motions of the first solvation shell are coupled to the back electron transfer. Overall, our molecular dynamics simulations identify that the observed coherent translational motions arise from hydrogen bonding changes between the solute and nearby water molecules upon photoexcitation, and have an amplitude of tenths of ångströms, 120–200 cm -1 frequency and ~100 fs relaxation time. This study provides an atomistic view of coherent solvent reorganization mediating ultrafast intramolecular electron transfer.

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

Sources of Uncertainty in Modeled Land Carbon Storage within and across Three MIPs: Diagnosis with Three New Techniques

Terrestrial carbon cycle models have incorporated increasingly more processes as a means to achieve more-realistic representations of ecosystem carbon cycling. Despite this, there are large across-model variations in the simulation and projection of carbon cycling. Several model intercomparison projects (MIPs), for example, the fifth phase of the Coupled Model Intercomparison Project (CMIP5) (historical simulations), Trends in Net Land–Atmosphere Carbon Exchange (TRENDY), and Multiscale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP), have sought to understand intermodel differences. Here in this study, the authors developed a suite of new techniques to conduct post-MIP analysis to gain insights into uncertainty sources across 25 models in the three MIPs. First, terrestrial carbon storage dynamics were characterized by a three-dimensional (3D) model output space with coordinates of carbon residence time, net primary productivity (NPP), and carbon storage potential. The latter represents the potential of an ecosystem to lose or gain carbon. This space can be used to measure how and why model output differs. Models with a nitrogen cycle generally exhibit lower annual NPP in comparison with other models, and mostly negative carbon storage potential. Second, a transient traceability framework was used to decompose any given carbon cycle model into traceable components and identify the sources of model differences. The carbon residence time (or NPP) was traced to baseline carbon residence time (or baseline NPP related to the maximum carbon input), environmental scalars, and climate forcing. Third, by applying a variance decomposition method, the authors show that the intermodel differences in carbon storage can be mainly attributed to the baseline carbon residence time and baseline NPP (>90% in the three MIPs). The three techniques developed in this study offer a novel approach to gain more insight from existing MIPs and can point out directions for future MIPs. Finally, since this study is conducted at the global scale for an overview on intermodel differences, future studies should focus more on regional analysis to identify the sources of uncertainties and improve models at the specified mechanism level.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

TU-CD-207-03: Time Evolution of Texture Parameters of Subtracted Images Obtained by Contrast-Enhanced Digital Mammography (CEDM)

Purpose: To evaluate the time evolution of texture parameters, based on the gray level co-occurrence matrix (GLCM), in subtracted images of 17 patients (10 malignant and 7 benign) subjected to contrast-enhanced digital mammography (CEDM). The goal is to determine the sensitivity of texture to iodine uptake at the lesion, and its correlation (or lack of) with mean-pixel-value (MPV). Methods: Acquisition of clinical images followed a single-energy CEDM protocol using Rh/Rh/48 kV plus external 0.5 cm Al from a Senographe DS unit. Prior to the iodine-based contrast medium (CM) administration a mask image was acquired; four CM images were obtained 1, 2, 3, and 5 minutes after CM injection. Temporal series were obtained by logarithmic subtraction of registered CM minus mask images.Regions of interest (ROI) for the lesion were drawn by a radiologist and the texture was analyzed. GLCM was evaluated at a 3 pixel distance, 0° angle, and 64 gray-levels. Pixels identified as registration errors were excluded from the computation. 17 texture parameters were chosen, classified according to similarity into 7 groups, and analyzed. Results: In all cases the texture parameters within a group have similar dynamic behavior. Two texture groups (associated to cluster and sum mean) show a strong correlation with MPV; their average correlation coefficient (ACC) is r{sup 2}=0.90. Other two groups (contrast, homogeneity) remain constant with time, that is, a low-sensitivity to CM uptake. Three groups (regularity, lacunarity and diagonal moment) are sensitive to CM uptake but less correlated with MPV; their ACC is r{sup 2}=0.78. Conclusion: This analysis has shown that, at least groups associated to regularity, lacunarity and diagonal moment offer dynamical information additional to the mean pixel value due to the presence of CM at the lesion. The next step will be the analysis in terms of the lesion pathology. Authors thank PAPIIT-IN105813 for support. Consejo Nacional de Ciencia Y Tecnologia, PAPIIT-IN105813.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Discharge dynamics and plasma density recovery by on/off switches of additional gas

Measurement of the plasma density is investigated to study plasma dynamics by adding reactive gas (O{sub 2}) or rare gas (He) in Ar plasmas. When the O{sub 2} or He gas is added, plasma density is suddenly decreased, while the plasma density recovers slowly with gas off. It is found that the recovery time is strongly dependent on the gas flow rate, and it can be explained by effect of gas residence time. When the He gas is off in the Ar plasma, the plasma density is overshot compared to the case of the O{sub 2} gas pulsing due to enhanced ionizations by metastable atoms. Analysis and calculation for correlation between the plasma density dynamics and the gas pulsing are also presented in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

An integrated analytical and experimental study of contact acoustic nonlinearity at rough interfaces of fatigue cracks

This paper presents an analytical model that describes the relationship between interface stiffness and separation of rough interfaces such as fatigue cracks. The contact acoustic nonlinearity at the interface is simulated by a quasi-static model based on Hertzian contact theory. The model is validated using the results of dynamic acoustoelastic testing (DAET) with a Rayleigh wave probe on fatigue cracks in two aluminum alloy samples. One novel aspect of this work is that all the required geometrical parameters for the model are acquired directly from the aperture profile of real cracks extracted from their scanning electron microscopy (SEM) images. Also, the separation of the crack faces during dynamic perturbation is independently measured using a 3D laser Doppler vibrometer (LDV). In addition, using DAET allows for an unprecedented direct comparison between experimental and analytical results. Overall, this is because unlike conventional nonlinear ultrasonic tests that are based on measuring the amplitude of higher harmonics, DAET outputs directly the strain-dependency of transmission and time delay of ultrasonic waves propagating across the interface. The model is found in good qualitative agreement with the experimental results, although the predictions tend to underestimate the variation of transmission coefficient and time delay. We conduct a sensitivity analysis to investigate the influence of different assumptions and simplifications on model predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Methodology and Tool for the Physical Security Analysis of Micro and Advanced Reactors

This work proposes a dynamic evaluation methodology to relax the conservatism in physical security evaluation, by leveraging an ongoing work in the Light Water Reactor Sustainability pathway. This methodology is implemented in a dynamic risk assessment tool named Event Modeling Risk Assessment using Linked Diagrams (EMRALD). The work extends EMRALD’s capability to support a sandbox feature where analysts can easily create attack scenarios and modify advanced/small modular reactor (A/SMR) security and safety features using templates. This approach saves time and cost since the analysis does not require creating detailed computer-aided design models, as is commonly required in commercial force-on-force software tools. EMRALD is completely free to use at https://emraldapp.inl.gov. We have developed basic templates including physical barriers, intrusion sensors, physical areas, and safety actions, that can be downloaded from EMRALD’s GitHub site: https://github.com/idaholab/EMRALD. These templates use generic data commonly used for training purposes, which do not reflect any actual operating nuclear reactor. Users may adjust the data in the templates with their own dataset and/or create new templates in EMRALD. The proposed methodology combines security and safety by assessing sabotage effects up to the radiological consequence to the public instead of merely the core damage state. This practice follows the industry standard for advanced non-light-water reactors currently proposed for endorsement by the Nuclear Regulatory Commission. The combination of security and safety is expressed in an achievability-consequence chart. EMRALD can be used to generate data for this chart. A hypothetical case study using a representative sodium-cooled fast reactor (SFR) facility is presented in this report to demonstrate this methodology. This case study does not contain any actual nuclear plant information. This work will benefit A/SMR vendors and utilities to implement security by design during the reactor design iteration phase, such that they do not have to perform upgrades and retrofits to the reactor after it is installed to improve its physical protection system. The tool may also be used to analyze domestic or foreign reactor designs to support the International Nuclear Security Techniques for Advanced Reactors (INSTAR) bilateral missions. Future works are planned to implement the methodology on a reference SFR reactor and a reference high-temperature gas-cooled reactor to obtain insights and lessons-learned for the A/SMR community.

97 MATHEMATICS AND COMPUTING↗