Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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,117 records · Page 62

Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale (Final Report)

This project titled “Use of Stable Mercury Isotopes to Assess Mercury and Methylmercury Transformation and Transport across Critical Interfaces from the Molecular to the Watershed Scale” represents a collaborative effort between the University of Michigan (Jason Demers, PI) and Oak Ridge National Laboratory (Scott Brooks, co-I). Much has been learned about mercury (Hg) cycling in stream ecosystems, and East Fork Poplar Creek (EFPC) in particular, through decades of previous research. Nevertheless, some of the most fundamental questions regarding the sources of bioavailable Hg and its transformation to toxic methylmercury (MeHg) have remained unanswered. These fundamental questions include: (1) what are the sources and biogeochemical processes that lead to the input of dissolved Hg to stream water across critical subsurface interfaces within stream ecosystems, and EFPC in particular? and (2) what are the sources and biogeochemical processes that control the production and fate of bioaccumulative MeHg within stream ecosystems, and in EFPC in particular? To address these fundamental questions, our project aimed to couple laboratory experiments and field observations, both utilizing natural abundance Hg stable isotope techniques, to identify the processes responsible for generating mobile, bioavailable dissolved Hg from recalcitrant legacy sources within critical subsurface zones (e.g., streambed hyporheic zone, riparian floodplain subsurface). We used the isotopic signature of this bioavailable dissolved Hg to track its mobilization across these critical interfaces in order to link diffuse subsurface sources of dissolved Hg with increases in surface water dissolved Hg flux measured at the watershed scale. Additionally, our research aimed to determine the isotopic composition of MeHg within these same critical subsurface zones. We directly assessed the isotopic composition of MeHg within biota in order to gain insight into which subsurface sources of inorganic Hg and toxic MeHg are available for bioaccumulation within the EFPC ecosystem. Net fluxes of dissolved Hg along the flow path of EFPC were shown to vary spatially and temporally. In the Upper EFPC, within the Y12 boundary, stream water flux of dissolved Hg consistently decreased between the outfall and the downstream boundary of Y12 (57% ± 29%, 1SD). Within the Upper EFPC, an assessment of Hg isotopic composition suggested that losses were strongly reaction-driven, although isotopic diagnostics did not conform to any known processes. Downstream of Y12, in the upper reach of the Lower EFPC, dissolved Hg fluxes tended to increase during the dormant season (net gain of 11-120%), and decrease during the growing season (net loss of 23% +/- 18%, 1SD). In the downstream-most reach of Lower EFPC, dissolved Hg fluxes increased by 12-108% in 9 out of 10 monthly assessments. Overall, diffuse fluxes from the non-Y12 watershed accounted for 34% (+/- 17%, 1SD) of all dissolved Hg exported during base flow. Within Lower EFPC, an assessment of Hg isotopic composition was consistent with the contribution of diffuse Hg inputs from high-concentration hotspots within riparian floodplains and streambed hyporheic pore water. To investigate remobilization of recalcitrant Hg from legacy sediment sources, we developed procedures that coupled isotopic analysis with sequential extractions of streambed sediment. We found that the proportion of weakly-bound Hg within EFPC streambed sediment was relatively small, but could still account for a large proportion of the annual flux of dissolved Hg from EFPC. These sequential extractions also showed that this weakly-bound Hg fraction could be replenished from the much larger fraction of recalcitrant Hg in sediment. The isotopic composition of these weakly-bound and remobilized recalcitrant Hg fractions within the sediment was consistent with high-concentration dissolved Hg hotspots within hyporheic pore water. Thus, this research provided novel evidence that legacy mercury sources within streambed sediment could provide an ongoing contribution of dissolved Hg to surface waters. Finally, we developed new methods for the direct determination of the MeHg isotopic composition of organisms, which allowed a more direct evaluation of inorganic Hg and MeHg sources accumulating in the food web. We found that fish and aquatic invertebrates in both EFPC and a regional background site obtained inorganic Hg and MeHg from multiple isotopically distinct sources, including sediment, suspended particulates, and periphyton. Photodemethylation was found to be an important reaction influencing MeHg dynamics at both sites. However, the balance of microbial methylation and demethylation processes differed between the two streams, with fractionation resulting from methylation and demethylation processes being relatively in balance within the regional background site, whereas microbial methylation appeared to be dominant over microbial demethylation within the EFPC ecosystem. Broadly, the application of Hg isotopic analysis in this study led to numerous novel insights regarding the biogeochemical cycling of Hg in stream ecosystems. This research project promoted the development of two new approaches, including the coupling of sequential extractions with Hg isotopic analysis to assess remobilization of recalcitrant Hg within sediments, and a new method for the isotopic analysis of MeHg isolated from environmental samples. Both of these efforts represent advances in capacity for the field of mercury isotopic analysis and environmental assessment. Overall, this study demonstrates that the application of Hg stable isotope techniques continues to provide new insights into the biogeochemical cycling of Hg in complex aquatic environments, both within the EFPC and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 1—Theory

Summary Core-scale measurements are considered the ground truth that oil and gas or electric utility companies use to predict the migration of fluids such as oil, natural gas, carbon dioxide, or brine deep underground during their extraction or injection operations. To provide a greater understanding of petrophysical properties of low-permeability geologic formations such as shales and tight gas sandstones, this study introduces a novel core-analysis procedure. The technique follows conventional pressure-pulse-decay permeametry, where the pressure in an inlet chamber adjacent to a cylindrical core plug undergoes a rapid pressurization, the system is shut in, and the pressure reaches a new equilibrium. However, unlike a standard unidirectional pressure-pulse decay, the full-immersion pressure-pulse decay (Hannon 2019) applies a pressure disturbance to the entire outer surface area of the sample. This article covers the numerical simulator designed to model flow through an anisotropic porous sample in this scenario. The model assumes distinct but uniform permeabilities along the radial and axial directions of the cylindrical plug sample. When extracting a plug vertically (or perpendicular to bedding), the permeability along the radial direction associates with the horizontal permeability (i.e., parallel to bedding), whereas in the axial direction, flow occurs perpendicular to bedding (a vertical permeability). Investigations of these model outputs demonstrate an approximately 20-fold decrease in time to complete a full-immersion experiment compared with conventional pressure-pulse decay. Furthermore, the pressure-decay curves resulting from the full-immersion method have slightly different shapes than those resulting from other unidimensional transient methods. These differences begin to demonstrate that, under achievable experimental conditions, the analysis of pressure data from one full-immersion test could enable the simultaneous estimation of the apparent permeabilities parallel and perpendicular to bedding of a cylindrical sample in addition to its porosity. A follow-up article finalizes the proof of this capability with a parameter-estimation procedure and presents experimental verification through a proof-of-concept study. Described in greater detail in that article, the parameter estimator requires multiple forward simulations to analyze the data, which behooves minimizing the compute time for each simulation. By using an alternating direction implicit time-marching scheme and a structured but variably spaced grid, the numerical simulator built for this purpose provides a forward model output with suitable accuracy in approximately 0.5 seconds.

Energy & Fuels↗

Spatial Mode Correction of Single Photons Using Machine Learning

Abstract Spatial modes of light constitute valuable resources for a variety of quantum technologies ranging from quantum communication and quantum imaging to remote sensing. Nevertheless, their vulnerabilities to phase distortions, induced by random media, impose significant limitations on the realistic implementation of numerous quantum‐photonic technologies. Unfortunately, this problem is exacerbated at the single‐photon level. Over the last two decades, this challenging problem has been tackled through conventional schemes that utilize optical nonlinearities, quantum correlations, and adaptive optics. In this article, the self‐learning and self‐evolving features of artificial neural networks are exploited to correct the complex spatial profile of distorted Laguerre–Gaussian modes at the single‐photon level. Furthermore, the potential of this technique is used to improve the channel capacity of an optical communication protocol that relies on structured single photons. The results have important implications for real‐time turbulence correction of structured photons and single‐photon images.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamic Strength of Copper at High Pressures Using Pressure Shear Plate Experiments

Measurement of strength at high pressures is critical to understanding the behavior of materials subjected to extreme loading conditions. Recent developments have extended the pressure shear plate impact (PSPI) technique to the high-pressure regime. Such experiments provide a unique opportunity to extract the complete stress-strain behavior of materials at relatively high pressures. The modified technique includes a new fiber-optic heterodyne transverse velocity interferometer system and new analysis methods to account for the inelastic response of the anvil plates. In this study, PSPI experiments are conducted on oxygen-free high conductivity (OFHC) copper at pressures ranging from 10 to 43 GPa and at strain rates of ~10 5 s -1 . Complete stress-strain curves of copper are obtained using a hybrid methodology that relies on numerical simulations and particle velocity records gathered in these experiments. Experiments are conducted at similar pressures and strain rates using two types of anvils that undergo different levels of plastic strains: tool steel and tungsten carbide, to check the robustness of the hybrid method. The behavior of copper at these high pressures and strain rates is compared with previous literature data, and the effect of pressure and strain rate on the strength behavior of copper are discussed. Here, it was observed that the rate of change in yield strength change in pressure was 2.5 times the rate of change of shear modulus with pressure. The strength models that predict the response of copper at these extreme conditions are also examined and validated.

36 MATERIALS SCIENCE↗

Reactive flash sintering and characterization of bulk high entropy nitrides

Over the past decade, numerous high-entropy ceramics have been synthesized, often displaying attractive properties. However, the study on facile preparation of bulk high entropy nitrides (HEN) are limited, despite its broad potential applications. This research demonstrates for the first time rapid fabrication (within ∼6 min) of bulk high-entropy nitrides, especially (Al 0.17 Nb 0.17 Ta 0.17 Ti 0.32 Zr 0.17 )N, from binary nitride powder mixtures using a highly efficient reactive flash sintering (RFS) technique. X-ray diffraction (XRD) shows the HENs from RFS are near single-phase solid solutions with a rock salt crystal structure, while in situ synchrotron study carried out during RFS captured in real time the formation of HEN, which was preserved upon cooling, suggesting thermodynamic stability of the HEN phase, even up to extreme pressure (∼35.6 GPa). Microscopic analyses using SEM, STEM, and EDS reveal decent uniformity for HEN with no obvious segregation of elements, even to submicron scale. Some properties of the obtained bulk HENs are consistent with expectations. For example, their hardness and bulk modulus are close to estimates based on rule-of-mixture (ROM) values from the constituent binary nitrides. Meanwhile, some other measured properties seem to show surprises. For example, the fracture toughness for the HENs (e.g., 7.81 ± 1.40 MPa•m 1/2 or higher) turns out to be more than double of the expected ROM estimates. The significantly improved fracture toughness is attributed to the observed nano-layered structure of the HENs, despite the HEN’s cubic crystal structure and high hardness. In addition, the oxidation resistance shows improvement up till ∼800°C, possibly due to Ta doping that suppress oxygen vacancy formation in the oxide shell, while the 5-metal HEN of (Al 0.17 Nb 0.17 Ta 0.17 Ti 0.32 Zr 0.17 )N displays superconductivity (T c of ∼5–7 K from magnetism and resistivity measurements, slightly lower than ROM estimate), despite insulating property of starting AlN. Furthermore, future study combining experimental investigation using larger samples to confirm the observed increase in fracture toughness and oxidation resistance, theoretical modeling at different length scale, and more detailed structural/chemical characterization, especially at the atomic scale, are all needed to fully understand the inter-relationships between composition, processing, structure, and novel properties for these HENs and the development of related new materials for different applications.

Flash sintering↗

Multi-fidelity physics-informed machine learning for probabilistic damage diagnosis

Machine learning (ML) models are gaining popularity in structural health monitoring (SHM) because of their ability to learn the complex relationship between damage and sensor data. However, the lack of sufficient experimental data for structures with different degrees of damage is a key problem in training ML models for SHM. This problem can be alleviated by using physics-based models to generate the required training data to build physics-informed ML (PIML) models for SHM. However, it takes significant computational effort to perform enough high-fidelity simulations of the diagnostic test. It is thus desirable to know whether the available computational resource budget should be expended on numerous low-fidelity physics simulations, or a small number of high-fidelity simulations, or their combination. In this paper, we investigate this aspect of generating adequate training data for PIML, by constructing multi-fidelity PIML models. We evaluate the performance of several PIML models, trained with different amounts of low-fidelity and high-fidelity data, in locating hidden cracks in concrete structures using a nonlinear dynamics-based diagnosis technique. Here, we find that high-fidelity physics simulations that do not cover the (test and damage) parameter space do not improve the performance of diagnostic PIML models built using data from many low-fidelity physics simulations.

42 ENGINEERING↗

Nanoscale wetting controls reactive Pd ensembles in synthesis of dilute PdAu alloy catalysts

The performance of bimetallic dilute alloy catalysts is largely determined by the size of minority metal ensembles on the nanoparticle surface. By analyzing the synthesis of catalysts comprising Pd 8 Au 92 nanoparticles supported on silica using surface-sensitive techniques, we report that whether Pd overgrowth occurs before or after Au nanoparticle deposition onto the support controls the surface Pd ensemble size and abundance. These differences in Pd ensembles influence catalytic reactivity in H 2 –D 2 isotope exchange and benzaldehyde hydrogenation, which, in correlation with theoretical calculations, is used to elucidate the active site(s) in each reaction. To clarify how the synthetic sequence controls the formation of Pd ensembles, we combine numerical wetting calculations and molecular dynamics simulations (with a machine-learned force field) to visualize Pd deposition and migration on the nanoparticle surface, respectively. Our results suggest that the nanoparticle–support interface restricts nanoparticle accessibility to Pd deposition, which consequently controls the Pd ensemble size, illustrating the critical role of nanoscale wetting phenomena during bimetallic catalyst preparation.

36 MATERIALS SCIENCE↗

Growth and Saturation of the Electron Drift Instability in a Crossed Field Plasma

The linear growth and nonlinear energy transfer of the electron drift instability (EDI) are experimentally measured in the plume of a low-temperature, Hall effect discharge. A frequency-based bispectral analysis technique applied to fast ion density fluctuation measurements shows a growth rate function that is qualitatively similar to predictions from the linear instability dispersion relation, but an order of magnitude smaller. Calculation of the nonlinear transfer function indicates multiple three-wave interactions between high-frequency resonances of the instability in addition to an inverse energy cascade toward lower-frequency modes. Furthermore, these results are discussed in the context of recent theoretical, numerical, and experimental efforts on the EDI in Hall effect discharges and how the EDI may impact anomalous cross field transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum Circuit Cutting for Classical Shadows

Classical shadow tomography is a sample-efficient technique for characterizing quantum systems and predicting many of their properties. Circuit cutting is a technique for dividing large quantum circuits into smaller fragments that can be executed more robustly using fewer quantum resources. We introduce a divide-and-conquer circuit cutting method for estimating the expectation values of observables using classical shadows. We derive a general formula for making predictions using the classical shadows of circuit fragments from arbitrarily cut circuits and provide the sample complexity analysis for the case when observables factorize across fragments. Then, we numerically show that our divide-and-conquer method outperforms traditional uncut shadow tomography when estimating high-weight observables that act non-trivially on many qubits and discuss the mechanisms for this advantage.

97 MATHEMATICS AND COMPUTING↗

A Provably Accurate Randomized Sampling Algorithm for Logistic Regression

In statistics and machine learning, logistic regression is a widely-used supervised learning technique primarily employed for binary classification tasks. When the number of observations greatly exceeds the number of predictor variables, we present a simple, randomized sampling-based algorithm for logistic regression problem that guarantees high-quality approximations to both the estimated probabilities and the overall discrepancy of the model. Our analysis builds upon two simple structural conditions that boil down to randomized matrix multiplication, a fundamental and well-understood primitive of randomized numerical linear algebra. We analyze the properties of estimated probabilities of logistic regression when leverage scores are used to sample observations, and prove that accurate approximations can be achieved with a sample whose size is much smaller than the total number of observations. To further validate our theoretical findings, we conduct comprehensive empirical evaluations. Overall, our work sheds light on the potential of using randomized sampling approaches to efficiently approximate the estimated probabilities in logistic regression, offering a practical and computationally efficient solution for large-scale datasets.

Chowdhury, Agniva↗

An exact inversion method for extracting orientation ordering by small-angle scattering

Here, we outline a nonparametric inversion strategy for determining the orientation distribution function (ODF) of sheared interacting rods using small-angle scattering techniques. With the presence of direct inter-rod interaction and fluid mechanical forces, the scattering spectra are no longer characterized by the azimuthal symmetry in the coordinates defined by the principal directions of simple shear conditions, which severely compounds the reconstruction of ODFs based on currently available methods developed for dilute systems. Using a real spherical harmonic expansion scheme, the real-space ODFs are uniquely determined from the anisotropic scattering spectra and their numerical accuracy is verified computationally. Our method can be generalized to extract ODFs of uniaxially anisotropic objects under different flow conditions in a properly transformed reference frame with suitable basis vectors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Applications of atomic layer deposition and chemical vapor deposition for perovskite solar cells

Metal halide perovskite solar cells (PSCs) have rapidly evolved over the past decade to become a photovoltaic technology on the cusp of commercialization. In the process, numerous fabrication strategies have been explored with the goal of simultaneously optimizing for device efficiency, stability, and scalability. Atomic layer deposition (ALD) and chemical vapor deposition (CVD) have proven to be effective tools for the fabrication of various components of PSCs. This review article examines the application of ALD and CVD for the deposition and modification of charge transport layers, passivation layers, absorber materials, encapsulants, and electrodes. It outlines the use of these vapor deposition techniques in state-of-the-art, multi-junction solar cell devices, and also contains a discussion of the stability of metal halide perovskite materials under ALD and CVD conditions based on in situ characterization reported in literature. Finally, this article concludes with insights into future ALD and CVD research directions that could be undertaken to further aid the deployment of PSCs in emerging solar photovoltaic markets.

14 SOLAR ENERGY↗

An extended moments model of quantum efficiency for metals and semiconductors

The complexity of photocathode designs and detector materials, and the need to model their performance for short pulse durations, response to high frequency photons, presence of coatings and/or thinness of the absorptive layer, necessitates modifications to three-step and moments models of photoemission that are used in simulation codes. Here, methods to include input from computationally intensive approaches, such as density functional theory to model optical properties and transfer matrix approaches to treat emission from the surface or transport past coatings, by means of parametric models are demonstrated. First, a technique to accurately represent optical behavior so as to model reflectivity and penetration depth is given. Second, modifications to bulk models arising from the usage of thin film architectures, and a means to rapidly calculate them, is provided. Third, a parameterization to model the impact of wells associated with coatings and surface layers on the transmission probably is given. In all cases, the methods are computationally efficient and designed to allow for including input from numerically intensive approaches that would otherwise be unavailable for simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anomalous compressible mode generation by global frame projections of pure Alfven mode

Alfven wave is the single most important physical phenomenon of magneto-hydrodynamic turbulence and has far-reaching impact to almost all studies related to astrophysical magnetic field. Yet the restoration of the Alfven wave fluctuations from a given magnetic field, aka the local Alfven wave problem, is never properly addressed in literature albeit its importance. Previous works model the Alfven wave fluctuation as the perturbation along a straight-line, constant magnetic field. However, Lazarian & Pogosyan (2012) suggested that the decomposition of Alfven wave along a straight line, aka. the global frame decomposition, has a factor of discrepancy to the true local Alfven wave fluctuation. Here, in this paper, we provide a geometric interpretation on how the local Alfven wave is related to the global frame through the use of vector frame formulation. We prove both analytically and numerically that the local frame Alfven wave is an orthogonal transformation of that of the global frame and related by the local Alfvenic Mach number. In other words, when we observe Alfven wave in the global frame of reference, some of the Alfven wave will be mistaken as compressible waves. The importance of frame choices have a far-reaching impact to the analytical studies of MHD turbulence. Combining the frame formalism and the new techniques we can have accurate measurement to some of the fundamental turbulence properties like the inclination angle of mean magnetic field relative to the line of sight.

79 ASTRONOMY AND ASTROPHYSICS↗

Multipartite entanglement in two-dimensional chiral topological liquids

The multipartite entanglement structure for the ground states of two-dimensional (2D) topological phases is an interesting albeit not well-understood question. Utilizing the bulk-boundary correspondence, the calculation of tripartite entanglement in 2D topological phases can be reduced to that of the vertex state, defined by the boundary conditions at the interfaces between spatial regions. In this paper, we use the conformal interface technique to calculate entanglement measures in the vertex state, which include area-law terms, corner contributions, and topological pieces, and a possible additional order-one contribution. Furthermore, this explains our previous observation of the Markov gap ℎ = $\frac{c}{3}$ ⁢ln⁡2 in the three-vertex state, and generalizes this result to the p-vertex state, general rational conformal field theories, and more choices of subsystems. Finally, we support our prediction by numerical evidence, finding precise agreement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multilevel Spectral Coarsening for Graph Laplacian Problems with Application to Reservoir Simulation

We extend previously developed two-level coarsening procedures for graph Laplacian problems written in a mixed saddle point form to the fully recursive multilevel case. The resulting hierarchy of discretizations gives rise to a hierarchy of upscaled models, in the sense that they provide approximation in the natural norms (in the mixed setting). This property enables us to utilize them in three applications: (i) as an accurate reduced model, (ii) as a tool in multilevel Monte Carlo simulations (in application to finite volume discretizations), and (iii) for providing a sequence of nonlinear operators in a full approximation scheme for solving nonlinear pressure equations discretized by the conservative two-point flux approximation. Finally, we illustrate the potential of the proposed multilevel technique in all three applications on a number of popular benchmark problems used in reservoir simulation.

multilevel Monte Carlo↗

Automated Process Planning for Embossing and Functionally Grading Materials via Site-Specific Control in Large-Format Metal-Based Additive Manufacturing

The potential for site-specific, process-parameter control is an attribute of additive manufacturing (AM) that makes it highly attractive as a manufacturing process. The research interest in the functionally grading material properties of numerous AM processes has been high for years. However, one of the issues that slows developmental progress in this area is process planning. It is not uncommon for manual programming methods and bespoke solutions to be utilized for site-specific control efforts. This article presents the development of slicing software that contains a fully automated process planning approach for enabling through-thickness, process-parameter control for a range of AM processes. The technique includes the use of parent and child geometries for controlling the locations of site-specific parameters, which are overlayed onto unmodified toolpaths, i.e., a vector-based planning approach is used in which additional information, such as melt pool size for large-scale metal AM processes, is assigned to the vectors. This technique has the potential for macro- and micro-structural modifications to printed objects. A proof-of-principle experiment is highlighted in which this technique was used to generate dynamic bead geometries that were deposited to induce a novel surface embossing effect, and additional software examples are presented that highlight software support for more complex objects.

36 MATERIALS SCIENCE↗