Lume-model
Lume-model defines data structures used in the LUME modeling tool set.
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.
Lume-model defines data structures used in the LUME modeling tool set.
Choosing an appropriate group structure for multigroup transport is far from an exact science. For some applications, one blindly uses a group structure developed years ago by forgotten methods. Furthermore, one sometimes uses the same group structure for a variety of problems, even if the group structure was originally developed with a certain application in mind. In this work, we create optimized group structures with simulated annealing for critical assembly test problems and apply a random forest regressor with bagging to choose the best group structure based on parameters of the different test problems. The optimized group structures were generated using a simulated annealing optimizer for several simple, spherical, and unreflected problems. The optimization was performed to minimize a cost function that included fission rate, absorption rate, leakage, and k{sub eff}. A random forest regressor was then trained on a set of International Criticality Safety Benchmark Evaluation Project inputs and used to select one of these six group structures. The trained machine learning model chose the best group structure 65% of the time, and one of the three best 89% of the time. Furthermore, it decreased the L2 error over all test problems by a factor of 25 when compared to the standard Los Alamos 70-group structure. In other words, the model chose group structure that were far more appropriate for the test problems than the traditional LANL group structure. (authors)
Proteins with intrinsic or unfolded state disorder comprise a new frontier in structural biology, requiring the characterization of diverse and dynamic structural ensembles. Here we introduce a comprehensive Bayesian framework, the Extended Experimental Inferential Structure Determination (X-EISD) method, which calculates the maximum log-likelihood of a disordered protein ensemble. X-EISD accounts for the uncertainties of a range of experimental data and back-calculation models from structures, including NMR chemical shifts, J-couplings, Nuclear Overhauser Effects (NOEs), paramagnetic relaxation enhancements (PREs), residual dipolar couplings (RDCs), hydrodynamic radii (R h ), single molecule fluorescence Förster resonance energy transfer (smFRET) and small angle X-ray scattering (SAXS). We apply X-EISD to the joint optimization against experimental data for the unfolded drkN SH3 domain and find that combining a local data type, such as chemical shifts or J-couplings, paired with long-ranged restraints such as NOEs, PREs or smFRET, yields structural ensembles in good agreement with all other data types if combined with representative IDP conformers.
Real-world datasets in chemical engineering and bioengineering processes—such as those from catalytic reactors, multiphase flows, polymerization reactors, bioreactors, and clinical trials—can often be unlabeled or disorganized, rendering the training of existing supervised learning models ineffective at learning the underlying dynamics. To salvage these datasets for decision-making, we first seek to obtain clarity from the cluttered data. Here, we present a framework for developing “structural” generative models, discovering emergent equations, and constructing efficient emulators from scrambled datasets by integrating unsupervised organizational learning techniques (Questionnaires) with advanced deep learning architectures (Deep Hidden Physics Models and Deep Operator Networks). Our approach is demonstrated on two illustrative model systems: (a) a 1D advection–diffusion partial differential equation representing a winding underground pipe and (b) an ensemble of Stuart–Landau oscillators, an agent-based system of coupled ordinary differential equations. In both cases, we successfully reconstruct meaningful spatial, temporal, and parameter embeddings from scrambled data, enabling good predictions of system dynamics. As a result, we highlight the framework’s potential for broader applications, enabling data-driven system identification in fields with inherently disorganized or hidden parameter spaces.
Here, the advent of high-throughput Density Functional Theory (DFT) calculations has supported the creation of large databases containing the quantitative output necessary for constructing theoretical phase diagrams and predicting physical properties. In this Article, we present a complementary resource, the Intermetallic Reactivity Database (IRD), focused on the chemical bonding features of solid-state structures and indicators of potential structural transformations. Each IRD entry augments common features, such as band structures and density of states (DOS) distributions, with chemically motivated information including DFT-Chemical Pressure (CP) schemes and visualizable representations of the atomic charges. Together, these data types enable the rationalization and prediction of potential structural phenomena encountered in intermetallic chemistry, as we illustrate with four examples: the origins of the Y 2 Ni 2 Mg structure in terms of CP features of its parent structures, the anticipation of intergrowth phases from the net atomic CPs collected in Al-containing binary phases, the correlation between trends in the CP schemes of CaCu 5 -type phases and experimentally observed structural variations, and finally, the development of theoretical methodology with the testing of a streamlined method generating DFT-CP schemes. Altogether, these examples highlight how the IRD supports the creation of models of structural chemistry that extend beyond the bounds of its entries.
A new, three-dimensional, analytical, steady-state wake model is presented that includes local flow acceleration near the rotor, improving the wake description compared to existing models. Wake structures such as the momentum deficit and regions of accelerated flow are concisely described with compound and normal Gaussian functions. Large-eddy simulations (LES) are used as training data to develop the model using two in-line turbines under various inflow conditions parameterized by hub-height wind speed and turbulence intensity. Mass conservation is considered by fixing two components of the wake velocity model and optimizing the third to best satisfy continuity; after which, the model performs comparably if not better than existing work with regards to both relative error and mass consistency. The final model demonstrates a high degree of flexibility, making use of empirical correlations to scale across different inflow conditions. The inclusion of these effects is capable of revealing unused opportunities for enhanced power generation by aligning wake trajectories with these regions of accelerated flow.
Topological solitons commonly appear as energy-minimizing field configurations, but examples of stable, spatially localized objects with coexisting solitonic structures and singular defects are rare. Here we use a nonpolar chiral liquid crystal system to show how twist domain walls can co-self-assemble with vortices to form spatially localized topological objects with spontaneous folding. These soliton–vortex assemblies, which we call ‘möbiusons’, have a topology of the molecular alignment field resembling that of the Möbius strip’s surface and package localized field excitations into folded structures within a confinement-frustrated uniform far-field background. Upon supplying energy in the form of electric pulses, möbiusons with different overall symmetries of structure exhibit folding-dependent rotational and translational motions, as well as topological cargo-carrying abilities that can be controlled by tuning the amplitude and frequency of the applied fields. We demonstrate on-demand transformations between various möbiusons and show examples of encoding information by manipulating folds in such structures. A model based on the energetics of solitons and vortices provides insights into the origins of the folding instability, whereas minimization of the Landau–de Gennes free energy closely reproduces details of their internal structure. Our findings may provide a route towards topology-enabled light-steering designs.
This study presents the executive disruption model (EDM) of tinnitus distress and subsequently validates it statistically using two independent datasets (the Construction Dataset: n = 96 and the Validation Dataset: n = 200). The conceptual EDM was first operationalised as a structural causal model (construction phase). Then multiple regression was used to examine the effect of executive functioning on tinnitus-related distress (validation phase), adjusting for the additional contributions of hearing threshold and psychological distress. For both datasets, executive functioning negatively predicted tinnitus distress score by a similar amount (the Construction Dataset: β = −3.50, p = 0.13 and the Validation Dataset: β = −3.71, p = 0.02). Theoretical implications and applications of the EDM are subsequently discussed; these include the predictive nature of executive functioning in the development of distressing tinnitus, and the clinical utility of the EDM.
This poster presents a novel method for validation and verification that focuses on identifying causal relationships between data elements, moving beyond traditional statistical and machine learning approaches. These methods employ causal discovery techniques to reveal the underlying mechanisms of data generation. The research utilizes structural causal models and directed acyclic graphs to depict causal relationships. This approach assists in achieving alignment between simulation models and reality.
Biochar is an engineered carbon-rich substance used for soil improvement, environmental management, and other diverse applications. To date, the understanding of how biomass affects biochar microstructure has been limited due to the complexity of analysis involved in tracing the changes in the physical structure of biomass as it undergoes thermochemical conversion. In this study, we used synchrotron x-ray micro-tomography to visualize changes in the internal structure of biochar from diverse feedstock (miscanthus straw pellets, wheat straw pellets, oilseed rape straw pellets, and rice husk) during pyrolysis by collecting a sequence of 3D scans at 50 °C intervals during progressive heating from 50 °C to 800 °C. The results show a strong dependence of biochar porosity on feedstock as well as pyrolysis temperature, with observed porosity in the range of 7.41–60.56%. Our results show that the porosity, total surface area, pore volume, and equivalent diameter of the largest pore increases with increasing pyrolysis temperature up to about 550 °C. The most dramatic development of pore structure occurred in the temperature range of 350–450 °C. This understanding is pivotal for optimizing biochar’s properties for specific applications in soil improvement, environmental management, and beyond. By elucidating the nuanced variations in biochar’s physical characteristics across different production temperatures and feedstocks, this research advances the practical application of biochar, offering significant benefits in agricultural, environmental, and engineering contexts.
Abstract Global warming modulates soil respiration (RS) via microbial decomposition, which is seasonally dependent. Yet, the magnitude and direction of this modulation remain unclear, partly owing to the lack of knowledge on how microorganisms respond to seasonal changes. Here, we investigated the temporal dynamics of soil microbial communities over 12 consecutive months under experimental warming in a tallgrass prairie ecosystem. The interplay between warming and time altered (P < 0.05) the taxonomic and functional compositions of microbial communities. During the cool months (January to February and October to December), warming induced a soil microbiome with a higher genomic potential for carbon decomposition, community-level ribosomal RNA operon (rrn) copy numbers, and microbial metabolic quotients, suggesting that warming stimulated fast-growing microorganisms that enhanced carbon decomposition. Modeling analyses further showed that warming reduced the temperature sensitivity of microbial carbon use efficiency (CUE) by 28.7% when monthly average temperature was low, resulting in lower microbial CUE and higher heterotrophic respiration (Rh) potentials. Structural equation modeling showed that warming modulated both Rh and RS directly by altering soil temperature and indirectly by influencing microbial community traits, soil moisture, nitrate content, soil pH, and gross primary productivity. The modulation of Rh by warming was more pronounced in cooler months compared to warmer ones. Together, our findings reveal distinct warming-induced effects on microbial functional traits in cool months, challenging the norm of soil sampling only in the peak growing season, and advancing our mechanistic understanding of the seasonal pattern of RS and Rh sensitivity to warming.
Low-density polymer foams of varying sizes, shapes, and densities are of specific interest to the inertial confinement fusion (ICF) program and related high-energy density plasma physics research. Historically, these foams are comprised of polystyrene or other low atomic number materials and have densities in the 30 to 300 mg/cm 3 range. However, at the lower end of this density range, these traditional polymer foams become fragile and difficult to cast and machine into the geometries needed. Recently, the need by experimentalists for materials with densities below 30 mg/cm 3 has increased. To address these needs, we are developing three-dimensional (3-D) printing techniques to create high-precision, low-density, and repeatable complex lattice structures. Using two-photon polymerization 3-D printing, we recently developed the first 5 mg/cm3 low-density lattice structure having an annular hemispherical shape. These microscale to mesoscale structures were modeled and designed using the nTopology software, specifically utilizing the "Voronoi volume lattice" and "random points in body" option blocks. All printing operations were performed using the Nanoscribe Photonic Professional GT instrument. Characterization of these 3-D structures was conducted using various microscopic and X-ray tomographic imaging techniques. Furthermore, overall printed part sizes ranged from 1 to 5 mm in diameter and were composed of lattice ligaments having thicknesses in the 3- to 5-µm range. These structures have been incorporated into ICF targets recently shot on both the University of Rochester’s Laboratory of Laser Energetics Omega laser and the National Ignition Facility.
The pore structure and density of lignocellulosic feedstocks dictate intraparticle transport phenomena and thereby play an important role in thermochemical conversion processes such as fast pyrolysis for biofuel and biochemical production. Variations in microstructure are inherent from different biomass species and can be introduced by preprocessing techniques such as cutting and pelletization. Morphological changes also occur during conversion and lead to vastly different pore structures and behavior during pyrolysis, which impact required conversion times and product distributions. The current work presents a comprehensive comparison of fast pyrolysis of neat and pelletized pine feedstocks, which includes single-particle experiments, modeling, and 3D imaging by X-ray computed tomography (XCT). The particle-scale model included anisotropic heat and mass transport in a shrinking particle with pyrolysis reactions based on the CRECK mechanism with boundary conditions informed by reactor-scale simulations of the single-particle reactor. The models were validated by measurements of the temperature and mass loss from single-particle pyrolysis experiments of neat and pelletized pine. Quantitative analysis of XCT geometries revealed that pyrolytic conversion yielded chars with increased porosity and permeability compared to the unpyrolyzed materials, along with decreased tortuosity and anisotropy. Pelletization of the pine feedstock resulted in a much denser, less permeable material, which converted slower and produced more residual char after pyrolysis compared to neat pine. The results from particle modeling revealed that accounting for the dynamic and anisotropic heat and mass transport caused by differences in pore structure is critical to achieving agreement with experimental results. Overall, this study highlights the dramatic differences in conversion behavior imparted by pelletization and the importance of capturing microstructural attributes in computational models to guide the design and optimization of pyrolysis processes for specific biomass feedstocks.
Since the discovery of its electrochemical activity, Li-excess disordered rock-salt (DRX) cathode material has received worldwide attention as it sets up a new way to exploit oxygen redox beyond the conventional layered structure with late-3d transition metals. However, the intricate structure-function relationship in the disordered lattice of the DRX material fogs the researcher's lens on the underlying redox mechanisms. Here, in this study, we employ a synergistic approach combining neutron total scattering with reverse Monte Carlo modeling and density functional theory calculations to unravel the landscape of oxygen redox reactions in DRX. Redox activities are evaluated in diverse oxygen clusters (OLi x TM 6-x ) and the spatial distribution of these clusters in the model DRX structure (Li 1.16 Ti 0.37 Ni 0.37 Nb 0.1 O 2 and Li 1.2 Ti 0.35 Ni 0.35 Nb 0.1 O 1.8 F 0.2 ) is explicitly determined. The results unveil that by regulating the short-range ordering between cations, fluorine atoms can effectively decouple the location of Li extraction and electron depletion. Such disentanglement between the Li reservoir and electron reservoir in the DRX lattice could play a pivotal role in protecting the oxidized oxygen and preserving the lattice framework during cycling. Through a tentatively designed non-fluorinated DRX oxide realizing similar Li-electron decoupling, an obvious enhancement of the cycling capability can be achieved without compromising the capacity release.
RNA-Puzzles is a collective endeavor dedicated to the advancement and improvement of RNA three-dimensional structure prediction. With agreement from structural biologists, RNA structures are predicted by modeling groups before publication of the experimental structures. We report a large-scale set of predictions by 18 groups for 23 RNA-Puzzles: 4 RNA elements, 2 Aptamers, 4 Viral elements, 5 Ribozymes and 8 Riboswitches. We describe automatic assessment protocols for comparisons between prediction and experiment. Our analyses reveal some critical steps to be overcome to achieve good accuracy in modeling RNA structures: identification of helix-forming pairs and of non-Watson–Crick modules, correct coaxial stacking between helices and avoidance of entanglements. Three of the top four modeling groups in this round also ranked among the top four in the CASP15 contest.
Billions of organic molecules have been computationally generated, yet functional inorganic materials remain scarce due to limited data and structural complexity. Here, in this work, we introduce Structural Constraint Integration in a GENerative model (SCIGEN), a framework that enforces geometric constraints, such as honeycomb and kagome lattices, within diffusion-based generative models to discover stable quantum materials candidates. SCIGEN enables conditional sampling from the original distribution, preserving output validity while guiding structural motifs. This approach generates ten million inorganic compounds with Archimedean and Lieb lattices, over 10% of which pass multistage stability screening. High-throughput density functional theory calculations on 26,000 candidates shows over 95% convergence and 53% structural stability. A graph neural network classifier detects magnetic ordering in 41% of relaxed structures. Furthermore, we synthesize and characterize two predicted materials, TiPd 0.22 Bi 0.88 and Ti 0.5 Pd 1.5 Sb, which display paramagnetic and diamagnetic behaviour, respectively. Our results indicate that SCIGEN provides a scalable path for generating quantum materials guided by lattice geometry.
A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high-fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not predicted by potential flow engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.
A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not captured in engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.