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 469 records · Page 26

Can general circulation models (GCMs) represent cloud liquid water path adjustments to aerosol–cloud interactions?

Abstract. General circulation models (GCMs), unlike other lines of evidence, indicate that anthropogenic aerosols cause a global-mean increase in cloud liquid water path (ℒ) and thus a negative adjustment to radiative forcing of the climate by aerosol–cloud interactions. In part 1 of this series of papers, we showed that this is true even in models that reproduce the negative correlation observed in present-day internal variability in ℒ and cloud droplet number concentration (Nd). We studied several possible confounding mechanisms that could explain the noncausal cloud–aerosol correlations in GCMs and that possibly contaminate observational estimates of radiative adjustments. Here, we perform single-column and full-atmosphere GCM experiments to investigate the causal model-physics mechanisms underlying the model radiative adjustment estimate. We find that both aerosol–cloud interaction mechanisms thought to be operating in real clouds – precipitation suppression and entrainment evaporation enhancement – are active in GCMs and behave qualitatively in agreement with physical process understanding. However, the modeled entrainment enhancement has a negligible global-mean effect. This raises the question of whether the GCM estimate is incorrect due to parametric or base-state representation errors or whether the process understanding gleaned from a limited set of canonical cloud cases is insufficiently representative of the diversity of clouds in the real climate. Regardless, even at limited resolution, the GCM physics appears able to parameterize the small-scale microphysics–turbulence interplay responsible for the entrainment enhancement mechanism. We suggest ways to resolve tension between current and future (storm-resolving) global modeling systems and other lines of evidence in synthesis climate projections.

54 ENVIRONMENTAL SCIENCES↗

Incorporating an Interactive Fire Plume-Rise Model in the DOE's Energy Exascale Earth System Model Version 1 (E3SMv1) and Examining Aerosol Radiative Effect

The vertical distribution of biomass burning aerosol (BBA) is important in regulating their impacts on weather and climate. The plume-rise process affects the injection height of BBA and interacts with the air parcel lifting and cloud processes. However, these processes are not represented in most global climate models. In this study, we replaced the fixed vertical profiles of monthly BBA emissions in the Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) with an interactive fire plume-rise model. The vertical distribution of BBA emissions was calculated as a function of ambient thermodynamic conditions from the host E3SMv1, with distributions of fire sizes and sensible heat fluxes derived from the observations. The maximum fire radiative power (FRP) technique was used to determine the fire size. Scaling-FRP technique is used to calculate the wildfire heat release. Daily BBA emission, superimposed with a fire diurnal cycle retrieved from the satellite observation, was included in model simulations. The model shows improved agreement with satellite retrievals and in situ measurement during the National Oceanic and Atmospheric Administration Wildfire Experiment for Cloud chemistry, Aerosol absorption, and Nitrogen campaign. The model-observation comparison demonstrates the importance of the plume-rise model and fire diurnal cycle assumption in determining the BBA fields. We also find that E3SMv1 with new features produces a larger carbonaceous aerosol burden, leading to 0.13 W m –2 warming at the top of atmosphere compared to the default E3SMv1. This highlights the importance of accurately representing the BBA injection height and indicates a no-linear nature in the BBA-induced radiative effect.

54 ENVIRONMENTAL SCIENCES↗

Understanding the causes of satellite–model discrepancies in aerosol–cloud interactions using near-LES simulations of marine boundary layer clouds

Aerosol–cloud interactions (ACI) remain the largest source of uncertainty in model estimates of anthropogenic radiative forcing, primarily because of deficiencies in representing aerosol–cloud microphysical processes that lead to inconsistent cloud liquid water path (LWP) responses to aerosol perturbations between observations and models. To investigate this discrepancy, we conducted a series of large-eddy-scale simulations driven by realistic meteorology over the eastern North Atlantic, and evaluated LWP susceptibility, precipitation processes, and boundary layer thermodynamics using satellite and ground-based observations. Simulated LWP responses show a strong dependence on cloud state. Non-precipitating thin clouds exhibit a modest LWP decrease with increasing cloud droplet number concentration (N d ), consistent in sign but weaker in magnitude than satellite estimates, reflecting enhanced turbulent mixing and evaporation. The largest model-observation discrepancy occurs in non-precipitating thick clouds, where simulated LWP susceptibilities are strongly positive (+0.32) while observations indicate large negative values (−0.69). This discrepancy stems from excessive precipitation driven by underestimated entrainment, overly active accretion, and overly broad drop-size distributions in polluted conditions. While our high-resolution setup mitigates the excessive drizzling common in coarser models and captures key regime transitions, these biases persist – highlighting that improved parameterizations of cloud-top processes, precipitation, and aerosol effects are needed beyond simply increasing model resolution. Additionally, misrepresented moisture inversions in reanalysis introduce a moist bias in cloud-top relative humidity, further amplifying positive LWP susceptibility. Our results also suggest that large negative N d –LWP relationships in observations may reflect internal cloud processes rather than true ACI effects.

Aerosol-cloud interaction↗

Interaction-expansion inchworm Monte Carlo solver for lattice and impurity models

Multiorbital quantum impurity models with general interaction and hybridization terms appear in a wide range of applications, including embedding, quantum transport, and nanoscience. However, most quantum impurity solvers are restricted to a few impurity orbitals, discretized baths, diagonal hybridizations, or density-density interactions. We generalize the inchworm quantum Monte Carlo method to the interaction expansion, and we explore its application to typical single- and multiorbital problems encountered in investigations of impurity and lattice models. Our implementation generically outperforms bare and bold-line quantum Monte Carlo algorithms in the interaction expansion. For the systems studied here, our implementation remains inferior to the more specialized hybridization expansion and auxiliary field algorithms. So, the problem of convergence to unphysical fixed points, which hampers so-called bold-line methods, is not encountered in inchworm Monte Carlo.

36 MATERIALS SCIENCE↗

A Thermodynamic Model for Nd(III)–Sulfate Interaction at High Ionic Strengths and Elevated Temperatures: Applications to Rare Earth Element Extraction

Neodymium (Nd), a rare earth element (REE), is critical to numerous industries. Neodymium can be extracted from ore concentrates, waste materials, or recycled materials such as recycled Nd-Fe-B permanent magnets. In a standard process, concentrated sulfuric acid (H 2 SO 4 ) is used as an extraction/leaching agent. Therefore, knowledge of Nd(III)–sulfate interaction at high ionic strengths is important for optimization of the extraction process. In addition, sulfate is also a major species in natural surface waters and present in nuclear waste streams. Nd(III) has been used a chemical analog to trivalent actinides in nuclear waste research and development. Consequently, knowledge of Nd(III)-sulfate interactions is also impactful to the field of nuclear waste management. Here, in this study, we have developed a thermodynamic model that can describe the interaction of Nd(III) with sulfate to ionic strengths up to ~ 16.5 mol·kg –1 and to temperatures up to 100 °C. The model adopts the Pitzer formulation to describe activity coefficients of aqueous species. This model can be used to design and optimize a chemical process for REE recovery from ore concentrates, recycled materials, and acid mine drainage (AMD) and to understand the mobility of REEs and actinides in the environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phenomenological theory of the magnetic 90° helical state

We explore a phenomenological phase diagram for the magnetic helical state with 90° turn angle between neighboring spins in the external magnetic field. Such a state is formed by the Eu spin layers in the superconducting iron arsenide RbEuFe 4 As 4 . The peculiarity of this spin configuration is that it is not realized in the standard Heisenberg model with bilinear exchange interactions. A minimum model allowing for such a state requires the biquadratic nearest-neighbor interaction term. In addition, in tetragonal materials, the 90° helix state may be stabilized by the in-plane fourfold anisotropy term, which also fixes helix orientation with respect to the crystal lattice. Such a system has a very rich behavior in the external magnetic field. The magnetic field induces the metamagnetic transition to the double-periodic state with the moment angles (α α -α,-α) with respect to the field for the four subsequent spins. The transition field to this state from the deformed helix is determined by the strength of biquadratic interaction. The transition is second-order for small biquadratic coupling and becomes first-order when this coupling exceeds the critical value. On the other hand, the aligned state at high magnetic field becomes unstable with respect to the formation of an incommensurate fan state, which transforms into the double-periodic state with decreasing magnetic field. The range of this incommensurate state near the saturation field is proportional to square of the biquadratic coupling. In addition, when the magnetic field is applied along one of four the equilibrium moment directions, the deformed helix state experience the first-order rotation transition at the field determined by the fourfold anisotropy.

36 MATERIALS SCIENCE↗

Quantitative dissection of Agrobacterium T-DNA expression in single plant cells reveals density-dependent synergy and antagonism

Agrobacterium pathogenesis, which involves transferring T-DNA into plant cells, is the cornerstone of plant genetic engineering. As the applications that rely on Agrobacterium increase in sophistication, it becomes critical to achieve a quantitative and predictive understanding of T-DNA expression at the level of single plant cells. Here we examine if a classic Poisson model of interactions between pathogens and host cells holds true for Agrobacterium infecting Nicotiana benthamiana. Systematically challenging this model revealed antagonistic and synergistic density-dependent interactions between bacteria that do not require quorum sensing. Using various approaches, we studied the molecular basis of these interactions. To overcome the engineering constraints imposed by antagonism, we created a dual binary vector system termed ‘BiBi’, which can improve the efficiency of a reconstituted complex metabolic pathway in a predictive fashion. Our findings illustrate how combining theoretical models with quantitative experiments can reveal new principles of bacterial pathogenesis, impacting both fundamental and applied plant biology.

Alamos, Simon↗

Shift Happens: Building Robust AI Models with Domain Adaptation

Artificial Intelligence (AI) is revolutionizing physics research—from probing the large-scale structure of the Universe to modeling subatomic interactions and fundamental forces. Yet, a major challenge persists: AI models trained on simulations or old experiment / astronomical survey often perform poorly when applied to new data—exposing issues of dataset (domain) shift, model robustness, and uncertainty in predictions. This summer school session will introduce students to common challenges in applying AI across domains and present solutions based on domain adaptation—a set of techniques designed to improve model generalization under domain shift. We will cover foundational ideas, practical strategies, and current research frontiers in this area. Through examples in astrophysics, we'll explore how domain adaptation can help bridge the gap between synthetic and real-world data, improve trust in model outputs, and advance scientific discovery. The concepts discussed are broadly applicable across physics and other scientific disciplines, making this a valuable topic for anyone interested in building robust, transferable AI models for science.

Ciprijanovic, A. [Fermilab] (ORCID:000000031281719↗

AI in Astrophysics: Tackling Domain Shift, Model Robustness and Uncertainty

Artificial Intelligence (AI) is revolutionizing physics research from probing the large-scale structure of the Universe to modeling subatomic interactions and fundamental forces. Yet, a major challenge persists: AI models trained on simulations or old experiment / astronomical survey often perform poorly when applied to new data exposing issues of dataset (domain) shift, model robustness, and uncertainty in predictions. This talk will introduce common challenges in applying AI across domains and present solutions based on domain adaptation a set of techniques designed to improve model generalization under domain shift. We will cover foundational ideas, practical strategies, and current research frontiers in this area. Through examples in astrophysics, we'll explore how domain adaptation can help bridge the gap between synthetic and real-world data, improve trust in model outputs, and advance scientific discovery.

Ciprijanvoic, Aleksandra [Fermilab] (ORCID:0000000↗

Extended particle absorber for efficient modeling of intense laser-solid interactions

Here, an extended thermal particle boundary condition is devised to more efficiently and accurately model laser–plasma interactions in overdense plasmas. Particle-in-cell simulations of such interactions require many particles per cell, and a large region of background plasma is often necessary to correctly mimic a semi-infinite plasma and avoid electron refluxing from a truncated plasma. For long-pulse lasers of many picoseconds, such constraints can become prohibitively expensive. Here, an extended particle boundary condition (absorber) is designed that instantaneously stops and re-emits energetic particles streaming toward the simulation boundary over a defined region, allowing sufficient time and space for a suitably cool return current to develop in the background plasma. Tunable parameters of the absorber are explained, and simulations using the absorber with a 3-ps laser are shown to accurately reproduce those of a causally separated boundary while requiring only 20% the number of particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Topological and magnetic phase transitions in the bilayer Kitaev-Ising model

We investigate the phase diagram of a bilayer Kitaev honeycomb model with Ising interlayer interactions, deriving effective models via perturbation theory and performing Majorana mean-field theory calculations. We show that a diverse array of magnetic and topological phase transitions occur, depending on the direction of the interlayer Ising interaction and the relative sign of Kitaev interactions. When two layers have the same sign of the Kitaev interaction, a first-order transition from a Kitaev spin liquid to a magnetically ordered state takes place. The magnetic order points along the Ising axis and it is (anti)ferromagnetic for (anti)ferromagnetic Kitaev interactions. However, when two layers have opposite signs of the Kitaev interaction, we observe a notable weakening of magnetic ordering tendencies and the Kitaev spin liquid survives up to a remarkably larger interlayer exchange. Our mean-field analysis suggests the emergence of an intermediate gapped ℤ 2 spin-liquid state, which eventually becomes unstable upon vison condensation. The confined phase is described by a highly frustrated 120° compass model. We furthermore use perturbation theory to study the model with the Ising axis pointing along the $\hat{z}$ axis or lying in the 𝑥⁢𝑦 plane. In both cases, our analysis reveals the formation of one-dimensional Ising chains, which remain decoupled in perturbation theory, resulting in a subextensive ground-state degeneracy. Our results highlight the interplay between topological order and magnetic ordering tendencies in bilayer quantum spin liquids.

Kitaev model↗

EnergyPlus-MCP: A model-context-protocol server for ai-driven building energy modeling

Traditional building energy modeling with the EnergyPlus building performance simulation engine requires domain expertise, programming skills, and intensive manual efforts limiting its effective adoption. This paper introduces EnergyPlus-MCP, the first open-source Model Context Protocol (MCP) server specifically designed for EnergyPlus simulation workflows, establishing a new foundational infrastructure for AI-driven building energy modeling. The MCP server implements a layered architecture with 35 specialized tools spanning model management, editing and analysis, HVAC and other systems configuration inspection, and simulation execution, enabling Large Language Models to interact with EnergyPlus through conversational interfaces. The server addresses critical workflow barriers by automating model validation, streamlining energy efficiency measures modification, and providing intelligent output management with interactive visualization. Through practical demonstrations using a multi-zone building retrofit analysis, we show how the EnergyPlus-MCP server significantly reduces manual efforts while maintaining full simulation rigor. By providing accessible natural language interfaces to sophisticated building energy analysis, this approach enables scalable deployment of simulation expertise across public and private organizations, educational institutions, and research teams, fundamentally transforming traditional building energy modeling practices.

AI↗

Non-standard interactions in radiative neutrino mass models

Models of radiative Majorana neutrino masses require new scalars and/or fermions to induce lepton-number-violating interactions. We show that these new particles also generate observable neutrino non-standard interactions (NSI) with matter. We classify radiative models as type-I or II, with type-I models containing at least one Standard Model (SM) particle inside the loop diagram generating neutrino mass, and type-II models having no SM particle inside the loop. While type-II radiative models do not generate NSI at tree-level, popular models which fall under the type-I category are shown, somewhat surprisingly, to generate observable NSI at tree-level, while being consistent with direct and indirect constraints from colliders, electroweak precision data and charged-lepton flavor violation (cLFV). We survey such models where neutrino masses arise at one, two and three loops. In the prototypical Zee model which generates neutrino masses via one-loop diagrams involving charged scalars, we find that diagonal NSI can be as large as (8%, 3.8%, 9.3%) for (ε$_{ee}$, ε$_{μμ}$, ε$_{ττ}$), while off-diagonal NSI can be at most (10$^{–3}$%, 0.56%, 0.34%) for (ε$_{eμ}$, ε$_{eτ}$, ε$_{μτ}$). In one-loop neutrino mass models using leptoquarks (LQs), (ε$_{μμ}$, ε$_{ττ}$) can be as large as (21.6%, 51.7%), while ε$_{ee}$ and (ε$_{eμ}$, ε$_{eτ}$, ε$_{μτ}$) can at most be 0.6%. Other two- and three-loop LQ models are found to give NSI of similar strength. The most stringent constraints on the diagonal NSI are found to come from neutrino oscillation and scattering experiments, while the off-diagonal NSI are mostly constrained by low-energy processes, such as atomic parity violation and cLFV. We also comment on the future sensitivity of these radiative models in long-baseline neutrino experiments, such as DUNE. While our analysis is focused on radiative neutrino mass models, it essentially covers all NSI possibilities with heavy mediators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Developing an Interactive OpenMP Book with Large Language Models

Abstract. This paper presents an approach to authoring a textbook titled Interactive OpenMP Programming with the assistance of Large Language Models (LLMs). The writing process utilized state-of-the-art LLMs, including Gemini Pro 1.5, Claude 3, and ChatGPT-4, to generate the initial structure and outline of the book, as well as the initial content for specific chapters. This content included detailed descriptions of individual OpenMP constructs and practical programming examples. The outline and content have then undergone extensive manual revisions to meet our book goals. In this paper, we report our findings about the capabilities and limitations of these LLMs. We address critical questions concerning the necessity of textbook resources and the effectiveness of LLMs in creating fundamental and practical programming content. Our findings suggest that while LLMs offer significant advantages in generating textbook content, they require careful integration with traditional educational methodologies to ensure depth, accuracy, and pedagogical effectiveness. The Interactive OpenMP Programming book is developed with the framework of Jupyter Book, enabling the execution of code within the book from the web browser, providing instant feedback and a dynamic learning experience that stands in contrast to traditional educational resources. The book represents a significant step towards modernizing programming education, offering insights into practical strategies for generating the textbook through advanced AI tools.

Large Language Model · OpenMP · Interactive Book ·↗

Towards a unified model of neutrino-nucleus interactions

The achievement of the goals of the ongoing and future accelerator-based neutrino experiments – notably the determination of CP violating phase in the lepton sector – will require the development of advanced models of neutrino-nucleus interactions. In this review, we summarise the present status of experimental studies of neutrino-nucleus scattering, and discuss the developments and perspectives of the theoretical approach based on factorisation of the nuclear cross section, which has recently emerged as a promising framework for the description of the variety of processes contributing to the detected signals.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neutrinoless double beta decay via light neutralinos in R-parity violating supersymmetry

We perform a study of neutrinoless double beta (0νββ) decay mediated by the lightest neutralino of arbitrary mass in the minimal supersymmetric Standard Model (MSSM) under the presence of R-parity violating trilinear interactions. In this scenario, the exchange of the lightest neutralino can result in 0νββ decay of either long-range or short-range behaviour, depending on the neutralino mass. Using nuclear matrix elements calculated in the interacting boson model, we use an interpolation between the long- and short-range behaviours with an approximate formula. The non-observation of 0νββ decay is then used to place constraints on the supersymmetry parameter space, compatible with constraints from collider experiments. We compare these constraints to bounds from pion decays, CKM unitarity and big bang nucleosynthesis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Molecular Dynamics Investigation of Lipid-Specific Interactions with a Fusion Peptide

The HIV-1 fusion peptide, which is a short hydrophobic peptide from the gp41 coat glycoprotein that participates in the infection of a cell, interacts with model lipid bilayer membranes in a concentration-dependent manner. The interaction of the peptide with the bilayer also strongly depends on the lipid composition. Here, molecular dynamics simulations were performed to investigate lipid-specific interactions that arise shortly after the binding of a less-fusogenic variant of the HIV-1 fusion peptide to a lipid bilayer composed of a mixture of dimyristoyl phosphatidylcholine and dimyristoyl phosphatidylglycerol. The impact of peptide concentration was also studied. An improved understanding was gained of the lipid-specific interactions experienced by the FP. New insight was also gained into how the peptide concentration changes these interactions.

60 APPLIED LIFE SCIENCES↗

Modeling of vapor-liquid interactions in condensing ejectors

Ejectors are compact mechanical devices that utilize the expansion of a high-pressure primary fluid to entrain and compress a low-pressure secondary fluid by means of momentum transfer between the two streams of fluid. Condensing ejectors feature both momentum and heat transfer through the interaction between the vapor and liquid streams. The goal of this paper is to develop and validate a one-dimensional slug flow model to simulate the vapor–liquid interactions in two types of condensing ejectors – one with primary liquid and secondary vapor flows (Type I), and the other with primary vapor and secondary liquid flows (Type II). Control volume analysis of the mass, momentum, and energy balance in each phase and across the interface was conducted for the ejectors. The slug flow models for both types of ejectors are validated with published results in the literature. The friction coefficient on the inner wall of the ejectors and the interfacial heat transfer coefficient are identified as controlling parameters for the ejector performance in terms of pressure and temperature distributions along the axis. The detailed parametric study shows that the liquid inlet velocity and mixing tube diameter have a significant impact on the performance of the Type I ejectors, and the performance of the Type II ejectors is mainly controlled by the secondary liquid flow inlet temperature.

42 ENGINEERING↗