Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability study”

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 1,621 records · Page 90

Impacts of variable renewable energy on wholesale markets and generating assets in the United States: A review of expectations and evidence

Here, we synthesize available literature, data, and analysis on the degree to which growth in variable renewable energy (VRE) has impacted or might in the future impact bulk power system assets, pricing, and costs in the United States. Most studies of future scenarios indicate that VRE reduces wholesale energy prices and capacity factors of thermal generators. Traditional baseload generators are more exposed to these changing market conditions than low-capital cost and more flexible intermediate and peak-load generators. From analysis of historical data we find that VRE is already influencing the bulk power market through changes in temporal and geographic patterns areas with higher levels of VRE. The most significant observed impacts have concentrated in areas with significant VRE and/or nuclear generation along with limited transmission, with negative pricing also often occurring during periods with lower system-wide load. So far, however, VRE, has had a relatively modest impact on historical average annual wholesale prices across entire market regions, at least in comparison to other drivers. The reduction of natural gas prices is the primary contributor to the decline in wholesale prices since 2008. Similarly, VRE impacts on thermal plant retirements have been limited and there is little relationship between the location of recent retirements and VRE penetration levels. Finally, although impacts on wholesale prices have been modest so far, impacts of VRE on the electricity market will be more significant under higher VRE penetrations.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Modeling hydrogen markets: Energy system model development status and decarbonization scenario results

Hydrogen can be used as an energy carrier and chemical feedstock to reduce greenhouse gas emissions, especially in difficult-to-decarbonize markets such as medium- and heavy-duty vehicles, aviation and maritime, iron and steel, and the production of fuels and chemicals. Significant literature has been accumulated on engineering-based assessments of various hydrogen technologies, and real-world projects are validating technology performance at larger scales and for low-carbon supply chains. While energy system models continue to be updated to track this progress, many are currently limited in their representation of hydrogen, and as a group they tend to generate highly variable results under decarbonization constraints. Here, the present work provides insights into the development status and decarbonization scenario results of 15 energy system models participating in study 37 of the Stanford Energy Modeling Forum (EMF37), focusing on the U.S. energy system. The models and scenario results vary widely in multiple respects: hydrogen technology representation, scope and type of hydrogen end-use markets, relative optimism of hydrogen technology input assumptions, and market uptake results reported for 2050 under various decarbonization assumptions. Most models report hydrogen market uptake increasing with decarbonization constraints, though some models report high carbon prices being required to achieve these increases and some find hydrogen does not compete well when assuming optimistic assumptions for all advanced decarbonization technologies. Across various scenarios, hydrogen market success tends to have an inverse relationship to success with direct air capture (DAC) and carbon capture and storage (CCS) technologies. While most model-scenario combinations predict modest hydrogen uptake by 2050 – <10 million metric tons (MMT) – aggregating the top 10 % of market uptake results across sectors suggests an upper range demand potential of 42–223 MMT. The high degree of variability across both modeling methods and market uptake results suggests that increased harmonization of both input assumptions and subsector competition scope would lead to more consistent results across energy system models. The wide variability in results indicates strongly divergent conclusions on the role of hydrogen in a decarbonized energy future.

08 HYDROGEN↗

Leveraging Structured Biological Knowledge for Counterfactual Inference: A Case Study of Viral Pathogenesis

Counterfactual inference is a useful tool for comparing outcomes of interventions on complex systems. It requires us to represent the system in form of a structural causal model, complete with a causal diagram, probabilistic assumptions on exogenous variables, and functional assignments. Specifying such models can be extremely difficult in practice. The process requires substantial domain expertise, and does not scale easily to large systems, multiple systems, or novel system modifications. At the same time, many application domains, such as molecular biology, are rich in structured causal knowledge that is qualitative in nature. This manuscript proposes a general approach for querying a causal knowledge graph with a causal question and converting the qualitative result into a quantitative structural causal model that can learn from data to answer the question. Here, we demonstrate the feasibility, accuracy and versatility of this approach using two case studies in systems biology. The first demonstrates the appropriateness of the underlying assumptions and the accuracy of the results. The second demonstrates the versatility of the approach by querying a knowledge base for the molecular determinants of a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-induced cytokine storm and performing counterfactual inference to predict the causal effect of medical countermeasures for severely ill COVID-19 patients.

60 APPLIED LIFE SCIENCES↗

The influence of variability in meridional over turning on global ocean circulation

The Atlantic Meridional Overturning circulation (AMOC) is an important global scale circulation and changes in AMOC can induce significant regional and global climate impacts. In this work, we study the stability of AMOC and its influence on global ocean circulation and the surface climate though analyzing a set of sensitivity experiments using the Community Climate System Model version 1 (CESM1). Results show that a collapsed AMOC can induce changes in global ocean circulation, such as reduced (or reversed) Bering Strait transport, weakened Indonesia throughflow and Agulhas current, but strengthened Drake Passage transport. It also changes the global wind pattern and surface temperature, such as a seesaw-like surface temperature change between Northern and Southern Hemispheres, a weakening of Indian-Australia summer monsoon, a southward shift of the Southern ocean westerlies, etc. We also found that AMOC and the Pacific deep meridional overturning circulation (or PMOC) do not form a natural seesaw under modern day climate and geography. A collapsed AMOC (active PMOC) is not stable under modern conditions if there is no additional freshwater (salt) input in the subpolar North Atlantic (Pacific), suggesting that the modern mean state of AMOC (PMOC) does not depend on local haline forcing although its variability and changes do.

54 ENVIRONMENTAL SCIENCES↗

Multiscale characterization and representation of variability in ceramic matrix composites

Low density, high strength, and high creep and oxidation resistance properties of ceramic matrix composites (CMCs) make them an ideal choice for use in extreme environments in space and military applications. This paper presents a detailed characterization study of structural and manufacturing flaws in Carbon fiber Silicon-Carbide-Nitride matrix (C/SiNC) CMCs at different length-scales. Energy-dispersive spectroscopy (EDS) is used for the chemical characterization of the material’s elemental constituents. High-resolution multiscale graphs obtained from scanning electron microscope (SEM) and confocal laser scanning microscope (LSM) are used to characterize the distribution and morphology of defects at different length scales. This is followed by the classification and quantification of the common manufacturing defects. An image processing algorithm based on the image segmentation process is developed to quantify the variability of various scale-dependent architectural parameters. Finally, a three-dimensional stochastic representative volume element (SRVE) generation algorithm is developed to provide precise representations of material textures at multiple length scales. The developed algorithm accurately accounts for material features and flaws based on a range of multiscale structural and defects characterization results.

36 MATERIALS SCIENCE↗

Probabilistic Forecast of All-Sky Solar Radiation Using Enhanced WRF-Solar: Preprint

This study presents enhancements of the Weather Research and Forecasting model with solar extensions (WRF-Solar) to provide probabilistic forecasts of solar radiation. Our approach builds ensemble WRF-Solar runs by introducing stochastic perturbations of variables that produce the largest uncertainties in predicting surface irradiance and clouds. The key variables are identified using tangent linear sensitivity analysis of six physics packages responsible for all-sky irradiance variability. An optimal strategy to stochastically perturb the selected variables is developed and applied to WRF-Solar to generate ensemble members for day-ahead solar prediction. The National Solar Radiation Database (NSRDB) is used to validate the ensemble forecast at arbitrary locations on the model grid. Preliminary results indicate that the proposed technique can potentially produce WRF-Solar ensembles providing reliable information of solar prediction uncertainty. This study describes the implemented methodology and initial results as well as future research to improve ensemble-based probabilistic forecasts with WRF-Solar.

all-sky↗

Generic and ML Workloads in an HPC Datacenter: Node Energy, Job Failures, and Node-Job Analysis

HPC datacenters offer a backbone to the modern digital society. Increasingly, they run Machine Learning (ML) jobs next to generic, compute-intensive workloads, supporting science, business, and other decision-making processes. However, understanding how ML jobs impact the operation of HPC datacenters, relative to generic jobs, remains desirable but understudied. In this work, we leverage long-term operational data, collected from a national-scale production HPC datacenter, and statistically compare how ML and generic jobs can impact the performance, failures, resource utilization, and energy consumption of HPC datacenters. Our study provides key insights, e.g., ML-related power usage causes GPU nodes to run into temperature limitations, median/mean runtime and failure rates are higher for ML jobs than for generic jobs, both ML and generic jobs exhibit highly variable arrival processes and resource demands, significant amounts of energy are spent on unsuccessfully terminating jobs, and concurrent jobs tend to terminate in the same state. We open-source our cleaned-up data traces on Zenodo (https://doi. org/10.5281/zenodo.13685426), and provide our analysis toolkit as software hosted on GitHub (https://github.com/atlarge-research/2024-icpads-hpc-workload-characterization). This study offers multiple benefits for data center administrators, who can improve operational efficiency, and for researchers, who can further improve system designs, scheduling techniques, etc.

crossanalysis↗

Decarbonization of the power sector with CCS: Case study in two regions in the U.S. and lessons for Latin America

This work consists of estimates of potential changes in the total systems cost (TSC) of two RTOs of the United States (U.S.), MISO-N and SPP RTO West, under the traditional decarbonization pathway of replacement of fossil-power plants with variable renewable energy (VRE) and the less traditional pathway of retrofitting fossil fuel units with carbon capture and storage (CCS). Although the power mixes of MISO-N and SPP RTO West are particular to these regions, results can apply to other regions, including Latin American countries that are planning to decarbonize their power sectors. This case study serves to highlight lessons on differences in technology costs between these two pathways, as well as the cost associated with decarbonization rates close to 100%.

Pena-Cabra, Ivonne A.↗

Chronological changes in soil biogeochemical properties of the glacier foreland of Midtre Lovénbreen, Svalbard, attributed to soil-forming factors

We report that Glacier forelands provide an excellent opportunity to investigate vegetation succession and soil development along the chronosequence; however, there are few studies on soil biogeochemical changes from environmental factors, aside from time. This study aimed to investigate soil development and biogeochemical changes in the glacier foreland of Midtre Lovenbreen, Svalbard, by considering various factors, including time. Eighteen vegetation and soil variables were measured at 38 different sampling sites of varying soil age, depth, and glacio-fluvial activity. Soil organic matter (SOM) was quantitatively measured, and the compositional changes in SOM were determined following size-density fractionation. In the topsoil, the soil organic carbon (SOC) and total nitrogen (N) content was found to increase along the soil chronosequence and were highly correlated with vegetation-associated variables. These findings suggest that plant-derived material was the main driver of the light fraction of SOM accumulation in the topsoil. The heavy fractions of SOM were composed of microbially transformed organic compounds, eventually contributing to SOM stabilization within short 90-yr deglaciation under harsh climatic conditions. In addition to time, the soil vertical profiles showed that other environmental parameters, also affected the soil biogeochemical properties. The high total phosphorous (P) content and electrical conductivity in the topsoil were attributed to unweathered subglacial materials and a considerable amount of inorganic ions from subglacial meltwater. The high P and magnesium content in the subsoil were attributed to parent materials, while the high sodium and potassium content in the surface soil were a result of sea-salt deposition. Glacio-fluvial runoff hampered ecosystem development by inhibiting vegetation development and SOM accumulation. This study emphasizes the importance of considering various soil-forming factors, including parent/subglacial materials, aeolian deposition, and glacio-fluvial runoff, as well as soil age, to obtain a comprehensive understanding of the ecosystem development in glacier forelands.

glacio-fluvial runoff↗

Impact of environmental variables on the degradation of photovoltaic components and perspectives for the reliability assessment methodology

Backsheet cracking has been a major issue observed in the field; however, standardized qualification tests, such as IEC61215, are inadequate to reliably identify such failures of PV modules due to the lack of the critical weathering factors applied sequentially or in combination, such as those found in the service environments. To address this problem, in this work we investigated the effects of various environmental variables on the degradation and failure behaviors of the polyamide-based backsheet in PV modules retrieved from five different locations, encompassing a variety of climates, including humid subtropical, hot-summer Mediterranean, tropical savanna climate and hot arid. The correlations between the degradation indicators and the weathering variables were further demonstrated by principle components analysis (PCA). We found strong relationships between: carbonyl formation and reflected solar radiation; hydroxyl formation and module temperature; yellowness and NO2 concentration, while no simple correlation could be found between a specific weathering factor and cracking. By introducing additional stress factors to the aged polyamide-based backsheet films with the novel 'fragmentation test', we successfully reproduced the field cracking behaviour. This study has demonstrated that different degradation modes of PV components respond differently to the environmental stresses encountered in service. Thereby, any accelerated laboratory test based on a single set condition or lacking key environmental variables would be inadequate to assess the long-term performance of PV modules and components. A new reliability-based methodology is proposed to quantitatively link laboratory testing with field results for the service life prediction of PV materials.

14 SOLAR ENERGY↗

Discrete versus continuous: Enhancing battery optimization in capacity expansion models

This study compares two battery modeling approaches for capacity expansion models: discrete-duration and continuous-duration formulations. In the discrete approach, battery duration is fixed, and power capacity is optimized. In the continuous approach, both power and energy capacities are decision variables, allowing storage duration to be optimized endogenously. Although both discrete-duration and continuous-duration battery formulations are used in long-term power system planning models, the literature has provided limited direct, systematic comparisons of their implications within a common modeling framework. To address this gap, this study implements both approaches in the Regional Energy Deployment System (ReEDS TM ) capacity expansion model using two resource adequacy methods, across a range of future system conditions, and with varying battery cost projections. Results show continuous-duration and high-resolution discrete approaches produce similar capacity expansion outcomes. The continuous formulation achieves faster runtimes compared to discrete-duration runs with many discrete-duration options. However, the discrete-duration approach allows users to choose to have limited fidelity for storage duration options, which in some cases can outperform the continuous formulation. The continuous formulation has the lowest overall system costs, indicating its ability to fine-tune storage duration to better meet specific system needs. This study's findings provide a side-by-side evaluation of discrete and continuous battery modeling approaches and offer guidance for improving the representation of real-world systems, flexibility, and computational efficiency for representing energy storage in long-term power system planning models.

25 ENERGY STORAGE↗

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗

Performance Testing of an Integrated Magnetic Power Take-Off

A wave energy converter (WEC) and the power take-off (PTO) generator system can be represented by using a mass-spring-dampener model. By incorporating a negative stiffness spring within the PTO, the overall stiffness of the PTO can be lowered allowing the impedance of the PTO to be more closely matched with the WEC. By designing for impedance matching the WEC can greatly enhance its power generation capability. This project has involved the design, fabrication and testing of a new type of linear-stroke length and rotary stroke length adjustable negative stiffness magnetic springs for use within a wave energy converter (WEC). The magnetic springs were studied by using 3-D finite element analysis with the objective of creating high energy density and a long linear stroke length. After the construction and testing of both a proof-of-principle adjustable linear and rotary magnetic spring prototype the rotary (torsional) magnetic spring was selected for scaling up analysis. The selected scaled-up proof-of-principle magnetic spring had a peak torque of 850 Nm with a ±45 degree stroke length. By translating the inner rotor relative to the outer rotor, the stiffness could be adjusted to be either negative or positive stiffness. During this project a magnetic lead screw was also studied, and it was shown that by combining the magnetic lead screw with a linearly translating magnetic spring a very long rotary stroke length could be attained. However, as the magnetic lead screw increased complexity and reduced overall energy density relative to a rotary (torsional) spring, this design approach was not further pursed. Dynamometer testing of the variable stiffness magnetic springs was first completed by Portland State University and following this the magnetic spring performance was independently verified by Sandia National Laboratory (Sandia). A WEC analysis when using a variable stiffness magnetic spring was completed by using the WecOptTool. WecOptTool is an open-source WEC optimization software developed by Sandia that supports efficient power take-off (PTO) and control optimization. Wave condition data from the Oregon PacWave test-site was used in this study. The analysis showed that a WEC with a tunable stiffness value could consistently achieve about 80% of its maximum theoretical power production. The use of a tunable stiffness WEC, rather than zero-stiffness or a constant stiffness WEC was also shown to lead to a smaller maximum PTO force requirement. The variable stiffness magnetic spring was integrated into an experimental WEC developed by Sandia, called a Wave-Bot. Sandia successfully completed water-tank testing of the Wave-Bot at the Navy’s Carderock, Maryland, wave-basin test site. The wave-basin testing helped to experimentally demonstrate the operating capabilities and increased power generation capability of a WEC containing an integrated variable stiffness magnetic spring. A WEC capacity factor analysis was also completed. The capacity factor was defined as the ratio of annual average WEC generator power to the maximum (RMS) generated power. This capacity factor provided a means of identifying the ratio of potential revenue to cost. It was calculated that if a tunable variable stiffness magnetic spring has RMS constraints it could operate with a capacity factor above 30%. This is comparable to a wind turbine’s capacity factor.

16 TIDAL AND WAVE POWER↗

The Role of Electrolyte Composition in Enabling Li Metal‐Iron Fluoride Full‐Cell Batteries

Abstract FeF 3 conversion cathodes, paired with Li metal, are promising for use in next‐generation secondary batteries and offer a remarkable theoretical energy density of 1947 Wh kg −1 compared to 690 Wh kg −1 for LiNi 0.5 Mn 1.5 O 4 ; however, many successful studies on FeF 3 cathodes are performed in cells with a large (>90‐fold) excess of Li that disguises the effects of tested variables on the anode and decreases the practical energy density of the battery. Herein, it is demonstrated that for full‐cell compatibility, the electrolyte must produce both a protective solid‐electrolyte interphase and cathode‐electrolyte interphase and that an electrolyte composed of 1:1.3:3 (m/m) LiFSI, 1,2‐dimethoxyethane, and 1,1,2,2‐tetrafluoroethyl‐2,2,3,3‐tetrafluoropropyl ether fulfills both these requirements. This work demonstrates the importance of verifying electrode level solutions on the full‐cell level when developing new battery chemistries and represents the first full cell demonstration of a Li/FeF 3 cell, with both limited Li and high capacity FeF 3 utilization.

25 ENERGY STORAGE↗

The effect of polymer grafting on the mechanical properties of PEG-grafted cellulose nanocrystals in poly(lactic acid)

Poly(lactic acid) (PLA) is a commercially available bio-based polymer that is a potential alternative to many commodity petrochemical-based polymers. However, PLA's thermomechanical properties limit its use in many applications. Incorporating polymer-grafted cellulose nanocrystals (CNCs) is one potential route to improving these mechanical properties. One key challenge in using these polymer-grafted nanoparticles is to understand which variables associated with polymer grafting are most important for improving composite properties. In this work, poly(ethylene glycol)-grafted CNCs are used to study the effects of polymer grafting density and molecular weight on the properties of PLA composites. All CNC nanofillers are found to reinforce PLA above the glass transition temperature, but non-grafted CNCs and CNCs grafted with short PEG chains (<2 kg mol –1 ) are found to cause significant embrittlement, generally resulting in less than 3% elongation-at-break. By grafting higher molecular weight PEG (10 kg mol –1 ) onto the CNCs at a grafting density where the polymer chains are predicted to be in the semi-dilute polymer brush conformation (~0.1 chains nm –2 ), embrittlement can be avoided.

36 MATERIALS SCIENCE↗

Enhancing methane production of anaerobic digested sewage sludge by advanced wet oxidation & steam explosion pretreatment

Anaerobic digestion (AD) of sewage sludge is associated with lowered carbon conversion efficiency due to the recalcitrance of microbial cells in the Waste Activated Sludge (WAS) from the aeration basins with 50 % of the dry weight still being present after AD. The effect of Advanced Wet Oxidation & Steam Explosion (AWOEx) pretreatment on the methane yield at variable temperatures (165–200 °C), residence time (15–45 min) and oxygen dosage (1%–10% based on VS concentration) was tested in this study. Central composite design (JMP® Pro) was used to design 19 runs of AWOEx pretreatment of dewatered digested sewage sludge (DADSS) at different operating conditions. After pretreatment the DADSS was fed to 19 bioreactors operated in semi-continues mode at 37 °C. The highest average methane yields of 183 mL/g VS and 170 mL/g VS were found for the bioreactors receiving DADSS pretreated at 165 °C with a retention time of 15 min and 10 % O2 (Condition 1), and, 182.5 °C with a retention time of 15 min and 5.5 % O2 (Condition 2). This corresponds to an increase in the methane production of 156.2 % and 140.5 % compared to the methane production from the DADSS without pretreatment.

09 BIOMASS FUELS↗

Rapid macrovoid characterization in membranes prepared via nonsolvent-induced phase separation: A comparison between 2D and 3D techniques

Optimizing the performance of asymmetric membranes prepared via nonsolvent-induced phase separation (NIPS) requires a quantitative understanding of how processing variables influence membrane morphology. Presently, the most useful structural quantification techniques require 3D visualization of the membrane structure and are best suited for studies seeking detailed information on small datasets. This study proposes and validates a rapid and accurate technique for quantifying macroporosity (i.e., D m ), a simple descriptor of sublayer macrovoid content in asymmetric membranes D m . values measured from segmented cross-sectional imaging performed via X-ray computed tomography (XCT) and scanning electron microscopy (SEM) are presented and compared for three asymmetric membranes prepared from commercial polymers. Importantly, analyses of 3D XCT membrane reconstructions reveal that D m is described by a single, centralized mean, which demonstrates that macrovoid content is spatially homogenous. Thus, D m can be approximated from limited sampling of the 2D cross-sectional membrane structure via SEM. A proposed 2D SEM sampling method provides D m estimates within ±6% of corresponding 3D XCT values with 30 independent measurements for the three membranes. Further sensitivity is achieved using complementary descriptors such as macrovoid count density (i.e., C m ). This technique is thus a useful tool for characterizing macroporosity from a broad selection of membrane samples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ink-substrate interactions during 3D printing revealed by time-resolved coherent X-ray scattering

Additive printing techniques are regarded as revolutionary and versatile methods of advanced device manufacturing, stemming from the possibility to pattern materials on a custom-based approach and the potential to create novel microstructures and achieve new functionalities. Despite these advantages, the inherent anisotropy of the printing process is a source of property gradients within the printed materials, often associated with variable and/or poor performance. Up to date, the evolutionary pathways associated with printing have largely remained unaddressed, mainly owing to the difficulty to study the transformations induced in the material during processing. Time-resolved coherent X-ray scattering techniques, such as X-ray photon correlation spectroscopy, enable the in situ study of transient nanoscale and mesoscale states in a large variety of materials, including amorphous ones, by directly accessing the most relevant timescales and length scales of their nanoscale and mesoscale dynamics, self-assembly, and mesostructure evolution. In this work, we conduct in operando studies of continuous-flow direct writing with colloidal inks, focusing on how the ink formulation and ink-substrate interactions affect the processes that determine the macroscopic properties of the printed materials. We find fundamental differences in the ink structural relaxations emerging from the primary colloid properties (monodisperse versus aggregated colloids) and the substrates surface energy and mechanical properties. Our work helps to reveal and quantify the basic science governing the evolution of 3D-printed materials during processing, ultimately improving engineering criteria for the design of printable materials.

36 MATERIALS SCIENCE↗