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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Design of Multi-Stage Solvent Extraction Process for Separation of Rare Earth Elements

Flowsheet design and stage determination for the separation of rare earth elements (REEs) using solvent extraction (SX) is a challenging task because of the chemical similarity of the REEs. Low separation factors between the elements and complex equilibrium chemistry provide unique challenges to designing an efficient flowsheet for the separation of elements. The multi-stage nature of the SX process adds further complexity, making the assessment of products for a proposed design and stage combination difficult. Therefore, to develop a SX flowsheet, it is essential to quantify the performance for various design and separation conditions. This paper attempts to address the challenge by utilizing an equilibrium and process modeling approach. Results from a bench-scale study performed on a 10 g/L rare earth salt mixture were used in studying the extraction/stripping behavior and developing equilibrium models. DEHPA with TBP as a phase modifier was used as an extractant, while hydrochloric acid was utilized as a stripping agent. The results obtained were used in developing extraction/stripping models, which were integrated into a process framework of a SX train in a Matlab/Simulink environment. The models were programmed as a function block routine and used for developing a flowsheet, which was simulated for differing separation and design conditions. To identify optimum stage combinations, a particle swarm optimization (PSO) routine was developed and implemented for each SX train. Recovery and purity of elements of interest were used as objective function criteria. The stage combination leading to the minimization of the objective function was used to identify the optimum stage combination for a series of SX trains to attempt a balance of purity and recovery. The models and optimization method were implemented to separate a feed mixture containing REEs, which indicated that 99.52 and 85.41 percent purity is achievable for Yttrium and Lanthanum separation using 8-12-3 and 10-3-5 stage combination for loading, scrubbing, and striping. The model also indicated difficult separability between neodymium, praseodymium, and cerium.

Srivastava, Vaibhav (ORCID:0000000212645987)↗
DOE OSTI2025

Chain Flexibility and Structure of a Polyimide Copolymer: Revisiting the Freely Rotating Chain Model

Poly(4–4′-oxydiphenylene-pyromellitimide)-based polyimides─trade name Kapton─have wide-ranging engineering applications owing to their thermal and mechanical stability, but little is known about underlying chain-level characteristics. While theoretical models have conceptualized Kapton as inflexible polycyclic rods separated by freely rotating diphenyl ether hinge groups, the model’s core predictions remain untested and subtleties of the relaxation behavior are missed, which atomistic modeling can resolve. To these ends, we examine all-atom Kapton structures in crystalline and glassy amorphous configurations using a DFT-validated class II force field. Constructing amorphous configurations is challenging, as the fused-ring-containing backbone has slow relaxation dynamics and scaling suggestive of entanglements even in oligomers. In conclusion, we find larger backbone rearrangements of the linear polycyclic segments about ether groups that are consistent with the rod-hinge picture on the monomer scale, whereas a ring rotation analysis suggests partially flexible rod-like segments and involves multiple facile rotational relaxation modes.

Liesen, Nicholas T. [Lawrence Livermore National L↗
DOE OSTI2021

$\mathcal{g}$ modes of neutron stars with hadron-to-quark crossover transitions

Here we perform the first study of the principal core g-mode oscillation in hybrid stars containing quark matter, utilizing a crossover model for the hadron-to-quark transition inspired by lattice QCD. The ensuing results are compared with our recent findings of g-mode frequencies in hybrid stars with a first-order phase transition using Gibbs constructions. We find that models using Gibbs construction yield g-mode amplitudes and the associated gravitational energy radiated that dominate over those of the chosen crossover model owing to the distinct behaviors of the equilibrium and adiabatic sound speeds in the various models. Based on our results, we conclude that were g modes to be detected in upgraded LIGO and Virgo detectors it would indicate a first-order phase transition akin to a Gibbs construction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Model Driven Deception for Defense of Operational Technology Environments

Due to the strong integration of real-world physics, OT deception platforms must operate differently than traditional IT deceptions. For instance, turning off a valve will be detected downstream by other sensors because the flow will reduce and stop. Additionally, controllers and applications leverage data from sensors to send control commands to each other. A believable deception must be integrated with the system to project the effects of events. An attack will likely attempt to control the physical process in a negative manner. To make the attacker believe they are achieving their objective, it must predict the effects of these actions, to a reasonable degree. Our approach to simulating a model to generate realistic decoy behavior is explored including description of two approaches: a physics model-based approach and a data driven approach. The performance of two machine learning techniques are investigated in their ability to learn a good enough model of the physics of the system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Pavement condition and climatic data in southeast Texas: A dataset for evaluating flood impacts on pavement performance

Effective pavement maintenance is essential for economic stability, optimal network performance, and roadway safety. Achieving this requires thorough evaluation of pavement conditions, including structural integrity, surface roughness, and distress characteristics. Pavement performance indicators play a critical role in influencing vehicle safety and ride quality. Recent advances have emphasized the use of data-driven modeling to anticipate pavement behavior, with the goal of optimizing resource allocation and refining Maintenance and Rehabilitation (M&R) strategies through accurate condition assessment. A foundational requirement for these modeling efforts is the availability of standardized, high-quality datasets that can support robust and reproducible infrastructure analysis. This data article presents a comprehensive dataset assembled to facilitate pavement performance prediction, with a geographic focus on Southeast Texas, particularly the flood-vulnerable area of Beaumont. The dataset encompasses pavement and traffic attributes, meteorological records, flood simulation outputs, ground deformation measurements, and topographic indices, enabling detailed examination of both load-associated and non-load-associated degradation mechanisms. Data preprocessing was performed using ArcGIS Pro, Microsoft Excel, and Python to ensure consistency and usability in data-driven modeling applications, including machine learning workflows. Key contributions of this dataset include its utility in analyzing the climatic and environmental factors affecting pavement conditions, identifying critical predictive features, and enabling in-depth correlation analysis across diverse variables. By filling existing gaps in input variable selection resources, this dataset supports the development of predictive tools for estimating future maintenance demand and enhancing the resilience of pavement networks in flood-impacted areas. The resource highlights the importance of standardized datasets for advancing pavement management practices and provides a robust foundation for ongoing infrastructure performance modeling.

42 ENGINEERING↗
DOE OSTI2020

A CALPHAD-informed approach to modeling constituent redistribution in Zr-based metallic fuels using BISON

Here, a CALPHAD-informed (Computer Coupling of Phase Diagrams and Thermochemistry) constituent redistribution model was developed for Zr-based metallic fuels and incorporated into the BISON fuel performance code. Three uncertain model parameters associated with β and γ phase kinetics were calibrated using integral test data from U-Zr fuel elements irradiated in Experimental Breeder Reactor II. The calibrated constituent redistribution model was shown to predict the behavior of U-Zr fuels with excellent accuracy. Model predictions for U-Pu-Zr fuels were physically reasonable but less accurate. Reduction of uncertainties in the ternary phase transition temperatures and collection of kinetic data for the ζ phase are expected to improve the model’s ternary predictions. Finally, the new model was coupled to existing thermomechanics models in BISON to simulate irradiation of an entire U-Zr fuel element, demonstrating its ability to accurately predict the behavior of U-Zr fuels with realistic geometries and mesh resolutions at the engineering scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Scattering neutrinos, spin models, and permutations

We consider a class of Heisenberg all-to-all coupled spin models inspired by neutrino interactions in a supernova with N degrees of freedom. These models are characterized by a coupling matrix that is relatively simple in the sense that there are only a few, relative to N , nontrivial eigenvalues, in distinction to the classic Heisenberg spin-glass models, leading to distinct behavior in both the high-temperature and low-temperature regimes. When the momenta of the neutrinos are uniform and random in directions, we can calculate the large- N partition function for the quantum Heisenberg model. In particular, the high-temperature partition function predicts a non-Gaussian density of states, providing interesting counterexamples showing the limits of general theorems on the density of states for quantum spin models. We can repeat the same argument for classical Heisenberg models, also known as rotor models, and we find the high-temperature expansion is completely controlled by the eigenvalues of the coupling matrix, and again predicts non-Gaussian behavior for the density of states as long as the number of eigenvalues does not scale linearly with N . Indeed, we derive the amusing fact that these partition functions are essentially the generating function for counting permutations in the high-temperature regime. Finally, for the case relevant to neutrinos in a supernova, we identify the low-temperature phase as a unique state with the direction of the momenta of the neutrino dictating its coherent state in flavor-space, a state we dub the “flavor-momentum-locked” state. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Impact analysis of DERs on bulk power system stability through the parameterization of aggregated DER_a model for real feeders

With an ever increasing percentage of distributed energy resources (DERs) connected behind the meter in the distribution system, it is becoming increasingly important to equip transmission planners with the visibility of DER dynamic performance in distribution system. Not having visibility of the disconnection of DERs with the occurrence of transmission events, could result in an erroneous view of the stability of bulk power system. Here, in this paper, a parameterized aggregated model (DER_a) serves as a representation of the distribution-level dynamics of real residential feeders, which is used for analysis of bulk power system stability. The parameters are obtained by executing dynamic Monte Carlo simulations. Faults are then induced at the substation level causing the DERs to trip which subsequently enables the parameterization of the low and high voltage breakpoints (v l0 , v l1 , v h0 , and v h1 ) of the DER_a model’s partial voltage trip block. These parameters are then utilized to study the effectiveness of the DER_a model to represent the behavior of the aggregated DERs’ response and their impact on the bulk power system. The case study shows the ability of the positive-sequence DER_a model to provide an accurate estimation of DERs that are susceptible to trip due to 3-$\phi$ and 1-$\phi$ faults of transmission.

14 SOLAR ENERGY↗
DOE OSTI2023

Comparative determinant quantum Monte Carlo study of the acoustic and optical variants of the Su-Schrieffer-Heeger model

We compare the acoustic Su–Schrieffer–Heeger (SSH) model with two of its optical variants where the phonons are defined on either on the sites or bonds of the system. First, we discuss how to make fair comparisons between these models in any dimension by ensuring their dimensionless coupling λ and relevant phonon energies are the same. We then use determinant quantum Monte Carlo to perform nonperturbative and sign-problem-free simulations of all three models on one-dimensional chains at and away from half-filling. By comparing the results obtained from each model, we demonstrate that the optical and acoustic models produce near identical results within error bars for suitably chosen phonon energies and λ at half-filling. In contrast, the bond model has quantitatively different behavior due to its coupling to the q = 0 phonon mode. These differences also manifest in the total length of the chain, which shrinks for the bond model but not for the acoustic and optical models when λ ≠ 0. Furthermore, our results have important implications for quantum Monte Carlo modeling of SSH-like interactions, where these models are sometimes regarded as being interchangeable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Discrete element model for powder grain interactions under high compressive stress

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. Finally, as validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Influence of alloy solidification path on melt pool behavior in additive manufacturing

Numerical models used to study transport phenomena in laser-powder bed fusion processes often rely on assumptions and simplifications to reduce their computational expense. One common simplification is in the description of latent heat evolution during the solid-liquid phase change (i.e., the solidification pathway), justified by the fact that the mushy zone thickness is similar to the numerical grid spacing used for continuum transport models. The lack of resolution of transport phenomena in the mushy zone motivates the use of computationally convenient solidification paths such as linear or sigmoidal relationships over pathways derived from fundamental solidification theory such as equilibrium or Scheil models. In the present work, an uncertainty quantification (UQ) framework is used to analyze the influence of solidification pathway selection on the solidification dynamics and melt pool geometries in laser based additive manufacturing (AM) of IN625. Results show the solidification pathway has a quantifiable influence on the cooling rate at the liquidus isotherm, mushy zone thickness, and solidification time. Due to similarities in the latent heat evolution at the beginning of solidification, the equilibrium and Scheil models predict similar cooling rates near the liquidus isotherm, however the wider freezing range of Scheil leads to a wider mushy zone compared to equilibrium. The non-physical latent heat release profiles of sigmoidal and linear paths lead to significant overpredictions of cooling rates at the liquidus isotherm compared to equilibrium and Scheil. Finally, these results indicate that careful consideration should be given to the choice of solidification pathway to ensure reliable model predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Comparison of PV Resource Modeling for Sizing Microgrid Components

Microgrid systems are being deployed with increased frequency to meet critical building loads during unplanned power outages. Determining the appropriate type and size of components that make up a microgrid can greatly affect its ability to meet this need. This paper presents a new tool to appropriately size PV, battery, and generator capacities to meet all site loads, given a resilience goal. It uses a statistical model to simulate microgrid behavior under a large range of solar resource conditions, providing resilience planners with confidence that critical loads will be met even under extreme weather conditions. We compare this model with one that uses typical weather data for sizing microgrids, and quantify the resilience effects of using these different methods. Our method produces more conservative generator capacities and fuel requirements in an outage, predicting up to 30% more fuel required in the most extreme case. Microgrid systems that are designed using more conservative estimates are more likely to continue to serve load during emergency outages under a large range of conditions.

Newman, Sarah F.↗
DOE OSTI2023

Effect of Na Concentration on Cs Distribution with Crystalline Silicotitanate in Tank Waste Simulants

Crystalline silicotitanate (CST) ion exchange media is currently utilized in the tank side cesium removal (TSCR) system on the Hanford site to remove Cs-137 from the tank waste supernate. As the main dose contributor to the liquid supernate in Hanford tank waste, it is necessary to remove the Cs-137 to expedite low activity waste processing and immobilization. Ongoing batch contact experiments with CST have been done to tease out the impact of group II metals, anions, potassium, and temperature on Cs removal to better understand bounding conditions for operations as well as aid in development of a wholesome isotherm model to predict Cs behavior in the tank waste. This paper discusses the impact of Na concentration on Cs exchange and provides insight into the challenges associated with activity coefficient estimations for modeling Cs distribution in tank waste matrices.

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

Type 2 diabetes impairs annulus fibrosus fiber deformation and rotation under disc compression in the University of California Davis type 2 diabetes mellitus (UCD-T2DM) rat model

Abstract Understanding the biomechanical behavior of the intervertebral disc is crucial for studying disease mechanisms and developing tissue engineering strategies for managing disc degeneration. We used synchrotron small-angle X-ray scattering to investigate how changes to collagen behavior contribute to alterations in the disc’s ability to resist compression. Coccygeal motion segments from 6-month-old lean Sprague-Dawley rats ( n=7) and diabetic obese University of California Davis type 2 diabetes mellitus (UCD-T2DM) rats ( n=6, diabetic for 68±7 days) were compressed during simultaneous synchrotron scanning to measure collagen strain at the nanoscale (beamline 7.3.3 of the Advanced Light Source). After compression, the annulus fibrosus was assayed for nonenzymatic cross-links. In discs from lean rats, resistance to compression involved two main energy-dissipation mechanisms at the nanoscale: (1) rotation of the two groups of collagen fibrils forming the annulus fibrosus and (2) straightening (uncrimping) and stretching of the collagen fibrils. In discs from diabetic rats, both mechanisms were significantly impaired. Specifically, diabetes reduced fibril rotation by 31% and reduced collagen fibril strain by 30% (compared to lean discs). The stiffening of collagen fibrils in the discs from diabetic rats was consistent with a 31% higher concentration of nonenzymatic cross-links and with evidence of earlier onset plastic deformations such as fibril sliding and fibril–matrix delamination. These findings suggest that fibril reorientation, stretching, and straightening are key deformation mechanisms that facilitate whole-disc compression, and that type 2 diabetes impairs these efficient and low-energy elastic deformation mechanisms, thereby altering whole-disc behavior and inducing the earlier onset of plastic deformation.

Rosenberg, James L.↗
DOE OSTI2021

Response of U.S. West Coast Mountain Snowpack to Local Sea Surface Temperature Perturbations: Insights from Numerical Modeling and Machine Learning

Sea surface temperature (SST) significantly modulates the precipitation and temperature over land, with important consequences on land surface processes such as snowpack. Compared to the impact of remote SST, the effect of nearshore/local SST is less well understood. In this study, the impact of local SST on the mountain snowpack of the U.S. West Coast is investigated using two 6-km regional climate simulations driven by the same lateral boundary conditions but with time-varying versus time-invariant and warmer local SSTs during 2003–15. Results show that local SST warming leads to warmer winter with more precipitation over the mountains. Meanwhile, the removal of SST temporal variability results in reduced temperature variability but increased precipitation variability. As a result, winter snow accumulation decreases by ~200 mm per season in the Cascade Mountains in the north but increases by ~100 mm per season in the Sierra Nevada in the south. Such a dipole response results from the competing effects of precipitation and temperature change at different elevations and are amplified by the enhanced atmospheric river moisture transport. To further delineate the relative contributions of different meteorological factors to the snowpack response, two neural network models were developed to predict the snow behaviors at seasonal and monthly scales. These models reveal the dominant influence of the total amount and the average temperature of precipitation on the snowpack response. Furthermore, these findings highlight the sensitivity of mountain snowpack to local SST in the western United States and underscore the importance of local SST and atmospheric rives to accurate snowpack estimations for water management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Scaling Behavior of a Turbulent Kinetic Energy Closure Scheme for the Stably Stratified Atmosphere: A Steady-State Analysis

Abstract We present a turbulent kinetic energy (TKE) closure scheme for the stably stratified atmosphere in which the mixing lengths for momentum and heat are not parameterized in the same manner. The key difference is that, while the mixing length for heat tends toward the stability independent mixing length for momentum in neutrally stratified conditions, it tends toward one based on the Brunt–Väisälä time scale and square root of the TKE in the limit of large stability. This enables a unique steady-state solution for TKE to be obtained, which we demonstrate would otherwise be impossible if the mixing lengths were the same. Despite the model’s relative simplicity, it is shown to perform reasonably well with observational data from the 1999 Cooperative Atmosphere–Surface Exchange Study (CASES-99) using commonly employed model constants. Analyzing the scaling behavior of the nondimensional velocity and potential temperature gradients, or of the stability (correction) functions, reveals that for large stability the present model scales in the same manner as the first-order operational scheme of Viterbo et al. Alternatively, it appears as a blend of two cases of the TKE closure scheme of Baas et al. Critically, because a unique steady-state TKE can be obtained, the present model avoids the nonphysical behavior identified in one of the cases of Baas et al.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Ion exchange behavior of astatine and bismuth

Astatine and bismuth sorption on several ion exchange resins from nitric acid media has been studied. This work covers commercially available resins Dowex 50 × 4, MP thiol, Dowex 1 × 8, and TEVA. One of the main advantages of using ion exchangers for astatine separation and purification is the absence of any organic media in the effluent in comparison with extraction chromatography resins. The behavior of the above-mentioned metals was investigated in up to 4 M HNO 3 solutions. The determined distribution coefficients are greater than 1 for all the resins studied, reaching 10 4 for MP thiol and TEVA resins. The interaction of nitric acid with an anion exchanger in Cl-form results in an exchange reaction, where AtO + chloro complexes can form. To understand the astatine behavior under these conditions, stability constants of AtOCl and AtOCl 2 – complexes have been reevaluated and discrepancy in the literature values has been eliminated. For each resin a thermodynamic model has been developed to suggest a possible mechanism of astatine sorption. Literature and new experimental data on bismuth sorption by studying resins in nitric acid media have been reviewed and a mathematical model to describe its behavior has been suggested. Furthermore, a ratio of corresponding fit functions of At and Bi assigned to the same resin and acidity has been used to estimate the separation factors of these elements.

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

Effect of particle size on the capture of uranium oxide colloidal particles from aqueous suspensions via high-gradient magnetic filtration

The effectiveness of High Gradient Magnetic Filtration (HGMF) in capturing uranium oxide particles from suspensions was investigated in this study. Two sets of experiments were performed to evaluate the importance of size on the capture of uranium oxide particles. The first considered two batches sieved into size bins of< 5, 5–10, 10–15, and 15–20 µm, while the second was performed using two suspensions with diameters smaller than 1.0 µm and between 1.0 and 1.5 µm. Iron oxide experiments, with particles between 0.3 and 0.8 µm, were performed for calibration purposes. In all experiments, a surfactant (Triton-X100 or sodium dodecyl sulfate) was used to prevent particle aggregation and limit the influence of non-magnetic capture mechanisms. A magnetic field of approximately 1.1 Tesla was generated using a water cooled electromagnet. HGMF was performed using tubular filters packed with ferromagnetic stainless-steel wool. Of the initial four uranium oxide particle sizes, magnetic capture was only observed for particles with a diameter of less than 5 µm, while larger particles experienced no magnetic and minimal total capture. For particles with diameters smaller than 1.0 µm and between 1.0 and 1.5 µm, capture efficiencies increased by 39 ± 9% and 34 ± 6% respectively, solely due to the magnetic field. Although the magnetic force is proportional to particle diameter, the capture efficiency decreased as diameter increased. So these results suggest that Brownian diffusion, which is influential for micron sized particles and increases with decreasing particle size, is acting in conjunction with the magnetic force to influence the efficacy of HGMF for uranium oxide. This important finding underscores the effectiveness of Brownian diffusion in increasing the rate of collision between particles and collector fibers. A stochastic trajectory model was developed to incorporate the influence of Brownian motion on particle behavior and filter removal efficiency. Modeling results are discussed and compared for uranium and iron oxide particles.

42 ENGINEERING↗