On rapid compaction of granular materials: Combining experiments with in-situ imaging and mesoscale modeling
Not provided.
SEARCH · Engineering Papers
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.
Not provided.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The chief objective of manufacturing process improvement efforts is to significantly minimize process resources such as time, cost, waste, and consumed energy while improving product quality and process productivity. This paper presents a novel physics-informed optimization approach based on artificial intelligence (AI) to generate digital process twins (DPTs). The utility of the DPT approach is demonstrated in the case of finish machining of aerospace components made from gamma titanium aluminide alloy (γ-TiAl). This particular component has been plagued with persistent quality defects, including surface and sub-surface cracks, which adversely affect resource efficiency. Previous process improvement efforts have been restricted to anecdotal post-mortem investigation and empirical modeling, which fail to address the fundamental issue of how and when cracks occur during cutting. In this work, the integration of in-situ process characterization with modular physics-based models is presented, and machine learning algorithms are used to create a DPT capable of reducing environmental and energy impacts while significantly increasing yield and profitability. Based on the preliminary results presented here, we report an improvement in the overall embodied energy efficiency of over 84%, 93% in process queuing time, 2% in scrap cost, and 93% in queuing cost has been realized for γ-TiAl machining using our novel approach.
Metallic Zn is a preferred anode material for rechargeable aqueous batteries towards a smart grid and renewable energy storage. Importantly, understanding how the metal nucleates and grows at the aqueous Zn anode is a critical and challenging step to achieve full reversibility of Zn battery chemistry, especially under fast-charging conditions. Here, by combining in situ optical imaging and theoretical modeling, we uncover the critical parameters governing the electrodeposition stability of the metallic Zn electrode, that is, the competition among crystallographic thermodynamics, kinetics, and Zn 2+ -ion diffusion. Moreover, steady-state Zn metal plating/ stripping with Coulombic efficiency above 99 % is achieved at 10-100 mA cm -2 in a reasonably high concentration (3 M) ZnSO 4 electrolyte. Significantly, a long-term cycling-stable Zn metal electrode is realized with a depth of discharge of 66.7% under 50 mA cm -2 in both Zn || Zn symmetrical cells and MnO 2 || Zn full cells.
Crystal plasticity finite element method (CPFEM) models are widely used to simulate the deformation behaviour of polycrystalline materials, but their calibration is often limited by their high computational cost and the non-convexity of the optimisation landscape. Here, this study develops a multi-objective surrogate-assisted calibration workflow that couples a multi-objective genetic algorithm (MOGA) with an adaptively trained deep neural network (DNN) surrogate model to efficiently identify CPFEM parameters from experimental data. The workflow is demonstrated on three crystal plasticity (CP) formulations of increasing complexity — Voce hardening (VH), two-coefficient latent hardening (LH2), and six-coefficient latent hardening (LH6) — using in situ electron backscatter diffraction (EBSD) measurements of Alloy 617 under uniaxial tensile loading. The CPFEM models are calibrated against the experimentally observed stress–strain response and reorientation trajectories of eight grains, then validated against eight additional trajectories and overall texture evolution. Across the CP formulations, the macroscopic response was reproduced reliably, while differences emerged in the robustness and accuracy of the grain-scale predictions. Including grain reorientation trajectories in the multi-objective calibration improved texture evolution predictions and filtered out physically inconsistent parameter sets that can arise from calibrating against only the stress–strain data. The workflow also demonstrates good transferability of calibrated parameters from a low- to a high-fidelity microstructural model. These results provide practical guidance for integrating in situ microstructural data into CPFEM through efficient, repeatable, and physically meaningful multi-objective calibration.
Since nuclear energy is crucial in the decarbonization of the energy supply, one hurdle to remove is the handling of high-level radioactive waste (HLW). Disposal of HLW in a deep geological repository has long been deemed a viable permanent option. In the design of a deep geological repository, compacted bentonite is the most commonly proposed buffer material. Predicting the long-term chemical evolution in bentonite, which is important for the safety assessment of a repository, has been challenging because of the complex coupled processes. Models for large-scale tests and predictions based on such models have been some of the best practices for such purposes. An 18-year-long in situ test with two dismantling events provided a unique set of chemical data that allowed for studying chemical changes in bentonite. In this paper, we first developed coupled thermal, hydrological, mechanical, and chemical (THMC) models to interpret the geochemical data collected in the in situ test and then extended the THMC model to 200 years to make long-term prediction of the geochemical evolution of bentonite. The interpretive coupled THMC model shows that the geochemical profiles were strongly affected by THM processes such as evaporation/condensation, porosity change caused by swelling, permeability change, and the shape of concentration profiles for major cations were largely controlled by transport processes, but concentration levels were regulated by chemical reactions, and the profiles of some species such as pH, bicarbonate, and sulfate were dominated by these reactions. The long-term THMC model showed that heating prolongs the time that bentonite becomes fully saturated in the area close to the heater/canister; however, once the bentonite becomes fully saturated, high concentrations of ions in bentonite near the heater, which was observed in the field test, will disappear; illitization continues for 50 years but will not proceed further.
The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).
Abstract. Growth in adoption of distributed wind turbines for energy generation is significantly impacted by challenges associated with siting and accurate estimation of the wind resource. Small turbines, at hub heights of 40 m or less, are greatly impacted by terrestrial obstacles such as built structures and vegetation that can cause complex wake effects. While some progress in high-fidelity complex fluid dynamics (CFD) models has increased the potential accuracy for modelling the impacts of obstacles on turbulent wind flow, these models are too computationally expensive for practical siting and resource assessment applications. To understand the efficacy of available models in situ, this study evaluates classic and commonly used methods alongside new state-of-the-art lower-order models derived from CFD simulations and machine learning approaches. This evaluation is conducted using a subset of an extensive original dataset of measurements from more than 300 operational wind turbines in the northern Netherlands. The results show that data-driven methods (e.g. machine learning and statistical modelling) are most effective at predicting production at real sites with an average error in annual energy production of 2.5 %. When sufficient data may not be available de novo to support these data-driven approaches, models derived from high-fidelity simulations show promise and reliably outperform classic methods. On average these models have 6.3 %–11.5 % error compared with 26 % for classic methods and 27 % baseline error for reanalysis data without obstacle correction. While more performant on average, these methods are also sensitive to the quality of obstacle descriptions and reanalysis inputs.
The Raft River geothermal field is the site of an innovative Department of Energy Enhanced Geothermal System (EGS) project to determine the viability of using combined thermal and hydraulic stimulation techniques to improve energy production. Well RRG-9 is currently undergoing a stimulation program using injectate from the US Geothermal Raft River Power Plant and cold water from a cooling tower make-up water well. The stimulation began on 13 June 2013 with injection from the power plant at a temperature of about 39 °C and a pressure of 275 psig. Next, two positive displacement plunger type pumps were used to increase the injection pressure and flow rate for about one month. The highest rate achieved was 258 gpm at a pressure of 741 psig. During this time, fluid from the cooler water well was injected for about 2 weeks at various pressures. Then, the pumps were removed and plant injection resumed on 25 September. Plant injection will continue until the spring of 2014, when a high pressure hydraulic stimulation will be conducted. A series of seismic monitoring stations deployed around the well are providing data on seismic events occurring at the site. Over the past year, 51 microseismic events have been recorded, all less than Magnitude 1. During injection, several diagnostic tests were conducted to gain a better understanding of the well and reservoir. A step-rate test was performed on 22 August to measure the in-situ stress and aid in modeling in-situ fractures. A tracer was injected into the well on 9 September. No tracer was detected in adjacent production wells after several months. A second borehole televiewer survey was conducted for comparison to pre-stimulation images. Here, a third borehole televiewer survey is planned after the high pressure stimulation. Injection test data is evaluated in real time. A modified Hall plot analysis indicates the effective permeability is increasing. The injectivity index supports the results of the modified Hall plot analysis. As the thermal stimulation has continued, the injectivity index has consistently followed an upward trend from 0.1 gpm/psi to 0.53 gpm/psi.
The ideal elastic limit is the upper bound of the achievable strength and elastic strain of solids. However, the elastic strains that bulk materials can sustain are usually below 2%, due to the localization of inelastic deformations at the lattice scale. In this study, we achieved >5% elastic strain in bulk quantity of metallic glass, by exploiting the more uniform and smaller-magnitude atomic-scale lattice strains of martensitic transformation as a loading medium in a bulk metallic nanocomposite. The self-limiting nature of martensitic transformation helps to prevent lattice strain transfer that leads to the localization of deformation and damage. This lattice strain egalitarian strategy enables bulk metallic materials in kilogram-quantity to achieve near-ideal elastic limit. This concept is verified in a model in situ bulk amorphous (TiNiFe)-nanocrystalline (TiNi(Fe)) composite, in which the TiNiFe amorphous matrix exhibits a maximum tensile elastic strain of similar to 5.9%, which approaches its theoretical elastic limit. As a result, the model bulk composite possesses a large recoverable strain of similar to 7%, a maximum tensile strength of above 2 GPa, and a large elastic resilience of similar to 79.4 MJ/m 3 . The recoverable strain and elastic resilience are unmatched by known high strength bulk metallic materials. This design concept opens new opportunities for the development of high-performance bulk materials and elastic strain engineering of the physiochemical properties of glasses.
The response of high explosives to shock loading is traditionally measured with a steady loading pressure. In many accident scenarios involving fragment impact, however, a loading duration that is shorter than the build up to detonation may occur. Fragments passing through multiple materials before reaching a high explosive charge may produce loading that is comprised of more than one shock wave. Additionally, the build up to detonation in high explosive corner turning loads the explosive a short duration pressure pulse, since rarefactions can often rapidly overtake the reactive wave. For these reasons, we have studied the response of the insensitive high explosive (IHE) materials PBX 9502 and LX-17 to complex loadings of varied intensity and duration. We refer to a single loading of limited duration as a “thin pulse”, whereas more complex scenarios were studied with an impactor that produces a double shock in the explosive. The following report presents experimental data and analyses of thin pulse shock initiation and double shock experiments designed to guide development of models of Insensitive High Explosives (IHEs) under controlled one-dimensional conditions relevant to accident scenarios and corner turning. Thin pulse shock initiation data on PBX 9502 and LX-17 were obtained under varied pulse duration, pressed density, and temperature conditions in order to probe various parameters essential for the development of a physics-based Cheetah reactive flow hotspot model. In situ pressure gauges provide insight into the degree of reaction in the explosive that are not obtainable with optical PDV measurements or distance measurements such as run to detonation. Double shock data was obtained to inform a Composition Aware Cheetah model which can be applied to any TATB-based IHEs. This model supports efforts to find a new IHE formulation and potentially incorporate new binders into IHE formulations. Simulations of each experiment are included to demonstrate the utility of these focused experiments to developing models of HE behavior. One-dimensional gas gun experiments are essential for characterizing shocked HE behavior and informing HE models.
The decommissioning of hardened nuclear facilities provides a unique challenge in balancing current worker dose with exposure to future receptors and media. In situ decommissioning (ISD) is a cost-effective and safe option for the closure of such facilities but requires a detailed understanding of the threat that radioactive and hazardous constituents remaining in the facilities pose to groundwater and surrounding surface water. The Savannah River Site (SRS) extensively employs contaminant migration models within the decommissioning process to develop site-specific removal, grouting, and monitoring strategies in support of safe and effective final end states of nuclear facilities. This work discusses the use of contaminant migration modeling in the ISD process at SRS, including specific examples from the closure of P-Reactor and R-Reactor, Building 235-F, F-Canyon and F/H Laboratory Complex. Migration models developed during the closure of P-Reactor and R-Reactor were divided into four source areas: the reactor vessel, the process area, the disassembly basin, and the purification wing. Each source area was assigned a specific inventory and migration pathway and was then subject to varying hypothetical removal, capping, and grouting schemes. Modeling indicated that groundwater Maximum Concentration Limits (MCLs) might be exceeded if no action was taken for eleven and ten constituents of concern (COCs) at P-Reactor and R-Reactor, respectively. These exceedances could occur in as few as 200 to 500 years. Alternatively, the migration modeling demonstrated that the selected ISD actions reduced contaminant mobility which allowed for significant radioactive decay and resulted in fewer predicted exceedances of groundwater MCLs (five COCs for P-Reactor; eight COCs for R-Reactor). In response to the modeling results for P-Reactor and R-Reactor, effectiveness monitoring programs were developed to target contaminants, specific to each reactor, that may migrate to groundwater. Contaminant migration modeling also revealed that roof collapse was a large factor in the release of COCs to the environment, giving rise to roof improvements, and an inspection and vegetation control program to ensure roof stability over time. Contaminant migration modeling is also aiding in the closure planning for hardened facilities in F Area. Building 235-F housed the Actinide Billet Line, which produced Np-237 billets for irradiation in SRS reactors, and the Plutonium Fuel Form (PuFF) facility that produced Pu-238 heat sources for the space program. As a result of these missions, areas within Building 235-F contain considerable residual amounts of both Pu-238 and Np-237. Contaminant migration modeling of Building 235-F was originally performed in 2012 to identify the feasibility of ISD and the amount of radioactive material removal required to prevent the exceedance of groundwater MCLs. The original model indicated that a 60% reduction in the PuFF facility Pu-238 inventory could keep groundwater concentrations below standards, while Np-237 did not pose a threat to groundwater. However, updates to the model with an emphasis on source impact pathways revealed that, due to the orientation of the source areas relative to groundwater flow, no amount of reasonable removal of Pu-238 would keep groundwater concentrations below MCLs and that Np-237 could be a large contributor to localized MCL exceedances. With this insight, the refined 2019 model is being used to assist in the development of grouting plans specific to each facility source area, where bentonite may be utilized to slow the migration of Pu-238 and its daughter products from the PuFF facility and a reducing grout may decrease Np-237 transport by ensuring the nuclide remains in the less mobile +IV oxidation state. The beginning phases of contaminant migration modeling are underway for F-Canyon and associated facilities using lessons learned from the Reactors and Building 235-F. A contaminant migration pathway, similar to the pathway used for the reactor vessels, is being developed for the hot and warm canyons. The F-Canyon inventories are being spatially refined to identify specific and localized source zones, comparable to Building 235-F, that may impact groundwater. Contaminant migration modeling provides vital input within the ISD process, from facility investigation to post closure effectiveness monitoring, making it a valuable tool in the development of safe and effective end states for hardened nuclear facilities at SRS. (authors)
We present a method for creating spacecraft-like data which can be used to train Machine Learning (ML) models to detect and classify structures in in situ spacecraft data. First, we use the Grad-Shafranov equation to numerically solve for several magnetohydrostatic equilibria which are variations on a known analytic equilibrium. These equilibria are then used as the initial conditions for Particle-In-Cell simulations in which the structures of interest are observed and labeled. We then take one-dimensional slices through the simulations to replicate what a spacecraft collecting data from the simulation would observe. This sliced data then can be used as training data for the initial training of ML models intended for use on spacecraft data. We demonstrate the method applied to the problem of detecting small-scale plasmoids in the magnetotail, which is important for understanding complex magnetotail reconnection dynamics. The simple 1D classifier we train is able to detect more than 70% of the plasmoid points in the data set but also produces a large number of false positives. Our further work on this example problem is detailed, and further potential uses of the method are discussed.
A multiscale simulation approach was developed and employed to optimize the sheet surface conditions for higher interfacial temperature and joint strength in ultrasonic welding of magnesium alloy AZ31 and dual-phase steel DP590. First, a mesoscale model was used to study the relationship between friction coefficient and surface roughness, which can be modified by various engineering methods. Then a macroscopic process model was employed to study the effects of surface roughness on heat generation, indicating that a temperature increase can be achieved with rougher surfaces on two sides of both DP590 and AZ31 sheets. Samples prepared by sanding and filing, as well as grinding, were first characterized for surface roughness and then welded under ultrasonic vibration. An infrared camera was used to measure temperatures in situ for model validation. Overall, lap shear test results for the welded joint showed that the joint strength can be improved by 10~25% using filing and round grinding methods as a result of the enhanced heat generation and mechanical interlocking on the interface.
Direct non-oxidative methane (CH 4 ) conversion to value-added hydrogen (H 2 ) and C 2 products remains hindered by fundamental catalytic scaling constraints and rapid surface deactivation at elevated temperatures. Plasma-enabled catalysis offers a promising route to overcome the thermodynamic and kinetic barriers of direct non-oxidative methane upgrading at mild conditions, yet control over C–C product selectivity and catalyst stability remains elusive. Here, we establish a unified mechanistic framework including Langmuir–Hinshelwood (L–H) and Langmuir–Rideal (L–R) mechanisms that disentangles the roles of plasma excitation (including vibrationally activated species and radicals), surface temperature (T sur ), and catalyst binding energy in steering CH 4 conversion to H 2 and C 2 hydrocarbons. Through a combination of density functional theory (DFT) informed microkinetic modeling, in situ and ex situ surface characterization, and product quantification under dielectric barrier discharge conditions, we show that vibrationally excited CH 4 lowers activation barriers selectively for dissociative chemisorption, enabling surface activation across a wide range of transition metal catalysts at low thermal energy input. We find that once CH 4 is dissociatively chemisorbed, the branching between C 2 H 2 , C 2 H 4 , and C 2 H 6 is governed by surface properties (carbon binding energy, T sur , etc), regardless of plasma excitation. The DFT informed microkinetic model decouples the effects of molecular activation from surface properties and indentifies operating windows that maximize target yields while suppressing carbon accumulation and subsequent catalytic inactivation. Experiments on polycrystalline Cu/Al 2 O 3 , Ni/Al 2 O 3 , and Pt/Al 2 O 3 validate these predictions, revealing catalyst-dependent branching toward ethane or ethylene and distinct deactivation profiles. We unify these trends into a generalized three-dimensional plasma-thermal-catalytic design space, from which reduced descriptors such as T vib /T sur in the limit of vibrationally excited L–H pathways emerge as predictive metrics. These results enable rational tuning of methane conversion pathways and unlock selective C 2 formation using earth-abundant metals under mild plasma conditions.