Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “load 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.

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

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sample glue layer investigation and mitigation for laser induced prompt impulse experiments

Understanding longer timescale material reactions under dynamic stress loading is critical for applications in materials engineering, shock physics, and planetary science. Prompt impulse experiments generate lower pressures since the ablator—the material first removed by the laser—is thicker and farther from the diagnostic plane, capturing aggregate material responses from the initial shock wave, rarefaction waves, and later time effects. This complexity demands thorough material characterization and simulation support. Since traditional sample construction is specific to supported shock experiments, designing prompt impulse experiments requires reconsideration around target design and sample engineering. Here, we present sample preparation techniques, experimental investigations, and theoretical simulations to investigate glue layer impacts, aiming to standardize samples for consistent data at lower laser fluences. We find that glue layers <30 μm have a minimal impact on peak velocity and pulse shape. The peak velocity scales linearly with glue layer thickness until a glue layer of 75 μm. For glue layers >75 μm, the peak velocity no longer scales with thickness; however, the pulse shape continues to degrade as described by simulations.

Lasers↗
DOE OSTI2024

Rapid detection of rare events from in situ X-ray diffraction data using machine learning

High-energy X-ray diffraction methods can non-destructively map the 3D microstructure and associated attributes of metallic polycrystalline engineering materials in their bulk form. These methods are often combined with external stimuli such as thermo-mechanical loading to take snapshots of the evolving microstructure and attributes over time. However, the extreme data volumes and the high costs of traditional data acquisition and reduction approaches pose a barrier to quickly extracting actionable insights and improving the temporal resolution of these snapshots. This article presents a fully automated technique capable of rapidly detecting the onset of plasticity in high-energy X-ray microscopy data. The technique is computationally faster by at least 50 times than the traditional approaches and works for data sets that are up to nine times sparser than a full data set. This new technique leverages self-supervised image representation learning and clustering to transform massive data sets into compact, semantic-rich representations of visually salient characteristics ( e.g. peak shapes). These characteristics can rapidly indicate anomalous events, such as changes in diffraction peak shapes. It is anticipated that this technique will provide just-in-time actionable information to drive smarter experiments that effectively deploy multi-modal X-ray diffraction methods spanning many decades of length scales.

Zheng, Weijian↗
DOE OSTI2022

High-temperature superelasticity in NiTiPt and NiTiPd shape memory alloys

Alloying NiTi shape memory alloys with sufficient Pt and Pd at the expense of Ni has been shown to increase transformation temperatures well above those possible with the binary composition, enabling consideration of a new class of high-temperature mechanical sensors and actuators. To achieve shape memory behavior at temperatures greater than 200 °C, previous studies of NiTiPt and NiTiPd compositions have focused primarily on alloys with greater than 16 at.% Pt and 25 at.% Pd. While generally not the focus, none of these previous studies demonstrated superelastic behavior at high temperatures. In an effort to identify alloys that transform to austenite at approximately 200 °C and exhibit superelastic behavior with recoverable strains greater than 2%, a series of Ti-rich NiTiPt ternary alloys with 15.5–16.5 at.% Pt, and Ti-rich NiTiPd ternary alloys with 25–26 at.% Pd have been targeted. The polycrystalline alloys were processed via arc melting followed by homogenization heat treatment. Thermomechanical properties were measured using electro-discharge machined tension test coupons and a miniature load frame equipped with a clamshell furnace. Ultimately, superelasticity in uniaxial tension with clear plateau regions was demonstrated and characterized in these alloys at temperatures between 219 °C and 270 °C.

Massad, Jordan E.↗
DOE OSTI2024

Multi–Stage Optimization of Pore Size and Shape in Pore–Space–Partitioned Metal–Organic Frameworks for Highly Selective and Sensitive Benzene Capture

Compared to exploratory development of new structure types, pushing the limits of isoreticular synthesis on a high-performance MOF platform may have higher probability of achieving targeted properties. Multi-modular MOF platforms could offer even more opportunities by expanding the scope of isoreticular chemistry. However, navigating isoreticular chemistry towards best properties on a multi-modular platform is challenging due to multiple interconnected pathways. Here on the multi-modular pacs (partitioned acs) platform, we demonstrate accessibility to a new regime of pore geometry using two independently adjustable modules (framework-forming module 1 and pore-partitioning module 2). A series of new pacs materials have been made. Benzene/cyclohexane selectivity is tuned, progressively, from 4.5 to 15.6 to 195.4 and to 482.5 by pushing the boundary of the pacs platform towards the smallest modules known so far. The exceptional stability of these materials in retaining both porosity and single crystallinity enables single-crystal diffraction studies of different crystal forms (as-synthesized, activated, guest-loaded) that help reveal the mechanistic aspects of adsorption in pacs materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Modeling the Yaw Behavior of Tail Fins for Small Wind Turbines: November 22, 2021 - May 21, 2024

This report describes the aerodynamic modelling of the yaw behavior of tail fins for small wind turbines (SWTs). The analysis is based on unsteady slender body theory (USBT) and is formulated to be used with the proposed tail fin module to be added to OpenFAST. Linearized USBT, which has been applied in the past to SWT tail fins, is developed in three main ways. First, USBT is extended to high yaw angles and aspect ratios by modelling the associated nonlinear vortex dynamics and the chordwise load distribution respectively. Second, we consider the effect of time-varying wind speed. The extended theory is compared to recent unpublished measurements at the University of Perugia, Italy, for the yaw behavior of delta, elliptical, and rectangular tail fins without a rotor and nacelle. The fins were released from initial yaw angles of -40 degrees and -80 degrees ; the latter is of sufficient magnitude to show the importance of the nonlinear yaw dynamics. Further, the friction in the model tail fin bearings was measured and modelled. The constants in the aerodynamic models were taken from the literature but significant improvement in accuracy is achieved by using system identification techniques to optimize the values. This leads to the third main finding: extended USBT is accurate for all generic tail fin shapes that we considered and is, therefore, very suitable for inclusion in aeroelastic codes for SWTs. The final development is of a general nonlinear equation for yaw response that is simpler than the extended USBT and is also suitable for inclusion in OpenFAST. The two nonlinear models are used to highlight the geometric requirements for good tail fin design.

distributed wind↗
DOE OSTI2024

A crystal-plasticity-informed Gaussian Process Regression model to capture anisotropy in single crystal shape memory alloys

This work presents a machine learning (ML) framework that model the anisotropic actuation responses in a shape memory alloy. A Gaussian Process Regression (GPR) based ML model is trained on a set of different crystal orientations subjected to different actuation conditions. The training employed thermo-mechanical responses from a crystal-plasticity model that captures phase-transformation, stress-induced plasticity, and transformation-induced plasticity. Further, on training the GPR-ML model at fixed stress level for different orientations, it captured the thermo-mechanical responses accounting for the anisotropy, and predicted responses for new orientations with good accuracy. The GPR-ML model is able to capture the transformation temperature variations even when trained using multiple stress levels, and the transformation strain showed significant deviations. The developed GPR-ML model gave reasonable predictions for an unexplored sample set of orientations and loading conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
DOE OSTI2025

Structured light approaches in laser-based plasma diagnostics

There is a growing demand for plasma diagnostics suitable for industrial plasma reactors employed in semiconductor nanofabrication, especially relevant to microelectronics and quantum information systems. Such reactors typically have limited optical access and pose considerable diagnostic challenges, including intense background emission, significant thermal loads, and contamination of optical viewports. In this study, we outline research into structured light techniques (laser beams with tailored spatial, temporal, or phase characteristics) that effectively overcome these issues using laser-induced fluorescence (LIF) as an example. The focus of presented diagnostics is on ion kinetics analysis within an industrial plasma source, although this approach is broadly applicable to other plasma systems and diagnostic contexts. We present a confocal LIF implementation using an axicon-generated Bessel annular beam, achieving spatial resolutions of approximately 5 mm at a focal distance of 300 mm, with potential improvements to about 1 mm. This approach matches conventional orthogonal LIF performance but requires only one optical port. Wavelength-modulation LIF employs nonlinear laser wavelength tuning to measure spectral line derivatives, suppressing background emission and enhancing details of spectral line shape. Additionally, we present new results on applying vortex beams (laser beams carrying orbital angular momentum, OAM) for LIF measurements in an industrial plasma device. These measurements enable simultaneous axial and tangential velocity determination using a single laser beam and have been tested with xenon ion transition. Initial quantification of results was performed. Together, these structured-light approaches provide robust, background-resilient, multi-dimensional diagnostics for complex plasma environments.

Romadanov, Ivan [Princeton Plasma Physics Laborato↗
DOE OSTI2026

Habitat Affinity of Riverine Dissolved Organic Matter Linked to Molecular Traits

Fluxes of organic matter across Earth habitats shape ecosystem function and carbon storage. Yet, it is challenging to predict howthe molecules are preferentially accumulated in specific habitats, that is, the habitat affinities, due to their diverse character-istics. Here, we develop an indicator of compositional-level habitat affinity for dissolved organic matter (DOM) by quantifyingthe aggregated affinities of individual molecules that are over- or under-represented relative to a reference habitat. Applyingthis indicator to 93 paired riverine water-sediment sites, we find that 65.1% of molecules show non-significant habitat affinitiesfor water or sediment and are thermodynamically favorable and susceptible to microbial degradation. Consistently, the habitataffinities of DOM assemblages decrease with lower molecular weights and recalcitrance for waters, and with higher thermody-namic favorability for sediments. These affinities are both reduced by elevated temperatures and anthropogenic pressures, suchas nitrogen loading and land-use change. Our projection of habitat affinities of global riverine DOM shows that regions withintense anthropogenic disturbances or warm climates, such as East Asia, Western Europe, and the Amazon, tend to exhibit loweraffinities. These findings suggest enhanced molecular homogenization between water and sediment habitats is associated withenvironmental changes and underscore the importance of habitat affinities for predicting the potential trajectories of organiccarbon across ecosystems

Cui, Yifan↗
DOE OSTI2022

Facet Impact of CuMn 2 O 4 Spinel Nanocatalysts on Enhancement of the Oxygen Reduction Reaction in Alkaline Media

Surface-structure engineering represents an attractive strategy to optimize the energy conversion performance of nanocatalysts using their deliberately controlled exposed facets. To further exploit the potential of non-PGM-based spinel catalysts for alkaline oxygen reduction reaction (ORR), a cathodic fuel cell reaction, we hereby report a strategy of ORR improvement by controlling the crystallographic facets of ultra-small CuMn 2 O 4 spinel nanocatalysts through a developed colloidal synthesis approach. The synthesis of CuMn 2 O 4 nanocrystals with morphological control relies on the design and selection of the Cu-/Mn-precursors with striking discrepancies in reaction kinetics. Following carbon loading and an annealing post-treatment of the as-synthesized nanocatalysts, the exclusively {101} facet-exposed CuMn 2 O 4 spinel nano-octahedra exhibit improved electrocatalytic activity toward ORR in 1 M KOH, when compared to their spherical counterparts, exhibiting a mass activity (MA) of 37.6 A/g at 0.85 V. After 10,000 cycles of ORR durability test, the nano-octahedra still retain an MA of 24.5 A/g, which is twice that of the CuMn 2 O 4 spinel nanospheres. Structural characterizations after durability testing indicate that the MA decay is likely associated with a decrease in the Mn 3+ fraction and the emergence of Cu + on CuMn 2 O 4 nano-octahedral surfaces. Finally, as a paradigm, this synthesis approach could be extended to other Mn-based spinel nanocatalysts with precise shape control, enabling us to understand and establish the relationship between the surface lattice/valence state and electrocatalytic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of a high current density, high temperature superconducting cable for pulsed magnets

Abstract A low AC loss Rare Earth Barium Copper Oxide (REBCO) cable, based on the VIPER cable technology has been developed by Commonwealth Fusion Systems for use in high field, REBCO based tokamaks. The new cable is composed of partitioned and transposed copper ‘petals’ shaped to fit together in a circular pattern with each petal containing a REBCO tape stack and insulated from each other to reduce AC losses. A stainless steel jacket adds mechanical robustness—also serving as a vessel for solder impregnation—while a tube runs through the middle for cooling purposes. Additionally, fiber optic sensors are placed under the tape stacks for quench detection. To qualify this design, a series of experiments were conducted as part of the SPARC tokamak Central Solenoid Model Coil program—to retire the risks associated with full scale, fast ramping, high flux HTS Central Solenoid (CS) and Poloidal Field (PF) coils for tokamak fusion power plants and net energy demonstrators. These risk study and risk reduction experiments include (1) AC loss measurement and model validation in the range of ~5 T/s, (2) an IxB electromagnetic loading of over 850 kN/m at the cable level and up to 300 kN/m at the stack level, (3) a transverse compression resilience of over 350 MPa, (4) manufacturability at tokamak relevant speeds and scales, (5) cable to cable joint performance, (6) fiber optic based quench detection speed, accuracy, and feasibility, and (7) overall winding pack integration and magnet assembly. The result is a cable technology, now referred to as PIT VIPER, with AC losses that measure fifteen times lower (at ~5 T/s) than its predecessor technology; a 2% or lower degradation of critical current (Ic) at high IxB electromagnetic loads; no detectable Ic degradation up to 570 MPa of transverse compression on the cable unit cell; end to end magnet manufacturing, consistently producing Ic values within 7% of the model prediction; cable to cable joint resistances at 20 K on the order of ~15 nΩ; and fast, functional quench detection capabilities that do not involve voltage taps. This cable technology will be tested comprehensively in a Central Solenoid Model Coil to prove its readiness for compact, high field tokamak operation.

Sanabria, Charlie (ORCID:0000000150175309)↗
DOE OSTI2023

Large-eddy simulations of turbulent flows in arrays of helical- and straight-bladed vertical-axis wind turbines

Effects of helical-shaped blades on the flow characteristics and power production of finite-length wind farms composed of vertical-axis wind turbines (VAWTs) are studied numerically using large-eddy simulation (LES). Two helical-bladed VAWTs (with opposite blade twist angles) are studied against one straight-bladed VAWT in different array configurations with coarse, intermediate, and tight spacings. Statistical analysis of the LES data shows that the helical-bladed VAWTs can improve the mean power production in the fully developed region of the array by about 4.94%–7.33% compared with the corresponding straight-bladed VAWT cases. The helical-bladed VAWTs also cover the azimuth angle more smoothly during the rotation, resulting in about 47.6%–60.1% reduction in the temporal fluctuation of the VAWT power output. Using the helical-bladed VAWTs also reduces the fatigue load on the structure by significantly reducing the spanwise bending moment (relative to the bottom base), which may improve the longevity of the VAWT system to reduce the long-term maintenance cost.

17 WIND ENERGY↗
DOE OSTI2024

Additive manufacturing of LiCoO 2 electrodes via vat photopolymerization for lithium ion batteries

Additive manufacturing has the potential to revolutionize the fabrication of lithium-ion batteries for a diversity of applications including in portable, biomedical, aerospace, and the transportation fields. Standard commercial batteries consist of stacked layers of various components (current collectors, cathode, anode, separator and electrolyte) in a two-dimensional manner. By leveraging the latest advances in additive manufacturing and computer-aided design, an improved geometric and electrochemical configuration of these batteries can maximize energy efficiency while allowing design optimization to reduce dead space for a given application. In this work, a composite UV photosensitive resin was prepared and used as feedstock in a vat photopolymerization system. The resin was loaded with LiCoO 2 acting as electrochemically active material for the cathode of a lithium-ion battery, and was further improved with the addition of conductivity-enhancing carbonaceous additives. Challenges to additive manufacturing arise from the opacity and high viscosity of the composite nature of these electrochemically-active resins, which cause light refraction during selective UV curing. Subsequently, items were printed and subjected to a thermal post-processing step to obtain an adequate compromise between electrochemical performance and mechanical integrity. Both sintered and green state 3D printed cathodes were assembled into half-cell lithium-ion batteries using lithium metal as a reference and counter electrode. Electrochemical cycling of these batteries yielded satisfactory results approaching commercial LiCoO 2 cathodes’ performance, with the potential advantages of additive manufacturing – high surface area anode–cathode configurations for power performance as well as shape conformability.

25 ENERGY STORAGE↗
DOE OSTI2026

Multipoint Aerostructural Optimization of Wind Turbine Rotors Using a Coupled Blade‐Resolved Aerostructural Solver

Physics‐based design optimization workflows thread the needle between computational cost limitations and simulation complexity, often compromising between modeling detail and the range of operating design conditions. Multipoint aerostructural optimization of wind turbine rotors has so far been confined to low‐fidelity analyses or to high‐fidelity studies with simplified structural models, leaving the most complex design trade‐offs unexplored. We close this gap by performing the first tightly coupled gradient‐based multipoint aerostructural rotor optimization using 3D aerodynamic and structural solvers with discrete coupled adjoints. The optimizer simultaneously varies blade planform, airfoil shapes, and structural thickness through more than 270 design variables, minimizing a weighted combination of rotor mass and power across multiple wind speeds. Applied to a modified DTU 10‐MW benchmark under conservative structural and aerodynamic constraints, our multipoint optimization reduces rotor mass by up to 36% and increases power by 12%–15% across the main operating conditions; biasing the objective toward power yields power gains up to 18% and a 17% mass reduction. For a nominal wind distribution, 3‐point rotor designs accounting for low RPM and high thrust conditions capture dominant trade‐offs and outperform single‐point designs. Adding two off‐design points changes individual‐condition power by less than 3% but leaves the weighted average within 0.5%, and the mass‐power bias has a stronger effect on the final design than the operating‐point weighting itself. Our framework extends naturally to richer load cases and site‐specific wind distributions, providing a basis for high‐fidelity multipoint design earlier in industrial workflows.

17 WIND ENERGY↗
DOE OSTI2023

Mechanics of Dynamic and Deformable DNA Nanostructures

In DNA nanotechnology, DNA molecules are designed, engineered, and assembled into arbitrary-shaped architectures with predesigned functions. Static DNA assemblies often have delicate designs with structural rigidity to overcome thermal fluctuations. Dynamic structures reconfigure in response to external cues, which have been explored to create functional nanodevices for environmental sensing and other applications. However, the precise control of reconfiguration dynamics has been a challenge due partly to flexible single-stranded DNA connections between moving parts. Deformable structures are special dynamic constructs with deformation on double-stranded parts and single-stranded hinges during transformation. These structures often have better control in programmed deformation. However, related deformability and mechanics including transformation mechanisms are not well understood or documented. In this review, we summarize the development of dynamic and deformable DNA nanostructures from a mechanical perspective. We present deformation mechanisms such as single-stranded DNA hinges with lock-and-release pairs, jack edges, helicity modulation, and external loading. Theoretical and computational models are discussed for understanding their associated deformations and mechanics. We elucidate the pros and cons of each model and recommend design processes based on the models. The design guidelines should be useful for those who have limited knowledge in mechanics as well as expert DNA designers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Mechanical Comparison of Short Models of Nb$_3$ Sn Low-$\beta$ Quadrupole for the Hi-Lumi LHC

MQXF is the Nb 3 Sn Low-β quadrupole magnet that the HL-LHC project is planning to install in the LHC interaction regions in 2026 to increase the LHC integrated luminosity. The magnet will be fabricated in two different lengths: 4.2 m for MQXFA, built in the US by the Accelerator Upgrade Project (AUP), and 7.15 m for MQXFB, fabricated by CERN. In order to qualify the magnet design and characterize its performance with different conductors, cable geometries and pre-load configurations, five short model magnets, called MQXFS, were fabricated, assembled and tested. We compare the mechanical behavior of short model magnets using experimental data and new numerical models that take into account the measured coil sizes as a function of position.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Ion Channel Laser Based on Direct Laser Acceleration of a Shaped Beam Driver (Final Technical Report)

The main focus of our research was to understand the synergies between an electron bunch and a laser pulse when the two co-propagate through the plasma. We have clearly demonstrated using theoretical and computational modeling that the propagation distance of both the bunch and the laser pulse could be extended. The ability of the combined bunch/pulse system to extend the propagation distance and the size of the plasma bubble enables highly-efficient sources of relativistic electrons. As pointed out by the Plasma Decadal Study, such electron sources can be used for generating extremely bright x-rays for a variety of applications can serve as probes and diagnostics for other plasma experiments: high energy density (HED) sciences, inertial confinement fusion (ICF), and potentially Fusion Materials and Technology (FM&T). Development of advanced diagnostics of plasma-based accelerators is yet another key area identified by recent reports. Even broader security and medical applications of compact accelerator-based radiation sources, such as very high energy electron (VHEE) sources for FLASH radiobiology, have been identified by a recent multi-agency panel. A number of key technical issues were considered and successfully resolved during the course of the grant. Those include: beam loading effect of both the driver and witness bunches on the wake, laser channeling by the bunch, and the direct laser acceleration (DLA) of the driver bunch by the laser pulse. To demonstrate clear synergy, we were able to ascertain that the total energy gain of a witness bunch in the wake of the combined bunch/laser complex is large than the sum of the energy gains in the wake of the laser pulse alone, and the driver bunch alone. We have also demonstrated that not only the energy gain is improved, but the energy spread is not sacrificed. To carry out these simulations, we have to develop a range of in-house computational tools, including a fully-3D code WAND-PIC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Competition between roughness and strength for scale-dependent surfaces

Rocks famously have scale-dependent strength, yet the actual dependence is notoriously hard to measure or incorporate into any theoretical framework. Natural rough surfaces present an opportunity to solve the problem. Surfaces sliding in shear evolve as protrusions collide. These asperities can deform or break, thus creating a new surface shape. In particular, natural surfaces have roughness at all scales as well as scale-dependent strength. Based on a scaling analysis, we have previously suggested that the scale-dependent aspect ratio of steady-state surfaces should be proportional to the scale-dependent shear strain at yield. If true, scale-dependent strength could easily be inferred from natural surfaces. Thus, moving beyond the scaling argument to a rigorous treatment of scale-dependent strength for multiscale rough surfaces in shear is important. However, analytic frameworks for analyzing multiscale problems are challenging, as conventional continuum mechanics typically involves a single value for a material property across scales. Here, in this work, we build on the formalism of Persson (2001) that presents a method to compute contact area for rough surfaces with a prescribed topographic spectrum using a stochastic differential equation. The Persson formalism allows for plastic yield under normal loading of otherwise elastic materials and leaves open the possibility of scale-dependent yield stress. In this study, we pursue this route to develop a theory and numerical results for the yielding of a rough, elastoplastic surface with scale-dependent yield stress. Here, we examine surfaces for which the power spectrum of the topography 𝐶 and yield stress 𝑌 follow power laws as a function of scale 𝜆, such that 𝐶∼𝜆 −𝑚 and 𝑌∼𝜆 −𝑛 , respectively. In this formal treatment of the problem, we focus on surfaces in contact and the resulting yield and do not impose shear. Numerical solutions show that the deviation from the elastic scaling solution is bounded as expected by the prior 1D heuristic scaling argument that anticipates the Hurst exponent as 1−𝑛. We also show that the plasticity is expected to erode the contacts if 𝑚 is lower than 𝑛−3, which corresponds to a Hurst exponent lower than 1−𝑛/2. This result is rigorously sound for 2D, i.e., realistic surfaces, and quantitatively different than the prior scaling argument. The theory now permits a correspondingly quantitative approach to interpreting natural surfaces.

elasticity↗