SEARCH · Engineering Papers

Results for “smoothing technique”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Computation of Direct Sensitivities of Spatial Multibody Systems With Joint Friction

Abstract Friction exists in most mechanical systems and may have a major influence on the dynamic performance of the system. The incorporation of friction in dynamic systems has been a subject of active research for several years owing to its high nonlinearity and its dependence on several parameters. Consequently, optimization of dynamic systems with friction becomes a challenging task. Gradient-based optimization of dynamical systems is a prominent technique for optimal design and requires the computation of model sensitivities with respect to the design parameters. The novel contribution of this paper is the derivation of the analytical methodology for the computation of direct sensitivities for smooth multibody systems with joint friction using the Lagrangian index-1 formulation. System dynamics have been computed using two different friction models; the Brown and McPhee, and the Gonthier et al. model. The methodology proposed to obtain model sensitivities has also been validated using the complex finite difference method. A case study has been conducted on a spatial multibody system to observe the effect of friction on the dynamics and model sensitivities, compare sensitivities with respect to different parameters and demonstrate the numerical and validation aspects. Since design parameters can have very different magnitudes and units, the sensitivities have been scaled with the parameters for comparison. Finally, a discussion has been presented on the interpretation of the case study results. Due to the incorporation of joint friction, ‘jumps’ or discontinuities are observed in the model sensitivities akin to those observed for hybrid dynamical systems.

Engineering↗

DOE OSTI2023

Influence of Realistic, Cyclic Atmospheric Cycles on the Pitting Corrosion of Austenitic Stainless Steels

Pitting corrosion was evaluated on stainless steels 304H, 304, and 316L the surfaces of which had ASTM seawater printed on them as a function of surface roughness after exposure to an exemplar realistic atmospheric diurnal cycle for up to one year. Methods to evaluate pitting damage included optical imaging, scanning electron microscopy imaging, profilometry analysis, and polarization scans. The developed cyclic exposure environment did not significantly influence pitting morphology nor depth in comparison to prior static exposure environments. Cross-hatching was observed in a majority of pits for all material compositions with the roughest surface finish (#4 finish) and in all surface finishes for the 304H composition. Evidence is provided that cross-hatched pit morphologies are caused by slip bands produced during the grinding process for the #4 finish or by material processing. Additionally, micro-cracking was observed in pits formed on samples with the #4 surface finish and was greatly reduced or absent for pits formed on samples with smooth surface finishes. This suggests that both a low RH leading to an MgCl 2 -dominated environment and a rough surface containing significant residual stress are necessary for micro-cracking. Finally, the use of various characterization techniques and cross sectioning was employed to both qualitatively and quantitatively assess pitting damage across all SS compositions and surface finishes.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Functional Data Analysis in Wearable Body Sensor Networks

Improving response time of indirect room-size calorimeters is still an outstanding problem in metabolic research. Accurate estimates of instantaneous rates of gaseous exchange require numerical differentiation of measured gaseousgas concentrations. We propose a new method to estimate the instantaneous gaseousgas exchange rates in indirect calorimetry. In contrast to the previously developed techniques, the method addresses the problem of differentiation of gaseous concentrations as an ill-posed problem. By applying the method of regularization, the problem of differentiation is converted into a well-posed problem resulting in smooth and consistent gaseous exchange rates. The validity of the method is tested on a large dataset of calorimeter experiments which included 313 human experiments along with 231 alcohol combustion experiments. It is demonstrated that the method is able to reliably differentiate between the “unphysiological” process of alcohol combustion and physiological variations produced by human metabolism. The method also allowed unraveling the previously unreported relative kinetics of O2 consumption and Respiratory Quotient (RQ) in humans. It was found that the kinetics of oxidative fuel selection lags behind the energy expenditure in humans exhibiting some sort of oxidative inertia. The time lag varies from 2-3 min up to 30 min, depending on particular individual. No such lag was found in alcohol combustion experiments. In addition to the relative kinetics of substrate oxidation, two statistical indexes reflecting variability of minute-by-minute RQ were estimated. The indexes were the RQ’s standard deviation and RQ’s first-order derivative. Both indexes showed statistically significant difference between human experiments and alcohol combustion experiments. We conclude that the proposed method can consistently extract physiologically-relevant information from noisy calorimetry data and the aforesaid information can provide additional insights into the mechanism of metabolic fuel selection in humans.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Functional Data Analysis in Wearable Body Sensor Networks

60 - APPLIED LIFE SCIENCES↗

DOE OSTI2022

Experimental and Computational Investigation of Integrated Internal and Film Cooling Designs Incorporating a Thermal Barrier Coating

Abstract Few studies in the open literature have studied the effect of thermal barrier coatings (TBC) when used in combination with shaped hole film cooling and enhanced internal cooling techniques. The current study presents Reynolds-averaged Navier–Stokes (RANS) conjugate heat transfer simulations that identify trends in cooling design performance as well as experimental measurements of overall effectiveness using a flat-plate matched-Biot number model with a simulated TBC layer of 0.42D thickness, where D is the film cooling hole diameter. Coolant is fed to the film cooling holes in a co-flow configuration, and the results of both the smooth and rib-turbulated channels are compared. At a constant coolant flow rate, enhanced internal cooling was found to provide a 44% increase in spatially-averaged overall effectiveness, $\bar{\bar{\phi}}$, without a TBC. The results show that the addition of a TBC can raise $\bar{\bar{\phi}}$ on a film-cooled component surface by 47%. The optimum velocity ratio was found to decrease with the addition of enhanced cooling techniques and a TBC as the film provided minimal benefit at the expense of reduced internal cooling. While the computational results closely identified trends in overall system performance without a TBC, the model over-predicted effectiveness on the metal–TBC interface. The results of this study will inform turbine component design as material science advances increase the reliability of the TBC.

Engineering↗

DOE OSTI2020

Heat Transfer Coefficients of Additively Manufactured Tubes with Internal Pin Fins for Supercritical Carbon Dioxide Cycle Recuperators

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in additively manufactured tubes with internal pin fins at the US DOE’s National Energy Technology Laboratory in Morgantown, WV. The measurement procedures were validated by conducting benchmark tests with smooth stainless-steel tube and comparing the results with published correlations for Nusselt number (Nu) and friction factor. Over Reynolds numbers (ReD) ranging from 5 × 10 4 to 2.5 × 10 5 , measured Nu was within 5% of the Dittus-Boelter correlation and measured friction factors were within 5% of the McAdams correlation for smooth tube flow. The candidate pin fin patterned pipes were additively manufactured (AM) at the Oak Ridge National Laboratory. The pins were circular or elliptical in cross-section and were printed at a 30⁰ angle relative to the inner wall (to meet AM constraints). The pin arrangement was helical to promote enhanced heat transfer due to swirl flow. Pin length to diameter aspect ratio was 1.33, 2, and 8, while the pin diameter to tube diameter ratio was 0.188, 0.125, and 0.063. Tests were performed for ReD varying from 6.9 × 10 4 to 2.2 × 10 5 and at conditions equivalent to the low pressure (LP) outlet (8.69 MPa, 361 K) and the high pressure (HP) inlet (20.7 MPa, 350 K) of the low temperature recuperator (LTR) in an indirect sCO 2 cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients. For the best performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.19), the local heat transfer coefficient increased by 136% relative to the Dittus-Boelter correlation at the LTR low pressure outlet and 194% at the LTR high pressure inlet. These correspond to a 282% and a 271% increase in the product of the heat transfer coefficient and surface area (adjusted for fin efficiency) product, respectively. Large pressure drops across the test articles were observed. For Tube Design A, the average friction factor, across the range of ReD considered, was significantly larger than the McAdams correlation at both the LTR LP outlet and the LTR HP inlet. A thermal performance factor was utilized to express the ratio of material required to build a finned heat exchanger relative to a finless heat exchanger with the same heat duty and pumping power. Tube Design A was estimated to decrease the required heat exchanger material by 13%.

Searle, Matthew↗

DOE OSTI2021

Production of desensitized, ultrafine PETN powder

Pentaerythritol tetranitrate (PETN) is a widely studied high explosive (HE), most commonly used in detonator applications. In this work, we use a spray drying technique to manufacture a new “nano-powder” morphology of PETN. The spray dried PETN is several orders of magnitude smaller in particle size than traditionally prepared PETN powders. The spray dried PETN has a mostly spherical and smooth morphology, while traditional crystalline PETN typically has sharp edges and crystal facets. Small-scale sensitivity tests including drop-weight impact, friction, and electrostatic discharge (ESD) indicate the spray dried PETN is less sensitive than traditionally used forms of PETN powders. Furthermore, we also observed no changes in chemical properties (melt temperature and onset of decomposition) after spray drying, and the material remains in the tetragonal, crystalline phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Unlocking hidden information in sparse small-angle neutron scattering measurements

Hypothesis Small-Angle Neutron Scattering (SANS) is a powerful technique for studying soft matter systems such as colloids, polymers, and lyotropic phases, providing nanoscale structural insights. However, its effectiveness is limited by low neutron flux, leading to long acquisition times and noisy data. Here, we hypothesize that Bayesian statistical inference using Gaussian Process Regression (GPR) can reconstruct high-fidelity scattering data from sparse measurements by leveraging intensity smoothness and continuity. Experiments and Simulations The method was benchmarked computationally and validated through SANS experiments on various soft matter systems, including wormlike micelles, colloidal suspensions, polymeric structures, and lyotropic phases. GPR-based inference was applied to both experimental and synthetic data to evaluate its effectiveness in noise reduction and intensity reconstruction. Findings GPR significantly enhances SANS data quality and therefore reducing measurement times by up to two orders of magnitude. This cost-effective approach maximizes experimental efficiency, enabling high-throughput studies and real-time monitoring of dynamic systems. It is particularly beneficial for weakly scattering and time-sensitive studies. Beyond SANS, this framework applies to other low-SNR techniques, including laboratory-based small-angle X-ray scattering and various dynamical scattering methods. Furthermore, it offers transformative potential for compact neutron sources, enhancing their viability for structural analysis in resource-limited settings.

Small angle neutron scattering↗

DOE OSTI2022

Nontrivial Fermi surface topology of the kagome superconductor CsV 3 Sb 5 probed by de Haas–van Alphen oscillations

Here, we have investigated the normal state Fermi-surface properties of the kagome superconductor CsV 3 Sb 5 using torque magnetometry with applied fields (H) up to 35 T. The torque signal shows clear de Haas-van Alphen (dHvA) oscillations above 15 T. The oscillations are smooth and consist of seven distinct frequencies with values from similar to ~18 T to 2135 T. The presence of higher frequencies in CsV 3 Sb 5 is further confirmed by carrying out additional measurements using the tunnel diode oscillator technique. All frequencies measured at different tilt angles (θ) of the field direction with respect to the c axis show a 1/cosθ dependence, implying that the Fermi surfaces corresponding to these frequencies are two dimensional (2D). The absence of dHvA oscillations at θ = 90° further supports the presence of 2D Fermi surfaces. The Berry phase (φ) determined from the Landau level fan diagram for all frequencies is similar to ~ 0.4. This value is close to the theoretical value of φ = 0.5 for a nontrivial system, which strongly supports the nontrivial topology of the Fermi surfaces of these frequencies. Several quantities characterizing the Fermi surface are calculated employing the Lifshitz-Kosevich theory. These findings are crucial for exploring the interplay between nontrivial band topology, charge-density wave, and unconventional superconductivity of CsV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Neural simulation-based inference approach for characterizing the Galactic Center γ-ray excess

The nature of the Fermi γ-ray Galactic Center Excess (GCE) has remained a persistent mystery for over a decade. Although the excess is broadly compatible with emission expected due to dark matter annihilation, an explanation in terms of a population of unresolved astrophysical point sources e.g., millisecond pulsars, remains viable. The effort to uncover the origin of the GCE is hampered in particular by an incomplete understanding of diffuse emission of Galactic origin. This can lead to spurious features that make it difficult to robustly differentiate smooth emission, as expected for a dark matter origin, from more “clumpy” emission expected from a population of relatively bright, unresolved point sources. Here, we use recent advancements in the field of simulation-based inference, in particular density estimation techniques using normalizing flows, in order to characterize the contribution of modeled components, including unresolved point source populations, to the GCE. Compared to traditional techniques based on the statistical distribution of photon counts, our machine learning-based method is able to utilize more of the information contained in a given model of the Galactic Center emission, and in particular can perform posterior parameter estimation while accounting for pixel-to-pixel spatial correlations in the γ-ray map. This makes the method demonstrably more resilient to certain forms of model misspecification. On application to Fermi data, the method generically attributes a smaller fraction of the GCE flux to unresolved point sources when compared to traditional approaches. We nevertheless infer such a contribution to make up a non-negligible fraction of the GCE across all analysis variations considered, with at least $38^{+9}_{-19}$% of the excess attributed to unresolved point sources in our baseline analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Geranium Photodiodes for Hard Xray Detection

This report summarizes the investigations engaged at UC Davis under subcontract B632083, as related to the LLNL-LDRD Project “High-Speed X-Ray Imager Arrays.” The first section entails work done within the Electrical and Computer Engineering (ECE) Department. The purpose of the work done so far is to understand photodiode design and which parameters are conducive to the design goals, which are: I: Maximize the detector quantum efficiency. The device must have a quantum efficiency of roughly 5% at 80 keV. II: Minimize temporal response. Collection time is limited to 1ns. All charge must be collected in that time. III: Minimize dark current. In this document, various photodiode parameters and designs are explored, and their impact on these goals are tabulated. In addition to these, various instruments that have been constructed for either testing these parameters or helping teammates with their experiments. The second section entails work done with the Materials Science Engineering (MSE) Department. In accordance with the proposed development of the germanium ROIC backside imagers, a critical step in the micromachining of these devices was the removal of the substrate on which the PIN junctions were epitaxially grown. A common technique for the micromachining of silicon wafers is the use of chemical etching alongside an “etch-stop” layer. In this technique, a chemical etchant will remove material until it reaches the etch-stop layer, at which the etch rate will drop significantly. This high selectivity (ratio of etch rates) allows any unevenness in the initial etching to be effectively smoothed away in the event that part of the wafer surface reaches the etch-stop layer sooner than another. An initial literature review found a paper by Divan et al. in which the authors used ion implantation to dope germanium wafers and achieve selectivity on the order of 100, sufficient for a reliable etch-stop.

42 ENGINEERING↗

DOE OSTI2025

Reliable and Efficient Machine Learning (Final Technical Report)

Modern scientific experiments generate massive amounts of data at a pace much faster than humans can manually analyze. While machine learning has revolutionized commercial data analysis (such as recommending movies or recognizing faces), applying these tools to complex scientific discovery is challenging because scientific answers must be precise, interpretable, and adhere to physical laws. The research under this project aims to develop new mathematical tools and computer algorithms specifically designed for scientific applications. Major progress has been made in automatically cleaning and deconstructing messy experimental data, analyzing the visual information of physical phenomena, determining the underlying physical variables, and providing rig orous mathematical analysis of interesting algorithms and concepts widely used in machine learning. This project addressed the critical gap between our ability to generate massive scientific data and our ability to extract interpretable information from it. We established mathematical foundations for Scientific Machine Learning (SciML) aimed at effective data analytics and automated discovery. Our work focused on three core objectives: (1) developing reliable feature extraction methods for dynamic high-dimensional data, (2) establishing mathematical foundations for discovering dynamics via neural networks, and (3) creating rigorous optimization techniques for these models. Key outcomes come from two fronts. On the practical side, they include the development of algorithms that significantly enhance the extraction of signals from field data, as well as the capability to handle situations that exhibit smooth variations or physical stretching due to temperature changes. They also include the creation of an automated framework for discovering fundamental state variables from raw experimental data, demonstrating the ability to identify intrinsic physical dimensions without prior knowledge of the governing laws. On the theoretical front, the research results in theoretical advances in Optimal Transport, a widely used notion in SciML, specifically regarding functions with fixed-size nodal sets, provide sharp bounds relevant to uncertainty quantification. Meanwhile, the outcomes also include the establishment of convergence theories for nonlocal gradient descent methods, enabling robust optimization with noisy data in high-dimensional settings commonly encountered in scientific modeling. The project also helps creating opportunities to train the next generation of researchers, equipping them with the necessary technical skills for today’s workplace and preparing them for future advances.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Observing signals of spectral features in the cosmic-ray positrons and electrons from Milky Way pulsars

The Alpha Magnetic Spectrometer (AMS-02) has provided unprecedented precision measurements of the electron and positron cosmic-ray fluxes and the positron fraction spectrum. At the higher energies, sources as energetic local pulsars, may contribute to both cosmic-ray species. The discreteness of the source population, can result in features both on the positron fraction measurement and in the respective electron and positron spectra. For the latter, those would coincide in energy and would contrast predictions of smooth spectra as from particle dark matter. In this work, using a library of pulsar population models for the local part of the Milky Way, we perform a power-spectrum analysis on the cosmic-ray positron fraction. We also develop a technique to cross-correlate the electron and positron fluxes. We show that both such analyses, can be used to search statistically for the presence of spectral wiggles in the cosmic-ray data. For a significant fraction of our pulsar simulations, those techniques are already sensitive enough to give a signal for the presence of those features above the regular noise, with forthcoming observations making them even more sensitive. Lastly, by cross-correlating the AMS-02 electron and positron spectra, we find an intriguing first hint for a positive correlation between them, of the kind expected by a population of local pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Toward Ideal Metal–Organic Framework Thin-Film Growth via Automated Layer-by-Layer Deposition: Examples Based on Perylene Diimide Linkers

Photoactive metal-organic framework (MOF) thin films offer an opportunity for translating the advantages of periodic, crystalline, and tailorable light-harvesting materials directly into devices such as those for photoelectrochemical solar energy conversion. In this study, we report the fabrication of light-absorbing perylene-diimide-containing pillared-paddlewheel MOF thin films using an automated layer-by-layer (LbL) deposition technique. Our focus here is on optimizing the growth of representative chromophoric MOFs as oriented films of uniform and predefined thickness. Growth was examined as a function of metal identity, pillaring ligand composition, and supporting-surface chemical functionality. Further, application of atomic force microscopy (AFM) and complementary techniques revealed that the surface-supported MOFs initially display island-type film growth (Volmer-Weber growth), resulting in comparatively rough films. Further growth is accompanied by the merging of islands, resulting in films that, depending on experimental details, can be remarkably smooth (i.e., roughness on the order of ±1 nm (one structural repeat unit in the pillaring direction)). These details include the use of 1,4-diazabicyclo[2.2.2]octane (DABCO) as a MOF pillar and ALD-grown zinc oxide as a film support (ALD = atomic layer deposition). Also helpful for mitigating island-type growth, at least in part, is the replacement of Zn 2+ by Cu 2+ as the metal component of the MOF. Notably, each of these adjustments entails replacing weaker chemical bonds with stronger ones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Controlling Preferred Grain Orientation in Vapor-Deposited Metal-Halide Perovskite Thin Films

Metal-halide perovskites are promising semiconductor materials for light-absorber layers in solar cells due to their optimized band gap, defect tolerance, and relative ease of synthesis. In this work, we demonstrate the viability of vapor transport deposition (VTD), a solvent-free processing method, to realize polycrystalline films of both lead iodide (PbI 2 ) and methylammonium lead iodide (MAPbI 3 ) perovskites showing preferential crystallographic orientation. Notable is the difference in perovskite film texturing observed for sequentially and co-deposited VTD films of MAPbI 3 , as well as the lack of strong preferential orientation in solution-processed MAPbI 3 films. While sequential deposition of the perovskite precursors PbI 2 and methylammonium iodide (MAI) leads to a MAPbI 3 film with a limited preferential grain orientation, co-deposited MAPbI 3 films mirror the orientation of excess PbI 2 present in the film. The latter is confirmed by the similar dependence of preferential orientation for PbI 2 and co-deposited MAPbI 3 on underlayer roughness, as well as the emergence of preferential orientation in MAPbI 3 only when excess PbI 2 is present. Both MAPbI 3 and PbI 2 films show a strong preferred lattice orientation on smooth substrates, while on rougher substrates, this orientation is disrupted. Furthermore, these results are interpreted as excess PbI 2 frustrating the diffusion of MAI in the film, which would otherwise disrupt the preferred orientation, pointing to a notable difference in lattice orientation arising from processing technique and conditions.

Deposition↗

DOE OSTI2022

3ω techniques for measurement of volumetric heat capacity and anisotropic thermal conductivity of a solution processable, hybrid organic/inorganic film, Te-PEDOT:PSS

We report that measuring the thermal properties of anisotropic films of hybrid materials poses a challenge to existing metrology techniques. We have developed a new approach for measuring the volumetric heat capacity and anisotropic thermal conductivity of these systems using the 3ω method. While there exist many avenues for measuring the thermal properties of thin films, most carry with them difficult requirements such as smooth surfaces or advanced lithography. Here, we present measurements of a film's in-plane and cross-plane conductance and its volumetric heat capacity using relatively simple sample configurations, each requiring a single heater. For the measurement of volumetric heat capacity, we present a new model fitting method, relying on a standard film-on-substrate configuration. For the measurement of in-plane thermal conductance by 3ω, we have developed the use of an embedded micro-wire heater in suspended drop cast films, allowing for a 12 μm wide heater without the need for advanced lithography. We also expose the surprisingly significant effect of thermal radiation in the suspended film measurement and its associated error. Our measurements reveal a large anisotropy in the thermal conductivity of our test material, Te-PEDOT:PSS, of k in-plane / k cross-plane = 19, consistent with the nanoscale morphology of the material.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A Minkowski difference-based advancing front packing technique for generating convex noncircular particles in complex domains

In this work, a Minkowski difference-based advancing front approach is proposed to generate convex and non-circular particles in a predefined computational domain. Two specific algorithms are developed to handle the contact conformity of generated particles with the boundaries of the computational domain. The first, called the open form, is used to handle the smooth contact of generated particles with (external) boundaries, while the other, called the closed form, is proposed to handle the internal boundaries of a computational domain with a complex cavity. The Gilbert-Johnson-Keerthi (GJK) method is used to efficiently solve the contact detection between the newly generated particle at the front and existing particles. Furthermore, the problem of one-sided particle lifting, which can cause some defects in the packing structure in existing advancing front methods during packing generation, is highlighted and an effective solution is developed. Several examples of increasing complexity are used to demonstrate the efficiency and applicability of the proposed packing generation approach. The numerical results show that the generated packing is not only more uniform, but also achieves a higher packing density than existing advancing front methods.

42 ENGINEERING↗