Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Compact model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 649 records · Page 36

Bubble Combustion

A method of energy production that is capable of low pollutant emissions is fundamental to one of the four pillars of NASA s Aeronautics Blueprint: Revolutionary Vehicles. Bubble combustion, a new engine technology currently being developed at Glenn Research Center promises to provide low emissions combustion in support of NASA s vision under the Emissions Element because it generates power, while minimizing the production of carbon dioxide (CO2) and nitrous oxides (NOx), both known to be Greenhouse gases. and allows the use of alternative fuels such as corn oil, low-grade fuels, and even used motor oil. Bubble combustion is analogous to the inverse of spray combustion: the difference between bubble and spray combustion is that spray combustion is spraying a liquid in to a gas to form droplets, whereas bubble combustion involves injecting a gas into a liquid to form gaseous bubbles. In bubble combustion, the process for the ignition of the bubbles takes place on a time scale of less than a nanosecond and begins with acoustic waves perturbing each bubble. This perturbation causes the local pressure to drop below the vapor pressure of the liquid thus producing cavitation in which the bubble diameter grows, and upon reversal of the oscillating pressure field, the bubble then collapses rapidly with the aid of the high surface tension forces acting on the wall of the bubble. The rapid and violent collapse causes the temperatures inside the bubbles to soar as a result of adiabatic heating. As the temperatures rise, the gaseous contents of the bubble ignite with the bubble itself serving as its own combustion chamber. After ignition, this is the time in the bubble s life cycle where power is generated, and CO2, and NOx among other species, are produced. However, the pollutants CO2 and NOx are absorbed into the surrounding liquid. The importance of bubble combustion is that it generates power using a simple and compact device. We conducted a parametric study using CAVCHEM, a computational model developed at Glenn, that simulates the cavitational collapse of a single bubble in a liquid (water) and the subsequent combustion of the gaseous contents inside the bubble. The model solves the time-dependent, compressible Navier-Stokes equations in one-dimension with finite-rate chemical kinetics using the CHEMKIN package. Specifically, parameters such as frequency, pressure, bubble radius, and the equivalence ratio were varied while examining their effect on the maximum temperature, radius, and chemical species. These studies indicate that the radius of the bubble is perhaps the most critical parameter governing bubble combustion dynamics and its efficiency. Based on the results of the parametric studies, we plan on conducting experiments to study the effect of ultrasonic perturbations on the bubble generation process with respect to the bubble radius and size distribution.

Corrigan, Jackie↗

NASA Tech Briefs, March 2010

Topics covered include: Software Tool Integrating Data Flow Diagrams and Petri Nets; Adaptive Nulling for Interferometric Detection of Planets; Reducing the Volume of NASA Earth-Science Data; Reception of Multiple Telemetry Signals via One Dish Antenna; Space-Qualified Traveling-Wave Tube; Smart Power Supply for Battery-Powered Systems; Parallel Processing of Broad-Band PPM Signals; Inexpensive Implementation of Many Strain Gauges; Constant-Differential-Pressure Two-Fluid Accumulator; Inflatable Tubular Structures Rigidized with Foams; Power Generator with Thermo-Differential Modules; Mechanical Extraction of Power From Ocean Currents and Tides; Nitrous Oxide/Paraffin Hybrid Rocket Engines; Optimized Li-Ion Electrolytes Containing Fluorinated Ester Co-Solvents; Probabilistic Multi-Factor Interaction Model for Complex Material Behavior; Foldable Instrumented Bits for Ultrasonic/Sonic Penetrators; Compact Rare Earth Emitter Hollow Cathode; High-Precision Shape Control of In-Space Deployable Large Membrane/Thin-Shell Reflectors; Rapid Active Sampling Package; Miniature Lightweight Ion Pump; Cryogenic Transport of High-Pressure-System Recharge Gas; Water-Vapor Raman Lidar System Reaches Higher Altitude; Compact Ku-Band T/R Module for High-Resolution Radar Imaging of Cold Land Processes; Wide-Field-of-View, High-Resolution, Stereoscopic Imager; Electrical Capacitance Volume Tomography with High-Contrast Dielectrics; Wavefront Control and Image Restoration with Less Computing; Polarization Imaging Apparatus; Stereoscopic Machine-Vision System Using Projected Circles; Metal Vapor Arcing Risk Assessment Tool; Performance Bounds on Two Concatenated, Interleaved Codes; Parameterizing Coefficients of a POD-Based Dynamical System; Confidence-Based Feature Acquisition; Algorithm for Lossless Compression of Calibrated Hyperspectral Imagery; Universal Decoder for PPM of any Order; Algorithm for Stabilizing a POD-Based Dynamical System; Mission Reliability Estimation for Repairable Robot Teams; Processing AIRS Scientific Data Through Level 3; Web-Based Requesting and Scheduling Use of Facilities; AutoGen Version 5.0; Time-Tag Generation Script; PPM Receiver Implemented in Software; Tropospheric Emission Spectrometer Product File Readers; Reporting Differences Between Spacecraft Sequence Files; Coordinating "Execute" Data for ISS and Space Shuttle; Database for Safety-Oriented Tracking of Chemicals; Apparatus for Cold, Pressurized Biogeochemical Experiments; Growing B Lymphocytes in a Three-Dimensional Culture System; Tissue-like 3D Assemblies of Human Broncho-Epithelial Cells; Isolation of Resistance-Bearing Microorganisms; Oscillating Cell Culture Bioreactor; and Liquid Cooling/Warming Garment.

Source record↗

The Population of Galactic Centre Filaments – III. Candidate Radio and Stellar Sources

Recent MeerKAT radio continuum observations of the Galactic Centre at 20 cm show a large population of non-thermal radio filaments (NRFs) in the inner few hundred pc of the Galaxy. We have selected a sample of 57 radio sources, mainly compact objects, in the MeerKAT mosaic image that appear to be associated with NRFs. The selected sources are about four times the number of radio point sources associated with filaments than would be expected by random chance. Furthermore, an apparent correlation between bright IR stars and NRFs is inferred from their similar latitude distributions, suggesting that they both co-exist within the same region. To examine if compact radio sources are related to compact IR sources, we have used archival 2MASS, and Spitzer data to make spectral energy distribution of individual stellar sources coincident or close to radio sources. We provide a catalogue of radio and IR sources for future detailed observations to investigate a potential three-way physical association between NRFs, compact radio and IR stellar sources. This association is suggested by models in which NRFs are cometary tails produced by the interaction of a large-scale nuclear outflow with stellar wind bubbles in the Galactic Centre.

plasmas↗

Long-term 3D MHD simulations of black hole accretion discs formed in neutron star mergers

ABSTRACT We examine the long-term evolution of accretion tori around black hole (BH) remnants of compact object mergers involving at least one neutron star, to better understand their contribution to kilonovae and the synthesis of r-process elements. To this end, we modify the unsplit magnetohydrodynamic (MHD) solver in flash 4.5 to work in non-uniform three-dimensional spherical coordinates, enabling more efficient coverage of a large dynamic range in length scales while exploiting symmetries in the system. This modified code is used to perform BH accretion disc simulations that vary the initial magnetic field geometry and disc compactness, utilizing a physical equation of state, a neutrino leakage scheme for emission and absorption, and modelling the BH’s gravity with a pseudo-Newtonian potential. Simulations run for long enough to achieve a radiatively inefficient state in the disc. We find robust mass ejection with both poloidal and toroidal initial field geometries, and suppressed outflow at high disc compactness. With the included physics, we obtain bimodal velocity distributions that trace back to mass ejection by magnetic stresses at early times, and to thermal processes in the radiatively inefficient state at late times. The electron fraction distribution of the disc outflow is broad in all models, and the ejecta geometry follows a characteristic hourglass shape. We test the effect of removing neutrino absorption or nuclear recombination with axisymmetric models, finding ${\sim} 50{{\ \rm per\ cent}}$ less mass ejection and more neutron-rich composition without neutrino absorption, and a subdominant contribution from nuclear recombination. Tests of the MHD and neutrino leakage implementations are included.

Fahlman, Steven (ORCID:0000000348759940)↗

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗

Improved Treatment of Snow Over Ice Sheets in the NASA GISS Climate Model: Towards Ice Sheet–Climate Coupling

Representing the interactions between ice sheets and climate is essential for more accurate prediction of climate change and sea level rise. Ice sheets interact with the overlying atmosphere via the accumulation of snow and its compaction into firn, then ice, as well as the melting of surface snow and ice and the creation of runoff water. Getting an adequate representation of heat transfer, compaction, and melting processes is essential for an accurate representation of snow on land ice in global climate models. We are implementing an improved snow model on top of land ice as part of an effort to couple the NASA GISS climate model with the PISM ice sheet model. The new snow model includes additional layers and processes that are not currently incorporated (e.g., liquid water retention, percolation and refreezing, and snow densification), and mass and energy transfer methods that are consistent with both static ice sheets (with implicit iceberg fluxes) and interactive ice sheets (with explicit dynamics). We are tuning the densification scheme of this snow model with temperature and density data from common FirnCover and SumUp observations at locations in the accumulation zone of Greenland, and we compare the resulting density profiles to other SumUp density profiles in Greenland and Antarctica. We will assess the impact of this new snow model in climate model simulations with a static ice sheet compared with the previous (simpler) 2-layer snow model. Finally, we aim to use the non-coupled simulations as a baseline to assess the impact of dynamic coupling with an interactive ice sheet model.

snow↗

SymbolNet: neural symbolic regression with adaptive dynamic pruning for compression

Abstract Compact symbolic expressions have been shown to be more efficient than neural network (NN) models in terms of resource consumption and inference speed when implemented on custom hardware such as field-programmable gate arrays (FPGAs), while maintaining comparable accuracy (Tsoi et al 2024 EPJ Web Conf. 295 09036). These capabilities are highly valuable in environments with stringent computational resource constraints, such as high-energy physics experiments at the CERN Large Hadron Collider. However, finding compact expressions for high-dimensional datasets remains challenging due to the inherent limitations of genetic programming (GP), the search algorithm of most symbolic regression (SR) methods. Contrary to GP, the NN approach to SR offers scalability to high-dimensional inputs and leverages gradient methods for faster equation searching. Common ways of constraining expression complexity often involve multistage pruning with fine-tuning, which can result in significant performance loss. In this work, we propose S y m b o l N e t , a NN approach to SR specifically designed as a model compression technique, aimed at enabling low-latency inference for high-dimensional inputs on custom hardware such as FPGAs. This framework allows dynamic pruning of model weights, input features, and mathematical operators in a single training process, where both training loss and expression complexity are optimized simultaneously. We introduce a sparsity regularization term for each pruning type, which can adaptively adjust its strength, leading to convergence at a target sparsity ratio. Unlike most existing SR methods that struggle with datasets containing more than O ( 10 ) inputs, we demonstrate the effectiveness of our model on the LHC jet tagging task (16 inputs), MNIST (784 inputs), and SVHN (3072 inputs).

Tsoi, Ho Fung (ORCID:0000000225502184)↗

Probing jet base emission of M87* with the 2021 Event Horizon Telescope observations

We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the significantly enhanced (u,v) coverage in the 2021 Event Horizon Telescope (EHT) observations. The integration of the 12−m Kitt Peak Telescope (USA) and NOEMA (France) stations into the array introduces two critical intermediate-length baselines to SMT (USA) and IRAM 30−m (Spain), providing sensitivity to emission structures at spatial scales of ∼250 μas and ∼2500 μas (∼ 0.02 pc and ∼ 0.02 pc). Without these new baselines, previous EHT observations of the source in 2017 and 2018 lacked the capability to constrain emission on large scales, where a “missing flux” of order ∼1 Jy is expected to reside. To probe these scales, we analyzed closure phases–robust against station-based gain calibration errors–and model the jet base emission using a simple Gaussian component offset from the compact ring emission at spatial separations > 100 μas. Our analysis revealed a Gaussian feature centered at (ΔRA ≈ 320 μ as, ΔDec. ≈ 60 μ as), projected separation of ≈ 5500 AU, with an estimated flux density of only ∼60 mJy, implying that most of the missing flux identified in previous EHT studies had to originate from different, larger scales. Brighter emission at the relevant spatial scales is firmly ruled out, and the data do not favor more complex models. This component aligns with the inferred position of the large-scale jet and is therefore physically consistent with the emission of the jet base. While our findings point to detectable jet base emission at 230 GHz, the limited coverage provided by only two intermediate baselines limits our ability to robustly reconstruct its morphology. Consequently, we treated the recovered Gaussian as an upper limit on the jet base flux density. Future EHT observations with expanded intermediate baseline coverage will be essential to constrain the structure and nature of this component with higher precision.

accretion, accretion disks↗

Kinetic frustration by limited bond availability controls the LAT protein condensation phase transition on membranes

LAT is a membrane-linked scaffold protein that undergoes a phase transition to form a two-dimensional protein condensate on the membrane during T cell activation. Governed by tyrosine phosphorylation, LAT recruits various proteins that ultimately enable condensation through a percolation network of discrete and selective protein-protein interactions. Here, we describe detailed kinetic measurements of the phase transition, along with coarse-grained model simulations, that reveal that LAT condensation is kinetically frustrated by the availability of bonds to form the network. Unlike typical miscibility transitions in which compact domains may coexist at equilibrium, the LAT condensates are dynamically arrested in extended states, kinetically trapped out of equilibrium. Modeling identifies the structural basis for this kinetic arrest as the formation of spindle arrangements, favored by limited multivalent binding interactions along the flexible, intrinsically disordered LAT protein. These results reveal how local factors controlling the kinetics of LAT condensation enable formation of different, stable condensates, which may ultimately coexist within the cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Magnetic canopies and models of the solar chromosphere

Magnetic canopies are the topologies formed by magnetic field lines as they spread from compact, nearly vertical concentrations of flux in the low photosphere into the large-scale organization of the corona and heliosphere. Analytical techniques for inferring the base-heights of canopies from magnetographic data are reviewed together with observational evidence that much of the sun is covered by canopies which lie two or more pressure scale heights below the level which is traditionally inferred from thin flux tube models. Implications of these results for modeling the structure and energetics of the chromosphere are discussed, and it is argued that future models should be based upon both spectroscopic and polarimetric data. Recent improvements in magnetostatic model atmospheres are reviewed, and new observational data, including C IV Doppler-grams from the Ultraviolet Spectrometer and Polarimeter on the Solar Maximum Mission, are considered. Directions for future research in MHD modeling of canopies, in simulating spectrographic and polarimetric data from such models, and in observational programs are discussed.

Jones, H. P.↗

NREL’s Wind Turbine Drivetrain Condition Monitoring and Wind Plant Operation and Maintenance Research During the 2010s: A US Land-Based Perspective

The wind industry has seen tremendous growth during the past two decades, with the global cumulative installation capacity reaching more than 650 gigawatts by the end of 2019. Despite performance and reliability improvements of utility-scale wind turbines over the years, the industry still experiences premature component failures, leading to increased operation and maintenance (O&M) costs. Among various turbine components, gearboxes—and, more broadly, drivetrains—have shown to be costly to maintain throughout the design life of a wind turbine. The problem of premature component failure is industry wide. As early as 2007, the US Department of Energy (DOE) started to address this challenge through the National Renewable Energy Laboratory’s (NREL’s) reliability initiative that first focused on gearboxes and more recently expanded to entire drivetrains. The wind turbine drivetrain condition monitoring and wind plant O&M research that is the subject of this paper is part of the NREL initiative and includes a few research and development (R&D) activities conducted during the 2010s. These activities included technology evaluation during the first few years; novel monitoring technique investigation (specifically, compact filter analysis) during the middle years; and data and physics domain modeling for fault detection and prediction in recent years. A high-level summary of these activities is provided in this paper along with some key observations from each activity. Most of the work discussed has been published and can be referred to for more information. They reflect the expected evolution of wind turbine condition monitoring and O&M in the US market—primarily, a land-based perspective. In addition, we have identified several R&D opportunities that can be picked up by the research community to help industry advance in related areas, making wind power more cost competitive in the future.

17 WIND ENERGY↗

Hyperspectral Thermal Fabry-Perot Modeling

Fabry-Perot interferometer are simple elegant, tunable filters that can be used to make compact hyperspectral thermal imaging system. To foster the development of these sensors, software tools for the design and simulation of tunable Fabry-Perot infrared imagers have been developed. The tools are provided at three levels: basic, design, and system. Basic tools describe a nearly ideal Fabry-Perot filter with perfectly flat and parallel mirrors in collimated space. Design tools that take into account non-ideal behavior such as mirror and collimation defects calculate free spectral range, finesse, and spectral width of the interferometer. System tools help analyze an integration of the Fabry-Perot filter into a camera system. They include spectral convolution, first-order optical layout, and an estimation of signal-to-noise ratio. The complete set of tools allows for simulations of system operation and performance with various illumination sources. Spectral images generated in such simulations were used to examine applicability of Fabry-Perot system in remote sensing of atmospheric gases including detection of environmental pollutants and hazardous gases. Different operating conditions and system configurations are presented.

Ryan, Robert↗

A numerical study of the meso-structure variability in the compaction process of prepreg platelet molded composites

Compaction deformations in a thermoplastic composite system molded from chopped prepreg platelets were analyzed to obtain statistical distributions of local distortions that define the final consolidated meso-structure. This study provides an insight into and allows quantification of the intrinsic geometric variability of the composite meso-structure that controls the variability of the effective composite mechanical properties. A two-phase viscoelastic finite element model was utilized to treat the behavior of molten prepreg platelets, which were assumed to be incompressible and inextensible in the fiber direction during processing. A two-step simulation procedure with discretely modeled platelets and platelet-to-platelet contacts is proposed. The analysis begins with a drop simulation of the uncontrolled material deposition followed by a 3D finite element analysis of the preform compaction and bulk material flow. The proposed model provides a method for generating a composite meso-structure ab initio starting from a statistical sampling of the initial conditions for platelet deposition. The predicted platelet thickness distribution, intra-platelet fiber waviness and out-of-plane platelet undulations were found to be in good agreement with experimental measurements. An estimate of the volume of resin pockets was inferred from the compaction simulation.

36 MATERIALS SCIENCE↗

Inverse Compton X-ray emission from the superluminal quasar 3C 345

In quasars with strong radio cores, the inverse-Compton process is believed to be the dominant source X-ray emission. For objects with parsec-scale radio jets, simple models have predicted that components in the jet emerging from the quasar nucleus generate the observed X-ray emission. We have tested this hypothesis in detail for the quasar 3C 345 using a ROSAT X-ray observation in 1990 July, together with quasi-simultaneous very long base interferometry (VLBI) imaging of the parsec-scale jet at five frequencies. The ROSAT spectrum is well fitted by a power law with index alpha = -0.96 +/- -0.13, consistent with models in which the X-ray emission results from inverse-Compton scattering of radio radiation from high-energy electrons in compact components. We show that the radio properties of brightest `knot' in the jet (`C5') can be fitted with a homogeneous sphere model whose parameters require bulk relativistic motion of the emitting material; otherwise the predicted model whose parameters require bulk relativistic motion of the emitting material; otherwise the predicted inverse-Compton X-ray emission exceeds the observed flux. If C5 is the origin of the X-ray emission, then it has a Doppler factor delta = 7.5((sup +3 sub -2)). If the nucleus or other components contribute to the X-ray emission, then this becomes a firm lower limit to delta. The inhomogeneous jet model of Koenigl is a good fit both to the barely resolved (less than 1 pc) flat-spectrum nucleus in the radio, and also to the ROSAT X-ray spectrum. The synchrotron and inverse-Compton emitting fluid moves down a narrow cone (opening angle 2 phi approximately 5 deg) nucleus relativistically, with delta approximately 4.6. Doppler factors for the nucleus and C5, derived from our ROSAT observation, provide evidence for bulk relativistic motion in the jet. By combining these constraints with well-known superluminal motion of jet components, we can deduce geometry. For epoch 1990.5 we infer the Lorentz factor gamma = 7.5 ((sup +1.0 sub -1.5)) and angle to the line of sight theta = 8((sup +2 deg sub -3 deg)) for H(sub 0) = 100 km/s/Mpc. These values are the most reliable yet derived using this method, because of the near-simultaneity of our X-ray and VLBI observations and the quality of the multifrequency of VLBI images and component radio spectra.

Unwin, S. C.↗

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗

Correlations of Power-law Spectral and QPO Features In Black Hole Candidate Sources

Recent studies have shown that strong correlations are observed between low frequency QPO s and the spectral power law index for a number of black hole candidate sources (BHCs), when these sources exhibit quasi-steady hard x-ray emission states. The dominant long standing interpretation of QPO's is that they are produced in and are the signature of the thermal accretion disk. Paradoxically, strong QPO's are present even in the cases where the thermal component is negligible. We present a model which identifies the origin of the QPO's and relates them directly to the properties of a compact coronal region which is bounded by the adjustment from Kepleriaa to sub-Kelperian inflow into the BH, and is primarily responsible for the observed power law spectrum. The model also predicts the relationship between high and low frequency QPO's and shows how BH's can be unique identified from observations of the soft states of NS's and BHC's.

Fiorito, Ralph↗

Collaboration for Advanced Modeling of Particle Accelerators

The pverarching purpose is to accelerate and expand the scope of discoveries from high energy physics (HEP) particle accelerators by enabling the design of accelerators that are significantly more compact and cheaper to build and run. This will be realized through (i) developing high-performance computing (HPC) accelerator and beam modeling capabilities to design the full range of systems required (ii) developing community simulation ecosystems that seamlessly integrate accelerator elements to facilitate the design and control of next-generation accelerators.

43 PARTICLE ACCELERATORS↗

A soft departure from jamming: the compaction of deformable granular matter under high pressures

Here, the high-pressure compaction of three-dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follows the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

36 MATERIALS SCIENCE↗