Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix power”

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 199 records · Page 11

Analysis of Various Polysulfone-nanodiamond Mixed-Matrix Membranes

Historically, coal combustion systems have contributed to the extensive network of energy production across the world. Oxycombustion has been shown to be more energy efficient and can facilitate carbon capture from fossil energy utilization. The challenge for small modular deployment of coal combustion or gasification systems is that current O2 enrichment technologies do not scale downward cost-effectively. Fortunately, membrane processes have been developed to offer a more viable polymer-based technology for this application. The central focus is to deliver a high-performance polymer membrane capable of enriching O2 to 90-95% in oxycombustion systems for small modular coal-based power or gasification plants. To accomplish this, various properties of polysulfone-nanodiamond mixed-matrix membrane are examined to determine the appropriate conditions for applicable use. The polymers are infused with nanoparticles to increase the mechanical functionality of the structure which in turn alters the permeability and selectivity of the membrane. Furthermore, the thermal history can also impact the functionality and characteristics of these polymer systems. This study primarily discusses the structure and function of these assorted membranes in terms of permeability and thermomechanical characteristics.

01 COAL, LIGNITE, AND PEAT↗

Interfacial defect reduction enhances universal power law response in Mo–SiNx granular metals

Granular metals (GMs), consisting of metal nanoparticles separated by an insulating matrix, frequently serve as a platform for fundamental electron transport studies. However, few technologically mature devices incorporating GMs have been realized, in large part because intrinsic defects (e.g., electron trapping sites and metal/insulator interfacial defects) frequently impede electron transport, particularly in GMs that do not contain noble metals. Here, we demonstrate that such defects can be minimized in molybdenum–silicon nitride (Mo–SiN x ) GMs via optimization of the sputter deposition atmosphere. For Mo–SiN x GMs deposited in a mixed Ar/N 2 environment, x-ray photoemission spectroscopy shows a 40%–60% reduction of interfacial Mo-silicide defects compared to Mo–SiN x GMs sputtered in a pure Ar environment. Electron transport measurements confirm the reduced defect density; the dc conductivity improved (decreased) by 10 4 –10 5 and the activation energy for variable-range hopping increased 10×. Since GMs are disordered materials, the GM nanostructure should, theoretically, support a universal power law (UPL) response; in practice, that response is generally overwhelmed by resistive (defective) transport. Here, the defect-minimized Mo–SiN x GMs display a superlinear UPL response, which we quantify as the ratio of the conductivity at 1 MHz to that at dc, Δσ ω ⁠. Remarkably, these GMs display a Δσ ω up to 10 7 , a three-orders-of-magnitude improved response than previously reported for GMs. By enabling high-performance electric transport with a non-noble metal GM, this work represents an important step toward both new fundamental UPL research and scalable, mature GM device applications.

36 MATERIALS SCIENCE↗

HYBRID Modeling Validation and Verification Status Matrix

The HYBRID modeling repository is a premier resource for integrated energy systems modeling. HYBRID models have been developed since 2015 to describe the physical operation of tightly coupled thermal systems including power generators, thermal transport systems, thermal storage, thermal-to-electric conversion systems, and other thermal applications. Due to the increased size of the repository, a concise summary matrix of available models is desired. This matrix will consolidate not only the list of available models but also indicate original information sources, publications that have model examples, and level of validation and verification that exists for the models. Validation and verification (V&V) levels begin from simplified algebraic relationship and advance to dynamic data validation. Moving forward, as new models are contributed to HYBRID, their V&V level will be included, updating this matrix.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Low Cost Glass-Ceramic Matrix Composite Heat Exchanger

As part of ARPA-E’s High Intensity Thermal Exchange through Materials and Manufacturing Processes (HITEMMP) program, this project sought to develop novel heat exchanger (HX) capabilities to enable efficient and power dense power generation cycles. This class of HX comes under the category of ceramic/composite materials with the higher temperature goal in the program of ≥1100 °C inlet temperature operation. The enabling capability of this effort is the use of glass-ceramic matrix composite (GCMC) material which provides the high temperature durability of a ceramic, the flaw tolerance of a composite, a significantly faster and lower cost manufacturing process than conventional matrix CMCs and very low porosity levels < 0.5%. For thin-walled HX structures and the need to minimize leakage, the low porosity differentiator is particularly important. RTRC has prior experience with this material system and in the current project advanced the component design and manufacturing methods into new territory to produce features required for effective heat exchange under high pressures. In this approach, silicon carbide fiber is fabricated into a fiber preform using various textile processes. Graphite tooling is used both during the build-up of the fiber preform (interior tooling) and after the fiber preform has been completed (exterior tooling). This tooling assembly is heated to high temperature in an environment that has been evacuated and backfilled with inert gas. A reservoir of specialty glass is present and once the desired temperature has been reached to achieve the desired glass viscosity, an actuator distributes the glass throughout the fiber preform using passageways which are part of the tooling design in a process known as glass transfer molding. After the tooling has been removed, the composite is heat treated to convert the amorphous glass to a crystalline ceramic, providing improved properties. The project was divided into three phases focusing on the following: 1) 10 kW HX design and coupon-level tube sheet fabrication, 2) 10 kW HX fabrication, 3) 50 kW HX fabrication. During Budget Period 1 (BP1), additional risks were encountered and the need for additional funds was agreed upon by ARPA-E program leadership. Due to a variety of factors, the contract modification required nominally 18 months to execute at which time the HITEMMP program was effectively concluding. Because of this and the time that would be required to perform BP2 tasks, it was decided to conclude the project at the end of BP1. During the design of the 10 kW HX, manufacturing constraints were learned and incorporated, leading to a revised configuration for the fiber preform and HX. Heat exchange and pressure drop predictions also played a role in modifying the original design concept to be a higher aspect ratio shell-and-tube HX, simplifying the manufacturing process and improving the heat exchanger performance. Good gravimetric and volumetric thermal power densities of 11.2 kW/kg and 10,200 kW/m3 for the entire HX were projected that involved thermo-structural Finite Element Analysis to determine the structural mass needed for the high operation pressures of 250 bar cold inlet and 80 bar hot inlet. Fiber preforms using textile processes were produced for multiple headered tube sheets. Additional challenges were encountered during the glass transfer molding step for which solutions were identified, but programmatics did not allow them to be implemented in BP1. While complete HX test articles were not fabricated, the benefits of this GCMC material for a variety of high temperature applications remain.

30 DIRECT ENERGY CONVERSION↗

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the Resistive DCON code [A.H. Glasser, Z.R. Wang and J.-K. Park, Physics of Plasmas, \textbf{23}, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Since the solution of ideal outer region is critical to the asymptotic matching and is challenging to be solved in the toroidal geometry due to the singular power series solution at the resonant surfaces, systematic verification of the outer region $\Delta^\prime$ matrix is made by reproducing the well known analytical $\Delta^{\prime}$ result in [H.P. Furth, P.H. Rutherford and H. Selberg, The Physics of Fluids, \textbf{16}, 1054-1063 (1073)] as well as by making a quantitative benchmark with the PEST3 code [A. Pletzer and R.L. Dewar, J. Plasma Physics, \textbf{45}, 427-451 (1991)]. Finally, the reconstructed numerical solution of resistive plasma response from the toroidal matching matrix is presented. Comparing with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

MHD simulation↗

Wave function methods for canonical ensemble thermal averages in correlated many-fermion systems

We present a wave function representation for the canonical ensemble thermal density matrix by projecting the thermofield double state against the desired number of particles. Furthermore, the resulting canonical thermal state obeys an imaginary-time evolution equation. Starting with the mean-field approximation, where the canonical thermal state becomes an antisymmetrized geminal power (AGP) wave function, we explore two different schemes to add correlation: by number-projecting a correlated grand-canonical thermal state and by adding correlation to the number-projected mean-field state. As benchmark examples, we use number-projected configuration interaction and an AGP-based perturbation theory to study the hydrogen molecule in a minimal basis and the six-site Hubbard model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

The method of solving the linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the resistive DCON code [A. H. Glasser, Z. R. Wang, and J.-K. Park, Phys. Plasmas 23, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Additionally, the solution of the ideal outer region is critical to the asymptotic matching and is challenging in toroidal geometry due to the singular power series solution at the resonant surfaces. Thus, a systematic verification of the outer region Δ' matrix is made by reproducing the well-known analytical Δ' result in Furth et al. [Phys. Fluids 16, 1054–1063 (1073)] and by making a quantitative benchmark with the PEST3 code [A. Pletzer and R. L. Dewar, J. Plasma Phys. 45, 427–451 (1991)]. Finally, the reconstructed numerical solution of the resistive plasma response from the toroidal matching matrix is presented. Compared with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dissecting Tensor Cores via Microbenchmarks: Latency, Throughput and Numeric Behaviors

Tensor Cores have been an important unit to accelerate Fused Matrix Multiplication Accumulation (MMA) in all NVIDIA GPUs since Volta Architecture. To program Tensor Cores, users have to use either legacy wmma APIs or current mma APIs. Legacy wmma APIs are more easy-to-use but can only exploit limited features and power of Tensor Cores. Specifically, wmma APIs support fewer operand shapes and can not leverage the new sparse matrix multiplication feature of the newest Ampere Tensor Cores. However, the performance of current programming interface has not been well explored. Furthermore, the computation numeric behaviors of low-precision floating points (TF32, BF16, and FP16) supported by the newest Ampere Tensor Cores are also mysterious. In this paper, we explore the throughput and latency of current programming APIs. Further, we intuitively study the numeric behaviors of Tensor Cores MMA and profile the intermediate operations including multiplication, addition of inner product, and accumulation. All codes used in this work can be found in https://github.com/sunlex0717/DissectingTensorCores.

97 MATHEMATICS AND COMPUTING↗

Electromagnetic Transient (EMT) Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (T&D Models)

Simulation of high-fidelity models of extreme fast charging (XFC) systems and large-area power grids with many XFCs can be time consuming in traditional simulators. Traditional simulators use a single method of discretization for all the components that results in imposing a large computational burden of inverting a large matrix as well as increased computations related to single method of discretization (that is typically a trapezoidal method). To overcome the problem of simulating large-area power grids with many XFCs, in this paper, advanced numerical simulation algorithms are applied for the first time together to reduce the dimension of matrix inversion. Here, the algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. These algorithms apply multiple discretization algorithms rather than a single discretization algorithm that further reduces the computational burden. The approaches mentioned here have resulted in speed-up of up to 18x in the simulation of a single distribution system with 15 XFCs and of up to 271x in the simulation of a transmission-distribution system with 300 XFCs in multiple distribution feeders with respect to conventional simulators (like power systems computer aided design [PSCAD]).

42 ENGINEERING↗

A universal variational quantum eigensolver for non-Hermitian systems

Abstract Many quantum algorithms are developed to evaluate eigenvalues for Hermitian matrices. However, few practical approach exists for the eigenanalysis of non-Hermintian ones, such as arising from modern power systems. The main difficulty lies in the fact that, as the eigenvector matrix of a general matrix can be non-unitary, solving a general eigenvalue problem is inherently incompatible with existing unitary-gate-based quantum methods. To fill this gap, this paper introduces a Variational Quantum Universal Eigensolver (VQUE), which is deployable on noisy intermediate scale quantum computers. Our new contributions include: (1) The first universal variational quantum algorithm capable of evaluating the eigenvalues of non-Hermitian matrices—Inspired by Schur’s triangularization theory, VQUE unitarizes the eigenvalue problem to a procedure of searching unitary transformation matrices via quantum devices; (2) A Quantum Process Snapshot technique is devised to make VQUE maintain the potential quantum advantage inherited from the original variational quantum eigensolver—With additional $$O(log_{2}{N})$$ O ( l o g 2 N ) quantum gates, this method efficiently identifies whether a unitary operator is triangular with respect to a given basis; (3) Successful deployment and validation of VQUE on a real noisy quantum computer, which demonstrates the algorithm’s feasibility. We also undertake a comprehensive parametric study to validate VQUE’s scalability, generality, and performance in realistic applications.

97 MATHEMATICS AND COMPUTING↗

Advanced Moderation Module for High-Temperature Micro-Reactor Applications

Development and deployment of micro-reactors that provide competitive efficiency, prevailing compactness, and inherent safety is an immediate focus of the U.S. nuclear industry. For thermal neutron micro-reactors operating at elevated temperature for optimized efficiency, such as micro molten-salt reactors (MSRs), heat-pipe reactors, and very-high temperature reactors (VHTRs), high-performance moderator based on metal hydride can enhance the neutron economy and therefore achieve reduced weight and enhanced portability. As unclad metal hydride inevitably decomposes at elevated temperature, an enclosure is required for hydride moderator to deliver desired performance at elevated temperatures. Conventional moderator enclosure solutions based on high-temperature alloys introduce extraneous neutron penalty into the micro-reactor, affecting the neutronic benefits provided by the hydride moderator. Additionally, the compatibility between the high-temperature alloy enclosure and the micro-reactor matrix is also a potential issue. In this report, we disclose an Advanced Moderator Module (AMM) concept enabled by an innovative enclosure solution combining the advantages of refractory metals, ceramic matrix composites (CMCs), and advanced coating technology to serve as hydrogen permeation barrier up to very-high temperatures. A schematic description of the AMM structure is illustrated in Figure A1. The AMM contains a moderating material core made of metal hydride with high thermal stability, such as YH 2-x . The hydride core is enclosed by an H 2 barrier layer coated on a ductile refractory metal liner to minimize hydrogen loss during high-temperature operation. A ceramic matrix composite (CMC) cladding is adopted to provide further structural strength, especially during power transients. Between the CMC cladding and metal liner, an extra diffusion barrier coating is inserted to suppress the chemical interaction at elevated temperatures. Hence, based on a series of innovative material solutions, the AMM is capable of containing the metal hydride core at elevated temperature (>900oC) inside coated and lined CMC envelop with negligible hydrogen loss throughout the microreactor lifetime. The benefit of the AMM technology was assessed based on a comprehensive multi-stage reactor-physics analysis. Advanced moderation based on hydride metals (such as YH2-x) enables reaching optimum moderation with higher fuel content than traditional VHTR technology, which is required to design compact micro-reactor cores targeting long-life operation. The AMM enclosure technology provides lower thermal neutron absorption rates than traditional solutions (based on stainless steel or Mo-based alloy such as TZM), which reduces the fissile enrichment requirements by 6-8% on a TRISO-fueled design based on the EMPIRE core. Finally, combining the hydride moderator with neutron transparent enclosure solutions provides significant potential to boost neutronics performance of micro-reactors in terms of increased core lifetime or reduced size and weight by 30-50% on a micro-reactor based on the Holos Quad technology. The progress and plans for on-going development and demonstration efforts are also discussed in this report. The current demonstration of AMM is focused on the hydrogen diffusion barrier demonstration under thermal cycling, while radiation tolerance demonstration is planned. The next step of the demonstration plan will be the assembly of miniature AMM for high-temperature testing, which can be used as a prototype for future pilot scale demonstration.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Using angular two-point correlations to self-calibrate the photometric redshift distributions of DECaLS DR9

Calibrating the redshift distributions of photometric galaxy samples is essential in weak lensing studies. The self-calibration method combines angular auto- and cross-correlations between galaxies in multiple photometric redshift (photo-z) bins to reconstruct the scattering rates matrix between redshift bins. In this paper, we test a recently proposed self-calibration algorithm using the DECaLS Data Release 9 and investigate to what extent the scattering rates are determined. We first mitigate the spurious angular correlations due to imaging systematics by a machine learning based method. We then improve the algorithm for χ 2 minimization and error estimation. Finally, we solve for the scattering matrices, carry out a series of consistency tests, and find reasonable agreements: (1) finer photo-z bins return a high-resolution scattering matrix, and it is broadly consistent with the low-resolution matrix from wider bins; (2) the scattering matrix from the Northern Galactic Cap is almost identical to that from Southern Galactic Cap; (3) the scattering matrices are in reasonable agreement with those constructed from the power spectrum and the weighted spectroscopic subsample. We also evaluate the impact of cosmic magnification. Although it changes little the diagonal elements of the scattering matrix, it affects the off-diagonals significantly. The scattering matrix also shows some dependence on scale cut of input correlations, which may be related to a known numerical degeneracy between certain scattering pairs. In conclusion, this work demonstrates the feasibility of the self-calibration method in real data and provides a practical alternative to calibrate the redshift distributions of photometric samples.

79 ASTRONOMY AND ASTROPHYSICS↗

Load following power generation and power storage using REP and PEM technology

A hydrogen generation system for generating hydrogen and electrical power includes a power supply, a reformer-electrolyzer-purifier (REP) assembly including at least one fuel cell including an anode and a cathode separated by an electrolyte matrix, at least one low temperature fuel cell, and a hydrogen storage. The at least one fuel cell is configured to receive a reverse voltage supplied by the power supply and generate hydrogen-containing gas in the anode of the at least one fuel cell. The at least one low temperature fuel cell is configured to receive the hydrogen-containing gas output from the REP assembly. The at least one low temperature fuel cell is configured to selectably operate in a power generation mode in which the hydrogen-containing gas is used to generate electrical power and a power storage mode in which the hydrogen-containing gas is pressurized and stored in the hydrogen storage.

Jahnke, Fred C.↗

Preliminary Observations on the Hydrogen Redistribution Feedback in YH-Moderated Monolithic Microreactors.

Nuclear microreactors promise to open new markets for nuclear industry. This is due to their potential competitive cost in non-traditional market segments, e.g., mines, forward military basis, extraterrestrial surfaces and remote areas, and inherently safe characteristics that make them deployable where other power sources are not available or difficult to exploit. Transition metal dihydrides have been considered among the most promising candidates for moderating nuclear microreactors. In particular, yttrium hydride (YH\textsubscript{x}) has been selected for high temperature applications due to its high thermal stability and relatively high hydrogen retention at temperatures exceeding $870^\circ C$. One of the main issues associated with the use of hydrides is that, when exposed to temperature, stress, or concentration gradients, the hydrogen contained in the metallic matrix tends to redistribute and leak from the moderating elements, potentially leading to reactivity losses and power swings. The purpose of this paper is to gain a better physical understanding of the neutronic feedback associated with the hydrogen redistribution in YH\textsubscript{x}, by using Griffin and Bison. This feedback is inherently multiphysics, since the hydrogen distribution is strongly dependent on the moderator temperature spatial distribution. In particular, we want to understand the sign (+/-) of the neutronic feedback, its order of magnitude, and its underlying physical causes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Estimation Matrix Calibration of PMU Data-driven State Estimation Using Neural Network

Linear state estimation (LSE) is a phasor measurement unit (PMU) data-based power system state estimation that incorporates a linear measurement model in rectangular coordinates. Due to the high computational efficiency and high observational time-resolution, LSE can act as a supplementary state estimation in a wide-area monitoring system (WAMS). The performance of LSE is relatively sensitive to noises in measurements. Therefore, the estimation accuracy relies heavily on the accuracy of the estimation matrix, which is directly influenced by the measurement weight matrix. This paper proposes two novel calibration method of the estimation matrix using neural networks. One is based on the minimum absolute network loss (ANL), and the other is based on the minimum average squared network loss (ASNL). Both methods are tested and compared with LSE algorithms on the IEEE 14-bus system

neural network↗

Solving a class of infinite-dimensional tensor eigenvalue problems by translational invariant tensor ring approximations

Here, we examine a method for solving an infinite-dimensional tensor eigenvalue problem Hx = λx, where the infinite-dimensional symmetric matrix H exhibits a translational invariant structure. We provide a formulation of this type of problem from a numerical linear algebra point of view and describe how a power method applied to e -Ht is used to obtain an approximation to the desired eigenvector. This infinite-dimensional eigenvector is represented in a compact way by a translational invariant infinite Tensor Ring (iTR). Low rank approximation is used to keep the cost of subsequent power iterations bounded while preserving the iTR structure of the approximate eigenvector. We show how the averaged Rayleigh quotient of an iTR eigenvector approximation can be efficiently computed and introduce a projected residual to monitor its convergence. In the numerical examples, we illustrate that the norm of this projected iTR residual can also be used to automatically modify the time step to ensure accurate and rapid convergence of the power method.

97 MATHEMATICS AND COMPUTING↗

Fabrication and Characterization of Single-Particle Compacts

The programmatic objective of the advanced low enriched uranium (aLEU) fuel Project within the Office of Nonproliferation Research and Development (NA-22) is to advance the development of nuclear reactor designs, fuel materials, and fabrication technologies capable of meeting a number of challenging customer requirements, including reactor stability, extended lifetime, and power density without refueling, while using only low-enriched uranium (i.e., less than 20% enrichment of 235 U). The tandem requirements of delivering sufficient power density and maintaining fuel performance over the course of a multidecade fuel lifetime tightly constrain the field of potential fuel options. Monolithic fuel forms (e.g., fuel pellets) provide maximized density of fissile material, but they may face challenges in maintaining their structure and performance over the long reactor lifetimes envisioned in this program. Dispersion fuel forms with an inert matrix provide enhanced long-term stability but sacrifice fissile material density, and they may not provide sufficient power density over the reactor lifetime with the limit of low enrichment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

On the melt pool dynamic of voxel-controlled metal matrix composites via hybrid additive manufacturing: Laser powder bed fusion and ink-jetting

In this study, the effect of the addition of reinforcement nanoparticles to the 316L matrix by adopting ex-situ and in-situ method (drop on demand jetting) to produce 316L/Al 2 O 3 nanocomposite was investigated. In the ex-situ method, the Al 2 O 3 nanoparticles (NPs) were lightly mixed with 316L powder and processed by laser powder bed fusion. In the in-situ method, an ethanol-based ink containing Al 13 nanoclusters (NCs) was added to 316L powder and then processed by laser. Both ex-situ and in-situ method produced nanocomposites with Al-Si-Mn-O-enriched precipitations within the 316L matrix. The addition of NPs/NCs to the 316L matrix, altered the geometrical characteristic of the single-track melt pools. At the same laser power, with increasing the amount of Al 2 O 3 NPs and Al 13 NCs the melt pool deepened due to reduced thermal conductivity and prolonged liquid presence. Further, as a result, 316L/1 wt% Al 13 NCs deposited single track showed larger grains in comparison to 316L single track. At a high laser power of 150W, the Marangoni flow and the buoyancy force caused the nanoparticles to agglomerate and float to the top surface of tracks; therefore, the wt% fraction of precipitation was drastically reduced due to the loss of Al. The 316L/Al 2 O 3 NPs and 316L/Al 13 NCs exhibited the microhardness of 285 ± 13 HV and 293 ± 7 HV, respectively, higher than the deposited 316L single track, 265 ± 15 HV. Lastly, a hybrid LPBF+ink-jet printer was adopted to selectively change the composition of different zones by adding Al 13 NCs ink to 316L and producing a voxel-controlled metal matrix composite.

316L↗