Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Behavior modeling”

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 1,045 records · Page 58

Oscillatory redox behavior in oxides: Cyclic surface reconstruction and reactivity modulation via the Mars–van Krevelen mechanism

The breaking of translational symmetry at oxide surfaces gives rise to coordinatively unsaturated cations/anions and surface restructuring—key factors that govern surface reactivity. Using direct in situ environmental transmission electron microscopy (TEM) observations along with atomistic modeling, we report oscillatory redox behavior in CuO under H 2 , where cyclic surface reconstruction and reactivity modulation occur via the Mars–van Krevelen (MvK) mechanism. We observe self-switching between oxygen-rich and oxygen-deficient surface reconstructions, alternately activating and deactivating the surface for H 2 O formation. During periods of chemical inactivity, the oxygen-deficient surface undergoes slow reoxidation via lattice oxygen diffusing from subsurface and bulk reservoirs, restoring the active oxygen-rich surface termination. The inherent disparity in chemical activity among undercoordinated surface ions, along with sluggish subsurface-to-surface oxygen replenishment, drives this oscillatory redox cycle, modulating H2-induced loss of lattice oxygen at the surface and its delayed replenishment from the subsurface. This creates spatiotemporally separated redox steps at the oxide surface. The phenomena and atomistic insights presented here have significant implications for manipulating the surface reactivity of oxides by tuning the separation of these redox steps.

36 MATERIALS SCIENCE↗

DEM analysis of the effect of interface transition zone on dynamic splitting tensile behavior of high-strength concrete based on multi-phase model

Highlights: • Crushable aggregates with the real geometry shape are generated by using clump-cluster method. • Flat-Joint Model is adopted as the constitutive law. • The tensile strength is almost not affected by ITZ strength when its strength exceeds 0.8 times that of mortar. • A modified formula that can reflect the functional relationship between the ITZ strength and the TDIF is proposed. The tensile strength is one of the main parameters in the design of concrete structures and has an important influence on the cracking performance and durability of concrete materials. As the weak region in concrete, the strength of the interface transition zone (ITZ) directly affects the dynamic mechanical response of concrete. Meanwhile, considering that the tensile performance of concrete is far weaker than its compressive performance, the study of the strength of ITZ is of great significance to the dynamic tensile fracture characteristics and mechanical properties of concrete. In this paper, a method combining 3D scanning technology and Clump-Cluster was proposed to generate a real geometry shape model of crushable aggregate. On this basis, the 3D discrete element model (DEM) of high-strength concrete was constructed using the Flat-Joint Model (FJM). The effects of the ITZ strength on the failure mode, number of microcracks, tensile strength and strain rate effect of high-strength concrete under strain rates from 0.01/s to 100/s were studied by splitting tensile test. The simulation results indicate that the number of microcracks generated decreases with the increase of ITZ strength, and the proportion of ITZ cracks decreases obviously. When the ITZ strength exceeds the critical value of 0.8 times that of mortar, the tensile strength is almost not affected by ITZ strength. In addition, based on the numerical simulation results, a modified formula that can reflect the functional relationship between the ITZ strength and the dynamic increasing factor of the concrete tensile strength (TDIF) is proposed, and the rationality of the TDIF modified formula is verified.

36 MATERIALS SCIENCE↗

Shaping the Water-Harvesting Behavior of Metal–Organic Frameworks Aided by Fine-Tuned GPT Models

Here, we construct a data set of metal–organic framework (MOF) linkers and employ a fine-tuned GPT assistant to propose MOF linker designs by mutating and modifying the existing linker structures. This strategy allows the GPT model to learn the intricate language of chemistry in molecular representations, thereby achieving an enhanced accuracy in generating linker structures compared with its base models. Aiming to highlight the significance of linker design strategies in advancing the discovery of water-harvesting MOFs, we conducted a systematic MOF variant expansion upon state-of-the-art MOF-303 utilizing a multidimensional approach that integrates linker extension with multivariate tuning strategies. We synthesized a series of isoreticular aluminum MOFs, termed Long-Arm MOFs (LAMOF-1 to LAMOF-10), featuring linkers that bear various combinations of heteroatoms in their five-membered ring moiety, replacing pyrazole with either thiophene, furan, or thiazole rings or a combination of two. Beyond their consistent and robust architecture, as demonstrated by permanent porosity and thermal stability, the LAMOF series offers a generalizable synthesis strategy. Importantly, these 10 LAMOFs establish new benchmarks for water uptake (up to 0.64 g g –1 ) and operational humidity ranges (between 13 and 53%), thereby expanding the diversity of water-harvesting MOFs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling strategy for progressive failure prediction in lithium-ion batteries under mechanical abuse

The prediction for the internal failure of lithium-ion batteries (LIBs) under external mechanical abuse loading remains a challenge for safe design. This paper systematically studies the modeling approach for progressive failure simulation and short-circuit prediction. Hemispherical indentation experiments are conducted on a battery specimen containing 5 representative sandwich (RS) layers to examine the deformation and progressive failure behavior of battery components. Three different models, viz, a high-fidelity detailed model, an intermediate homogenized model and a fully homogenized model, are developed and employed to simulate the progressive failure in a multi-layer battery specimen. The simulation results obtained using the three models all correlate well with experimental phenomena, capturing the overall stress-strain response and ultimate failure of the battery structure. The numerical results can help us understand the failure behavior and the mechanism for mechanical abuse-induced short circuit in LIBs. The feasibility of the three models for use in progressive failure prediction is compared by analyzing the capability of the models in predicting the critical failure events, the computational efficiency of the models, and the feasibility to implement coupled mechanical-electrochemical-thermal simulations. The results of this work provide useful insights on practical choices for the modeling strategy and safety design of LIBs under mechanical abuse conditions.

25 ENERGY STORAGE↗

Sensitivity Analyses of Natural Convection in the AP1000 Passive Containment Cooling System Following LBLOCA Using CFD

The AP1000 power plant has implemented multiple layers of passive safety systems to enhance reactor safety and ensure system integrity during postulated scenarios, such as loss of coolant accident (LOCA) and main steam line break (MSLB). Among these, the passive containment cooling system (PCCS) is a safety-related system designed to prevent the containment from exceeding the design limits of both pressure and temperature following a postulated design basis accident, by transferring heat from the steel containment vessel to the atmosphere. During LOCA, natural air convection plays a vital role in removing heat from the steel containment vessel to the atmosphere. This is integrated with the condensate film from the Passive Containment Cooling Water Storage Tank (PCCWST). Air natural convection is also the only heat removal mechanism after the dry-out of the PCCWST and during the loss of shutdown decay heat removal (LOSDHR) event. Therefore, it is essential to investigate the thermal hydraulics behavior of the containment under transient conditions to ensure its safety and integrity. This paper presents a simplified CFD/ANSYS FLUENT model developed to analyze the impact of different parameters, such as air relative humidity, air temperature, and steel containment temperature, on the natural convection capability to remove the decay heat following LOCA. The model aims to examine the thermal behavior of the containment and provides a comprehensive understanding of the system's natural convection capability during postulated accidents. The results obtained from the model can be used to improve the safety and integrity of the containment system and provide valuable insights for future research in the field.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Polyubiquitin ligand-induced phase transitions are optimized by spacing between ubiquitin units

Biomolecular condensates form via multivalent interactions among key macromolecules and are regulated through ligand binding and/or posttranslational modifications. One such modification is ubiquitination, the covalent addition of ubiquitin (Ub) or polyubiquitin chains to target macromolecules. Specific interactions between polyubiquitin chains and partner proteins, including hHR23B, NEMO, and UBQLN2, regulate condensate assembly or disassembly. Here, we used a library of designed polyubiquitin hubs and UBQLN2 as model systems for determining the driving forces of ligand-mediated phase transitions. Perturbations to either the UBQLN2-binding surface of Ub or the spacing between Ub units reduce the ability of hubs to modulate UBQLN2 phase behavior. By developing an analytical model based on polyphasic linkage principles that accurately described the effects of different hubs on UBQLN2 phase separation, we determined that introduction of Ub to UBQLN2 condensates incurs a significant inclusion energetic penalty. This penalty antagonizes the ability of polyUb hubs to scaffold multiple UBQLN2 molecules and cooperatively amplify phase separation. The extent to which polyubiquitin hubs promote UBQLN2 phase separation is encoded in the spacings between Ub units. This spacing is modulated by chains of different linkages and designed chains of different architectures, thus illustrating how the ubiquitin code regulates functionality via the emergent properties of the condensate. The spacing in naturally occurring linear polyubiquitin chains is already optimized to promote phase separation with UBQLN2. We expect our findings to extend to other condensates, emphasizing the importance of ligand properties, including concentration, valency, affinity, and spacing between binding sites in studies and designs of condensates.

97 MATHEMATICS AND COMPUTING↗

Effects of local order parameter dependent transport coefficient in diblock copolymers under applied electric fields

In this work, we present an approach for constructing thermodynamically consistent time-dependent models relevant to thin films of diblock copolymers in applied electric fields. The approach is based on the principles of linear irreversible thermodynamics, and, in this work, it is applied to study the effects of electric fields on thin films of incompressible diblock copolymers. Enforcement of local incompressibility constraint at all times leads to a local order parameter dependent transport coefficient in the model for the diblock copolymers. The dependence of the transport coefficient on the local order parameter is used to relate it with the diffusion constant of Rouse chains and leads to sensitivity of the model to initial conditions. In addition, transient behavior is found to be affected when compared with an ad hoc model assuming a constant transport coefficient. Numerical results such as electric field induced alignment of lamellae domains due to the field are found to be in qualitative agreement with experiments. This approach opens up a systematic way of developing kinetic models for simulating effects of electrolytes added to thin films containing diblock copolymers in the presence of applied electric fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Model Based Validation of Intelligent Powertrain Strategies for Connected and Automated Vehicles

Systems incorporating Vehicle to Everything (V2X) and conventional cellular based communication in vehicles can significantly help improve energy consumption via a combination of intelligent powertrain control strategies, smarter routing algorithms and driving in such a way as to minimize fuel economy and the emission of carbon dioxide, known as "eco-driving." In projects led by the Southwest Research Institute (SwRI), large-scale traffic simulations are created to model real-world scenarios with dynamic behavior that is reactive to imposed changes. Coupled with high fidelity powertrain models, the closed loop framework enables research and development of such Connected and Automated Vehicle (CAV) enabled technologies at scale. This paper will discuss a traffic system simulation environment that was built based on the High Street urban corridor in Columbus, Ohio. Eco-driving strategies were tested at scale on a variety of powertrain platforms – internal combustion engines, hybrid electric and fully electric vehicles. Furthermore, the paper will focus on hybrid electric powertrain modeling along with details on how the powertrain model was leveraged to develop a sophisticated clustering scheme to help down-select speed traces from large scale simulation studies for validation on vehicle dynamometer. Nominal energy consumption improvement around 12% was observed with good match between simulation studies and vehicle testing.

33 ADVANCED PROPULSION SYSTEMS↗

Hydrogen Leak Modeling for Development of Smart Distributed Monitoring Under Unintended Releases

Hydrogen is a versatile and clean energy carrier that can be produced from various renewable sources such as wind, solar, and hydropower. Hydrogen has the potential to play a crucial role in decarbonizing industrial processes that are currently reliant on fossil fuels and provide long-duration and/or seasonal energy storage to enable electricity decarbonization. Hydrogen can also be used as a fuel for fuel cell vehicles, providing a zero-emission alternative to traditional internal combustion engines. DOE launched the Hydrogen Energy Earthshot (Hydrogen Shot) in June 2021 to reduce the cost of clean hydrogen by 80% to $1 per 1 kilogram in 1 decade ("1 1 1"). While promising, Hydrogen is highly-flammable, and in the presence of oxygen, it can form explosive mixtures. . Therefore, understanding leak scenarios is essential to evaluate and mitigate the safety risks associated with potential hydrogen leaks. An increased understanding of leak behavior, and having tools to model leaks, can help assess how hydrogen would disperse in different environments, influencing emergency response plans and safety measures, and identify potential issues with materials and design systems that can withstand the challenges posed by hydrogen. Recently, researchers have attempted to study hydrogen leaks for development of risk management strategies. However, the focus has been on closed or semi-closed spaces like storage rooms, vehicles, garages, and fueling stations - all promising locations for future hydrogen infrastructure. In this presentation, the modeling environment extends the span of research further by modeling hydrogen leak in an outdoor, open space. We will present the key challenges with modeling hydrogen leaks in an uncontrollable environment, how they were handled, and how modeling results informed sensor selection and placement. A Hydrogen research facility at the National Renewable Energy Laboratory (NREL) was used as a case study to model hydrogen leaks. In the future, Hydrogen wide area detection methodologies will be developed and tested at this site to monitor for unintended and operational hydrogen releases. The data generated from modeling will be used to develop a predictive model to detect hydrogen leak location based on concentration measured by sensors in this open space. Furthermore, the facility was also chosen because controlled hydrogen releases can be performed. A computational fluid dynamics (CFD) based modeling approach was taken to model hydrogen leak. The full-scale hydrogen facility was modeled with a large ambient domain. The electrolyzer at the facility can produce a controlled release rate of 27 kg-H2/hr. Site-specific atmospheric and weather condition data such as wind direction, wind speed at various altitudes, and temperature were used as inputs to the model. To capture the variability of weather conditions, a subset of the weather conditions experienced during daytime hours without precipitation over the course of three months was generated; using established data clustering techniques, a total of 100 condition sets were chosen. The results show statistical distributions and ranges of hydrogen concentrations at locations throughout the domain. These distributions are compared to experimental data from a constant mass flow, controlled hydrogen release at the facility. The stochastic wind conditions of the release make direct validation difficult, therefore, statistical comparison approaches were used. Wind conditions are found to significantly impact the release behavior, including direction and concentration. Sensor selection and placement is proposed for the facility and is now based on release behavior predicted for the facility given its weather patterns; this is much more informed than without the modeling results. The methodology and analysis procedure can be translated to other facilities using modified geometries and site-specific weather conditions. Hydrogen holds great promise as a renewable energy fuel, but ensuring safety in its production, storage, and use is paramount. Studying potential leak scenarios in an open space will help develop sensors to detect hydrogen on a large spectrum of concentration and eventually build a smart distributed monitoring system.

CFD↗

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗

Extending the Gutowski model to kidney-bean and elliptically shaped fibers

Efforts by manufacturers to produce more cost-effective carbon fibers have resulted in fibers with irregular cross-section often referred to as kidney-bean shaped fibers. In this research, compaction experiments were performed with a modified laser light section method to evaluate the compaction behavior of kidney-bean shaped carbon fibers. The kidney-bean shaped fibers followed a different compaction behavior compared to the Gutowski model for circular fibers. Additionally, these fibers required an order of magnitude larger force to compact than circular fibers to achieve similar fiber volume fraction, which has implications in infiltration and consolidation efficiency for composites manufacturing. A shape correction factor based on the fiber cross-sectional aspect ratio was proposed to extend the Gutowski model to fibers with irregular cross-sectional shapes. The modified Gutowski model provided an appropriate order of magnitude fit for the kidney-bean fibers. Furthermore, this modification to the Gutowski model recovered the original solution for circular fibers (cross-sectional aspect ratio = 1).

Materials Science↗

Implanted satellite transmitters affect sea duck movement patterns at short and long timescales

Abstract Studies of the effects of transmitters on wildlife often focus on survival. However, sublethal behavioral changes resulting from radio-marking have the potential to affect inferences from telemetry data and may vary based on individual and environmental characteristics. We used a long-term, multi-species tracking study of sea ducks to assess behavioral patterns at multiple temporal scales following implantation of intracoelomic satellite transmitters. We applied state-space models to assess short-term behavioral patterns in 476 individuals with implanted satellite transmitters, as well as comparing breeding site attendance and migratory phenology across multiple years after capture. In the short term, our results suggest an increase in dispersive behavior immediately following capture and transmitter implantation; however, behavior returned to seasonally average patterns within ~5 days after release. Over multiple years, we found that breeding site attendance by both males and females was depressed during the first breeding season after radio-marking relative to subsequent years, with larger relative decreases in breeding site attendance among males than females. We also found that spring and breeding migrations occurred later in the first year after radio-marking than in subsequent years. Across all behavioral effects, the severity of behavioral change often varied by species, sex, age, and capture season. We conclude that, although individuals appear to adjust relatively quickly (i.e. within 1 week) to implanted satellite transmitters, changes in breeding phenology may occur over the longer term and should be considered when analyzing and reporting telemetry data.

Lamb, Juliet S.↗

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗

Meso-scale modeling of UO2 nuclear fuel to high burnup

To improve the economics of light water reactors for commercial nuclear energy generation, utility operators are seeking to obtain regulatory approval to run UO2 fuel to higher levels of burnup. One potential impediment to obtaining this approval is the phenomenon of fuel fragmentation, relocation, and dispersal (FFRD). FFRD can result when fuel experiences a rapid temperature transient, such as that occurring during a Loss Of Coolant Accident (LOCA). FFRD has historically been most associated with the rim region in UO2 fuel pellets, where the phenomenon of fragmentation is also referred to as pulverization due to the small size of the fragments. More recent evidence suggests that the so-called “dark zone” (due to its appearance in micrographs) that can be observed in the mid-radial regions of high burnup fuel is also susceptible to FFRD. Although empirical fuel performance models have been developed that can adequately predict pulverization in the rim region under typical LWR conditions, a scientific understanding of what underlies fuel restructuring and subsequent FFRD is lacking even in the rim region, and no models are currently available for the behavior the dark zone. To address these challenges, the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has employed a multi-scale modeling approach to improve scientific understanding and develop new fuel performance models. In this talk, I will focus on meso-scale efforts, which form a crucial link between atomic-scale and engineering-scale models. Phase-field modeling combined with cluster dynamics is used to predict the restructuring process in the rim region. Phase-field fracture modeling, informed by atomistic simulations, is used to predict the onset of pulverization in the rim region. Combining these techniques together allows the extent of rim pulverization to be predicted. The formation and evolution of the dark zone has also been simulated with the phase-field method, using an improved approach to vacancy source term parameterization. The work shows the important impact of microstructure on fuel performance.

fracture↗

Uncertainty Quantification for Component Modeling Using the Discrete-Direct Approach

Threaded fastener behavior can be an important aspect of complex component and system behavior, but there is no one-size-fits-all finite element analysis technique. Proper modeling of threaded fastener joints requires careful consideration of many details, from test setup and data acquisition to constitutive modeling and uncertainty quantification approaches. This report details analysis of a “mini-radax” bolted-joint exemplar where a Discrete-Direct uncertainty quantification approach is employed to evaluate margin of the component. The mini-radax geometry is tested to failure on a drop table, and single-coupon tests of individual fasteners serve as foundational data for the analysis. Analysis predictions complement the test data well and provide additional context for engineering decision-making.

42 ENGINEERING↗

A Nonstationary and Non-Gaussian Moving Average Model for Solar Irradiance

Historically, power has flowed from large power plants to customers. Increasing penetration of distributed energy resources such as solar power from rooftop photovoltaic has made the distribution network a two-way-street with power being generated at the customer level. The incorporation of renewables introduces additional uncertainty and variability into the power grid. Distribution network operation studies are being adapted to include renewables; however, such studies require high quality solar irradiance data that adequately reflect realistic meteorological variability. Data from satellite-based products are spatially complete, but temporally coarse, whereas solar irradiances exhibit high frequency variation at very fine timescales. We propose a new stochastic method for temporally downscaling global horizontal irradiance (GHI) to 1 min resolution, but we do not consider the spatial aspect due to limited availability of the in situ irradiance measurements. Solar irradiance's first and second-order structures vary diurnally and seasonally, and our model adapts to such nonstationarity. Empirical irradiance data exhibits highly non-Gaussian behavior; we develop a nonstationary and non-Gaussian moving average model that is shown to capture realistic solar variability at multiple timescales. We also propose a new estimation scheme based on Cholesky factors of empirical autocovariance matrices, bypassing difficult and inaccessible likelihood-based approaches. The model is demonstrated for a case study of three locations that are located in diverse climates through the United States. The model is compared against competitors from the literature and is shown to provide better uncertainty and variability quantification on testing data.

Cholesky factor↗