Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional 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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

EnergyPlus Model Context Protocol Server (EnergyPlus-MCP) v0.1

EnergyPlus-MCP is the first open-source Model Context Protocol server specifically designed for EnergyPlus building energy simulation. This innovative software enables AI assistants and other applications to interact programmatically with EnergyPlus through a standardized, secure interface, eliminating traditional technical barriers in building energy modeling. The software provides specialized tools across five functional domains: server management, model configuration and loading, comprehensive building component inspection, systematic model modification, and simulation execution with results visualization. Key features include automated HVAC system discovery and topology mapping, advanced schedule analysis, intelligent model validation, and interactive visualization capabilities. EnergyPlus-MCP's layered architecture ensures robust separation between protocol communication and domain expertise, enabling scalable deployment across organizations, educational institutions, and research teams. Unlike direct LLM approaches that suffer from inconsistent results and security gaps, EnergyPlus-MCP provides validated, reliable interactions while maintaining scientific rigor. This democratizes sophisticated building energy analysis, making EnergyPlus accessible to broader audiences through conversational interfaces and streamlined workflows.

Li, Han [Lawrence Berkeley National Laboratory (LB↗
DOE OSTI2020

Comparing event generator predictions and ab-initio calculations of $\nu$-$^{12}$C neutral current quasi-elastic scattering at 1 GeV

The measurement of neutrino oscillations and exotic physics searches are important parts of the physics program in the near future, with new state-of-the-art experiments planned within the next decade. Future and modern experiments in these fields will make use of nuclear targets. Event Generators (EGs) are software used in the analysis of neutrino oscillation experiments. EGs use to predict kinematic observables for a range of neutrino energies. These simulations may lack physics captured by more rigorous theoretical calculations. This work compares EG performance to nuclear theory calculations by comparing observables generated in the two frameworks. We provide a common set of definitions between theory and experiment and assess the physics contained in EG simulations. Neutral current quasi-elastic (NCQE) scattering events for neutrinos and anti-neutrinos on a $^{12}$C target are simulated with a specific EG, NEUT, used by the T2K experiment for its analysis. The simulated cross sections are compared to analytic calculations from nuclear theory within the factorization scheme. We compare the NEUT implementation of two different models on nuclear spectral functions: the Relativistic Fermi Gas (RFG) and the correlated basis spectral function (CBF) to analytic calculations of the same models in the factorization scheme. For both nuclear physics models, we compare the appearance of features in the distributions relevant to experimental analyses. Qualitatively, the shape of the simulated distribution is similar to the one obtained through theory calculations; however, there are some discrepancies between the theory calculations and the NEUT simulation. While the EG simulations and analytic calculations with the same model of nuclear dynamics show similar overall features, there are still differences between the two.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The effect of temperature on photosystem II efficiency across plant functional types and climate

Modeling terrestrial gross primary productivity (GPP) is central to predicting the global carbon cycle. Much interest has been focused on the environmentally induced dynamics of photosystem energy partitioning and how improvements in the description of such dynamics assist the prediction of light reactions of photosynthesis and therefore GPP. The maximum quantum yield of photosystem II (Φ PSIImax ) is a key parameter of the light reactions that influence the electron transport rate needed for supporting the biochemical reactions of photosynthesis. Φ PSIImax is generally treated as a constant in biochemical photosynthetic models even though a constant Φ PSIImax is expected only for non-stressed plants. We synthesized reported Φ PSIImax values from pulse-amplitude-modulated fluorometry measurements in response to variable temperatures across the globe. We found that Φ PSIImax is strongly affected by prevailing temperature regimes with declined values in both hot and cold conditions. To understand the spatiotemporal variability in Φ PSIImax , we analyzed the temperature effect on Φ PSIImax across plant functional type (PFT) and habitat climatology. The analysis showed that temperature's impact on Φ PSIImax is shaped more by climate than by PFT for plants with broad latitudinal distributions or in regions with extreme temperature variability. There is a trade-off between the temperature range within which Φ PSIImax remains maximal and the overall rate of decline of Φ PSIImax outside the temperature range such that species cannot be simultaneously tolerant and resilient to extreme temperatures. Our study points to a quantitative approach for improving electron transport and photosynthetic productivity modeling under changing climates at regional and global scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Distinguishing Electron and Hole Dynamics in Functionalized CdSe/CdS Core/Shell Quantum Dots Using Complementary Ultrafast Spectroscopies and Kinetic Modeling

The evolution of excitation energy and photogenerated charges in semiconductor quantum dots (QDs) functionalized with molecular acceptors can be probed on ultrafast time scales using techniques such as transient absorption (TA) spectroscopy. However, historical interpretations that the 1S(e)-1S 3/2 (h) transition in Cd-chalcogenide QDs is fully attributable to electrons may be misleading, and multiexponential models used to fit TA kinetics do not correspond directly to specific photophysical processes. Here, we present visible-wavelength and mid-IR TA and time-resolved photoluminescence measurements to inform a comprehensive kinetic model of the photoexcited CdSe/CdS core/shell QDs functionalized with passivating oleic acid (OA), hole-accepting ferrocene, or electron-accepting naphthalene bisimide (NBI). We show that ~30% of the 1S signal and 72% of the IR signal can originate from holes in well-passivated core/shell QDs. We also demonstrate evidence of electron trapping in OA-capped core/shell QDs, with additional electron transfer and hole trapping in the QDs functionalized with NBI. Electron (hole) trapping and detrapping occur in 450 ± 100 ps (430 ± 70 ps) and 340 ± 100 ps (1.1 ± 0.4 ns) respectively, while the time constant for electron transfer to NBI is ~1.8 ns. The comprehensive picture of photophysical processes provided by the complementary ultrafast techniques and kinetic modeling can accelerate both the fundamental science and application development of nanostructured and molecular systems.

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

Non-smooth Bayesian optimization in tuning scientific applications

Tuning algorithmic parameters to optimize the performance of large, complicated computational codes is an important problem involving finding the optima and identifying regimes defined by non-smooth boundaries in black-box functions. Within the Bayesian optimization framework, the Gaussian process surrogate model produces smooth mean functions, but functions in the tuning problem are often non-smooth, which is exacerbated by the fact that we usually have limited sequential samples from the black-box function. Here, motivated by these issues encountered in tuning, we propose a novel Gaussian process model called a clustered Gaussian process (cGP), where the components are dynamically updated by clustering. In our studies, the performance of cGP can be better than stationary GPs in nearly 90% of the experiments and better than non-stationary GPs in nearly 70% of the repeated experiments while requiring less computational cost. cGP provides a novel approach for dynamic GP, computes more efficiently than recursive partitioning, and discovers non-smoothness regimes. We provide extensive experiments including high-performance computing (HPC) and industrial simulation functions to show the effectiveness of our methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Crowdsourcing biocuration: The Community Assessment of Community Annotation with Ontologies (CACAO)

Experimental data about gene functions curated from the primary literature have enormous value for research scientists in understanding biology. Using the Gene Ontology (GO), manual curation by experts has provided an important resource for studying gene function, especially within model organisms. Unprecedented expansion of the scientific literature and validation of the predicted proteins have increased both data value and the challenges of keeping pace. Capturing literature-based functional annotations is limited by the ability of biocurators to handle the massive and rapidly growing scientific literature. Within the community-oriented wiki framework for GO annotation called the Gene Ontology Normal Usage Tracking System (GONUTS), we describe an approach to expand biocuration through crowdsourcing with undergraduates. This multiplies the number of high-quality annotations in international databases, enriches our coverage of the literature on normal gene function, and pushes the field in new directions. From an intercollegiate competition judged by experienced biocurators, Community Assessment of Community Annotation with Ontologies (CACAO), we have contributed nearly 5,000 literature-based annotations. Many of those annotations are to organisms not currently well-represented within GO. Over a 10-year history, our community contributors have spurred changes to the ontology not traditionally covered by professional biocurators. The CACAO principle of relying on community members to participate in and shape the future of biocuration in GO is a powerful and scalable model used to promote the scientific enterprise. It also provides undergraduate students with a unique and enriching introduction to critical reading of primary literature and acquisition of marketable skills.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Find the Gap: Black Hole Population Analysis with an Astrophysically Motivated Mass Function

We introduce a novel black hole mass function that realistically models the physics of pair-instability supernovae with a minimal number of parameters. Applying this to all events in the LIGO-Virgo Gravitational-Wave Transient Catalog 2 (GWTC-2), we detect a peak at M BHMG = 46$^{+17}_{-6}$ M ⊙ . Repeating the analysis without the black holes from the event GW190521, we find this feature at M BHMG =54± 6 M ⊙ . These results establish the edge of the anticipated "black hole mass gap" at a value compatible with the expectation from standard stellar structure theory. The mass gap manifests itself as a discontinuity in the mass function and is populated by a distinct, less-abundant population of higher-mass black holes. We find that the primary black hole population scales with power-law index -1.95 ± 0.51 (-1.97 ± 0.44) with (without) GW190521, consistent with models of star formation. Using Bayesian techniques, we establish that our mass function fits a new catalog of black hole masses approximately as well as pre-existing phenomenological mass functions. We also remark on the implications of these results for constraining or discovering new phenomena in nuclear and particle physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

In vivo visualization of butterfly scale cell morphogenesis in Vanessa cardui

Many organisms exhibit functional micro- and nanoscale materials with structural definition and performance that challenge synthetic fabrication techniques, yet we know little about the processes that enable their formation. Using butterfly scales as a model system for functional biomaterials, we establish a timeline of scale formation and quantify relevant structural parameters for developing painted lady butterflies. We overcome challenges of previous efforts by imaging structure formation directly in living organisms, which allows us to continuously observe the evolving wing tissue and the fine details of individual scale cells. Visualization of scale structure formation in live butterflies forms the basis for modeling the underlying biomechanical processes and opens avenues for their translation into advanced fabrication strategies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Precise Modeling of a Complex Solenoidal Magnetic Field Using a Combination of Analytic Functions and a PINN

We demonstrate an iterative approach to modeling a sparsely measured magnetic field in a large-bore solenoid. This approach uses a hybrid of traditional and machine learning techniques. The traditional technique is a linear least-squares fit using a series solution to Laplace's equation, while the machine learning technique involves the training of a physics-informed neural network (PINN) on the least-squares fit residuals. We use a newly defined activation function "DELTAsnake," a modification to the snake activation function proposed by Ziyin et al. that allows for stronger curvature and non-monotonicity. The combined model approximately obeys Maxwell's equations to a level sufficient for producing high quality physics simulations and analysis. Our approach is applied to a highly realistic calculation of the expected magnetic field in the Mu2e experiment's Detector Solenoid which includes a simple model for the expected statistical measurement uncertainties. Using ten toy measurement simulations, we demonstrate the capabilities of our model in comparison to the least-squares method alone; the least-squares method alone results in a reduced chi-squared statistic of ${2.15 \pm 0.01}$, while our approach improves the reduced chi-square to ${1.034 \pm 0.005}$. Furthermore, for an average toy simulation, we show that the range of the RMS of the three field component residuals reduces from ${0.07-0.37}$ Gauss to ${0.05-0.07}$ Gauss. We find that this novel method is robust against a realistic systematic uncertainty deriving from Hall probe calibration bias and can be used to significantly reduce the number of measurements required to achieve an accurate model.

Kampa, Cole [Caltech] (ORCID:0000000192972920)↗
DOE OSTI2020

Updates to the Regional Seismic Travel Time (RSTT) Model: 1. Tomography

Abstract A function of global monitoring of nuclear explosions is the development of Earth models for predicting seismic travel times for more accurate calculation of event locations. Most monitoring agencies rely on fast, distance-dependent one-dimensional (1D) Earth models to calculate seismic event locations quickly and in near real-time. RSTT (Regional Seismic Travel Time) is a seismic velocity model and computer software package that captures the major effects of three-dimensional crust and upper mantle structure on regional seismic travel times, while still allowing for fast prediction speed (milliseconds). We describe updates to the RSTT model using a refined data set of regional phases (i.e., Pn, Pg, Sn, Lg) using the Bayesloc relative relocation algorithm. The tomographic inversion shown here acts to refine the previous RSTT public model ( rstt201404um ) and displays significant features related to areas of global tectonic complexity as well as further reduction in arrival residual values. Validation of the updated RSTT model demonstrates significant reduction in median epicenter mislocation (15.3 km) using all regional phases compared to the iasp91 1D model (22.1 km) as well as to the current station correction approach used at the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (18.9 km).

58 GEOSCIENCES↗
DOE OSTI2025

sparse_bias

This is a python package used to fit an unknown function from data that potential contains systematic biases related to metadata. The model fits the unknown function and uses a Bayesian horseshoe prior model to impose sparsity on the bias terms. This code has been generalized from research code developed for AIACHNE into a package that should have more general application in a wider class of statistical models.

Walton, Noah↗
DOE OSTI2022

Calibration, modeling, parameterization, and verification of the instrument transfer function of an interferometric microscope

Interferometric microscopes are used to measure surface roughness, from which the computed power spectral density function can be used to extract bandwidth-limited values of the various surface properties, such as root-mean-square (rms) height and slope errors. Measurements with a microscope equipped with different objectives that have an overlapping spatial frequency range usually give different rms results over the common frequency bandwidth. This is a result of different instrument transfer functions (ITFs) that attenuate spatial frequencies by different amounts over the overlapping range. We report on the use of binary pseudo-random array (BPRA) standards to characterize the ITF of an interferometric microscope with the various objectives. We use a simple model of a 1D binary Results show that the spectrum for an undersampled array is a cosine function, rather than a straight line, constant white noise spectrum. We have an analytical model of the ITF that includes the effects of an obscured aperture and defocus, in addition to the usual parameters of numerical aperture, wavelength, and sampling period. In addition, the model includes the effect of aliasing of spatial frequency components beyond the Nyquist back into the sub-Nyquist region. We compare the model PSD predictions to the measurements performed with different objectives. Departures of the measured PSDs from the model predictions indicate that there are higher-order ITF corrections yet to be identified and included in the model.

Takacs, PZ↗
DOE OSTI2021

Bayesian model averaging for analysis of lattice field theory results

Statistical modeling is a key component in the extraction of physical results from lattice field theory calculations. Although the general models used are often strongly motivated by physics, many model variations can frequently be considered for the same lattice data. Model averaging, which amounts to a probability-weighted average over all model variations, can incorporate systematic errors associated with model choice without being overly conservative. We discuss the framework of model averaging from the perspective of Bayesian statistics, and give useful formulae and approximations for the particular case of least-squares fitting, commonly used in modeling lattice results. In addition, we frame the common problem of data subset selection (e.g. choice of minimum and maximum time separation for fitting a two-point correlation function) as a model selection problem and study model averaging as a straightforward alternative to manual selection of fit ranges. Numerical examples involving both mock and real lattice data are given.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Material Effects on Electron-Capture Decay in Cryogenic Sensors

Several current searches for physics beyond the standard model are based on measuring the electron-capture (EC) decay of radionuclides implanted into cryogenic high-resolution sensors. The sensitivity of these experiments has already reached the level where systematic effects related to atomic state energy changes from the host material are a limiting factor. One example is a neutrino mass study based on the nuclear EC decay of 7 Be to 7 Li inside cryogenic Ta-based sensors. Here, to understand the material effects at the required level, we use density-functional theory to model the electronic structure of lithium atoms in different atomic environments of the polycrystalline Ta absorber film. The calculations reveal that the Li 1s binding energies can vary by more than 2 eV due to insertion at different lattice sites, at grain boundaries, in disordered Ta, and in the vicinity of various impurities. However, the total range of Li 1s shifts does not exceed 4 eV, even for extreme amorphous disorder. Furthermore, when investigating the effects on the Li 2s levels, we find broadening of more than 5 eV due to hybridization with the Ta band structure. Material effects are shown to contribute significantly to peak broadening in Ta-based sensors that are used to search for physics beyond the standard model in the EC decay of 7 Be, but they do not explain the full extent of observed broadening. Understanding these in-medium effects will be required for current- and future-generation experiments that observe low-energy radiation from the EC decay of implanted isotopes to evaluate potential limitations on the measurement sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural regulation-induced Li-electron disentanglement for stabilized oxygen redox of Li-excess disordered rock-salt cathode materials

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.

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

Pressure Drop in Seven-Pin Wire-Wrapped Rod Bundle for the Sodium Cartridge Loop in Versatile Test Reactor

This work studies the hydrodynamics of the seven pin wire-wrapped rod bundle in the sodium cartridge loop for the Versatile Test Reactor (VTR) through scaled water experiments and computational fluid dynamics (CFD) simulations. Here, the scaling analysis is first performed to demonstrate the hydrodynamic similarity between water and sodium flows at the same Reynolds number Re. A separate-effects test facility is designed and constructed based on the scaling analysis. Detailed experimental data on the pressure drop covering a wide range of Re values (1165 to 27 689) are obtained, which are used to evaluate existing correlations for friction factor and to benchmark CFD simulations. The experimentally determined friction factors f Exp agree well with the Upgraded Cheng and Todreas Detailed Correlation and Pacio-Chen-Todreas Detailed Model within ±7% but are significantly underpredicted by Rehme's correlation by 25%. Various CFD near-wall treatment methods are tested using ANSYS Fluent and evaluated by experimental data. It is found that when the recommended wall y + values are met, most of the near-wall treatment methods can give accurate friction factor predictions. The resolved near-wall method (y + ~ 1) with the Shear Stress Transport k - ω turbulence model and the scalable wall functions (y + > 12) with the realizable k - ε turbulence model can predict f Exp within ±5% The standard wall functions (y + > 30) and nonequilibrium wall functions (y + > 30) with the realizable k - ε model can predict f Exp within ±10%.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Search for pair-produced scalar and vector leptoquarks decaying into third-generation quarks and first- or second-generation leptons in $pp$ collisions with the ATLAS detector

A search for pair-produced scalar and vector leptoquarks decaying into quarks and leptons of different generations is presented. It uses the full LHC Run 2 (2015–2018) data set of 139 fb -1 collected with the ATLAS detector in proton–proton collisions at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV. Scalar leptoquarks with charge -(1/3)$e$ as well as scalar and vector leptoquarks with charge +(2/3)$e$ are considered. All possible decays of the pair-produced leptoquarks into quarks of the third generation ($t, b$) and charged or neutral leptons of the first or second generation ($e, μ, ν$) with exactly one electron or muon in the final state are investigated. No significant deviations from the Standard Model expectation are observed. Upper limits on the production cross-section are provided for eight models as a function of the leptoquark mass and the branching ratio of the leptoquark into the charged or neutral lepton. In addition, lower limits on the leptoquark masses are derived for all models across a range of branching ratios. Two of these models have the goal of providing an explanation for the recent $B$-anomalies. In both models, a vector leptoquark decays into charged and neutral leptons of the second generation with a similar branching fraction. Lower limits of 1980 GeV and 1710 GeV are set on the leptoquark mass for these two models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗