Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulator study”

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 811 records · Page 45

Understanding changes in cloud simulations from E3SM version 1 to version 2

Abstract. This study documents clouds simulated by the Energy Exascale Earth System Model (E3SM) version 2 (E3SMv2) and attempts to understand what causes the model behavior change in clouds relative to E3SMv1. This is done by analyzing the last 30-year (1985–2014) data from the 165-year historical simulations using E3SMv1 and v2 and four sensitivity tests to isolate the impact of changes in model parameter choices in its turbulence, shallow convection, and cloud macrophysics parameterization (Cloud Layers Unified By Binormals, CLUBB); microphysical parameterization (MG2); and deep-convection scheme (ZM), as well as model physics changes in convective triggering. It is shown that E3SMv2 significantly improves the simulation of subtropical coastal stratocumulus clouds and clouds with optical depth larger than 3.6 over the stratocumulus-to-cumulus transition regimes, where the shortwave cloud radiative effect (SWCRE) is also improved, and the Southern Ocean (SO) while seeing an overall slight degradation in low clouds over other tropical and subtropical oceans. The better performance in E3SMv1 over those regions is partially due to error compensation between its simulated optically thin and intermediate low clouds for which E3SMv2 actually improves simulation of optically intermediate low clouds. Sensitivity tests indicate that the changes in low clouds are primarily due to the tuning done in CLUBB. The impact of the ZM tuning is mainly on optically intermediate and thick high clouds, contributing to an improved SWCRE and longwave cloud radiative effect (LWCRE). The impact of the MG2 tuning and the new convective trigger is primarily on the high latitudes and the SO. They have a relatively smaller impact on clouds than CLUBB tuning and ZM tuning do. This study offers additional insights into clouds simulated in E3SMv2 by utilizing multiple data sets and the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP) diagnostic tool as well as sensitivity tests. The improved understanding will benefit future E3SM developments.

58 GEOSCIENCES↗

A Comparison of Human Error Probabilities Collected from HuREX and SHEEP Frameworks

This paper discusses how different are the HEPs collected from HuREX and SHEEP frameworks. This study is a preceding research to infer full-scope HRA data based on the data collected from the SHEEP framework. For the comparison, we used HEPs in the HuREX database published in KAERI-TR-6649 [6]. The HEPs have been collected from actual licensed operators when manipulating MCR simulators of Westinghouse type and Optimized Power Reactor (OPR1000) type in South Korea. For the HEPs based on the SHEEP framework, these have been collected from actual licensed operators and students when using a simplified simulator, i.e., Rancor Microworld.

99 GENERAL AND MISCELLANEOUS↗

CO 2 Capture Characteristics of Hyperbranched Poly(alkylene imine): A Molecular Dynamics Simulation Approach

This study explores the CO 2 capture characteristics of hyperbranched poly(ethylenimine) (HB-PEI) and poly- (propyleneimine) (HB-PPI) through molecular dynamics simulations using density functional theory-calibrated force fields. Key features such as density, free volume, glass transition temperature, CO 2 /H 2 O distribution, and molecular diffusion are systematically investigated to elucidate structure−function relationships under dry and hydrated conditions. HB-PEI demonstrates a slightly higher density and lower free volume compared to HB-PPI yet shows superior CO 2 capture due to the high amine concentration. Glass transition analysis indicates a higher thermal mobility in HBPEI, enhancing the CO 2 diffusivity. Pair correlation and coordination analyses confirm a stronger affinity of CO 2 with primary and secondary amines, particularly in hydrated environments where water competes with CO 2 for binding sites. Despite its more compact structure, HB-PEI outperformed HB-PPI in CO 2 and H 2 O transport, as confirmed by higher diffusion coefficients across all hydration levels. These findings highlight a critical balance among polymer architecture, amine accessibility, and hydration in designing next-generation solid amine sorbents for efficient direct air capture applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lattice Symmetry-Guided Charge Transport in 2D Supramolecular Polymers Promotes Triplet Formation

Singlet-to-triplet intersystem crossing (ISC) in organic molecules is intimately connected with their geometries: by modifying the molecular shape, symmetry selection rules pertaining to spin-orbit coupling can be partially relieved, leading to extra matrix elements for increased ISC. As an analog to this molecular design concept, the study finds that the lattice symmetry of supramolecular polymers also defines their triplet formation efficiencies. A supramolecular polymer self-assembled from weakly interacting molecules is considered. Its 2D oblique unit cell effectively renders it as a coplanar array of 1D molecular columns weakly bound to each other. Using momentum-resolved photoluminescence imaging in combination with Monte Carlo simulations, the study found that photogenerated charge carriers in the supramolecular polymer predominantly recombine as spin-uncorrelated carrier pairs through inter-column charge transfer states. This lattice-defined recombination pathway leads to a substantial triplet formation efficiency (≈60%) in the supramolecular polymer. These findings suggest that lattice symmetry of micro-/macroscopic structures relying on intermolecular interactions can be strategized for controlled triplet formation.

36 MATERIALS SCIENCE↗

Cooperative Role of Mixed Solvent in the Evaporation-Induced Self-Assembly of Polypeptoid Nanocrystals

Peptoids, or polypeptoids, are biomimetic polymers that can self-assemble into nanocrystals for biomedical and biotechnological applications. Polypeptoid nanocrystals can be prepared by evaporation-induced self-assembly, but the roles of solvent components for this process have long been overlooked at the molecular level, leaving a tunable parameter for improving self-assembly protocols. This work utilized molecular dynamics simulations to study the effects of water and the commonly used tetrahydrofuran (THF) on the assembly of nanosheets from molecules of acetylated diblock polypeptoid, poly-(N-decylglycine)-b-poly-(N-2-(2-(2-methoxyethoxy)-ethoxy) ethylglycine), abbreviated as Ac-Ndc10-Nte10. To probe the stages of self-assembly, isolated molecules and preassembled nanofibers/nanosheets were simulated in pure THF, water, and their mixtures, respectively. The assembly energies show that the THF/water mixture has a greater tendency to form nanosheets than pure water. In a THF/water mixture, polypeptoids were found more uncoiled in isolated states, less compact in disordered agglomerates, and with reduced requirement for the Nte block to cover hydrophobic Ndc surfaces in the nanocrystals. Mixed solvent is vital to initiating self-assembly, as THF assists in the opening of coiled polypeptoid molecules, while water provides the thermodynamics to aggregate and ultimately form nanocrystals. To obtain wider nanosheets, it is recommended that some THF be maintained in the aqueous solvent before it becomes exhausted by evaporation. Near the nanosheet surface, the THF concentration is higher than that in bulk solution (3-4 times in 4 M THF/water). The strong adsorption of THF indicates the self-assembly in a de facto mixed solvent. These results are expected to guide the refinement of evaporation-induced self-assembly protocols for polypeptoid nanocrystals.

Luo, Xubo↗

Machine Learning-Driven Optimization of Building Enclosures for Moisture Durability and Thermal Performance

The design of moisture-durable building enclosures with low embodied carbon often involves an iterative process of selecting the materials for the specific exposure conditions to meet the performance requirements. While hygrothermal simulations are commonly used to evaluate moisture durability, they often require advanced expertise for proper implementation. Machine learning (ML) provides a promising alternative by streamlining the design process and minimizing the reliance on complex simulations. This study presents a machine learning-based approach for predicting moisture durability in residential wall assemblies. The ML model was trained to estimate the mold index and maximum moisture content of various layers under typical exposure conditions. The model achieved a high predictive accuracy, with a coefficient of determination (R²) exceeding 0.90 when compared to traditional hygrothermal simulations on materials that were not part of training the ML model. Building on these results, the ML model was developed into a practical tool for optimizing wall assembly designs. This tool allows users to automatically optimize material selections based on energy, moisture, and carbon performance criteria. By incorporating multi-objective optimization, the tool identifies configurations that minimize embodied carbon while maintaining moisture safety and code-compliant thermal performance. Additionally, it provides insights into how material choices influence assembly durability, energy efficiency, and carbon reduction. The tool will be implemented in the Building Science Advisor (BSA) to enhance its performance and provide more granularity on the results. This research highlights the potential for ML-driven tools to simplify the design of high-performance building enclosures, offering architects and engineers a faster, more efficient way to balance critical performance factors.

Salonvaara, Mikael [ORNL] (ORCID:0000000318991554)↗

Evolution of Thyroglobulin Loop Kinetics in EpCAM

Epithelial cell-activating molecule (EpCAM) is an important cancer biomarker and therapeutic target given its elevated expression in epithelial cancers. EpCAM is a type I transmembrane protein that forms cis-dimers along the thyroglobulin type-1A-like domain (TYD) in the extracellular region. The thyroglobulin loop (TY loop) within the TYD is structurally dynamic in the monomer state of human EpCAM, binding reversibly to a TYD site. However, it is not known if this flexibility is prevalent across different species. Here, we conduct over 17 μs of all-atom molecular dynamics simulations to study EpCAM TY loop kinetics of five different species, including human, mouse, chicken, frog, and fish. We find that the TY loop remains dynamic across evolution. In addition to the TYD binding site, we discover a second binding site for the TY loop in the C-terminal domain (CTD). Calculations of the dissociation rate constants from the simulation trajectories suggest a differential binding pattern of fish EpCAM and other organisms. Whereas fish TY loop has comparable binding for both TYD and CTD sites, the TY loops of other species preferably bind the TYD site. A hybrid construct of fish EpCAM with human TY loop restores the TYD binding preference, suggesting robust effects of the TY loop sequence on its dynamic behavior. Our findings provide insights into the structural dynamics of EpCAM and its implication in physiological functions.

59 BASIC BIOLOGICAL SCIENCES↗

Classical Simulation of Boson Sampling Based on Graph Structure

Boson sampling is a fundamentally and practically important task that can be used to demonstrate quantum supremacy using noisy intermediate-scale quantum devices. In this Letter, we present classical sampling algorithms for single-photon and Gaussian input states that take advantage of a graph structure of a linear-optical circuit. The algorithms’ complexity grows as so-called treewidth, which is closely related to the connectivity of a given linear-optical circuit. Using the algorithms, we study approximated simulations for local Haar-random linear-optical circuits. For equally spaced initial sources, we show that, when the circuit depth is less than the quadratic in the lattice spacing, the efficient simulation is possible with an exponentially small error. Notably, right after this depth, photons start to interfere each other and the algorithms’ complexity becomes subexponential in the number of sources, implying that there is a sharp transition of its complexity. Finally, when a circuit is sufficiently deep enough for photons to typically propagate to all modes, the complexity becomes exponential as generic sampling algorithms. We numerically implement a likelihood test with a recent Gaussian boson sampling experiment and show that the treewidth-based algorithm with a limited treewidth renders a larger likelihood than the experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

What do far-infrared spectra of solitary water in “water-in-solvent” systems reveal about water’s solvation and dynamics?

Classical molecular dynamics simulations of water in ionic and dipolar solvents were used to interpret the far-infrared (FIR) rotation/libration spectra of “solitary water” in terms of water’s rotational dynamics and interactions with solvents. Seven solvents represented by nonpolarizable all-atom force fields and a series of idealized variable-charge solvents were used to span the range of solvent polarities (hydrogen bonding) studied experimentally. Simulated spectra capture the solvent dependence observed, as well as the relationship between the frequencies of water libration (νL) and OH stretching bands (νOH). In more strongly interacting solvents, simulated νL are ∼20% higher than those of experiment. In all solvents, the simulated spectra are composites of rotational motions about the two axes perpendicular to water’s dipole moment, and the different frequencies of these two motions are responsible for the breadth of the libration band and the bimodal shape observed in halide ionic liquids. Simulations overestimate the separation of these two components in most solvents. The character of water rotational motions changes markedly with solvent polarity, from quasi-free rotation in nonpolar and weakly polar solvents to highly constrained libration in strongly hydrogen bonding environments. The changeover to librational motions dominating the spectrum occurs between solvents such as benzene (νL ∼ 250 cm−1) and acetonitrile (νL ∼ 400 cm−1). For solvents in the latter category, the mean frequency of the experimental FIR band provides a direct measure of mean-squared torques and, therefore, force constants associated with interactions constraining water’s librational motion.

Chemistry↗

Modeling Adsorption in Silica Pores via Minkowski Functionals and Molecular Electrostatic Moments

Capillary condensation phenomena are important in various technological and environmental processes. Using molecular simulations, we study the confined phase behavior of fluids relevant to carbon sequestration and shale gas production. As a first step toward translating information from the molecular to the pore scale, we express the thermodynamic potential and excess adsorption of methane, nitrogen, carbon dioxide, and water in terms of the pore’s geometric properties via Minkowski functionals. This mathematical reconstruction agrees very well with molecular simulations data. Our results show that the fluid molecular electrostatic moments are positively correlated with the number of adsorption layers in the pore. Moreover, stronger electrostatic moments lead to adsorption at lower pressures. These findings can be applied to improve pore-scale thermodynamic and transport models.

42 ENGINEERING↗

Grid impact analysis using controller-hardware-in-the-loop for high-power vehicle charging stations

A controller-hardware-in-the-loop (CHIL) architecture for the evaluation of grid impacts arising from high-power vehicle charging stations is presented in this paper. Unlike simulation-based studies, the proposed method can be used to capture the interactions of the grid and the charging load along with charger controllers in real time. The proposed method can be used to evaluate the impact of charging load on the grid in terms of voltage variations and line congestion. The proposed CHIL platform allows for de-risking the vehicle charging station deployment by simulating the complex interactions among all the components of a vehicle charging station - i.e., the grid, vehicle, and charger controller - in a realistic manner before using the charging station in a grid. Further, the proposed CHIL approach can be used to evaluate the voltage regulation causalities of the vehicle charging station. Experimental results are presented in the paper to illustrate the applicability of the proposed method in a laboratory environment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Ensuring Transient Stability With Guaranteed Region of Attraction in DC Microgrids

DC microgrids have promising applications in renewable integration due to their better energy efficiency when connecting DC components. However, they might be unstable since many loads in a DC microgrid are regulated as constant power loads (CPLs) that have a destabilizing negative impedance effect. As a result, the state trajectory displacement caused by abrupt load changes or contingencies can easily lead to instability. Many existing works have been devoted to studying the region of attraction (ROA) of a DC microgrid, in which the system is guaranteed to be asymptotically stable. Nevertheless, existing work either focuses on using numerical methods for ROA approximations that generally have no performance guarantees or cannot ensure a desired ROA for a general DC microgrid. To close this gap, this paper develops an innovative control synthesis algorithm to make a general DC microgrid have a theoretically guaranteed ROA, for example, to cover the entirety of its operating range regarding state trajectories. Here, we first study the nonlinear dynamics of a DC microgrid to derive a novel transient stability condition to rigorously certify whether a given operating range is a subset of the ROA; then, we formulate a control synthesis optimization problem to guarantee the condition’s satisfaction. This condition is a linear constraint, and the optimization problem resembles an optimal power flow problem and has a good computational behavior. Simulation case studies verify the validity of the proposed work.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Large-Scale Adoption of Self-Synchronized Universal Droop Controller-Based Inverters to Enable Ancillary Services for Different Modes of Distribution System Operations

This work proposes the large-scale adoption of self-synchronized universal droop controller (SUDC)-based inverters to enable ancillary services for different modes of distribution system operations. The IEEE 123 bus system was modeled on a real-time simulator to study the performance of large-scale adoption of SUDC inverters in a distribution system. The resulting data collected shows that the voltage and the frequency were regulated within ranges, such as less than 5% for voltage and less than 0.5% for frequency, under different load variations and grid operations. Also, the black start was achieved within 0.4 s without any voltage overshoot. Through the simulation and validation on a small microgrid and the IEEE 123 bus distribution system, it can be concluded that the SUDC was successfully adopted to regulate the voltage and the frequency within the given ranges, and black start achieved within 1 s without voltage overshoot for different modes of distribution system operations.

14 SOLAR ENERGY↗

Simulation of Multipaction in LB650 Cavity with Weld Profile

This study aims to simulate multipaction in LB650 cavity with weld profile, a first step to replicating an actual cavity, whereas past studies have ignored the same. We attempt to align the simulation data with the experimental data and report the subsequent findings.

Ninan, Z.↗

Numerical simulation of gas production from natural gas hydrate deposits with multi-branch wells: Influence of reservoir properties

Vast amounts of natural gas hydrate are buried in subseafloor sediments without impermeable boundaries, which is recognized as an essential energy source for the future. Previously, the multi-branch well was proposed to enhance the recovery efficiency of natural gas hydrate, and the gas production rate has been dramatically improved comparing with the vertical well. However, the multi-branch well shows a terrible performance in gas production duration. Here, as a continuation of the previous study, numerical simulations were conducted to investigate the influence of hydrate reservoir properties on the gas production potential. Results indicate that it is hard to extract hydrate commercially for hydrate accumulations without impermeable boundaries. A high initial hydrate saturation leads to a long gas production duration but a low gas production rate. An increase in the intrinsic permeability of isotropic reservoirs would shorten the gas production duration and result in a low gas recovery ratio. Permeability anisotropy shows a noticeable effect on enhancing the gas recovery ratio and the gas production duration due to the improved pressure propagation pattern. Therefore, in the upcoming field tests, reservoir reconstructions that enhance permeability anisotropy are strongly suggested to obtain better outcomes.

03 NATURAL GAS↗

Optimal Gas-Electric Energy System Decarbonization Planning

As energy utilities implement climate change mitigation policies, system planners require strategies for achieving affordable emissions reductions. Coordinated planning of electric power and natural gas systems will allow synergistic investments to address cross-sector operational constraints, competing uses for net-zero emissions fuels, and shifts in energy demands across energy carriers. In this study, we develop a novel optimization program that finds the cost-minimizing mix of infrastructure expansion or reduction across gas and electric systems to satisfy sector-specific emissions constraints. Alongside energy supply resources, our framework allows for central-planning of end-use equipment stocks to allow switching between gas and electric appliances upon failure or premature replacement. The proposed model is used to simulate case study scenarios for a benchmark 24-pipe gas network coupled to a 24-node power system test network. We find that electrification of greater than 80% of core gas demands is a component of the least-cost solution for modeled energy systems. Despite this substitution, the gas system is maintained to service difficult-to-electrify customers and to deliver net-zero emissions gas to electricity generators in times of peak electricity demand. Restricting electrification of gas appliances increases reliance on power- to-gas technologies and increases annual costs by 15% in 2040. Neglecting constraints on pipeline blending of hydrogen can produce a misleading result that relies on hydrogen blend fractions of greater than 50%. In all cases, we find the average costs of delivered gas increase nearly 5-fold across the decarbonization transition, highlighting the importance of future work investigating cost-allocation strategies for ensuring an equitable energy transition.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Aitken Mode Aerosols Buffer Decoupled Mid‐Latitude Boundary Layer Clouds Against Precipitation Depletion

Aerosol-cloud-precipitation interactions are a leading source of uncertainty in estimating climate sensitivity. Remote marine boundary layers where accumulation mode (~100–400 nm diameter) aerosol concentrations are relatively low are very susceptible to aerosol changes. These regions also experience heightened Aitken mode aerosol (~10–100 nm) concentrations associated with ocean biology. Aitken aerosols may significantly influence cloud properties and evolution by replenishing cloud condensation nuclei and droplet number lost through precipitation (i.e., Aitken buffering). We use a large-eddy simulation with an Aitken-mode enabled microphysics scheme to examine the role of Aitken buffering in a mid-latitude decoupled boundary layer cloud regime observed on 15 July 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic flight campaign: cumulus rising into stratocumulus under elevated Aitken concentrations (~100–200 mg -1 ). In situ measurements are used to constrain and evaluate this case study. Our simulation accurately captures observed aerosol-cloud-precipitation interactions and reveals time-evolving processes driving regime development and evolution. Aitken activation into the accumulation mode in the cumulus layer provides a reservoir for turbulence and convection to carry accumulation aerosols into the drizzling stratocumulus layer above. Further Aitken activation occurs aloft in the stratocumulus layer. Together, these activation events buffer this cloud regime against precipitation removal, reducing cloud break-up and associated increases in heterogeneity. We examine cloud evolution sensitivity to initial aerosol conditions. With halved accumulation number, Aitken aerosols restore accumulation concentrations, maintain droplet number similar to original values, and prevent cloud break-up. Without Aitken aerosols, precipitation-driven cloud break-up occurs rapidly. In this regime, Aitken buffering sustains brighter, more homogeneous clouds for longer.

54 ENVIRONMENTAL SCIENCES↗

Control Design and EMT Simulation Toward the Restoration of Faulted Wind-Dominant Grids

Present power engineers have been challenged to employ electromagnetic transient (EMT) simulations to study the restoration of wind-dominant grids. This paper addresses this hurdle by developing strategies to: (i) reliably control the dc-link voltage of grid-forming wind turbines with dc-coupled batteries, (ii) tune controllers to stably interconnect wind power plants into a grid under recovery, and (iii) autonomously make restoration decisions if transmission faults occur during restoration. These contributions are assessed via EMT simulations of wind-dominant versions of the WSCC 9-bus grid and a 22-bus power system. Furthermore, this paper is significant to address recommendations by the North American Electric Reliability Corporation.

17 WIND ENERGY↗