Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shape analysis”

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 775 records · Page 43

Energy dependence of $\textit{p}$ + 232 Th fission mass distributions: Mass-asymmetric standard I and standard II modes, and multichance fission

Background: The predominant mass-asymmetric fission of actinide nuclides occurs mainly through the so-called standard I and standard II modes. Though understood to be caused by shape-dependent shell structures encountered between the fission barrier deformation and scission, the most relevant shell gaps are still not firmly established. The standard I mode had been associated with the spherical doubly magic 132 Sn, and thus the $\textit{Z}$ = 50 proton shell, but recently it has been proposed that standard I and standard II are associated with quadrupole and octupole deformed gaps at $\textit{Z}$ = 52 and 56, respectively. Purpose: In this work, we investigate how the relative probabilities of the standard I and standard II modes vary with excitation energy near threshold, probing where the two modes bifurcate. Methods: The Australian National University Heavy Ion Accelerator Facility and CUBE fission spectrometer have been used to measure fission mass distributions for the $\textit{p}$ + 232 Th reaction (forming 233 Pa) at closely spaced bombarding energy intervals from 6.5 to 28 MeV. Results: A model-independent analysis of the energy dependence of the shape of the mass-asymmetric peak shows a strong dependence of the standard I and standard II relative probability on excitation energy near threshold. Conclusions: The results are consistent with the standard II mode having a lower fission barrier than standard I in 233 Pa, with the latter increasing continually in relative probability above its barrier energy. It is concluded that multichance fission, in particular last chance fission, plays a strong role in determining the observed energy dependence of all fission modes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effects of Ar-40 and Fe-56 ions on retinal photoreceptor cells of the rabbit: Implications for manned missions to Mars

Losses of photoreceptor cells (rods) from the retinas of New Zealand white (NZW) rabbits were detectable within 2 years after localized acute irradiation of optic and proximal tissues with greater than or equal to 7 Gy of 530 MeV u(exp -1) Ar-40 ions or greater than or equal to 2 Gy of 465 MeV u(exp -1) Fe-56 ions in the Bragg plateau region of energy deposition. Those limits were determined only from an analysis of variance of dose groups because the shapes of the dose response curves at early post-irradiation times are not known, a concern being addressed by experiments in progress. Losses of photoreceptor cells for the period 0.5-2.5 years post-irradiation, determined by provisional linear regression analysis, were approxiamtely 1.7% Gy(exp -1) and 2.5% Gy(exp.-1) for Ar-40 and Fe-56 ions, respectively.

Williams, G. R.↗

Effects of 40Ar and 56Fe ions on retinal photoreceptor cells of the rabbit: implications for manned missions to Mars

Losses of photoreceptor cells (rods) from the retinas of New Zealand white (NZW) rabbits were detectable within 2 years after localized acute irradiation of optic and proximal tissues with > or = 7 Gy of 530 MeV u-1 40Ar ions or > or = 2 Gy of 465 MeV u-1 56Fe ions in the Bragg plateau region of energy deposition. Those limits were determined only from an analysis of variance of dose groups because the shapes of the dose response curves at early post-irradiation times are not known, a concern being addressed by experiments in progress. Losses of photoreceptor cells for the period 0.5-2.5 years post-irradiation, determined by provisional linear regression analysis, were approximately 1.7% Gy-1 and 2.5% Gy-1 for 40Ar and 56Fe ions, respectively.

NASA Discipline Radiation Health↗

Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption in the ASCENDS 2011 Airborne Campaign

We have previously demonstrated an efficient pulsed, wavelength-resolved IPDA lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission. Our team participated in the 2010 ASCENDS airborne campaigns we flew airborne version of the CO2 and O2 lidar on the NASA DC-8. The CO2 lidar measures the atmospheric backscatter profiles and shape of the 1572.33 nm absorption line using 250 mW average laser power, 30 wavelength samples per scan and 300 scans per second. Most flights had 5-6 altitude steps to greater than 12 km, and clear CO2 line shapes were observed at all altitudes. Our post-flight analysis estimated the Iidar range and pulse energies at each wavelength every second. We then solved for the best-fit CO2 absorption line shape, and calculated the Differential Optical Depth (DOD) at the line peak. We compared these to CO2 DODs calculated from spectroscopy based on HITRAN 2008 and the conditions from airborne in-situ readings. Analysis of the 2010 measurements over the Pacific Ocean and Lamont OK shows the expected -linear change of the peak DOD with altitude. For measurements at altitudes greater than 6 km the random errors were approximately 0.3 ppm for 80 sec averaging times. After the 2010 flights we improved the airborne lidar's scan uniformity, calibration and receiver sensitivity. Our team participated in the seven ASCENDS science flights during late July and August 2011. These flights were made over a wide variety of surface and cloud conditions near the US, including over the central valley of California, over several mountain ranges, over both broken and solid stratus cloud deck over the Pacific Ocean, snow patches on mountain tops, over thin and broken clouds above the US Southwest and Iowa, and over forests near the WLEF tower in Wisconsin. Analyses show the retrievals of lidar range and CO2 column absorption, as wen as estimates of CO2 mixing ratio worked well when measuring over topography with rapidly changing height and reflectivity, through thin clouds and to stratus cloud tops. For regions where the CO2 concentration was relatively constant, the measured CO2 absorption profile (averaged for 50 sec) matched the predicted profile to better than 1% RMS error for all flight altitudes. For 1 & 10 second averaging, the scatter in the retrievals was limited by signal shot noise (i.e. the sigual photon count). Analysis to date shows the decrease in CO2 due to vegetation when flying easterward over the Great Plains as well as the increase in CO2 concentration in the vicinity ofthe coal-fired power plant in New Mexico. Examples of these and other results will be presented.

Abshire, James B.↗

Exploring the Influence of Impervious Surface Density and Shape on Urban Heat Islands in the Northeast USA Using MODIS and Landsat

Impervious surface area (ISA) from the National Land Cover Database (NLCD) 2001 and land surface temperature (LST) from MODIS averaged over three annual cycles (2003-2005) are used in a spatial analysis to assess the urban heat island (UHI) signature and its relationship to settlement size and shape, development intensity distribution, and land cover composition for 42 urban settlements embedded in forest biomes in the Northeastern United States. Development intensity zones, based on percent ISA, are defined for each urban area emanating outward from the urban core to nearby rural areas and are used to stratify land surface temperature. The stratification is further constrained by biome type and elevation to insure objective intercomparisons between urban zones within an urban settlement and between settlements. Stratification based on ISA allows the definition of hierarchically ordered urban zones that are consistent across urban settlements and scales. In addition to the surrounding ecological context, we find that the settlement size and shape as well as the development intensity distribution significantly influence the amplitude of summer daytime UHI. Within the Northeastern US temperate broadleaf mixed forest, UHI magnitude is positively related to the logarithm of the urban area size. Our study indicates that for similar urban area sizes, the development intensity distribution is one of the major drivers of UHI. In addition to urban area size and development intensity distribution, this analysis shows that both the shape of the urban area and the land cover composition in the surrounding rural area play an important role in modulating the UHI magnitude in different urban settlements. Our results indicate that remotely sensed urban area size and shape as well as the development intensity distribution influence UHI amplitude across regional scales.

Zhang, Ping↗

Enhanced Modeling of First-Order Plant Equations of Motion for Aeroelastic and Aeroservoelastic Applications

A methodology is described for generating first-order plant equations of motion for aeroelastic and aeroservoelastic applications. The description begins with the process of generating data files representing specialized mode-shapes, such as rigid-body and control surface modes, using both PATRAN and NASTRAN analysis. NASTRAN executes the 146 solution sequence using numerous Direct Matrix Abstraction Program (DMAP) calls to import the mode-shape files and to perform the aeroelastic response analysis. The aeroelastic response analysis calculates and extracts structural frequencies, generalized masses, frequency-dependent generalized aerodynamic force (GAF) coefficients, sensor deflections and load coefficients data as text-formatted data files. The data files are then re-sequenced and re-formatted using a custom written FORTRAN program. The text-formatted data files are stored and coefficients for s-plane equations are fitted to the frequency-dependent GAF coefficients using two Interactions of Structures, Aerodynamics and Controls (ISAC) programs. With tabular files from stored data created by ISAC, MATLAB generates the first-order aeroservoelastic plant equations of motion. These equations include control-surface actuator, turbulence, sensor and load modeling. Altitude varying root-locus plot and PSD plot results for a model of the F-18 aircraft are presented to demonstrate the capability.

Pototzky, Anthony S.↗

Pure quark and gluon observables in collinear drop

We construct a class of pure quark and gluon observables by using the collinear drop grooming technique. The construction is based on linear combinations of multiple cumulative distributions of the jet mass in collinear drop, whose specific weights are fully predicted perturbatively. This yields observables which obtain their values purely from quarks (or purely from gluons) in a wide region of phase space. We demonstrate this by showing that these observables are effective in two phase space regions, one dominated by perturbative resummation and one dominated by nonperturbative effects. The nonperturbative effects are included using shape functions which only appear as a common factor in the linear combinations constructed. We test this construction using a numerical analysis with next-to-leading logarithmic resummation and various shape function models, as well as analyzing these observables with PYTHIA and VINCIA. Choices for the collinear drop parameters are optimized for experimental use.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

New class of asymmetric shapes of rotating liquid drops

Shapes and stability of surface-tension-endowed drops rotating rigidly at fixed angular momentum are calculated by finite-element analysis. A new family of asymmetric two-lobed drop shapes is discovered that branches from, and rejoins, the Pik-Pichak family of symmetric two-lobed shapes. The computations are verified for axisymmetric and symmetric two-lobed drop shape by comparison with previous approximations.

Brown, R. A.↗

Design of fuselage shapes for natural laminar flow

Recent technological advances in airplane construction techniques and materials allow for the production of aerodynamic surfaces without significant waviness and roughness, permitting long runs of natural laminar flow (NLF). The present research effort seeks to refine and validate computational design tools for use in the design of axisymmetric and nonaxisymmetric natural-laminar-flow bodies. The principal task of the investigation involves fuselage body shaping using a computational design procedure. Analytical methods were refined and exploratory calculations conducted to predict laminar boundary-layer on selected body shapes. Using a low-order surface-singularity aerodynamic analysis program, pressure distribution, boundary-layer development, transition location and drag coefficient have been obtained for a number of body shapes including a representative business-aircraft fuselage. Extensive runs of laminar flow were predicted in regions of favorable pressure gradient on smooth body surfaces. A computational design procedure was developed to obtain a body shape with minimum drag having large extent of NLF.

Dodbele, S. S.↗

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λc, that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λc-based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

Amanda Stricklan↗

Meshfree simulation and experimental validation of extreme thermomechanical conditions in friction stir extrusion

Abstract Friction stir extrusion (FSE) is a novel solid-phase processing technique that consolidates and extrudes metal powders, flakes, chips, or billets into high-performance parts by plastic deformation, which has the potential to save substantial processing time and energy. Currently, most studies on FSE are experimental and only a few numerical models have been developed to explain and predict the complex physics of the process. In this work, a meshfree simulation framework based on smoothed particle hydrodynamics (SPH) was developed for FSE. Unlike traditional grid-based methods, SPH is a Lagrangian particle-based method that can handle severe material deformations, capture moving interfaces and surfaces, and monitor the field variable histories explicitly without complicated tracking schemes. These aspects of SPH make it attractive for the FSE process, where in situ evolution of field variables is difficult to observe experimentally. To this end, a 3-D, fully thermomechanically coupled SPH model was developed to simulate the FSE of aluminum wires. The developed model was thoroughly validated by comparing the numerically predicted material flow, strain, temperature history, and extrusion force with experimental results for a certain set of process parameters. The validated SPH model can serve as an effective tool to predict and better understand the extreme thermomechanical conditions during the FSE process.

99 GENERAL AND MISCELLANEOUS↗

Quasi-Static Analysis of LaRC THUNDER Actuators

An analytic approach is developed to predict the shape and displacement with voltage in the quasi-static limit of LaRC Thunder Actuators. The problem is treated with classical lamination theory and Von Karman non-linear analysis. In the case of classical lamination theory exact analytic solutions are found. It is shown that classical lamination theory is insufficient to describe the physical situation for large actuators but is sufficient for very small actuators. Numerical results are presented for the non-linear analysis and compared with experimental measurements. Snap-through behavior, bifurcation, and stability are presented and discussed.

Campbell, Joel F.↗

Analysis of SMA Hybrid Composite Structures using Commercial Codes

A thermomechanical model for shape memory alloy (SMA) actuators and SMA hybrid composite (SMAHC) structures has been recently implemented in the commercial finite element codes MSC.Nastran and ABAQUS. The model may be easily implemented in any code that has the capability for analysis of laminated composite structures with temperature dependent material properties. The model is also relatively easy to use and requires input of only fundamental engineering properties. A brief description of the model is presented, followed by discussion of implementation and usage in the commercial codes. Results are presented from static and dynamic analysis of SMAHC beams of two types; a beam clamped at each end and a cantilevered beam. Nonlinear static (post-buckling) and random response analyses are demonstrated for the first specimen. Static deflection (shape) control is demonstrated for the cantilevered beam. Approaches for modeling SMAHC material systems with embedded SMA in ribbon and small round wire product forms are demonstrated and compared. The results from the commercial codes are compared to those from a research code as validation of the commercial implementations; excellent correlation is achieved in all cases.

Turner, Travis L.↗

Macroscopic Traffic Modeling Using Probe Vehicle Data: A Machine Learning Approach

Abstract The macroscopic fundamental diagram (MFD) captures an orderly relationship among traffic flow, density, and speed at the network level. It is a simple yet powerful tool for modeling traffic dynamics in large urban networks with broad application in traffic control and management. However, empirically derived MFDs in urban regions require high-resolution traffic data from the network. Having the network flow and vehicular density estimated at the (granular) census tract level using vehicle probe data, we apply machine learning methods to predict the MFDs across U.S. urban areas and capture the impacts of location-specific input features on the network flow–density relationships at a large scale. The results show that, among the four tested machine learning approaches (Random Forest, XGBoost, Support Vector Machine, and Neural Network), XGBoost delivers the best performance in predicting network traffic flow based on vehicular density and location attributes. Using interaction Shapley Additive explanation (SHAP) values and partial correlation analysis, we examine the factors influencing MFD shapes across different locations. Our empirical findings reveal that across U.S. urban areas, network topology, transportation infrastructure, and land use are primary factors shaping MFD curves, while demand and trip-related factors play a lesser role. Specifically, higher ranking roads, centrality, and development levels correlate positively with network capacity and critical density, whereas negative associations are observed for network connectivity, mixed-use development, and road roughness levels.

Jin, Ling↗

How Flood Hazards in a Warming Climate Could Be Amplified by Changes in Spatiotemporal Patterns and Mechanisms of Water Available for Runoff

Prior research on climate change impacts on flooding has primarily focused on changes in extreme rainfall magnitudes, often neglecting snow processes and spatiotemporal storm patterns, such as hyetograph shapes and areal reduction factors (ARFs). This study examines projected changes in extreme water available for runoff (W) events in two snow-dominated basins in the western United States: the Yakima River Basin (YRB) in Washington State and the Walker River Basin (WRB) spanning the California-Nevada border. We analyze changes in W magnitudes, mechanisms, hyetograph shapes, and ARFs, and study their compounded impacts on flood hazard. Our findings suggest increased extreme W magnitudes across a large portion of the basins, with steeper or flatter hyetographs, and higher ARF values under the future climate. These changes are driven by a shift from seasonal snowmelt to more rain-on-snow events at higher elevations and by increased rainfall at lower elevations. We then use a single event-based rainfall-runoff model to estimate flood hazard changes based on extreme W magnitudes, hyetograph shapes, ARFs, and their compounded impacts. Our analysis reveals that focusing solely on the magnitude of changes in extreme W can significantly underestimate future flood hazards and uncertainties. Ignoring future changes in spatiotemporal patterns can underestimate future flood hazards by 63% and underestimate the uncertainty in future flood events by 18% in the WRB. These results underscore the necessity of incorporating spatiotemporal dynamics into future flood hazard assessments to provide a more accurate evaluation of potential impacts.

54 ENVIRONMENTAL SCIENCES↗

Creating ground truth for nanocrystal morphology: a fully automated pipeline for unbiased transmission electron microscopy analysis

Control over colloidal nanocrystal morphology (size, size distribution, and shape) is important for tailoring the functionality of individual nanocrystals and their ensemble behavior. Despite this, traditional methods to quantify nanocrystal morphology are laborious. New developments in automated morphology classification will accelerate these analyses but the assessment of machine learning models is limited by human accuracy for ground truth, causing even unsupervised machine learning models to have inherent bias. Herein, we introduce synthetic image rendering to solve the ground truth problem of nanocrystal morphology classification. By simulating 2D images of nanocrystal shapes via a function of high-dimensional parameter space, we trained a convolutional neural network to link unique morphologies to their simulated parameters, defining nanocrystal morphology quantitatively rather than qualitatively. An automated pipeline then processes, quantitatively defines, and classifies nanocrystal morphology from experimental transmission electron microscopy (TEM) images. Using improved computer vision techniques, 42,650 nanocrystals were identified, assessed, and labeled with quantitative parameters, offering a 600-fold improvement in efficiency over best-practice manual measurements. Further, a classification algorithm was trained with a prediction accuracy of 99.5%, which can successfully analyze a range of concave, convex, and irregular nanocrystal shapes. The resulting pipeline was applied to differentiating two syntheses of nominally cuboidal CsPbBr 3 nanocrystals and uniquely classifying binary nickel sulfide nanocrystal phase based on morphology. This pipeline provides a simple, efficient, and unbiased method to quantify nanocrystal morphology and represents a practical route to construct large datasets with an absolute ground truth for training unbiased morphology-based machine learning algorithms.

77 NANOSCIENCE AND NANOTECHNOLOGY↗