SEARCH · Engineering Papers

Results for “shapes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2026

How nitrogen and oxygen shape SRF cavity performance

Nitrogen and oxygen-based surface treatments have revolutionized the performance of superconducting radiofrequency (SRF) cavities, enabling them to reach higher gradients and lower losses. However, the exact mechanisms by which these treatments improve cavity performance remain largely unknown. This work provides new insights into the role of nitrogen and oxygen in SRF cavity performance by using time-of-flight secondary ion mass spectrometry (TOF-SIMS) to precisely quantify the concentrations and depth profiles of these impurities within niobium cutouts. We correlate these impurity profiles with detailed cavity performance measurements, including surface resistance and quality factor, and compare our findings with predictions from BCS theory. The results demonstrate that while both nitrogen and oxygen enhance performance, ten times more oxygen is required to achieve the same reduction in BCS resistance as interstitial nitrogen. We present a potential model in which the observed variation arises from nitrogen's greater effectiveness in trapping hydrogen, thus reducing the formation of niobium hydrides and enhancing superconducting gap.

Hu, Hannah [Chicago U.]↗

DOE OSTI2024

Rapid Net Shape Fabrication of Complex Shapes

Explore the source record for details and available documents.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A Computational Tool Compatible with NEAMS Code Packages for Optimizing the Shape of Nuclear Reactor Components and of Whole Core Performance

We designed and implemented a shape optimization tool that functions with NEAMS codes, and that nuclear scientists and engineers can employ to optimize the shape of individual components and the whole core under the applicable single- or multi-physics model comprising the employed code(s). The shape-optimization tool enables varying the geometric shape itself as well as its dimensions to yield, potentially, new component designs that are not limited by the designer’s intuition and previous experience. In cases where the optimal-shape object is an individual component, we provide the capability for additional verification that the whole-core performance using the optimized component performs better, under the prescribed optimization criteria, than the initial design. Our shape-optimization tool couples to NEAMS codes via a flexible input- composer interface and enables the user to constrain the shape’s evolution to ensure the component’s manufacturability. Finally, we demonstrate our shape-optimization tool with single- and multi-physics NEAMS codes. This objective is motivated by the recent advances in manufacturing technology that, combined with rising interest in novel reactor concepts, are creating new opportunities for innovation in the design of individual components that affect the performance of the full reactor system. In particular, Additive Manufacturing (AM) enables mass production of highly precise, intricate and complex component shapes that are not feasible with traditional manufacturing techniques. To accomplish this goal we developed and implemented in MOOSE: (1) discrete shape optimization capability based on a state-space search that uses Artificial Intelligence strategies to find the optimal state/shape; (2) smooth shape optimization tool that employs PETSc’s toolkit for advanced optimization (TAO) to optimize node-displacement of the components’ model sidesets; (3) hierarchical core optimization workflow that recognizes the repeating patterns typical in a nuclear reactor and performs the optimization one level at a time with increasing length scale. Each of these tools is equipped with user-specified constraints to avoid optimal shapes that are not manufacturable. The developed shape optimization tool is verified and demonstrated on various nuclear reactor core components and models. The optimization process accounts for tightly coupled physics that govern the behavior of these target reactors, and exercises several NEAMS codes in a coupled multiphysics fashion. The impact of the delivered shape optimization tool will materialize in the optimal design, from the outset, of advanced reactors currently contemplated to regain the US’s leadership in nuclear energy R&D. Novel reactor concepts, e.g. Molten Salt Reactors, and sizes/capacities, e.g. micro- reactors, provide a unique opportunity to optimize performance from the early stages of development, before the investment in components’ production lines, validation experiments, and licensing regimes make future improvements in performance prohibitively expensive and force sub-optimal performance on the affected reactor concept in perpetuity. This benefit will be realized by the delivered shape optimization tool regardless of the applicable manufacturing process whether traditional or AM, thereby broadening the impact of this project on current and future reactor concepts and technologies

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing

Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2026

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation) [SWR-26-095]

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation): Multifidelity aerodynamic polar data generation for hydrofoil/tidal-turbine airfoil sections. Foilpolars ties together three pieces: *AeroSandbox supplies the baseline airfoil coordinates (UIUC database). *G2Aero parameterizes those shapes on a Grassmannian manifold (Karcher mean + PGA basis) and samples new perturbed shapes around that basis. *XFoil (panel method) and NeuralFoil (neural-network surrogate, shipped with AeroSandbox) each solve the resulting shapes for lift, drag, moment, and pressure at the swept angles of attack, Reynolds numbers, and n_crit values. Design optimization of foil shapes in a computationally efficient way requires polars data across many candidate shapes, not just a handful of baseline foils. However, high-fidelity CFD at that scale is too costly, and naive shape perturbation strays from realistic geometries. FOILPOLARS addresses this by loading baseline airfoils (via AeroSandbox) and mapping them onto a Grassmannian manifold (via G2Aero), computing a Karcher mean and principal geodesic analysis (PGA) basis. New shapes are sampled by perturbing PGA coefficients, keeping them close to the manifold of realistic foils. Each sampled shape is evaluated across a configurable sweep of angle of attack, Reynolds number, and critical amplification factor using two solvers: XFoil (panel method) and NeuralFoil (neural-network surrogate), producing a paired dataset of lift, drag, moment, pressure, convergence, and confidence, indexed alongside each shape's PGA coefficients and shared Grassmannian basis in a single xarray dataset. From this, FOILPOLARS produces convergence summaries and comparison plots per shape, Reynolds number, and n_crit. A command-line interface exposes each pipeline stage independently, supporting data-driven design, optimization, and machine-learning workflows for foils.

Sandhu, Rimple [National Laboratory of the Rockies↗

DOE OSTI2023

Learning Canonical Embeddings for Unsupervised Shape Correspondence With Locally Linear Transformations

We present a new approach to unsupervised shape correspondence learning between pairs of point clouds. We make the first attempt to adapt the classical locally linear embedding algorithm (LLE)-originally designed for nonlinear dimensionality reduction-for shape correspondence. The key idea is to find dense correspondences between shapes by first obtaining high-dimensional neighborhood-preserving embeddings of low-dimensional point clouds and subsequently aligning the source and target embeddings using locally linear transformations. We demonstrate that learning the embedding using a new LLE-inspired point cloud reconstruction objective results in accurate shape correspondences. More specifically, the approach comprises an end-to-end learnable framework of extracting high-dimensional neighborhood-preserving embeddings, estimating locally linear transformations in the embedding space, and reconstructing shapes via divergence measure-based alignment of probability density functions built over reconstructed and target shapes. Our approach enforces embeddings of shapes in correspondence to lie in the same universal/canonical embedding space, which eventually helps regularize the learning process and leads to a simple nearest neighbors approach between shape embeddings for finding reliable correspondences. Comprehensive experiments show that the new method makes noticeable improvements over state-of-the-art approaches on standard shape correspondence benchmark datasets covering both human and nonhuman shapes.

deformation↗

DOE OSTI2022

Shape memory polymers

New shape memory polymer compositions, methods for synthesizing new shape memory polymers, and apparatus comprising an actuator and a shape memory polymer wherein the shape memory polymer comprises at least a portion of the actuator. A shape memory polymer comprising a polymer composition which physically forms a network structure wherein the polymer composition has shape-memory behavior and can be formed into a permanent primary shape, re-formed into a stable secondary shape, and controllably actuated to recover the permanent primary shape. Polymers have optimal aliphatic network structures due to minimization of dangling chains by using monomers that are symmetrical and that have matching amine and hydroxl groups providing polymers and polymer foams with clarity, tight (narrow temperature range) single transitions, and high shape recovery and recovery force that are especially useful for implanting in the human body.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Shape memory polymers

Wilson, Thomas S.↗

DOE OSTI2022

Effect of post-extrusion heat treatment on mechanical property of aluminum alloy 2024 tube produced using Shear Assisted Processing and Extrusion (ShAPE)

The present work deals with the Shear Assisted Processing and Extrusion (ShAPE) of aluminum alloy 2024 tubes and mechanical property optimization via post extrusion heat treatments (T3510 and T8510). The ShAPE process enables extruding AA2024 at an extrusion speed of 7.4 m/min and extrusion temperature of 480 °C, which is otherwise beyond the extrudability limit of this alloy using conventional extrusion techniques. The ultimate tensile strength reaches 466 MPa for ShAPE+T3510 condition and 503 MPa for ShAPE+T8510 condition, both of which exceed the corresponding ASTM/ASM (minimum/typical) values for the studied tube geometry. The yield strength reaches 382 MPa for ShAPE+T3510 condition and 481 MPa for ShAPE+T8510 condition, both of which surpass any reported values in the literature or the highest industrial values. In addition, the ductility of AA2024 tube fabricated from ShAPE+T3510 or ShAPE+T8510 are evidently improved by the grain size refinement, second phase refinement, and enhanced uniformity in second phase distribution due to the severe plastic deformation of ShAPE process.

Shear Assisted Processing and Extrusions (ShAPE), ↗

DOE OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗

DOE OSTI2025

Revealing mechanisms of processing defect mitigation in laser powder bed fusion via shaped beams using high-speed X-ray imaging

The laser powder bed fusion (LPBF) process utilizing a focused Gaussian-shaped beam faces challenges, including pore formation, melt pool fluctuation and liquid spattering. While beam shaping technology has been explored as a potential approach for defect mitigation, the beam-matter interaction dynamics during melting with shaped beams remain unclear. Here, we report the direct observation of ring-shaped beam-matter interaction dynamics, including pore formation, melt pool fluctuation and liquid spattering, and unveil defect mitigation mechanisms in ring-shaped beam laser powder bed fusion process. Here, we find that, by spatially manipulating incident laser rays, the ring-shaped beam controls keyhole morphology, thereby managing the distribution of the reflected rays. This manipulation can effectively eliminate the formation of an unstable cavity at the keyhole tip, stabilizing the keyhole and mitigating keyhole pores. This enhanced keyhole stability effectively reduces the melt pool fluctuation, the formation of liquid breakup induced spatters and liquid droplet colliding induced large spatters in the laser powder bed fusion process. Additionally, the high-energy forefront of the ring-shaped beam effectively melts the powder bed, reducing agglomeration liquid spatter in the laser powder bed fusion process. The discovered defect mitigation mechanisms may guide the design of beam shaping strategies for simultaneously increasing the quality and productivity of metal additive manufacturing.

Beam shaping↗

DOE OSTI2022

Hard superellipse phases: particle shape anisotropy & curvature

In this work, we report computer simulations of two-dimensional convex hard superellipse particle phases vs. particle shape parameters including aspect ratio, corner curvature, and sidewall curvature. Shapes investigated include disks, ellipses, squares, rectangles, and rhombuses, as well as, shapes with non-uniform curvature including rounded squares, rounded rectangles, and rounded rhombuses. Using measures of orientational order, order parameters, and a novel stretched bond orientational order parameter, we systematically identify particle shape properties that determine liquid crystal and crystalline phases including their coarse boundaries and symmetry. We observe phases including isotropic, nematic, tetratic, plastic crystals, square crystals, and hexagonal crystals (including stretched variants). Our results catalog known benchmark shapes, but include new shapes that also interpolate between known shapes. Our results indicate design rules for particle shapes that determine two-dimensional liquid, liquid crystalline, and crystalline microstructures that can be realized via particle assembly.

42 ENGINEERING↗

DOE OSTI2026

Shape evolution in neutron-rich odd-even 105–109 Nb isotopes

Background: Neutron-rich nuclei around 𝑍 ≈ 40 are well known for exhibiting multiple shape transitions. Here, this region shows one of the sharpest shape transitions in the nuclear chart, evolving from a spherical vibrator at 𝑁 = 58 to a strongly deformed prolate shape at 𝑁 = 60. The largest deformations are observed for 38 Sr and 40 Zr . This abrupt shape transition disappears at 𝑍 = 36 and below, where a shape transition from spherical to oblate nuclei is predicted. On the other hand, for 𝑍 ≥ 42 and 𝑁 ≥ 60, the shape is known to evolve from axial to triaxial. While the even-𝑍 nuclei in this region have already been extensively studied, new insights can be gained from the studies of odd-𝑍 isotopes for a better understanding of the underlying mechanisms driving these phenomena. Purpose: The 41 Nb nuclei lie at the boundary between axially deformed Zr and triaxially deformed Mo nuclei. This work investigates the nuclear structure of very neutron-rich Nb nuclei up to 𝑁 = 68. The goal is to understand how the nuclear shape evolves as a function of isospin in this isotopic chain and provide new insights into the emergence of triaxial deformation. Methods: The structure of the neutron-rich Nb isotopes was investigated using state-of-the-art high-resolution 𝛾-ray spectroscopy of fission fragments produced via two different fission reactions. The use of 9 Be ⁢( 238 U, 𝑓) inverse kinematics, with a detection system comprising AGATA, EXOGAM, and VAMOS++, enabled the measurement of prompt and delayed 𝛾 rays from isotopically identified fission fragments, and 𝛾−𝛾−𝛾−𝛾 highfold data were obtained from a spontaneous fission source of 252 Cf using the Gammasphere array. Results: The level scheme of 105 Nb has been significantly extended, with the addition of two negative-parity bands observed for the first time. A new level scheme is proposed for 107 Nb, which is not in agreement with an earlier measurement, and new levels and transitions have been added to the very neutron-rich 109 Nb. The degree of triaxiality of the new bands is discussed on the basis of signature splitting analysis. The recently reported level scheme of 99 Nb has been revised. Conclusions: This systematic study on the Nb isotopic chain, compared to Zr and Mo, indicates that while the ground-state band exhibits a triaxial deformation, attributed to a proton hole coupled to a triaxially deformed Mo core, the negative-parity bands, based on isomeric bandheads, display an axially symmetric deformed structure, similar to that observed in the Zr isotopes, revealing the existence of a shape coexistence in the neutron-rich Nb nuclei.

Abushawish, M. [Université Claude-Bernard Lyon 1 (↗

DOE OSTI2023

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗

DOE OSTI2022

Permeability simulation of kidney-bean shaped carbon fibers

Cost saving efforts in carbon fiber manufacturing may result in irregularly shaped cross-sections, often as a kidney-bean shape. Here, in this research, flow simulations were performed to determine transverse permeability of these kidney-bean shaped fibers. Three different unit-cell packing arrangements for the kidney-bean shaped fibers and unit-cell orientations of 0°, + 90°, and – 90° were considered over the fiber volume fraction range 0.30 – 0.75. The kidney-bean shaped fibers had two distinct hexagonally packed unit-cell geometries because of the fiber’s asymmetry with different packing efficiency. The maximum fiber packing fraction for kidney-bean shaped fibers was 0.923, which was larger than the 0.907 for circular fibers. When compared to circular fibers with hexagonal packing, the kidney-bean shaped fibers had up to a 74% reduction in permeability. The ±90° orientations for KB-fibers with hexagonal packing had an infiltration time up to 3.86 longer than the hexagonal packed circular fibers for the same volume fraction.

36 MATERIALS SCIENCE↗

DOE OSTI2025

X-Ray Microtomography Measurements of Coevolving Particle Size and Shape

This study measures simultaneous changes in particle size and shape during grain crushing and proposes quantitative metrics to track their coevolving statistical distribution. Idealized granular materials with widely different initial particle shapes (i.e., spheres, quasi-ellipsoids, plates, and rods) are compressed oedometrically at vertical stress levels sufficient to induce particle breakage. X-ray microtomography and digital image analysis are then used to track changes in particle-scale characteristics. The observations indicate that more irregular particles tend to exhibit a higher degree of breakage compared to more spherical grains subjected to an equivalent stress level. Despite these differences, the measurements also indicate that all the tested particle sets approach a common attractor (expressed in terms of aspect ratio, flatness, and elongation), regardless of their initial morphology. A new shape evolution index is proposed to quantify such trends. With this new metric, it is shown that in all the tested materials, the particle shape evolves faster than size, especially for the more regular morphologies, thus implying that the shape tends to reach its attractor point before the grading reaches its ultimate particle size distribution. Furthermore, this finding underpins an ultimate stage of crushing-driven compression that occurs with particle size reduction but minimal shape changes (i.e., self-replicating shapes due to fractures).

Oedometric compression↗

DOE OSTI2021

Solid-Density Ion Temperature from Redshifted and Double-Peaked Stark Line Shapes

Heβ spectral line shapes are important for diagnosing temperature and density in many dense plasmas. This work presents Heβ line shapes measured with high spectral resolution from solid-density plasmas with minimized gradients. The line shapes show hallmark features of Stark broadening, including quantifiable redshifts and double-peaked structure with a significant dip between the peaks; these features are compared to models through a Markov chain Monte Carlo framework. Line shape theory using the dipole approximation can fit the width and peak separation of measured line shapes, but it cannot resolve an ambiguity between electron density n e and ion temperature T i , since both parameters influence the strength of quasistatic ion microfields. Here a line shape model employing a full Coulomb interaction for the electron broadening computes self-consistent line widths and redshifts through the monopole term; redshifts have different dependence on plasma parameters and thus resolve the n e – T i ambiguity. The measured line shapes indicate densities that are 80–100% of solid, identifying a regime of highly ionized but well-tamped plasma. This analysis also provides the first strong evidence that dense ions and electrons are not in thermal equilibrium, despite equilibration times much shorter than the duration of x-ray emission; cooler ions may arise from nonclassical thermalization rates or anomalous energy transport. Furthermore, the experimental platform and diagnostic technique constitute a promising new approach for studying ion-electron equilibration in dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Discrete element modeling of granular hopper flow of irregular-shaped deformable particles

Many natural and engineered granular materials have relatively deformable particles. Besides particle size and shape, particle deformability is another salient factor that significantly impacts the material’s flow behavior. Here, in this work, the flow of irregular-shaped deformable particles in a wedge-shaped hopper is investigated using discrete element simulations. A bonded-sphere model is developed to simultaneously capture irregular particle shapes and particle-wise deformations (e.g., compression, deflection, and distortion). Quantitative analysis of the effects of irregular shapes and particle deformations shows that the increase in particle stiffness tends to increase initial packing porosity and decrease the flow rate in the hopper. Rigid particles tend to have clogging issues, whereas deformable particles have less chance to, indicating particle deformation reduces the critical bridging width in the hopper flow. Detailed analysis of stress fields is also conducted to provide insights into the mechanism of particle flow and clogging. Stresses and discharge rates calculated from numerical simulations are compared and show good agreement with Walker’s theory and the extended Beverloo formula. Simulations with various particle shape combinations are also performed and show that the initial packing porosity decreases with an increasing percentage of fibers while the discharge rate has a complex dependency on particle shapes.

09 BIOMASS FUELS↗