Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy simulation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Phase space distribution functions and energy distributions of dark matter particles in haloes

ABSTRACT For a spherical dark matter halo with isotropic velocity distribution, the phase space distribution function (DF), the energy distribution, and the density profile form a set of self-consistent description of its equilibrium state, and knowing one is sufficient to determine the other two. The Navarro–Frenk–White density profile (NFW profile) is known to be a good approximation to the spherically averaged density distribution in simulated haloes. The DARKexp energy distribution is also known to compare well with the simulated energy distribution. We present a quantitative assessment of the NFW and DARKexp fits to the simulated DF and energy distribution for a wide range of haloes in a dark-matter-only simulation from the IllustrisTNG Project. As expected, we find that the NFW fits work well except at low energy when the density at small radii deviates from the NFW profile. Further, the NFW and DARKexp fits have comparable accuracy in the region where both fit well, but the DARKexp fits are better at low energy because they require matching of the central gravitational potential. We also find an approximate relation between the energy scale parametrizing the DARKexp energy distribution and that defined by the characteristic density and radius of the NFW profile. This relation may be linked to the relaxation process during halo formation.

Astronomy & Astrophysics↗
DOE OSTI2023

Hybrid renewable energy systems: the value of storage as a function of PV-wind variability

As shares of variable renewable energy (VRE) on the electric grid increase, sources of grid flexibility will become increasingly important for maintaining the reliability and affordability of electricity supply. Lithium-ion battery energy storage has been identified as an important and cost-effective source of flexibility, both by itself and when coupled with VRE technologies like solar photovoltaics (PV) and wind. In this study, we explored the current and future value of utility-scale hybrid energy systems comprising PV, wind, and lithium-ion battery technologies (PV-wind-battery systems). Using a price-taker model with simulated hourly energy and capacity prices, we simulated the revenue-maximizing dispatch of a range of PV-wind-battery configurations across Texas, from the present through 2050. Holding PV capacity and point-of-interconnection capacity constant, we modeled configurations with varying wind-to-PV capacity ratios and battery-to-PV capacity ratios. We found that coupling PV, wind, and battery technologies allows for more effective utilization of interconnection capacity by increasing capacity factors to 60%–80%+ and capacity credits to close to 100%, depending on battery capacity. We also compared the energy and capacity values of PV-wind and PV-wind-battery systems to the corresponding stability coefficient metric, which describes the location-and configuration-specific complementarity of PV and wind resources. Our results show that the stability coefficient effectively predicts the configuration-location combinations in which a smaller battery component can provide comparable economic performance in a PV-wind-battery system (compared to a PV-battery system). These PV-wind-battery hybrids can help integrate more VRE by providing smoother, more predictable generation and greater flexibility.

14 SOLAR ENERGY↗
DOE OSTI2025

Convergence and Quantum Advantage of Trotterized MERA for Strongly-Correlated Systems

Strongly-correlated quantum many-body systems are difficult to study and simulate classically. We recently proposed a variational quantum eigensolver (VQE) based on the multiscale entanglement renormalization ansatz (MERA) with tensors constrained to certain Trotter circuits. Here, we determine the scaling of computation costs for various critical spin chains which substantiates a polynomial quantum advantage in comparison to classical MERA simulations based on exact energy gradients or variational Monte Carlo. Algorithmic phase diagrams suggest an even greater separation for higher-dimensional systems. Hence, the Trotterized MERA VQE is a promising route for the efficient investigation of strongly-correlated quantum many-body systems on quantum computers. Furthermore, we show how the convergence can be substantially improved by building up the MERA layer by layer in the initialization stage and by scanning through the phase diagram during optimization. For the Trotter circuits being composed of single-qubit and two-qubit rotations, it is experimentally advantageous to have small rotation angles. We find that the average angle amplitude can be reduced considerably with negligible effect on the energy accuracy. Benchmark simulations suggest that the structure of the Trotter circuits for the TMERA tensors is not decisive; in particular, brick-wall circuits and parallel random-pair circuits yield very similar energy accuracies.

Miao, Qiang [Duke Quantum Center, Duke University,↗
DOE OSTI2025

Refining fast calorimeter simulations with a Schrödinger Bridge

Machine learning-based simulations, especially calorimeter simulations, are promising tools for approximating the precision of classical high energy physics simulations with a fraction of the generation time. Nearly all methods proposed so far learn neural networks that map a random variable with a known probability density, like a Gaussian, to realistic-looking events. In many cases, physics events are not close to Gaussian and so these neural networks have to learn a highly complex function. We study an alternative approach: Schrödinger bridge Quality Improvement via Refinement of Existing Lightweight Simulations (SQuIRELS). SQuIRELS leverages the power of diffusion-based neural networks and Schrödinger bridges to map between samples where the probability density is not known explicitly. We apply SQuIRELS to the task of refining a classical fast simulation to approximate a full classical simulation. On simulated calorimeter events, we find that SQuIRELS is able to reproduce highly non-trivial features of the full simulation with a fraction of the generation time.

Calorimeter methods↗
DOE OSTI2020

Evaluation of Uplight Energy Saving Actions Using ResStock TM (Summary of Results)

As part of a portfolio of product offerings, Uplight produces energy savings estimates to inform a utility’s residential customers of the expected benefits that would arise from adopting energy efficiency measures, which are termed “energy saving actions. ”Uplight’s current practice is to use a proprietary building simulation model to estimate energy and cost benefits from these actions, as discussed in Maguire et al. Uplight requested that the National Renewable Energy Laboratory (NREL) provide an independent estimate of the energy savings for a defined list of energy saving actions. NREL performed this study using two established building energy simulation tools: EnergyPlus ® ( EnergyPlus ) and ResStock TM .(

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

LINAC Longitudinal Simulation and Measurement of Output Energy

The Linac at Fermilab, one of the oldest accelerators on-site, plays a crucial role in delivering 400 MeV proton beams to other accelerating facilities. The beam traveling through Linac needs to be carefully monitored for better accelerator operation. One of the key factors of efficient Linac operation is its longitudinal stability. The stability can minimize Linac loss and maximize beam current. For this purpose, I examined the beam's longitudinal profile, which is crucial for proper beam transfer between its two sections. After careful monitoring of the beam through Linac, we need to ensure that it leaves with desired parameters. The beam coming out from the Linac is subsequently injected into the Booster, a synchrotron that relies on dipole magnets to bend the beam at injection. If the output energy of the beam deviates from the expected value, the beam will bend at an incorrect radius, leading to potential beam loss. Therefore, precise regulation of the output energy is critical for the efficient operation of the accelerator. To address this, I have conducted a thorough examination of the Linac longitudinal simulations and data from Beam Position Monitors (BPMs). My objective is to measure and stabilize the output energy and verify its correspondence with the predicted simulation results. All these analyses are important for efficient operation of the accelerator.

Safaryan, Milena↗
DOE OSTI2024

Data Driven Commercial Building Energy Code Compliance and Technology Inventory for New York City

Building Performance Standards (BPS) are gaining national traction. A BPS will require new processes in the design, construction, and operation of buildings that take the occupants into account and enable predictive analysis to ensure compliance with current and future GHG emissions caps. In New York City, most buildings over 25,000 square feet will be regulated by a BPS starting in 2024, regardless of whether it is new construction permitted under current energy codes or an existing building. This research is one of the first to begin the evaluation of a long-term series of building policies in the context of an open data ecosystem, in cooperation with city agencies. Existing building policies enacted in NYC have ranged from building energy benchmarking and labeling to energy audits to the regulation of GHG emission in buildings. Through the development of a dataset related to building technologies and energy consumption, this project can help to evaluate if meaningful conclusions can be drawn for the data that has been largely self-reported in compliance with city regulations. This project will also provide lessons learned from a deep dive into these types of datasets to provide best practices for municipalities or states seeking to embark on policies like those enacted in NYC. In addition, a Building Automation System (BAS) Stretch Standard of Care (SSOC) for owners, designers, and building operators will enable the measurement and predictive analysis of energy consumption and GHG emissions at the plant, system, or component level, in anticipation of regulated GHG limits on buildings based on energy use. The SSOC is expected to be suitable for use on a national level. The primary feature of an SSOC is a standardized format for a set of BAS points that can be used to control and to gather data from individual plants, systems, or components that are related to building energy consumption. This project examined how measurements compare to prescriptive or simulation-based energy code targets, finding little correlation between predictive 8760-hour energy modeling and actual energy consumption for a small sample (n=27) of buildings constructed after 2015. Other analysis found that, while large multifamily housing (MFH) buildings showed a general trend similar to predicted reductions in energy use from the implementation of model commercial energy codes, this trend was not evident in the office, K-12 school, and hotel use groups in NYC. No upward or downward trends in energy consumption were found when buildings were grouped by size. Energy audit data were analyzed and it appears that there is bias by audit company on measures recommended to clients. Further research should be performed to cross-analyze this with other attributes, such as building size, vintage, and number of stories. Analysis found that for 281 buildings that were permitted and completed after 2015 and had submitted benchmarking data in 2022, between 81% and 96% (by use group) were found to be in compliance with the 2024 to 2029 NYC BPS emission caps, and between 55% and 89% were in compliance with the 2030-2034 caps. This work is beneficial to the public in helping policymakers and building stakeholders better understand the wide-ranging implications of a BPS.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Neutral Pion Electroproduction and development of a Neutral Particle Spectrometer

Protons and neutrons, i.e. nucleons, are the basic building blocks of the matter in the visible universe. The strong force binds the nucleons to form nuclei. The electromagnetic force forms the atoms by binding the electrons with the nuclei. The electromagnetic interaction is well understood by Quantum Electrodynamics (QED), which shows the most precise predictability amongst all the theories in physics. In QED, charges interact with each other by exchanging photons. The Quantum Chromodynamics (QCD) describes the strong interaction. Its degrees of freedom are quarks and gluons, the fundamental constituents of the nucleons. The quarks interact with each other by exchanging gluons. However, unlike QED, the gluons interact amongst themselves. This feature of self bindings of the gluons confines the quarks and gluons in the nucleons/hadrons, never to be seen as free. In order to study some of the features of QCD, such as confinement or the structure of the nucleon, one usually needs to rely on experiments. Electromagnetic probes, governed by the well-understood QED, are excellent tools to probe the nucleon. In general, different scales, e.g. electron beam energies, probe different regions of the nucleon. At low energy, of the order of a few GeV, the electron probes the nucleon in the valence quark region. As its energy increases, the electron probes the sea quark and gluon regions. The study of the nucleon structure in all these regions is needed to fully understand QCD. Form factors and parton distribution functions measured from elastic scattering and deep inelastic scattering of leptons off nucleons have provided a partial view of the internal structure of the nucleon. In the mid-1990s, Generalized Parton Distributions (GPDs) were developed. These new objects are a generalization of the form factors and parton distribution functions, but contain richer information on the nucleon internal structure. GPDs are accessible experimentally by deep exclusive reactions. Deeply virtual Compton scattering (DVCS) and deeply virtual meson production (DVMP) are some examples. The first dedicated DVCS/DVMP experiment took place in 2004 in Hall A at The Thomas Jefferson National Accelerator Facility, i.e. Jefferson Lab, in Virginia, U.S.A. A new DVCS/DVMP experiment, after the beam energy upgrade of Jefferson Lab, was carried out in Hall A in a wider kinematic range. Its data were taken from 2014 to 2016. In Hall C at Jefferson Lab, the next DVCS/DVMP experiment will take place. The Hall C experiment will further exploit the kinematic range with higher precision. A Neutral Particle Spectrometer (NPS) is in development to measure DVCS/DVMP events under high background conditions. Jefferson Lab will provide the highest precision data in the valence quark region for various exclusive reactions. The Electron-Ion Collider (EIC) is a future experimental facility currently planned to start operations around 2030 in the U.S.A. Its high energy and high luminosity will probe the sea quark and gluon regions providing answers to the outstanding questions of QCD, in particular in the region where matter is dominated by gluons. First of all, this document describes the data analysis and results of the Hall A neutral pion electroproduction off the proton, from the data taken in 2014-2016. Later, some of the developments towards the construction of the electromagnetic calorimeter of the NPS for the upcoming DVCS/DVMP experiment in Hall C are presented. Finally, one of the candidate materials for the EIC calorimeter, a glass scintillator, will be briefly introduced. I have participated to all these projects, in collaboration with many colleagues. I present in this thesis my contributions to each of these projects. My contributions to the neutral pion data analysis were focused on background subtractions on the calorimeter, acceptance calculations, and the estimation of the systematic uncertainty associated to the event selection cuts. Some necessary information on calibrations of the detectors and data analysis methods are also described. In the NPS project, I performed background dose calculations and energy and position resolution studies of the calorimeter, all using Monte Carlo simulations, with realistic geometries of the experimental apparatus. Characterization of the crystals of the calorimeter was also done. Additionally, I measured the radiation hardness of some glass scintillator in its early stage of development. In order to have a future reference when the glass calorimeter prototype will be tested, I simulated the energy resolution of the prototype.

Ko, Ho-San↗
DOE OSTI2023

Advances in the Co-Simulation of Detailed Electrical and Whole-Building Energy Performance

This article describes recent co-simulation advances for the simultaneous modeling of detailed building electrical distribution systems and whole-building energy performance. The co-simulation architecture combines the EnergyPlus® engine for whole-building energy modeling with a new Modelica library for building an electrical distribution system model that is based on harmonic power flow. This new library allows for a higher-fidelity modeling of electrical power flows and losses within buildings than is available with current building electrical modeling software. We demonstrate the feasibility of the architecture by modeling a simple, two-zone thermal chamber with internal power electronics converters and resistive loads, and we validate the model using experimental data. The proposed co-simulation capability significantly expands the capabilities of building electrical distribution system models in the context of whole-building energy modeling, thus enabling more complex analyses than would have been possible with individual building performance simulation tools that are used to date.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

RouteE: A Vehicle Energy Consumption Prediction Engine

The emergence of connected and automated vehicles and smart cities technologies create the opportunity for new mobility modes and routing decision tools, among many others. To achieve maximum mobility and minimum energy consumption, it is critical to understand the energy cost of decisions and optimize accordingly. The Route Energy prediction model (RouteE) enables accurate estimation of energy consumption for a variety of vehicle types over trips or sub-trips where detailed drive cycle data are unavailable. Applications include vehicle route selection, energy accounting and optimization in transportation simulation, and corridor energy analyses, among others. The software is a Python package that includes a variety of pre-trained models from the National Renewable Energy Laboratory (NREL). However, RouteE also enables users to train custom models using their own data sets, making it a robust and valuable tool for both fast calculations and rigorous, data-rich research efforts. The pre-trained RouteE models are established using NREL’s Future Automotive Systems Technology Simulator paired with approximately 1 million miles of drive cycle data from the Transportation Secure Data Center, resulting in energy consumption behavior estimates over a representative sample of driving conditions for the United States. Validations have been performed using on-road fuel consumption data for conventional and electrified vehicle powertrains. Transferring the results of the on-road validation to a larger set of real-world origin-destination pairs, it is estimated that implementing the present methodology in a green-routing application would accurately select the route that consumes the least fuel 90% of the time. The novel machine learning techniques used in RouteE make it a flexible and robust tool for a variety of transportation applications.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Small Hydropower Plant Response Improvement Using Energy Storage

Small hydropower plants contribute significantly to global power generation. However, due to limited storage, these can have low ramping capacity and poor load-frequency regulation. These can prevent small hydropower plants from being used in isolated grids as backup power to critical loads and rural areas. In this paper, a control architecture for frequency control is proposed that facilitates the use of energy storage to improve the response of standalone small hydropower plants. The frequency controller generates power commands using proportional control on frequency deviation and a user-defined power setpoint, and this incorporates with a current controller that facilitates power injection. The distinctive feature of the controller design is that it takes into account frequency thresholds and storage constraints to enable power injections, and it provides the capability to perform automated and manual recharge of the energy storage. Simulations using detailed nonlinear models demonstrate the improved performance of the hydropower with energy storage using the proposed controller. Automated and manual recharge and the impact of the energy storage constraints on the performance of the hydropower plant are demonstrated.

13 HYDRO ENERGY↗
DOE OSTI2021

A review of tidal energy—Resource, feedbacks, and environmental interactions

The ocean contains a variety of renewable energy resources, little of which has been exploited. Here, we review both tidal range and tidal stream energy, with a focus on the resource, feedbacks, and environmental interactions. The review covers a wide range of timescales of relevance to tidal energy, from fortnightly (spring-neap) and semi-diurnal variability, down to array, and device-scale turbulence. When simulating the regional tidal energy resource, and to assess environmental impacts, it is necessary to account for feedbacks between the tidal array and the resource itself. We critically review various methods for simulating energy extraction, from insights gained through theoretical studies of “tidal fences” in idealized channels, to realistic three-dimensional model studies with complex geometry and arrays of turbines represented by momentum sinks and additional turbulence due to the presence of rotors and support structures. We discuss how variability can be reduced by developing multiple (aggregated) sites with a consideration of the enhanced phase diversity offered by exploiting less energetic tidal currents. This leads to future research questions that have not yet been explored in depth at first-generation tidal sites in relatively sheltered channels (e.g., the interaction of waves with currents). Such enhanced understanding of real sea conditions, including the effects of wind and waves, leads to our other identified primary future research direction—reduced uncertainties in turbulence predictions, including the development of realistic models that simulate the interaction between ambient turbulence and the turbulence resulting from multiple wakes, and changes to system-wide hydrodynamics, water quality, and sedimentation.

Environmental impacts↗
DOE OSTI2020

Towards isolated attosecond electron bunches using ultrashort-pulse laser-solid interactions

We investigate MeV-level attosecond electron bunches from ultrashort-pulse laser-solid interactions through similarities between experimental and simulated electron energy spectra. We show measurements of the bunch duration and temporal structure from particle-in-cell simulations. The experimental observation of such bunches favors specular reflection direction when focusing the laser pulse onto a subwavelength boundary of thick overdense plasmas at grazing incidence. Particle-in-cell simulation further reveals that the attosecond duration is a result of ultra-thin (~tenth of a micron) gaps of zero electromagnetic energy density in the modulated reflected radiation, while the bunching (locally peaked electron concentration) comes from the highly-directional electron angular distribution acquired by the electrons in a grazing incidence setup. To isolate a single electron bunch, we perform simulations using 1-cycle laser pulses and analyze the effect of carrier-envelop phase with particle tracking. The duration of the electron bunch can be further decreased by increasing the laser intensity and the focal spot size, while its direction can be changed by tuning the preplasma density gradient.

High-harmonic generation↗
DOE OSTI2024

Numerical challenges for energy conservation in N -body simulations of collapsing self-interacting dark matter halos

Dark matter (DM) halos can be subject to gravothermal collapse if the DM is not collisionless, but engaged in strong self-interactions instead. When the scattering is able to efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. However, a comparison with observations requires accurate predictions. One method to obtain such predictions is to use N-body simulations. Collapsed halos are extreme systems that pose severe challenges when applying state-of-the-art codes to model self-interacting dark matter (SIDM). In this work, we investigate the root of such problems, with a focus on energy non-conservation. Moreover, we discuss possible strategies to avoid them. We ran N-body simulations, both with and without SIDM, of an isolated DM-only halo and we adjusted the numerical parameters to check the accuracy of the simulation. We find that not only the numerical scheme for SIDM can lead to energy non-conservation, but also the modelling of gravitational interaction and the time integration are problematic. The main issues we find are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; and (c) SIDM velocity kicks breaking the time symmetry. Tuning the parameters of the simulation to achieve a high level of accuracy allows us to conserve energy not only at early stages of the evolution, but also later on. However, the cost of the simulations becomes prohibitively large as a result. Some of the problems that make the simulations of the gravothermal collapse phase inaccurate can be overcome by choosing appropriate numerical schemes. However, other issues still pose a challenge. Our findings motivate further works on addressing the challenges in simulating strong DM self-interactions.

dark matter↗
DOE OSTI2024

Structure preserving and energy dissipative contact approaches for implicit dynamics

In this work, several structure preserving and energy dissipative contact approaches are proposed and evaluated. The time integration schemes considered are general with regard to the version of constraint type, but here the emphasis was on mortar contact. The proposed mortar contact approach conserves both linear and angular momentum for mortar contact in a novel way. The proposed time integration scheme can conserve energy or provide strict contact dissipation. In addition, the proposed scheme enforces both gap constraints and gap velocity constraints (i.e., persistency). Using a midstep time integrator often causes energy dissipation during initial impact, here this energy can be recovered in a novel way. Enforcing the gap velocity constraint mitigates the contact chatter of the contact pressure and nodes in many problems. Whereas some approaches enforce the persistency condition and gap constraints simultaneously during the solution of the equations of motion (EOM) requiring multipliers for both constraints included in the equation set, here the gap constraint is solved through the equations of motion and the persistency condition is satisfied in the time integration scheme by the velocity update after the equations of motion. It is shown that this approach is strictly dissipative in that a plastic contact condition can be achieved. Analogous to a coefficient of restitution for rigid bodies, any dissipated energy can then be returned upon release if energy conservation is desired. Structure preserving methods are good for long-time dynamics simulations and energy conserving and strictly dissipative methods can overcome stability issues associated with standard time integration algorithms such as the Newmark method.

42 ENGINEERING↗
DOE OSTI2021

Shallow-circuit variational quantum eigensolver based on symmetry-inspired Hilbert space partitioning for quantum chemical calculations

Development of resource-friendly quantum algorithms remains highly desirable for noisy intermediate-scale quantum computing. Based on the variational quantum eigensolver (VQE) with unitary coupled-cluster Ansatz, we demonstrate that partitioning of the Hilbert space made possible by the point-group symmetry of the molecular systems greatly reduces the number of variational operators by confining the variational search within a subspace. In addition, we found that instead of including all subterms for each excitation operator, a single-term representation suffices to reach required accuracy for various molecules tested, resulting in an additional shortening of the quantum circuit by a factor of 4–8. With these strategies, VQE calculations on a noise-free quantum simulator achieve energies within a few meVs of those obtained with the full unitary coupled-cluster Ansatz with single and double excitations for the H 4 -square, H 4 -chain, and H 6 -hexagon molecules, while the number of cnot gates, a measure of the quantum-circuit depth, is reduced by a factor of as large as 35. Furthermore, we introduced an efficient “score” parameter to rank the excitation operators, so that the operators causing larger energy reduction can be applied first. Using the H 4 square and H 4 chain as examples, We demonstrated on noisy quantum simulators that the first few variational operators can bring the energy within the chemical accuracy, while additional operators do not improve the energy since the accumulative noise outweighs the gain from the expansion of the variational Ansatz.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

LINAC Longitudinal Simulation and Measurement of Output Energy

The Fermilab Linac, a pivotal and historic accelerator at Fermilab, is crucial to the laboratory's operations, supplying a 400 MeV beam to various acceleration facilities. Due to daily variations in the Linac's output energy, precise monitoring and machine tuning are important to ensure the exiting energy meets the Booster's acceptance criteria. To address this, Beam Position Monitors (BPMs) are employed to assess and adjust the beam's energy, providing essential data on both the transverse position and longitudinal phase of the beam. We also aim to regulate the longitudinal phase profile to match our simulation predictions. By utilizing a Python-based simulation and beam data from three BPMs, we can determine the optimal phasing correction required for the final RF stage to achieve the desired energy. Currently, the calculations of longitudinal phase profile are based on simulations, so additional research is needed to verify if the predicted longitudinal phase distribution aligns with actual real-world data.

Safaryan, Milena↗
DOE OSTI2023

An interplay between a hydrogen atmosphere and dislocation characteristics in BCC Fe from time-averaged molecular dynamics

The interplay between hydrogen and dislocations (e.g., core and elastic energies, and dislocation–dislocation interactions) has implications on hydrogen embrittlement but is poorly understood. Continuum models of hydrogen enhanced local plasticity have not considered the effect of hydrogen on dislocation core energies. Energy minimization atomistic simulations can only resolve dislocation core energies in hydrogen-free systems because hydrogen motion is omitted so hydrogen atmosphere formation can’t occur. Additionally, previous studies focused more on face-centered-cubic than body-centered-cubic metals. Discrete dislocation dynamics studies of hydrogen–dislocation interactions assume isotropic elasticity, but the validity of this assumption isn’t understood. Here, we perform time-averaged molecular dynamics simulations to study the effect of hydrogen on dislocation energies in body-centered-cubic iron for several dislocation character angles. We see atmosphere formation and highly converged dislocation energies. We find that hydrogen reduces dislocation core energies but can increase or decrease elastic energies of isolated dislocations and dislocation–dislocation interaction energies depending on character angle. We also find that isotropic elasticity can be well fitted to dislocation energies obtained from simulations if the isotropic elastic constants are not constrained to their anisotropic counterparts. These results are relevant to ongoing efforts in understanding hydrogen embrittlement and provide a foundation for future work in this field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗