Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unit dynamics”

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 523 records · Page 29

U.S. national water and energy land dataset for integrated multisector dynamics research

Abstract Understanding resource demands and tradeoffs among energy, water, and land socioeconomic sectors requires an explicit consideration of spatial scale. However, incorporation of land dynamics within the energy-water nexus has been limited due inconsistent spatial units of observation from disparate data sources. Herein we describe the development of a National Water and Energy Land Dataset (NWELD) for the conterminous United States. NWELD is a 30-m, 86-layer rasterized dataset depicting the land use of mappable components of the United States energy sector life cycles (and related water used for energy), specifically the extraction, development, production, storage, distribution, and operation of eight renewable and non-renewable technologies. Through geospatial processing and programming, the final products were assembled using four different methodologies, each depending upon the nature and availability of raw data sources. For validation, NWELD provided a relatively accurate portrayal of the spatial extent of energy life cycles yet displayed low measures of association with mainstream land cover and land use datasets, indicating the provision of new land use information for the energy-water nexus.

58 GEOSCIENCES↗

Quintessential cosmological tensions

Several cosmological tensions have emerged in light of recent data, most notably in the inferences of the parameters and . We explore the possibility of alleviating both these tensions simultaneously by means of the Albrecht-Skordis “quintessence” potential. The field can reduce the size of the sound horizon while concurrently suppressing the power in matter density fluctuations before it comes to dominate the energy density budget today. Interestingly, this rich set of dynamics is governed entirely by one free parameter that is of in Planck units. We find that the inferred value of can be increased, while that of can be decreased, both by compared to the case. However, ultimately the model is disfavored by Planck and BAO data alone, compared to the standard model, with a . When including large scale structure and supernova data we find . We note that historically much attention has been focused on preserving the three angular scales , , and to their values. Our work presents an example of how, while doing so indeed maintains a relatively good fit to the cosmic microwave background data for an increased number of ultrarelativistic species, it is a priori insufficient in maintaining such a fit in more general model spaces.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling and Power-Hardware-in-the-Loop Validation of Synchronous Wind: An Inverterless Grid-Forming Wind Power Plant

Grid-forming (GFM) control of Type-3 and Type-4 wind turbine generators (WTGs) has attracted substantial attention in power systems research; however, the limited overcurrent capability of power electronics converters continues to deteriorate the grid strength of the evolving power systems. Synchronous wind, also known as a Type-5 WTG, offers a unique GFM solution to address grid integration and grid strength issues by keeping the grid largely synchronous at very high integration levels of renewable generation. A Type-5 WTG interfaces with the electric grid via a synchronous generator driven by a variable-speed hydraulic torque converter; hence, the wind rotor operates in variable-speed mode for maximum power generation, and the generator shaft remains synchronous to the grid. This paper develops and tests a high-fidelity model of a Type-5 WTG in a power-hardware-in-the-loop (PHIL) testing environment. The PHIL demonstration shows that a Type-5 WTG inherently behaves as a GFM unit and can obtain similar performance in terms of power responses, wind rotor dynamics, and efficiency compared to a Type-3 WTG in high-wind conditions. The developed model provides further insight into how Type-5 WTGs can benefit the smooth transition to power systems with high integration levels of inverter-based resources.

grid strength↗

MVP: a modular viromics pipeline to identify, filter, cluster, annotate, and bin viruses from metagenomes

While numerous computational frameworks and workflows are available for recovering prokaryote and eukaryote genomes from metagenome data, only a limited number of pipelines are designed specifically for viromics analysis. With many viromics tools developed in the last few years alone, it can be challenging for scientists with limited bioinformatics experience to easily recover, evaluate quality, annotate genes, dereplicate, assign taxonomy, and calculate relative abundance and coverage of viral genomes using state-of-the-art methods and standards. Here, we describe Modular Viromics Pipeline (MVP) v.1.0, a user-friendly pipeline written in Python and providing a simple framework to perform standard viromics analyses. MVP combines multiple tools to enable viral genome identification, characterization of genome quality, filtering, clustering, taxonomic and functional annotation, genome binning, and comprehensive summaries of results that can be used for downstream ecological analyses. Overall, MVP provides a standardized and reproducible pipeline for both extensive and robust characterization of viruses from large-scale sequencing data including metagenomes, metatranscriptomes, viromes, and isolate genomes. As a typical use case, we show how the entire MVP pipeline can be applied to a set of 20 metagenomes from wetland sediments using only 10 modules executed via command lines, leading to the identification of 11,656 viral contigs and 8,145 viral operational taxonomic units (vOTUs) displaying a clear beta-diversity pattern. Further, acting as a dynamic wrapper, MVP is designed to continuously incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field of viromics. MVP is available at https://gitlab.com/ccoclet/mvp and as versioned packages in PyPi and Conda.

59 BASIC BIOLOGICAL SCIENCES↗

A spatially regularized detector for emergent/re-emergent disease outbreaks

Early detection of outbreaks caused by emergent pathogens, using epidemiological surveillance data i.e., daily case counts, is difficult. This is because the data tend to be noisy during the early epoch of the outbreak. In contrast, the spread-rate of the disease tends to be well-behaved, as it depends only on the mixing patterns of the population and the characteristics of the pathogen, neither of which behave erratically in space-time. In this report, we explore whether the spread-rate can be used for epidemiological surveillance, conditional on case count data. Estimating the spread-rate from case count data allows us to exploit exogenous information, e.g., incubation period distributions etc., which can considerably smooth out any erratic temporal behavior. Further, epidemiological dynamics are spatially correlated, and if case counts are available for multiple areal units e.g., counties, these correlations could potentially be used to suppress noise in the early epoch data. These exogenous information and structure are not exploited by conventional syndromic surveillance detectors to extract a well-behaved latent variable for monitoring purposes. The technical challenge lies in the estimation of the spread-rate field jointly over a collection of areal units; further, the spread-rate varies over time. We develop a method based on mean-field variational inference to approximately estimate the spread-rate field, using a Gaussian Random Field Model for spatial regularization. The method is tested on the estimation of spread-rate in the thirty-three counties of New Mexico and detect the arrival of the Fall 2020 COVID-19 wave in September 2020. We find that the method is scalable, but underestimates the uncertainty in the estimated spread-rate field. We detect the arrival of the Fall 2020 wave a week ahead of conventional syndromic surveillance algorithms, but our simplistic detection algorithm, based on simple anomaly detection, suffers from a high false positive rate, similar to conventional detectors.

59 BASIC BIOLOGICAL SCIENCES↗

Isolation of Histone from Sorghum Leaf Tissue for Top Down Mass Spectrometry Profiling of Potential Epigenetic Markers

Histones belong to a family of highly conserved proteins in eukaryotes. They pack DNA into nucleosomes as functional units of chromatin. Post-translational modifications (PTMs) of histones, which are highly dynamic and can be added or removed by enzymes, play critical roles in regulating gene expression. In plants, epigenetic factors including histone PTMs are related to their adaptive responses to the environment. Understanding the molecular mechanisms of epigenetic control can bring unprecedented opportunities for engineering solutions to increase the resilience of crops to climate change. Herein, we describe a protocol to isolate the nuclei and purify histones from sorghum leaf tissue. The extracted histones can be analyzed as their intact forms by top-down mass spectrometry (MS) coupled to online reversed-phase (RP) liquid chromatography (LC). Combinations and stoichiometry of multiple PTMs on the same histone proteoform can be readily identified. In addition, histone tail clipping can be detected using the top-down LC-MS workflow thus yielding the global PTM profile of core histones (H4, H2A, H2B, H3). By comparing the PTM profile among samples corresponding to different conditions (e.g. drought vs. control), potential epigenetic marks can be discovered as targets for further characterization using approaches such as chromatin immunoprecipitation – sequencing (ChIP-seq).

59 BASIC BIOLOGICAL SCIENCES↗

Maximum respiration rates in hyporheic zone sediments are primarily constrained by organic carbon concentration and secondarily by organic matter chemistry

Abstract. River corridors are fundamental components of the Earth system, and their biogeochemistry can be heavily influenced by processes in subsurface zones immediately below the riverbed, referred to as the hyporheic zone. Within the hyporheic zone, organic matter (OM) fuels microbial respiration, and OM chemistry heavily influences aerobic and anaerobic biogeochemical processes. The link between OM chemistry and respiration has been hypothesized to be mediated by OM molecular diversity, whereby respiration is predicted to decrease with increasing diversity. Here we test the specific prediction that aerobic respiration rates will decrease with increases in the number of unique organic molecules (i.e., OM molecular richness, as a measure of diversity). We use publicly available data across the United States from crowdsourced samples taken by the Worldwide Hydrobiogeochemical Observation Network for Dynamic River Systems (WHONDRS) consortium. Our continental-scale analyses rejected the hypothesis of a direct limitation of respiration by OM molecular richness. In turn, we found that organic carbon (OC) concentration imposes a primary constraint over hyporheic zone respiration, with additional potential influences of OM richness. We specifically observed respiration rates to decrease nonlinearly with the ratio of OM richness to OC concentration. This relationship took the form of a constraint space with respiration rates in most systems falling below the constraint boundary. A similar, but slightly weaker, constraint boundary was observed when relating respiration rate to the inverse of OC concentration. These results indicate that maximum respiration rates may be governed primarily by OC concentration, with secondary influences from OM richness. Our results also show that other variables often suppress respiration rates below the maximum associated with the richness-to-concentration ratio. An important focus of future research will identify physical (e.g., sediment grain size), chemical (e.g., nutrient concentrations), and/or biological (e.g., microbial biomass) factors that suppress hyporheic zone respiration below the constraint boundaries observed here.

58 GEOSCIENCES↗

NETL Direct Air Capture Center

In 2022 Congress authorized the $25M National Energy Technology Laboratory Direct Air Capture Center. The facility will be specifically targeted at accelerating the commercialization of technologies beyond the conceptual stage which have not yet reached full pilot-scale (TRL 3 to 6). The ability to operate over a wide range of conditions will help the developers to understand how their technologies respond in different climates, from summer to winter and arid to tropical. NETL recognizes that DAC technology is still early in its evolution, and a wide variety of highly varied technologies are being explored. The Center will be designed with substantial flexibility to accommodate the rapidly evolving technological landscape. Testing systems at three scales will be included: lab-scale systems designed to examine the long-term stability of DAC materials, bench-scale module testing systems capable of probing flow dynamics, and small pilot-scale skid rooms able to test prototype DAC units under a wide variety of climate conditions. The NETL DAC Center will feature dedicated engineering, technician, scientific and logistical support for experimental system design, installation, operation of experiments, and interpretation of results for projects including but not limited to kinetics of absorption/desorption, thermal duty, and post-mortem elemental characterization to understand degradation. The NETL DAC Center projects will also have access to dedicated process modeling and analysis to evaluate the technoeconomic aspects of new technologies. NETL’s Systems Engineering and Analysis (SEA) team, the Carbon Capture Simulation Initiative (CCSI), and the Institute for the Design of Advanced Energy Systems (IDAES) are all potential collaborators with NETL’s testing partners from industry, academia, and other research institutions. NETL also has world class capabilities in device scale modeling (i.e., MFIX) that will be available to help design more efficient contactors and optimize unique internal configurations that can be realized by advanced manufacturing. The talk will cover the design and capabilities of the center, the planned availability, and feedback will be solicited from the audience on DAC testing needs.

Luebke, David↗

Accelerating commercialization of direct air capture technology

In 2022 Congress authorized the $25M National Energy Technology Laboratory Direct Air Capture Center. The facility will be specifically targeted at accelerating the commercialization of technologies beyond the conceptual stage which have not yet reached full pilot-scale (TRL 3 to 6). The ability to operate over a wide range of conditions will help the developers to understand how their technologies respond in different climates, from summer to winter and arid to tropical. NETL recognizes that DAC technology is still early in its evolution, and a wide variety of highly varied technologies are being explored. The Center will be designed with substantial flexibility to accommodate the rapidly evolving technological landscape. Testing systems at three scales will be included: lab-scale systems designed to examine the long-term stability of DAC materials, bench-scale module testing systems capable of probing flow dynamics, and small pilot-scale skid rooms able to test prototype DAC units under a wide variety of climate conditions. The NETL DAC Center will feature dedicated engineering, technician, scientific and logistical support for experimental system design, installation, operation of experiments, and interpretation of results for projects including but not limited to kinetics of absorption/desorption, thermal duty, and post-mortem elemental characterization to understand degradation. The NETL DAC Center projects will also have access to dedicated process modeling and analysis to evaluate the technoeconomic aspects of new technologies. NETL’s Systems Engineering and Analysis (SEA) team, the Carbon Capture Simulation Initiative (CCSI), and the Institute for the Design of Advanced Energy Systems (IDAES) are all potential collaborators with NETL’s testing partners from industry, academia, and other research institutions. NETL also has world class capabilities in device scale modeling (i.e., MFIX) that will be available to help design more efficient contactors and optimize unique internal configurations that can be realized by advanced manufacturing. The talk will cover the design and capabilities of the center, the planned availability, and feedback will be solicited from the audience on DAC testing needs.

Luebke, David↗

Modeling and Power-Hardware-in-the-Loop Validation of Synchronous Wind: An Inverterless Grid-Forming Wind Power Plant: Preprint

Grid-forming (GFM) control of Type-3 and Type-4 wind turbine generators has attracted substantial attention in power systems research; however, the limited over-current capability of power electronics converters continues to deteriorate the grid strength of the evolving power systems. Synchronous wind, also known as Type-5 wind turbine generator (WTG), offers a unique GFM solution to address grid integration and grid strength issues by keeping the grid largely synchronous at very high penetration levels of renewable generation. A Type-5 WTG interfaces to the electric grid via a synchronous generator (SG) driven by a variable-speed hydraulic torque converter; hence, the wind rotor operates in variable-speed mode for maximum power generation and the generator shaft remains synchronous to the grid. This paper developed and tested a high-fidelity model of Type-5 WTG under power-hardware-in-the-loop (PHIL) testing environment. The PHIL demonstration showed that a Type-5 WTGs inherently behaves as a GFM unit and can obtain similar performance in terms of power responses, wind rotor dynamics, and efficiency compared to Type-3 WTG in high wind conditions. The developed model also provides further insight on how Type-5 WTGs can benefit the smooth transition to power systems with high integration level of inverter-based resources.

grid strength↗

On the transferability of residence time distributions in two 10-km long river sections with similar hydromorphic units

Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface and their residence time distributions (RTDs) in the subsurface are important for managing the water quality and ecosystem health in dynamic river corridors. However, direct simulating high-spatial resolution HEFs and RTDs can be time-consuming, especially for watershed-scale modeling. Efficient surrogate models linking RTDs to hydromorphic units (HUs) can be alternatives for simulating RTDs in large-scale models. A common concern of these surrogate models, though, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this issue, this work evaluates the HEFs and resulting RTD-HU relationships for two 10-km long river corridors along the Columbia River leveraging a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework we previously developed. Applying such a framework at the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. Finally, our comparison shows that the similarity and transferability of the RTD-HU relationship is very low for the two investigated river sections, which suggests that devising a general algorithm to estimate RTDs based solely on surface water hydrodynamics and short-distance river channel topography data, as well as HU classification, might be nearly impossible.

54 ENVIRONMENTAL SCIENCES↗

Carbon-carbon interactions in warm dense titanium carbide

Encouraged by experimental reports of diamond precipitation in carbon-containing materials under warm dense matter conditions and in shock compression experiments, we examine the behavior of carbon in TiC using density functional theory-molecular dynamics simulations. Two polymorphs of TiC are considered, the ambient-pressure B1 ($Fm$$\overline{3}$$m$) structure in which C-C interactions are prevented by geometric constraints, and a high-pressure Cmcm structure in which C atoms condense into zigzag chains. Chemistry-inspired bonding analyses confirm the covalently bound nature of these chains and further illuminate important interatomic interactions in both structures. Upon melting of B1 TiC, new short-range C-C interactions develop in the liquid, while the pre-existing C-C interactions in Cmcm TiC persist in the liquid. Here, the resulting carbon networks in the melt provide a promising environment for eventual diamond nucleation.

Ab initio molecular dynamics↗

Data driven methods to recognize patterns in EIC weak-strong simulation

Beam-Beam simulations are currently being studied in preparation for future EIC experiments to study beam-beam effects and, in turn, maximize luminosity. Weak-strong methods are studied for single-particle dynamics during collision. 1 million macro-particles for 1 million turns are typically tracked, corresponding to only 10 seconds in the EIC. The goal of this study is to predict beam properties over the scale of hours. A potential solution focuses on using data-driven methods such as machine learning methods to analyze and extend the insights of the beam properties such as long-term nonlinear effects. This would aid in long-term predictions where results would be more efficiently acquired than a typical tracking simulation. Some limitations such as inaccurate predictions and spatial complexity are also discussed. These methods can then be applied to strong-strong simulations in the future studies.

Accelerator Physics↗

Imaging nanoscale carrier, thermal, and structural dynamics with time-resolved and ultrafast electron energy-loss spectroscopy

Time-resolved and ultrafast electron energy-loss spectroscopy (EELS) is an emerging technique for measuring photoexcited carriers, lattice dynamics, and near-fields across femtosecond to microsecond timescales. When performed in either a specialized scanning transmission electron microscope or ultrafast electron microscope (UEM), time-resolved and ultrafast EELS can directly image charge carriers, lattice vibrations, and heat dissipation following photoexcitation or applied bias. Yet, recent advances in theoretical calculations and electron optics are often required to realize the full potential of ultrafast EEL spectrum imaging. Here, in this review, we present a comprehensive overview of the recent progress in the theory and instrumentation of time-resolved and ultrafast EELS. We begin with an introduction to the technique, followed by a physical description of the loss function. We outline approaches for calculating and interpreting ground-state and transient EEL spectra spanning low-loss plasmons to core-level excitations analogous to x-ray absorption. We then survey the current state of time-resolved and ultrafast EELS techniques beyond photon-induced near-field electron microscopy, highlighting abilities to image carrier and thermal dynamics. Finally, we examine future directions enabled by emerging technologies, including electron beam monochromation, in situ and operando cells, laser-free UEM, and high-speed direct electron detectors. These advances position time-resolved and ultrafast EELS as a critical tool for uncovering nanoscale dynamic processes in quantum materials and solar energy conversion devices.

Computational methods↗

A Fast, Accurate Prediction for System-Wide Damage Due to Dynamic Wind Loading

The complex relationship between photovoltaic (PV) hardware configurations, overall system dynamics, and turbulent aerodynamic phenomena generates highly unsteady, non-uniform loads that can lead to damaging instabilities. These effects may result in glass breakage, cell cracking, and structural failures in frames and mounting systems, even under moderate wind conditions. Addressing industry concerns about premature system failures in field conditions deemed survivable, our research aims to develop a fast and accurate predictive model for system damage. This model integrates configurable hardware choices with advanced simulation tools to represent the overall system-specific dynamics effectively. Using this model, we predict responses under varying weather conditions and hardware setups, translating these predictions into pre-trained surrogate models capable of accurately identifying failure risks and rapidly testing new system hardening measures. In this presentation, we will showcase preliminary results in capturing system dynamics through our customizable library of PV hardware configurations. Additionally, we will highlight how these new tools build upon PVade's established wind load modeling capabilities and foster the development of advanced AI/ML surrogates for improving system robustness.

97 MATHEMATICS AND COMPUTING↗

Prediction of dynamic behavior of a system undergoing an exchange of components using a transmission simulator method

In systems with modular design, it is of interest to perform component exchanges to upgrade components and explore new system configurations. If the dynamics of the components are known, then dynamic substructuring can be used to obtain a prediction of system performance for the new configuration. The transmission simulator method is a substructuring technique that utilizes an additional fixture to apply an elastic boundary condition to a test article such that an appropriate dynamic response can be obtained for a given substructuring application, and the technique has been used widely for performing coupling and decoupling operations. Here, in this work, a modified transmission simulator method is proposed that allows a direct exchange of components to be performed. In a numerical demonstration, the proposed method is shown to accurately predict the dynamics of an assembly that has undergone a component exchange in which the exchanged components have different material properties and different geometric properties. This approach can be applied to dynamic qualification frameworks in which structures undergo modifications or upgrades rather than complete redesigns.

Component exchange↗

Impact of Heaters on Molten Salt Reactors Dynamics

Recently, there has been renewed interest in the Molten Salt Reactor (MSR) concept. This interest is mainly due to its advantages over Light Water Reactors (LWRs). Among these advantages are the flexibility in the fuel choice and the possibility to burn actinides [1]. Currently, two molten research reactors are under development. A 2 MW molten salt experimental reactor has been constructed and is planned for operation in China. A 1 MW Molten Salt Research Reactor is under construction in the United States [2]. Thus, Research and Development (R&D) programs are needed for the MSR technology. Reactor dynamics studies are crucial safety evaluation studies. The development of reactor dynamics tools for MSRs started during the Molten Salt Reactor Program (MSRP) at ORNL, focusing on the Molten Salt Reactor Experiment (MSRE) [3]. Although various dynamic simulation tools are available, R&D efforts to develop specialized tools for MSRs are still ongoing [4]. The primary motivation for this interest is the goal of commercializing MSRs. Developing a dynamics tool for MSRs is crucial for analyzing their safety and supporting their demonstration efforts. A crucial safety issue in MSRs is ensuring that the fuel salt in the primary loop remains molten. This concern becomes particularly important when the reactor operates at lower power levels. Most chloride and fluoride salts used as fuel have a melting point of approximately 450 °C [5]. Therefore, maintaining fuel salt temperatures above this melting point is essential to prevent it from freezing. Freezing of the fuel salt can lead to volume expansion, potentially causing damage to reactor components [6]. Moreover, it can create local blockages within fuel salt channels, reducing cooling efficiency and resulting in hot spots [7]. Various methods can prevent fuel salt from freezing at low power levels. These include reducing the fuel mass flow rate in the primary loop, decreasing heat exchange between the primary and secondary loops in the reactor, and using electric fuel salt heaters. Reducing the mass flow rate of fuel salt in the primary loop can increase the fuel salt transit time in the core, enhancing heat production and increasing temperature. However, it can potentially lead to challenges. These include increasing the potential for fuel salt deposition on the reactor channel walls [8], affecting the overall heat transfer processes. Thus, alternative methods should be used [9]. This research investigates the effects of fuel salt heaters on the dynamics and stability during MSR operation. The investigation includes studying the reactor stability with the fuel salt heaters on and off. In addition, the reactor response to transients is compared with and without the operation of fuel salt heaters. These transients include reactivity insertion, primary and secondary pump failure, and heat sink temperature increase. The results of the reactivity insertion transient are presented in the summary.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗