SEARCH · Engineering Papers

Results for “Simulation Results”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗

DOE OSTI2024

Orbit-following simulations of fast-ion transport and losses due to the Alfvén eigenmode burst in the Large Helical Device

Orbit-following simulations of fast-ion transport and losses with time-dependent electromagnetic perturbations are performed to clarify the roles of Alfvén eigenmodes (AEs) and the low-frequency magnetohydrodynamic (MHD) mode observed in the kinetic-MHD hybrid simulation of AE bursts in the Large Helical Device. Fast-ion pressure profile flattening in the kinetic-MHD hybrid simulation can be reproduced by an orbit-following simulation with only the primary single AE of the time-dependent amplitude following the kinetic-MHD hybrid simulation result, while orbit-following simulations with constant AE amplitude of average level during AE burst cannot reproduce the fast-ion pressure profile flattening observed. The effects of other modes are negligible on the fast-ion pressure profile flattening. The fast-ion losses in kinetic-MHD hybrid simulation can be reproduced by an orbit-following simulation with time-dependent amplitude when the low-frequency MHD mode is considered in addition to multiple AEs. This indicates the synergetic effect of multiple AEs and the low-frequency MHD mode on fast-ion losses.

Seki, Ryohsuke↗

DOE OSTI2022

Non-Conductor-Contact Surface Wave Reflectometry for Cable Insulation Damage Detection

Here, this paper presents a non-conductor-contact surface wave-based reflectometry technique to detect and locate missing or damaged insulation and/or jacket of electrical cables. The described sensor may be uniquely and practically implemented on live cables as a non-destructive evaluation (NDE) method since it may be applied to the outer insulation/jacket of a cable and requires no direct contact to the conductor. The innovation lies in a high-frequency surface wave launcher that can introduce broadband signals in the GHz frequencies that insure maximum electromagnetic (EM) energy localization within the insulation and/or jacket material. By further processing the reflection data in the frequency domain, the location and severity of insulation damage is detected. Simulated results of thin, thick, short, and long cables show the results expected from small, medium, and large damage regions. For example, the detection of 40 mm-long damage region at 40 m distance on a 50 m long cable (2.44 cm diameter) is demonstrated. Experimental results are presented that validate simulation findings.

Non Destructive Evaluation↗

DOE OSTI2025

Data Files for Runoff Evaluation in an Earth System Land Model for Permafrost Regions

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce.This data set includes all files that were produced and applied in the paper Runoff Evaluation in an Earth System Land Model for Permafrost Regions [Xiang et al. in review]. The paper is in review as of July 1 2025 in Geoscientific Model Development (GMD). In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-rich integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM’s parameterized representation of total runoff. This dataset contains 2 figure image files (*.png, *jpg) that describe the study site and methods, as well as folders (Figure*.zip) that contain the associated data files (*.csv, *.dat) and python code notebooks (*.ipynb) for figures 3-7 in the paper. Jupyter notebook (*.ipynb) files that produce the figure files using the associated data files will run within a python environment configured with Jupyter Lab or Notebook packages.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Preliminary Design of Reactor Cavity Cooling System for a Horizontal Compact HTGR

The Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed MIGHTR in 3 years and support its commercialization as a safe and low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. This report documents the preliminary design study of the RCCS for the HC-HTGR. It includes the establishment of the design requirements, a high-level design study by initial scoping calculations, and preliminary performance calculations of the HC-HTGR RCCS design. Design requirements for the HC-HTGR RCCS have been established to guide preliminary design activities and scoping performance calculations. Initial scoping calculations including estimation of the water inventory, estimation of HVAC thermal capability, and a parametric study on loop dimensions by standalone RCCS analysis. Based on scoping calculation results, a set of baseline dimensions of the HC-HTGR RCCS was derived. A water panel modeling approach was investigated to explore various potential design options for the water panel under consideration for the HC-HTGR RCCS using RELAP5-3D. A test case study was performed to assess the prediction capability of two modeling approaches. The results were compared with CFD simulations conducted in constant RPV temperature and heat flux boundary conditions. It confirms the capability of the RELAP5-3D modeling approach to include all important heat transfer mechanisms expected in the HC-HTGR RCCS operation conditions. Then, a reference RELAP5-3D model for the 1/8 th of a compartment of the preliminary design of the HC-HTGR RCCS was developed. A preliminary performance analysis was conducted to evaluate single-phase natural circulation performance with different top tank temperature values and panel conduction performance in various operation conditions. From single-phase natural circulation performance analysis, the system operation mode was investigated in normal operating and limiting design conditions. It showed operation mode in a subcooled state with a proper top tank water cooling system. Parasitic heat loss by both internal air flow and RCCS was estimated, showing it satisfies maintaining below target maximum heat loss of the HC-HTGR RCCS. From the panel conduction performance analysis, two candidate materials for the riser tube such as carbon steel and stainless steel were compared in the thermal performance of HC-HTGR RCCS. From a single water panel test compared with CFD simulation results, it was confirmed that the current capability of the RELAP5-3D modeling approach for the water panel predicts the thermal conduction of two different materials of the water panel. Then, system-level thermal performance analysis was performed for 1/8 th of the compartment of the preliminary HC-HTGR RCCS design. It was first observed that the current preliminary HC-HTGR RCCS design had minimal impact on the overall thermal performance of the water panel by changing pipe material from carbon steel to stainless steel. From Argonne’s effort on the ongoing water-based NSTF testing program, several considerations other than the thermal performance point of view were addressed to be considered in selecting pipe materials.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Large‐Eddy Simulations of a Convection Cloud Chamber: Sensitivity to Bin Microphysics and Advection

Abstract Bin microphysics schemes are useful tools for cloud simulations and are often considered to provide a benchmark for model intercomparison. However, they may experience issues with numerical diffusion, which are not well quantified, and the transport of hydrometeors depends on the choice of advection scheme, which can also change cloud simulation results. Here, an atmospheric large‐eddy simulation model is adapted to simulate a statistically steady‐state cloud in a convection cloud chamber under well‐constrained conditions. Two bin microphysics schemes, a spectral bin method and the method of moments, as well as several advection methods for the transport of the microphysical variables are employed for model intercomparison. Results show that different combinations of microphysics and advection schemes can lead to considerable differences in simulated cloud properties, such as cloud droplet number concentration. We find that simulations using the advection scheme that suffers more from numerical diffusion tends to have a smaller droplet number concentration and liquid water content, while simulation with the microphysics scheme that suffers more from numerical diffusion tends to have a broader size distribution and thus larger mean droplet sizes. Sensitivities of simulations to bin resolution, spatial resolution, and temporal resolution are also tested. We find that refining the microphysical bin resolution leads to a broader cloud droplet size distribution due to the advection of hydrometeors. Our results provide insight for using different advection and microphysics schemes in cloud chamber simulations, which might also help understand the uncertainties of the schemes used in atmospheric cloud simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

HERA Modeling and Simulation Exercise: BISON Results

A Modeling and Simulation (M&S) exercise is being performed for the High burnup Experiments for Reactivity initiated Accident (HERA) project under the Nuclear Energy Agency (NEA) Framework for Irradiation Experiments (FIDES) program. The goal of the M&S exercise is to improve M&S and experiment integration, facilitate community involvement in experiment design and interpretation, facilitate community collaboration, and aid in ensuring program data meet fuel performance code needs. The M&S exercise will compile and compare results from over 20 international organizations using 14 different fuel performance codes. This paper presents the results from the BISON fuel performance code generated by the Idaho National Laboratory participants.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Projecting climate change in South America using variable‐resolution Community Earth System Model: An application to Chile

Abstract We introduce variable‐resolution enabled Community Earth System Model (VR‐CESM) results simulating historical and future climate conditions at 28 km over South America and 14 km over the Andes. Three 30‐year simulations are performed: a historic (1985–2014), a near future (2030–2059), and an end‐century (2070–2099) simulation under the RCP8.5 scenario. Historic results compare favourably to several temperature and precipitation reanalysis products, though local biases are present, particularly during austral summer. Future simulations highlight broad warming patterns (+3–6°C by end‐century) and heterogeneous precipitation responses across South America that qualitatively agree with prior modelling efforts. Our results reveal that the interaction between temperature and precipitation changes produce shifts in several Köppen–Geiger climates. Notable changes include the near‐elimination of the Andean Tundra or Alpine climates, a 15% decrease in Tropical Rainforests and a Tropical Savannah expansion of 20%. To provide a regionally focused analysis of projected climate change and to illustrate the benefits of variable resolution modelling, we analyse changes in the magnitude and trend in seasonal and daily temperature and precipitation in Chile. We also examined several metrics [e.g., snow water equivalent (SWE), temperatures on wet days, and days below 0°C] to evaluate potential impacts of climate change on the Chilean cryosphere between the end‐of‐century and historic periods, finding wide‐ranging indications of cryospheric decline. These changes are interpreted through reductions in the timing (1–2.5 months earlier peak SWE) and magnitude (200–1,000 mm SWE decreases) of water stored as snow in the Andes, a 10–30% decrease in number of cool season wet days with temperatures below 1°C, and 50–200 fewer days (annually) with minimum temperatures below 0°C. Our aim in producing a high‐resolution dataset of climate projections from VR‐CESM is to support analyses of climate change throughout South America but especially in vulnerable montane regions and to provide additional results for comparison with previous, ongoing, and upcoming modelling efforts.

variable-resolution global climate models↗

DOE OSTI2021

Examining the Performance of Walsh-DSSS Against FBMC-SS in HF Channels

Abstract—Filter bank multicarrier spread spectrum (FBMCSS) has proven to be a robust and reliable waveform choice for communication over high frequency (HF) skywave links. However, the performance of this waveform has yet to be contextualized against typical robust HF waveforms, such as the Walsh-encoded waveform detailed in the MIL-STD-188-110D, Appendix D document. In this paper, we first outline the advantages of both the Walsh and FBMC-SS waveforms as well as present their developments. Simulation results are then presented for ideal, simulated HF, and HF with interference channel conditions. Lastly, skywave-HF results are presented for these two waveforms both with and without interference.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2025

Predicting long-term stress relaxation on Alloy 709 using the crystal plasticity finite element method

This report summarizes the results of physics-based, crystal plasticity simulations for the long-term stress relaxation behavior of Alloy 709. The purpose of the study was to provide insight into five key questions related to long-term behavior in high temperatures materials which are difficult or impossible to answer experimentally: (1) is there a threshold stress for long-term relaxation? (2) is there strain threshold for relaxation damage, below which significant damage does not accumulate? (3) does damage continue to accumulate as the material relaxes or will damage accumulation plateau under some loading conditions? (4) does stress relaxation loading inevitably lead to failure? and (5) which, if any, engineering models for relaxation damage accumulation reasonably match the simulation results? The report summarizes the numerical simulations used to address these five questions and provides at least partial answers to each question.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Impact of Cloud-Base Turbulence on CCN Activation: CCN Distribution

Abstract Following our previous investigation of the turbulence impact on cloud-base single-size CCN activation, this study considers a similar problem assuming CCN size distribution obtained from field measurements. The total CNN concentration is taken as either 200 cm −3 to represent clean conditions, or as 2000 cm −3 to represent polluted conditions. CCN is assumed to be sodium chloride. The CCN activation in the rising nonturbulent adiabatic parcel is contrasted with the activation within a rising adiabatic parcel filled with inertial-range homogeneous isotropic turbulence. The turbulent parcel of 64 3 m 3 and the turbulent kinetic energy dissipation rate of 10 −3 m −2 s −3 are used in most of the simulations. Results for a range of mean parcel ascent rates, between 0.125 and 8 m s −1 , are discussed. Overall, the adiabatic turbulent parcel simulations show results consistent with the adiabatic nonturbulent parcel, with higher activated CCN concentrations for stronger parcel ascent rates. The key difference is a blurriness of the separation between dry CCN size bins featuring activated and nonactivated (haze) CCN, especially for weak mean ascent rates. The blurriness comes from CCN getting activated and subsequently deactivated in the fluctuating supersaturation field, instead of all becoming cloud droplets above the cloud base. This leads to significantly larger spectral widths in turbulent parcel simulations compared to the nonturbulent parcel when activation is completed. Modeling results are discussed in the context of the impact of turbulent fluctuations on CCN activation documented in laboratory experiments using the Pi chamber.

Meteorology & Atmospheric Sciences↗

DOE OSTI2020

Modeling and analysis of heat emissions from buildings to ambient air

Heat emissions from buildings is a significant source of anthropogenic heat influencing the urban microclimate; however, they are usually oversimplified in urban climate and microclimate modeling. This study developed a bottom-up physics-based approach to calculate heat emissions from buildings to the ambient air and implemented the approach in EnergyPlus. In this work, a simple result verification was conducted by comparing the EnergyPlus simulated results against the spreadsheet calculations. Simulations covering 16 commercial building types, four climates, and two energy efficiency levels were conducted to understand and evaluate the building heat emissions and their temporal patterns as well as three major components: (1) building envelope (convective heat transfer to ambient air), (2) zones (air exfiltration and exhaust air), and (3) HVAC systems (relief air and heat rejection from condensers or cooling towers). The main findings are: (1) heat emissions are usually higher than the site energy use (about 2.5 times), and their dynamics should be considered; (2) building characteristics and their energy systems lead to differences in heat emission contributions from the three components, and their dynamics, for example, in the warehouse models, the envelope component accounts for 90.4%, while it is 12.7% for the large office models; (3) for most building typologies, the climate has a strong impact on heat emissions, for example, buildings with dominant heat emissions from the zone exhaust air and/or the HVAC reject heat, a general decrease in heat emissions in hotter climates is observed, while envelope-dominated buildings show the opposite; and (4) building technologies that reduce energy use in buildings may perform differently in reducing heat emissions. The developed heat emissions calculation method can be adopted in EnergyPlus and most other building energy modeling programs. It can provide dynamic building heat emissions as an input to urban climate computational fluid dynamics (CFD) models at a higher spatial and temporal resolution than is currently available, to improve the simulation accuracy of the urban microclimate and capture the urban heat island effect and urban overheating.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Optimization of Window Positions for Wind-Driven Natural Ventilation Performance

This paper optimizes opening positions on building facades to maximize the natural ventilation’s potential for ventilation and cooling purposes. The paper demonstrates how to apply computational fluid dynamics (CFD) simulation results to architectural design processes, and how the CFD-driven decisions impact ventilation and cooling: (1) background: A CFD helps predict the natural ventilation’s potential, the integration of CFD results into design decision-making has not been actively practiced; (2) methods: Pressure data on building facades were obtained from CFD simulations and mapped into the 3D modeling environment, which were then used to identify optimal positions of two openings of a zone. The effect of the selected opening positions was validated with building energy simulations; (3) results: The cross-comparison study of different window positions based on different geographical locations quantified the impact on natural ventilation effectiveness; and (4) conclusions: The optimized window position was shown to be effective, and some optimal solutions contradicted the typical cross-ventilation strategy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

METHODOLOGY AND APPLICATION OF PHYSICAL SECURITY EFFECTIVENESS BASED ON DYNAMIC FORCE-ON-FORCE MODELING

This paper describes ongoing work within the Light Water Reactor Sustainability (LWRS) Program at Idaho National Laboratory (INL) to optimize security and cost of nuclear power plants (NPPs). It reviews the conservatisms in conventional physical security posture and regulations. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). The dynamic assessment methodology leverages EMRALD to process results of force-on-force (FOF) simulations and crediting safety mitigation actions from probabilistic risk assessment (PRA) models as well as diverse and flexible coping strategies (FLEX) mitigation strategies. Timing information from these simulations are compared against the available time to perform mitigations obtained from Reactor Excursion and Leak Analysis Program (RELAP5) simulations. To illustrate the methodology, a station blackout (SBO) attack scenario was modeled in commercially available FOF simulation tools. The simulation results provide valuable insights into possible attack outcomes and as the probabilistic risk of a core damage event given these outcomes. Safety mitigation procedures were modeled in EMRALD, and were dependent on the attack outcomes by considering human operator uncertainties. RELAP5 simulations incorporating human and hardware uncertainties were performed to estimate the distribution of time-to-core damage. The results demonstrate that, even in the extreme case of a successful adversarial attack, plant mitigation strategies provide significantly high-likelihood of preventing radiological release. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the NPPs to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Analysis of buffer-IPyC separation in TRISO fuel particles

During High Temperature Gas-cooled Reactor (HTGR) operation, and due to the neutron irradiation in the reactor core, damage of the nuclear fuel coating layers occurs. The mechanism of damage formation in the TRISO fuel is explored by the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program, in which the debonding process between coating layers was also investigated. The purpose of this paper is to report simulation results for two models. Firstly, the debonding restricted model, where no gap formation between buffer and IPyC layers is permitted. Secondly, a debonding enabled model, where the gap between those layers is created. The simulations were performed with the Bison code. The inputs of the simulated models are based on data from the AGR-1 experiment. Further, the research included simulations on spherical and aspherical fuel types. The simulations match results obtained by the AGR-1 experiment, which as such shows that the Bison code is a good computational method for simulating the behavior of the gap between buffer and IPyC layers in TRISO fuel. Based on the irradiation experiments, and the Bison simulations, it was concluded that the most common scenario is a gap formation along the buffer-IPyC interface, while the least possible scenario is the situation where there is no gap formation at the buffer-IPyC junction. The computational results confirmed that the sphericity of the fuel influences the thickness of the gap that occurs at the buffer-IPyC junction, in a way that with increasing aspect ratio the gap thickness increases. In addition, the fuel sphericity does not influence the Weibull failure probability. The results obtained for spherical and aspherical fuel are nearly identical.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Sparse chronology strategy for integrating seasonal energy storage in capacity expansion models

Here, this study develops the sparse chronology method to enhance the representative period framework in capacity expansion models, enabling the effective integration of long-duration energy storage modeling. Traditional representative period methods cannot capture the state of charge of seasonal energy storage systems because they do not establish effective inter-day linkages to connect the state of charge between periods. The sparse chronology approach addresses this limitation by establishing inter-day linkages that allow state of charge to shift inter-seasonally. At the same time, it groups identical representative days into partitions, applying constraints sparsely and implicitly to reduce computational load further. Validation results demonstrate that this method successfully simulates long-duration energy storage patterns, achieving close alignment with a continuous yearly benchmark model, with seasonal trends and state of charge cycles clearly represented. The computational load analysis reveals that the sparse chronology method efficiently applies constraints on maximum and minimum state of charge limits within the representative day framework, eliminating the need for detailed constraints on each individual day. By partitioning representative days and constraining only the start and end of each partition, the method significantly decreases computational requirements. Simulation results show that sparse chronology closely approximates the continuous yearly method's accuracy, even with as few as 20 representative days, achieving correlation values with the benchmark of nearly 0.9 in state of charge plots. Furthermore, it maintains computational efficiency, requiring only 4 % of the solver time compared to the continuous yearly method with 20 representative days. This approach allows capacity expansion models to incorporate long-duration energy storage with high temporal, spatial, and technological resolution, enabling more detailed modeling for large-scale power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Boost of the Bio-memristor Performance for Artificial Electronic Synapses by Surface Reconstruction

Biomaterial-based memristors (bio-memristors) are often adopted to emulate biological synapse functions and applied to construct neural computing networks in brain-inspired chip systems. However, the randomness of conductive filament formation in bio-memristors inhibits their switching performance by causing the dispersion of the device-switching parameters. In this case, a facile porous silk fibroin (p-SF) memristor was obtained through a protein surface reconstruction strategy, in which the size of the hole can be adjusted by the density of hybrid nanoseeds. Here, the porous SF memristors exhibit greatly enhanced electrical characteristics, including uniform I-V cycles, centralized distribution of the switching voltages, and both high and low resistances, compared to devices without pores. The results of three-dimensional (3D) simulations based on classical density functional theory (cDFT) suggest that the reconstructed pores in the SF layers guide the formation and fracture of Ag filaments under an electric field and enhance the overall conductivity by separating Ag + ion and electron diffusion pathways. Ag + ions are predicted to preferentially diffuse through pores, whereas electrons diffuse through the SF network. Interestingly, the device conductance can be bidirectionally modulated gradually by positive and negative voltages, can faithfully simulate short-term and long-term plasticity, and can even realize the triplet-spike-timing-dependent plasticity (triplet-STDP) rule, which can be used for pattern recognition in biological systems. The simulation results reveal that a memristor network of this type has an accuracy of ~95.78% in memory learning and the capability of pattern learning. This work provides a facile technology route to improve the performance of bionic-material memristors.

36 MATERIALS SCIENCE↗