Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive model”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Results of a Geant4 benchmarking study for bio‐medical applications, performed with the G4‐Med system

Geant4, a Monte Carlo Simulation Toolkit extensively used in bio-medical physics, is in continuous evolution to include newest research findings to improve its accuracy and to respond to the evolving needs of a very diverse user community. In 2014, the G4-Med benchmarking system was born from the effort of the Geant4 Medical Simulation Benchmarking Group, to benchmark and monitor the evolution of Geant4 for medical physics applications. The G4-Med system was first described in our Medical Physics Special Report published in 2021. Results of the tests were reported for Geant4 10.5. Purpose In this work, we describe the evolution of the G4-Med benchmarking system. Methods The G4-Med benchmarking suite currently includes 23 tests, which benchmark Geant4 from the calculation of basic physical quantities to the simulation of more clinically relevant set-ups. New tests concern the benchmarking of Geant4-DNA physics and chemistry components for regression testing purposes, dosimetry for brachytherapy with a 125 I source, dosimetry for external x-ray and electron FLASH radiotherapy, experimental microdosimetry for proton therapy, and in vivo PET for carbon and oxygen beams. Regression testing has been performed between Geant4 10.5 and 11.1. Finally, a simple Geant4 simulation has been developed and used to compare Geant4 EM physics constructors and physics lists in terms of execution times. Results In summary, our EM tests show that the parameters of the multiple scattering in the Geant4 EM constructor G4EmStandardPhysics_option3 in Geant4 11.1, while improving the modeling of the electron backscattering in high atomic number targets, are not adequate for dosimetry for clinical x-ray and electron beams. Therefore, these parameters have been reverted back to those of Geant4 10.5 in Geant4 11.2.1. The x-ray radiotherapy test shows significant differences in the modeling of the bremsstrahlung process, especially between G4EmPenelopePhysics and the other constructors under study (G4EmLivermorePhysics, G4EmStandardPhysics_option3, and G4EmStandardPhysics_option4). These differences will be studied in an in-depth investigation within our Group. Improvement in Geant4 11.1 has been observed for the modeling of the proton and carbon ion Bragg peak with energies of clinical interest, thanks to the adoption of ICRU90 to calculate the low energy proton stopping powers in water and of the Linhard–Sorensen ion model, available in Geant4 since version 11.0. Nuclear fragmentation tests of interest for carbon ion therapy show differences between Geant4 10.5 and 11.1 in terms of fragment yields. In particular, a higher production of boron fragments is observed with Geant4 11.1, leading to a better agreement with reference data for this fragment. Conclusions Based on the overall results of our tests, we recommend to use G4EmStandardPhysics_option4 as EM constructor and QGSP_BIC_HP with G4EmStandardPhysics_option4, for hadrontherapy applications. The Geant4-DNA physics lists report differences in modeling electron interactions in water, however, the tests have a pure regression testing purpose so no recommendation can be formulated.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Cascade of phase transitions in a planar Dirac material

Here we investigate a model of interacting Dirac fermions in 2 + 1 dimensions with M flavors and N colors having the U(M)×SU(N ) symmetry. In the large-N limit, we find that the U(M) symmetry is spontaneously broken in a variety of ways. In the vacuum, when the parity-breaking flavor-singlet mass is varied, the ground state undergoes a sequence of M first-order phase transitions, experiencing M + 1 phases characterized by symmetry breaking U(M)→U(M – k)×U(k) with k ϵ {0, 1, 2, · · · , M}, bearing a close resemblance to the vacuum structure of three-dimensional QCD. At finite temperature and chemical potential, a rich phase diagram with first and second-order phase transitions and tricritical points is observed. Also exotic phases with spontaneous symmetry breaking of the form as U(3)→U(1) 3 , U(4)→U(2)×U(1) 2 , and U(5)→U(2)2×U(1) exist. For a large flavor-singlet mass, the increase of the chemical potential μ brings about M consecutive first-order transitions that separate the low-μ phase diagram with vanishing fermion density from the high-μ region with a high fermion density.

1/N expansion↗
DOE OSTI2025

Community Dynamics Drive Calcium Carbonate Production in an Enriched Consortium of Soil Microbes

Recently, there has been a focus on using soil microbes as a means to store carbon in the soil in the form of calcium carbonate, outcomes of which include soil stabilization and biocementation. The molecular processes involved in microbially induced calcium carbonate formation are known, but there is still a significant knowledge gap regarding how community interactions, emergent processes that are distinct from the roles of individual members, may drive the formation of carbonate. To answer these questions, we describe the development and application of a consortium of soil microbes consisting of one species each of the Rhodococcus , Microbacterium , and Curtobacterium genera and two species from the Bacillus genus. We term these five species cultivated together carbon storing consortium A (CSC-A). Growth assays show that only a subset of CSC-A members produces CaCO 3 with Rhodococcus producing the most CaCO 3 but the complete CSC-A produces significantly higher amounts of CaCO 3 compared to the sum total carbonate produced by all member species. The development of CSC-A shows that CaCO 3 production may be as much a community process as it is the contribution of individual species, requiring us to move beyond single species analysis to fully understand carbonate formation by microbial communities in nature. CSC-A will allow the scientific community to ask and answer key questions about the molecular interactions surrounding inorganic carbon formation in soil, an important knowledge gap that must be filled if we wish to stabilize soils and harness microbial processes for materials production.

calcium carbonate↗
DOE OSTI2023

High salt-induced PSI-supercomplex is associated with high CEF and attenuation of state transitions

Abstract While PSI-driven cyclic electron flow (CEF) and assembly of thylakoid supercomplexes have been described in model organisms like Chlamydomonas reinhardtii , open questions remain regarding their contributions to survival under long-term stress. The Antarctic halophyte, C. priscuii UWO241 (UWO241), possesses constitutive high CEF rates and a stable PSI-supercomplex as a consequence of adaptation to permanent low temperatures and high salinity. To understand whether CEF represents a broader acclimation strategy to short- and long-term stress, we compared high salt acclimation between the halotolerant UWO241, the salt-sensitive model, C. reinhardtii , and a moderately halotolerant Antarctic green alga, C. sp. ICE-MDV (ICE-MDV). CEF was activated under high salt and associated with increased non-photochemical quenching in all three Chlamydomonas species. Furthermore, high salt-acclimated cells of either strain formed a PSI-supercomplex, while state transition capacity was attenuated. How the CEF-associated PSI-supercomplex interferes with state transition response is not yet known. We present a model for interaction between PSI-supercomplex formation, state transitions, and the important role of CEF for survival during long-term exposure to high salt.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Discovery of multi-functional polyimides through high-throughput screening using explainable machine learning

Polyimides have been widely used in modern industries because of their excellent mechanical and thermal properties, e.g., high-temperature fuel cells, displays, and aerospace composites. However, it usually takes decades of experimental efforts to develop a successful product. Aiming to expedite the discovery of high-performance polyimides, we utilize computational methods of machine learning (ML) and molecular dynamics (MD) simulations. Our study provides compelling evidence for the effectiveness of a data-driven approach in discovering novel polyimides. We first build a comprehensive library of more than 8 million hypothetical polyimides based on the polycondensation of existing dianhydride and diamine/diisocyanate molecules. Then we establish multiple ML models for the thermal and mechanical properties of polyimides based on their experimentally reported values, including glass transition temperature, Young’s modulus, and tensile yield strength. The obtained ML models demonstrate excellent predictive performance in identifying the key chemical substructures influencing the thermal and mechanical properties of polyimides. The use of explainable machine learning describes the effect of chemical substructures on individual properties, from which human experts can understand the cause of the ML model decision. Applying the well-trained ML models, we obtain property predictions of the 8 million hypothetical polyimides. Then, we screen the whole hypothetical dataset and identify three (3) best-performing novel polyimides that have better-combined properties than existing ones through Pareto frontier analysis. For an easy query of the discovered high-performing polyimides, we also create an online platform https://polyimide-explorer.herokuapp.com/ that embeds the developed ML model with interactive visualization. Furthermore, we validate the ML predictions through all-atom MD simulations and examine their synthesizability. The MD simulations are in good agreement with the ML predictions and the three novel polyimides are predicted to be easy to synthesize via Schuffenhauer’s synthetic accessibility score. Following the proposed ML guidance, we successfully synthesized a novel polyimide and the experimentally obtained high glass transition/thermal decomposition temperature demonstrated its excellent thermal stability. Here our study demonstrates an efficient way to expedite the discovery of novel polymers using ML prediction and MD validation. The high-throughput screening of a large computational dataset can serve as a general approach for new material discovery in other polymeric material exploration problems, such as organic photovoltaics, polymer membranes, and dielectrics.

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

Damage progression and failure of SiC/SiC composite tubes under hard-contact radial expansion

The response of silicon carbide (SiC) fiber-reinforced SiC matrix (SiC/SiC) composite cladding to mechanical interaction with fissile fuel is a knowledge gap that must be overcome to design and assess SiC-based cladding systems for advanced nuclear applications. This study developed the relevant mechanical testing capability and identified the failure behavior and the critical microstructural features and processing defects. Sections of SiC composite tube were subjected to a modified expansion-due-to-compression (EDC) test in an X-ray computed tomography microscope: a polyurethane plug pressed surrogate Al 2 O 3 into the inner walls of the SiC/SiC composite tubes to achieve hard contact. A pure EDC test with just a polyurethane plug was also performed as a reference. Through the use of displacement fields, digital volume correlation revealed inhomogeneous deformation fields in the tubes, even for pure EDC, which was related to the inherent defects in the structure. Deep learning–aided segmentation and systematic data analysis revealed that the presence of inhomogeneous deformation applied by the hard contact was exaggerated by the presence of inner surface imperfections left behind from the matrix densification process. In conclusion, the findings provide insights into the applications, highlighting the necessity for improvements in inner surface roughness and the incorporation of localized contacts in pellet–cladding mechanical interaction computational models.

Composites↗
DOE OSTI2020

A fully implicit, scalable, conservative nonlinear relativistic Fokker–Planck 0D-2P solver for runaway electrons

Upon application of a sufficiently strong electric field, electrons break away from thermal equilibrium and approach relativistic speeds. These highly energetic ‘runaway’ electrons (~ MeV) play a significant role in tokamak disruption physics, and therefore their accurate understanding is essential to develop reliable mitigation strategies. As such, we have developed a fully implicit solver for the 0D-2P (i.e., including two momenta coordinates) relativistic nonlinear Fokker–Planck equation (rFP). As in earlier implicit rFP studies (NORSE, CQL3D), electron–ion interactions are modeled using the Lorentz operator, and synchrotron damping using the Abraham–Lorentz–Dirac reaction term. However, our implementation improves on these earlier studies by (1) ensuring exact conservation properties for electron collisions, (2) strictly preserving positivity, and (3) being scalable algorithmically and in parallel. Key to our proposed approach is an efficient multigrid preconditioner for the linearized rFP equation, a multigrid elliptic solver for the Braams–Karney potentials, and a novel adaptive technique to determine the associated boundary values. We verify the accuracy and efficiency of the proposed scheme with numerical results ranging from small electric-field electrical conductivity measurements to the accurate reproduction of runaway tail dynamics when strong electric fields are applied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

LeptonInjector and LeptonWeighter: A neutrino event generator and weighter for neutrino observatories

In this work, we present a high-energy neutrino event generator, called LeptonInjector, alongside an event weighter, called LeptonWeighter. Both are designed for large-volume Cherenkov neutrino telescopes such as IceCube. The neutrino event generator allows for quick and flexible simulation of neutrino events within and around the detector volume, and implements the leading Standard Model neutrino interaction processes relevant for neutrino observatories: neutrino-nucleon deep-inelastic scattering and neutrino-electron annihilation. In this paper, we discuss the event generation algorithm, the weighting algorithm, and the main functions of the publicly available code, with examples. Program summary: Program Titles: LeptonInjector and LeptonWeighter CPC Library link to program files: https://doi.org/10.17632/662gkpjfd9.1 Developer's repository links: https://github.com/icecube/LeptonInjector and https://github.com/icecube/LeptonWeighter Licensing provisions: GNU Lesser General Public License, version 3. Programming Language: C++11 External Routines: • Boost • HDF5 • nuflux (https://github.com/icecube/nuflux) • nuSQuIDS (https://github.com/arguelles/nuSQuIDS) • Photospline (https://github.com/icecube/photospline) • SuiteSparse (https://github.com/DrTimothyAldenDavis/SuiteSparse) Nature of problem: LeptonInjector: Generate neutrino interaction events of all possible topologies and energies throughout and around a detector volume. LeptonWeighter: Reweight Monte Carlo events, generated by a set of LeptonInjector Generators, to any desired physical neutrino flux or cross section. Solution method: LeptonInjector: Projected ranges of generated leptons and the extent of the detector, in terms of column depth, are used to inject events in and around the detector volume. Event kinematics follow distributions provided in cross section files. LeptonWeighter: Event generation probabilities are calculated for each Generator, which are then combined into a generation weight and used to calculate an overall event weight.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Dynamic probabilistic risk assessment for electric grid cybersecurity

Electric grid cybersecurity risk has become a significant concern of industries and governments. This paper proposes a dynamic probabilistic risk assessment method for electric grid cybersecurity risk analysis. The proposed method helps reduce the reliance on expert judgment, capture a broad range of components and system dynamics, and model the interactions between various contributing entities (e.g., attacker, operator). In addition, the scenarios with multiple events, such as the occurrence of both cyberattacks and failures of physical components, the occurrence of both cyberattacks and operators’ (in)correct reactions, are considered and analyzed. Further, for each cyberattack scenario, Monte Carlo simulations are used to obtain possible sequences of the system's evolution under study and then derive risk estimates. As an application of the proposed method, the risk assessment method serves as the basis of risk-informed defense resource allocation to improve electric grid cybersecurity. The proposed method is verified using the IEEE 14-bus system by evaluating different security resource allocations for selected cyberattack scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Vessel Rupture Thresholds for Vessel–Bubble Interactions Using an Earthworm Vasculature Model

Intravenous microbubble oscillation in the presence of ultrasound has the potential to yield a wide range of therapeutic benefits. However, the likelihood of vessel damage caused by mechanical effects has not been quantified as a function of the numerous important parameters in therapeutic ultrasound procedures. In this study, we examined the effects of microbubbles injected into the vasculature of the earthworm. It was found that the elastic properties of earthworm blood vessels are similar to those of arteries in older humans, and that earthworms are well suited to the large number of experiments necessary to investigate safety of procedures involving microbubble oscillation in sonicated vessels. Microbubbles were infused into earthworm vessels, and the rupture time during sonication was recorded as a function of ultrasound frequency, pulse repetition frequency and acoustic pressure. A modified mechanical index (MMI) was defined that successfully captured the trends in rupture probability and rupture time for the different parameter values, creating a database of vessel rupture thresholds. In the absence of bubbles, the product of MMI squared and rupture time was approximately constant, indicating a possible radiation-force effect. Here, the MMI was an effective correlating parameter in the presence of bubbles, though the mathematical dependence is not yet apparent. The results of the study are expected to be valuable in designing more refined studies in vertebrate models, as well as informing computational models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Multi-mode quasi-static excitation for systems with nonlinear joints

Finite element models can be used to model and predict the hysteresis and energy dissipation exhibited by nonlinear joints in structures. As a result of the nonlinearity, the frequency and damping of a mode is dependent on excitation amplitude, and when the modes remain uncoupled, quasi-static modal analysis has been shown to efficiently predict this behavior. However, in some cases the modes have been observed to couple such that the frequency and damping of one mode is dependent on the amplitude of other modes. To model the interactions between modes, one must integrate the dynamic equations in time, which is several orders of magnitude more expensive than quasi-static analysis. This work explores an alternative where quasi-static forces are applied in the shapes of two or more modes of vibration simultaneously, and the resulting load–displacement curves are used to deduce the effect of other modes on the effective frequency and damping of the mode in question. This methodology is demonstrated on a simple 2D cantilever beam structure with a single bolted joint which exhibits micro-slip nonlinearity over a range of vibration amplitudes. The predicted frequency and damping are compared with those extracted from a few expensive dynamic simulations of the structure, showing that the quasi-static approach produces reasonable albeit highly conservative bounds on the observed dynamics. This framework is also demonstrated on a 3D structure where dynamic simulations are infeasible.

42 ENGINEERING↗
DOE OSTI2022

Shape Control of Colloidal Semiconductor Nanocrystals through Thermodynamically Driven Aggregative Growth

The optoelectronic properties of colloidal semiconductor nanocrystals (NCs) can be manipulated by changing their geometric shapes. The precise synthetic control over particle morphologies, however, has remained elusive. Conventional growth techniques rely on the kinetic assembly of atomic units, where supersaturation and precipitation processes can lead to a broad distribution of particle shapes. In this paper, we demonstrate that replacing atomic precursors with small-size nanocrystals as building blocks for larger colloids offers an easier, more predictive control over nanoparticle shape evolution. The reported growth strategy is illustrated via shape-selective syntheses of CdSe and CdS NC cubes, spheres, rods, as well as unprecedented “donut” and ring-like structures. Growth of desired particle morphologies was achieved by choosing an appropriate reaction solvent that minimizes the surface free energy through a unique shape formation. This so-called "aggregative growth mechanism" is explained using a thermodynamic model for interacting viscous colloids. We expect that this novel growth strategy will offer a viable platform for shape–selective synthesis of many inorganic colloids, providing pathways to optimizable optoelectronic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Toward Engineering of Solution Microenvironments for the CO 2 Reduction Reaction: Unraveling pH and Voltage Effects from a Combined Density-Functional–Continuum Theory

Engineering the electrolyte microenvironment represents an attractive route to tuning the selectivity of electrocatalytic reactions beyond catalyst composition and morphology. However, harnessing the full potential of this approach requires understanding the interplay between voltage, electrolyte composition, and adsorbate binding within the electrical double layer, which is absent from the usual theoretical approaches. In this study, we apply a recently developed density functional theory (DFT)–continuum approach based on the effective screening medium method and reference interaction site model (ESM–RISM) to explore electrolyte effects with an enhanced description of the electrochemical interface. Furthermore, applying this method to the binding of CO adsorbates in potassium-containing electrolytes on copper, a problem of direct relevance to CO 2 electroreduction to value-added products, we show that the interdependence of voltage and pH leads to an unexpected change in adsorption site preference on Cu(001) terraces. Our findings highlight the often-overlooked importance of the electrical double-layer structure for predicting catalyst operation.

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

Mesoscale Modeling of Hydrogels Under Frictional Shear Stress

Hydrogels are three-dimensional networks of hydrophilic polymers often used as a simplified model of hydrated biological materials, from cartilaginous joints to the ocular tear film. However, the lubrication mechanisms of hydrogels remain poorly understood, partly due to their complex polymeric structure, which creates blurred interfaces during sliding that are challenging to study experimentally. In this study, we employ dissipative particle dynamics (DPD) to investigate the frictional behavior of a polymeric hydrogel network sliding against a solid wall in an explicit viscous solvent. This computational approach enables us to model hydrodynamic interactions and mesoscale polymer dynamics, capturing key aspects of hydrogel friction. Our simulations reveal that hydrogel friction is governed by the interplay between polymer relaxation and viscous shear, characterized by the Weissenberg number (Wi). At low Wi, friction coefficient remain nearly constant, dominated by polymer relaxation. However, at higher Wi, friction is dominated by viscous drag within a near-wall solvent layer, leading to a linear increase in friction coefficient with Wi. Furthermore, our results demonstrate an inverse relationship between the friction coefficient and the applied normal load, consistent with experimental observations. This work provides new insights into the fundamental tribological properties of hydrogels, shedding light on the micromechanics of hydrogel friction. Improving our understanding of hydrogel structure and dynamics under friction advances our knowledge of the mechanisms regulating biological lubrication in health and disease.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Segmental Dynamics and Vitrification in Associating Copolymer Melts: Role of Cluster Formation, Microdomains, and Cross-Linking

A microscopic statistical mechanical theory of the structure, self-assembly, and activated segmental relaxation is employed to study associating copolymer melts with high attractive sticker fractions, local clustering, and disordered microphase ordering. The stickers are dynamically pinned in a manner that does not affect equilibrium structure which mimics the much slower physical bond breaking process or postassembly cross-linking of sticky monomers. Local sticker clustering and microdomain spatial correlations significantly modify the activated relaxation of nonstickers and glass transition temperature, T g . A re-entrant glass-melting feature is predicted as sticker attraction strength is initially increased corresponding to a speed up of segmental relaxation, and hence reduction of T g relative to the cross-linked homopolymer network. A mechanistic analysis reveals three competing effects: a purely kinetic slowing down of nonstickers down due to cross-linking, disordering of the nonsticker local cage and weakening of effective forces they experience due to sticker physical clustering, and a longer range impact of microdomain scale correlations that results in nonmonotonic dynamical effects. At high enough attraction strength, a qualitative change emerges corresponding to a sticker fraction dependent elevation of T g , which eventually surpasses that of the cross-linked homopolymer network. Dynamically, the new physics arises from a complex evolution of the amplitude of the collective elastic field that dresses the large amplitude mobile segment hopping within a coupled local-nonlocal description of the alpha relaxation. Furthermore, the results are qualitatively consistent with recent experiments on associating PDMS and PPG telechelics of fixed sticker fraction but with chemically different end groups of variable attraction strengths. Finally, possible tests using simulation and the influence of material or model specific interaction potentials and other real world complications are discussed.

Copolymers↗