Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical model”

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

Stability study of aqueous foams under high-temperature and high-pressure conditions relevant to Enhanced Geothermal Systems (EGS)

In Enhanced Geothermal System (EGS) hydraulic fracturing is carried out by injecting cold water into deep Hot Dry Rocks (HDR) under carefully controlled conditions to create new or reopen existing fractures. Water-based fracturing fluids demonstrate some challenges including immense quantity of water usage, water sensitivity of the formations, water blocking, and lack of proppant carrying capacity and transportation. Thus, an alternative is to use foam-based fracturing fluid which offers potential advantage over conventional water-based fracturing fluid such as minimum water usage, reduced wellbore damage, high proppant carrying capacity, and less environmental damage. However, foams are complex mixture of gaseous phase and liquid phase which are thermodynamically unstable at downhole conditions, and their stability can decrease over time due to liquid drainage, bubble coarsening, and coalescence. This paper shows laboratory experiments executed to study foam stability at high temperature (200°C) and high pressure (6.9MPa) conditions which simulates the geothermal environment. Foam stability was characterized by half-life of foam, which is defined as the time taken by the foam to decreases by 50% of its original height due to drainage. In this paper, two types of gaseous phases, nitrogen (N 2 ) and carbon dioxide (CO 2 ) were investigated. Also, based on successful practice of foam-based fracturing fluid in oil and gas industries, four surfactants, including Alpha olefin sulfonate (AOS), Sodium dodecyl sulfonate (SDS), Tergitol™ (NP – 40), and Cetyltrimethylammonium chloride (CTAC) at optimum concentration of 1wt.% were tested for best stability performance. In addition, different stabilizing agents including guar gum, bentonite clay, crosslinker, silicon dioxide nanoparticles (SiO 2 ), graphene oxide (GO) were also studied. Experimental results showed that N 2 foams were more stable than CO 2 foams. It was observed that foam half-life decreased with the increase in temperature. Among all the surfactants, AOS foams showed the most promising thermal stability at high temperatures. Moreover, with the addition of stabilizing agents, foam's half-life was enhanced. Stabilizing agents such as crosslinker and GO dispersion showed the most stable foams with half-life recorded at 20min and 17min, respectively, at 200°C and 6.9MPa. Finally, pressure also showed a positive effect on foam stability; with increased pressure, foam half-life was increased. Based on the experimental data, analytical models for the effect of temperature and pressure were developed, considering foam degradation is a first-order kinetic reaction that linearly depends on the foam drainage mechanism. The effect of temperature on foam half-life was studied as an exponential decay model. In this model, foam half-life is a function of drainage rate constant (D A ) and activation energy (E a ) of the foam system. The effect of pressure on foam half-life was found to obey a power-law model where an increase in pressure showed an increase in foam half-life. Furthermore, a linear relation was studied for the effect of pressure on foam activation energy and drainage rate. Then the combined effects of temperature and pressure were studied, which yielded an analytical model to predict the foam stabilities in terms of half-life for different foam compositions. In conclusion, this research indicates that with an appropriate selection of surfactants and stabilizing agents, it is possible to obtain stable foams, which could replace conventional water fracturing fluid under EGS conditions.

15 GEOTHERMAL ENERGY↗

Closed-form theoretical model of the secondary drop size in partial coalescence—Capturing pertinent timescales and viscous forces

One of the most important outcomes of partial drop coalescence is the ratio of the secondary drop radius to the primary drop radius, known as the drop ratio, ri. This ratio is thought to be approximately constant and independent of physical parameters of the fluids involved. However, here we reveal that ambient fluid viscosity can alter the size of the secondary drop and the drop ratio consequently. Using scaling analysis, we derive a model that predicts the behavior of the drop ratio as a function of the Ohnesorge number, a dimensionless ratio of viscous to inertial forces. In addition, we present our experimental results of coalescing drops on a planar interface under the influence of surface tension gradients. A high-speed digital camera is used to observe the evolution of drops as they coalesce with a bulk liquid. We show that this process is influenced by the surface tension gradient between the drop and the bulk liquid. The ratio of the secondary drop to the primary drop in partial coalescence is smaller than the reported values for coalescence without a surface tension gradient. The analytical model derived through this study is based on a new modified Ohnesorge number that includes surface tension gradients. Our analytical model is compared against other models and the results illustrate good agreement with our experimental findings and experimental data in the literature.

42 ENGINEERING↗

Elucidating Cycling Rate-Dependent Electrochemical Strains in Sodium Iron Phosphate Cathodes for Na-ion Batteries

Battery electrodes materials undergo significant mechanical instabilities which affects their longevity and exert rate-limitations during the cycling process. In this study, we investigate the rate-dependent mechanical response of sodium iron phosphate (NaFePO4, NFP) cathodes during Na intercalation via galvanostatic cycling at different rates by employing digital image correlation, electrochemical methods, and mathematical model. The mechanical behaviour of the electrode shows strong dependance on the applied scan rate. At slower rates, electrode shows asymmetrical strain generation between anodic and cathodic cycles, which is attributed to the formation of cathode-electrolyte interface layers. The electrode undergoes smaller strain generation when cycled at slower rates when the same amount of Na ions is removed or inserted into the electrode. A mathematical model was developed to predict strain evolution in the composite electrode as well as the concentration profile of the Na ions in the electrode particles. Rate-dependent and time dependent factors on the strain generation in the electrode are attributed to the capacity-dependent intercalation strains, rate-dependent mismatch strains, and time-dependent irreversible strains. The combination of in situ strain measurements with the analytical model provided new insight into the electrochemically induced mechanical deformations in Na-ion cathode electrodes.

Ozdogru, Bertan↗

Back-Azimuth Estimation of Air-to-Ground Coupled Infrasound from Transverse Coherence Minimization

We present the transverse coherence minimization method (TCM)—an approach to estimate the back–azimuth of infrasound signals that are recorded on an infrasound microphone and a colocated three–component seismometer. Accurate back–azimuth information is important for a variety of monitoring efforts, but it is currently only available for infrasound arrays and for seismoacoustic sensor pairs separated by 10 s of meters. Our TCM method allows for the analysis of colocated sensor pairs, sensors located within a few meters of each other, which may extend the capabilities of existing seismoacoustic networks and supplement operating infrasound arrays. This approach minimizes the coherence of the transverse component of seismic displacement with the infrasound wave to estimate the infrasound back–azimuth. After developing an analytical model, we investigate seismoacoustic signals from the August 2012 Humming Roadrunner experiment and the 26 May 2021 eruption of Great Sitkin Volcano, Alaska, U.S.A., at the ranges of 6.5–185 km from the source. We discuss back–azimuth estimates and potential sources of deviation (1°–15°), such as local terrain effects or deviation from common analytical models. This practical method complements existing seismoacoustic tools and may be suitable for routine application to signals of interest.

58 GEOSCIENCES↗

Highly selective Si 3 N 4 /SiO 2 etching using an NF 3 /N 2 /O 2 /H 2 remote plasma. I. Plasma source and critical fluxes

Highly selective plasma etching of silicon nitride (Si 3 N 4 ), while not etching silicon dioxide (SiO 2 ), is a critical step in the fabrication of microelectronics devices. In many applications, this etching must be damage-free and isotropic, which then motivates the use of remote plasmas where the reactants interacting with the substrate are dominantly neutral species. In this paper and Paper II, mechanisms for highly selective Si 3 N 4 etching in remote plasmas are discussed based on results from experiments and simulations. It has been shown experimentally that high Si 3 N 4 /SiO 2 etch selectivity (≈380) can be achieved in the downstream effluent of an NF 3 /N 2 /O 2 /H 2 plasma. The authors found that H2 plays a principal role in the reaction mechanism as Si 3 N 4 /SiO 2 selectivity shows a sharp maximum as a function of the H2 flow rate. Based on this observation, and measured densities of F-atoms and H2 in the process chamber, a mechanism of selective Si 3 N 4 /SiO 2 etching is proposed in which HF molecules in vibrationally excited states accelerate etching reactions. A reaction mechanism for NF 3 /N 2 /O 2 /H 2 plasmas and its afterglow was developed to computationally determine the species densities and fluxes on the wafer level, validated by comparing with experimentally measured F-atom and H 2 densities. The calculated species densities and fluxes were used as input to an analytical model of Si 3 N 4 and SiO 2 etching based on the results of quantum chemistry simulations. This paper presents experimental results (etching data and species densities), the reaction mechanism for NF 3 /N 2 /O 2 /H 2 plasmas, and the results of simulations of gas phase chemistry. Quantum chemistry simulations of elementary etching reactions, description of the analytical model of Si 3 N 4 and SiO 2 etching, calculations of the etch rates, and Si 3 N 4 and SiO 2 selectivity with this model are presented in Paper II.

42 ENGINEERING↗

The Effect of Grain Boundary Facet Junctions on Segregation and Embrittlement

Junctions are discontinuities in flat grain boundaries that arise in all polycrystalline materials and are thought to play important roles in the response of a grain boundary network to thermal and mechanical loads. A key open question concerns the mechanisms by which solute segregation to junctions impacts properties of the grain boundary. Here, in this work, we investigate the influence of grain boundary facet junctions on solute embrittlement, and we present an analytical model that uses the hydrostatic stress field contributed by dislocations at multiple junctions to describe these effects. Specifically, we study junctions between {112} facets of various lengths in Au $\langle111\rangle$ Σ3 tilt grain boundaries. Copper and silver solutes are employed to determine if the effect of junctions on solute segregation and embrittlement is dependent on size relative to the host. Combined, atomistic simulation data and the analytical model show that Cu and Ag have opposite segregation responses to junctions due to the sign of the hydrostatic stress field induced by junctions. However, a positive shift in the embrittling potency is computed near junctions regardless of solute type or the stress state of the segregation site. Hence, for the conditions studied, junctions consistently shift the energetic landscape towards embrittlement.

36 MATERIALS SCIENCE↗

Static and dynamic properties of multi-ionic plasma mixtures

Complex plasma mixtures with three or more components are often encountered in astrophysics or in inertial confinement fusion (ICF) experiments. For mixtures containing species with large differences in atomic number Z, the modeling needs to consider at the same time the kinetic theory for low-Z elements combined with the theory of strongly coupled plasma for high-Z elements, as well as all the intermediate situations that can appear in multi-component systems. For such cases, we study the pair distribution functions, self-diffusions, mutual diffusion and viscosity for ternary mixtures at extreme conditions. These quantities can be produced from first principles using orbital free molecular dynamics at the computational expense of very intensive simulations to reach good statistics. Utilizing the first-principles results as reference data, we assess the merit of a global analytic model for transport coefficients, "Pseudo-Ions in Jellium" (PIJ), based on an iso-electronic assumption (iso-n e ). With a multi-component hypernetted-chain integral equation, we verify the quality of the iso-n e prescription for describing the static structure of the mixtures. This semi-analytical modeling compares well with the simulation results and allows one to consider plasma mixtures not accessible to simulations. Further, applications are given for the mix of materials in ICF experiments. A reduction of a multicomponent mixture to an effective binary mixture is also established in the hydrodynamic limit and compared with PIJ estimations for ICF relevant mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of discharge polarity on streamer breakdown criterion of ambient air in a non-uniform electric field

Streamer breakdown of atmospheric air with non-uniform dc electric field in a needle-to-plate electrode configuration is studied using a semi-analytic model and experimental measurements. A high voltage (either positive or negative) is applied to a hollow needle with 0.51 mm outer diameter and 0.25 mm inner diameter separated from a planar ground electrode by a gap distance of 0.1–1.4 cm. Breakdown voltages are recorded for both positive and negative discharge polarities. Empirical relations between the critical avalanche size for streamer breakdown and the gap distance are proposed. Using these empirical relations, a semi-analytic model based on Meek's criterion for streamer breakdown is developed to accurately predict the measured breakdown voltages. It is found that for pd > 380 Torr cm (or d > 0.5 cm at one atmosphere) streamer breakdown of ambient air occurs at a lower applied voltage for a positively biased needle compared to that with a negatively biased needle, referred as the polarity effect. For pd < 380 Torr cm breakdown is attained at a lower applied voltage with a negatively biased needle compared to that with a positively biased needle, and breakdown mode transits from the polarity effect to the so called inverted polarity effect.

42 ENGINEERING↗

The three-dimensional structure of black hole accretion flows within the plunging region

ABSTRACT We analyse, using new analytical models and numerical general relativistic magnetohydrodynamic simulations, the three-dimensional properties of accretion flows inside the plunging region of black hole spacetimes (i.e. at radii smaller than the innermost stable circular orbit). These simulations are of thick discs, with aspect ratios of order unity $h/r \sim 1$, and with a magnetic field geometry given by the standard low-magnetization ‘SANE’ configuration. This work represents the first step in a wider analysis of this highly relativistic region. We show that analytical expressions derived in the ‘thin disc’ limit describe the numerical results remarkably well, despite the large aspect ratio of the flow. We further demonstrate that accretion within this region is typically mediated via spiral arms, and that the geometric properties of these spiral structures can be understood with a simple analytical model. These results highlight how accretion within the plunging region is fundamentally two-dimensional in character, which may have a number of observational implications. We derive a modified theoretical description of the pressure within the plunging region which accounts for turbulent heating and may be of use to black hole image modelling.

Astronomy & Astrophysics↗

Resilience Analytics for Interdependent Power and Water Distribution Systems

The occurrence of large-scale outages in power distribution systems (PDSs) caused, for instance, by extreme weather events has raised concerns on the cascading impact of the power outages on the resilience of other critical infrastructures, including water distribution systems (WDSs), whose healthy operation highly relies on the availability of power from PDSs. Here, this paper proposes an analytics model for quantifying the interdependence between the resilience of power and water distribution systems. The proposed model first performs spatio-temporal outage analysis on the PDS and WDS to determine the PDS power serving capability and the resulting cascading impact on the WDS water serving capability during the occurrence of disruptions. The resilience interdependence is then determined by a proposed set of six metrics calculated using the power and water serving capability profiles, which trace and quantify the cascading impacts of power outages in the operation of WDSs during the degradation and recovery stages. The proposed analytics model is implemented on the IEEE 33-bus power distribution system supplying power for pumping stations of a 15-node test WDS, considering multiple cases of water storage capacity, single and multiple outage scenarios, and availability of distributed generation. The numerical results demonstrate that the proposed metrics track the temporal and spatial interdependence between the resilience of power and water distribution systems, providing a mechanism to identify weak spots, and assess options to enhance the resilience of the critical water infrastructure.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Collaborative Research: Explosive reconnection in relativistic magnetically-dominated plasmas (Final Report)

Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theoretical and Experimental Study of Current from Non-Disintegrable Suspended Particles at a Rotating Disk Electrode

Understanding the current response at an electrode from suspended solid particles in an electrolyte is crucial for developing materials to be used in semi-solid electrodes for energy storage applications. Here, an analytical model is proposed to predict and understand the current response from non-disintegrable solid particles at a rotating disk electrode. The current is shown to be limited by a combination of ion diffusion within the solid particle and the mean residence time of the particle at the rotating disk electrode. This results in a relationship between current and angular frequency of I ∝ ω 3 / 4 , instead of the classical I ∝ ω 1 / 2 predicted by Levich theory. Specifically, the current response of Li 4 Ti 5 O 12 (LTO) microparticles suspended in a non-aqueous electrolyte of lithium hexafluorophosphate (LiPF 6 ) in ethylene carbonate: diethyl carbonate (EC:DEC) was determined experimentally and compared favorably with predictions from the proposed analytical model using fitting parameters consistent with the experimental conditions.

25 ENERGY STORAGE↗

Machine Learning Full NMR Chemical Shift Tensors of Silicon Oxides with Equivariant Graph Neural Networks

The nuclear magnetic resonance (NMR) chemical shift tensor is a highly sensitive probe of the electronic structure of an atom and furthermore its local structure. Recently, machine learning has been applied to NMR in the prediction of isotropic chemical shifts from a structure. Current machine learning models, however, often ignore the full chemical shift tensor for the easier-to-predict isotropic chemical shift, effectively ignoring a multitude of structural information available in the NMR chemical shift tensor. Here we use an equivariant graph neural network (GNN) to predict full 29 Si chemical shift tensors in silicate materials. The equivariant GNN model predicts full tensors to a mean absolute error of 1.05 ppm and is able to accurately determine the magnitude, anisotropy, and tensor orientation in a diverse set of silicon oxide local structures. When compared with other models, the equivariant GNN model outperforms the state-of-the-art machine learning models by 53%. The equivariant GNN model also outperforms historic analytical models by 57% for isotropic chemical shift and 91% for anisotropy. The software is available as a simple-to-use open-source repository, allowing similar models to be created and trained with ease.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fusion gain from cylindrical liner-driven implosions of field reversed configurations

MagLIF experiments [M.R. Gomez et al., Phys. Plasmas 22, 056306 (2015)] on Z have demonstrated the basic principles of Magneto-Inertial Fusion (MIF) for wall confined plasmas. Other MIF schemes have been proposed based on the liner implosion of closed field magnetically confined plasmas such as Field Reversed Configurations (FRCs) [T. P. Intrator et al., Phys. Plasmas 15, 042505 (2008)]. We present a semi-analytical model of liner driven FRC implosions that predicts the fusion gain of such systems. The model predicts a fusion gain near unity for an FRC imploded by a liner driven by the Z Machine. We show that FRCs could be formed and imploded at the Z facility using the AutoMag liner concept [S. A. Slutz et al., Phys. Plasmas 24, 012704 (2017)]. Furthermore, an initial bias magnetic field can be supplied by the external magnets used in MagLIF experiments. The reverse field is then supplied by an AutoMag liner, which has helical conducting paths imbedded in an insulating substance. Experiments [Shipley et al., Phys. Plasmas 26, 052705 (2019)] have demonstrated that AutoMag can generate magnetic fields greater than 30 Tesla inside of the liner. We have performed 2D Radiation MHD simulations of the formation and implosion of an FRC, which are in good agreement with the analytical model. The FRC formation process could be studied on small pulsed power machines delivering about 1 MA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A unified thermionic and thermionic-field emission (TE–TFE) model for ideal Schottky reverse-bias leakage current

We present a unified thermionic emission (TE) and thermionic-field emission (TFE) model for the ideal reverse-bias leakage current in Schottky junctions. The unified TE–TFE analytical model advances upon previous analytical TFE models by Murphy–Good and Padovani–Stratton, which are the two most widely adopted models by the community, in two major aspects: (i) the applicability of the TFE expression therein is extended to near-zero surface electric fields by an error-function correction, allowing for the calculation of the total current by a nontrivial sum of TE (over-the-barrier current) and TFE (below-the-barrier current) contributions; therefore, an accurate description of the TE-to-TFE transition region is captured analytically for the first time; (ii) image-force lowering is considered with much-simpler correction terms. Comparisons with the reference numerical model show that the unified TE–TFE model has excellent accuracy, as well as a 10 000× reduction in computation time. The unified model is further tested against experimental data from Schottky barrier diodes based on Si, 4H-SiC, GaN, and β-Ga 2 O 3 , revealing accurate extractions of barrier heights and correct descriptions of the ideal reverse leakage characteristics. With the extended applicable range, improved accuracy, and computational efficiency, the unified TE–TFE model is highly valuable for the design and analysis of devices based on Schottky junctions, as well as for potential integration in technology computer-aided design (TCAD) tools.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Simulating image coaddition with the Nancy Grace Roman Space Telescope – II. Analysis of the simulated images and implications for weak lensing

One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom – an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image – on the Roman image simulations. In this paper, we measure the output noise power spectra, identify the sources of the major features in the power spectra, and show that simple analytic models that ignore sampling effects underestimate the power spectra of the coadded noise images. We compute the moments of both idealized injected stars and fully simulated stars in the coadded images, and their one- and two-point statistics. We show that the idealized injected stars have root-mean-square ellipticity errors (1–6) × 10 -4 per component depending on the band; the correlation functions are ≥2 orders of magnitude below requirements, indicating that the image combination step itself is using a small fraction of the overall Roman second moment error budget, although the fourth moments are larger and warrant further investigation. The stars in the simulated sky images, which include blending and chromaticity effects, have correlation functions near the requirement level (and below the requirement level in a wide-band image constructed by stacking all four filters). We evaluate the noise-induced biases in the ellipticities of injected stars, and explain the resulting trends with an analytical model. We conclude by enumerating the next steps in developing an image coaddition pipeline for Roman.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparative Studies of Optoelectrical Properties of Prominent PV Materials: Halide Perovskite, CdTe, and GaAs

We compare three representative high performance PV materials: halide perovskite MAPbI3, CdTe, and GaAs, in terms of photoluminescence (PL) efficiency, PL lineshape, carrier diffusion, and surface recombination and passivation, over multiple orders of photo-excitation density or carrier density appropriate for different applications. An analytic model is used to describe the excitation density dependence of PL intensity and extract the internal PL efficiency and multiple pertinent recombination parameters. A PL imaging technique is used to obtain carrier diffusion length without using a PL quencher, thus, free of unintended influence beyond pure diffusion. Our results show that perovskite samples tend to exhibit lower Shockley-Read-Hall (SRH) recombination rate in both bulk and surface, thus higher PL efficiency than the inorganic counterparts, particularly under low excitation density, even with no or preliminary surface passivation. PL lineshape and diffusion analysis indicate that there is considerable structural disordering in the perovskite materials, and thus photo-generated carriers are not in global thermal equilibrium, which in turn suppresses the nonradiative recombination. This study suggests that relatively low point-defect density, less detrimental surface recombination, and moderate structural disordering contribute to the high PV efficiency in the perovskite. This comparative photovoltaics study provides more insights into the fundamental material science and the search for optimal device designs by learning from different technologies.

14 SOLAR ENERGY↗

Probing the Mechanical Performance of Micro-architected Porous Structures Through In Situ Characterization and Analysis

Micropores play critical roles in both natural and man-made materials. Such pores take on a variety of shapes and sizes ranging from spherical to irregular sphere-like voids with diameters spanning from the nanometer to millimeter scales. When porous structures are mechanically loaded, the pores direct the stress around their free surfaces, altering the material’s mechanical response relative to fully dense materials. In this study, for the first time, we create micro-architected porous samples using nanolithography and investigate the role of pore morphology by conducting a series of in situ micropillar compression tests in scanning electron microscopy. The findings demonstrate that porosity is the primary factor influencing the mechanical response of these micro-architected materials, as often seen on a macroscopic level. Additionally, we observed that pore geometry had a significant impact on Young’s modulus, yield stress, and strain energy density as a secondary parameter. Then, the extracted Young’s modulus was compared to macroscopic empirical models and determined the analytical models sufficiently described the impact of porosity in the microscopic scale but failed to capture the impact of second-order parameters. Furthermore, these results suggest how porous materials can be tailored to achieve desired mechanical properties based on the engineering applications of interest.

36 MATERIALS SCIENCE↗