Search NASA⌕ Search

SEARCH · Search NASA

Results for “loading simulations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

OpenStudio®-MCP [SWR-26-035]

OpenStudio®-MCP is a Model Context Protocol (MCP) server that lets AI assistants perform building energy modeling through natural language. Rather than requiring users to learn the OpenStudio® SDK, EnergyPlus® scripting, or Ruby/Python automation, the server translates conversational requests into sequences of tool calls that create models, design HVAC systems, run simulations, and extract results — all within a single chat session. The server's 124 tools are organized into a skills architecture where each skill encapsulates a domain of building energy modeling (envelope, HVAC, loads, weather, simulation, results) behind typed, LLM-friendly interfaces. High-leverage operations like applying ASHRAE 90.1 baseline systems or generating standards-compliant typical buildings are exposed as single tool calls that internally wire dozens of OpenStudio® objects. Bundled measures from ComStock™ and Openstudio® -common-measures-gem are wrapped with dedicated tools and typed arguments rather than exposed through a generic measure interface, so AI models get consistent, error-resistant recipes without needing to discover measure arguments at runtime. A key design decision is structured results extraction: six SQL-based tools return surgical ~300–1,000 token responses (end-use breakdowns, envelope summaries, HVAC sizing, timeseries data) instead of requiring the AI to parse ~100K-token raw HTML reports, making iterative design exploration practical within context window limits. The codebase is designed as a reference implementation — explicit, well-commented, and modular — so that other simulation engines (EnergyPlus® standalone, TRNSYS, DOE-2) can use it as a template for building their own MCP servers.

Ball, Brian [National Laboratory of the Rockies (N↗
DOE OSTI2025

A VSC-HVDC-Assisted Black-Start Strategy in Bulk Power Systems a Case Study in San Diego

With the worldwide growth in deploying high-voltage direct current (HVDC) transmission systems, their ability to facilitate black-start (BS) restoration has been a research topic of interest. In this context, voltage source converter (VSC)-HVDC is regarded as a BS resource, and this paper proposes a VSC-HVDC-assisted parallel BS restoration strategy in bulk power systems. The proposed strategy consists of two stages: 1) determination of the VSC and generator startup sequence and 2) load restoration simulation. In the first stage, the entire blackout system is sectionalized into multiple subsystems. Each subsystem includes a VSC-HVDC station or traditional BS unit, it independently determines its generator startup timeline and the energization timelines for buses and lines. The second stage involves load restoration, conceptualized as a modified unit commitment problem, with the timelines established in the first stage work as critical inputs. The proposed BS restoration strategy is tested on the San Diego power system to simulate the 2011 Southwest blackout. The simulation results validate the effectiveness of using VSC-HVDC links as a BS resource which not only speeds up the restoration process but also reduces both energy and economic losses.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Modeling of Induction Motors and Variable Speed Drives for Multi-Domain System Simulations Using Modelica and the OpenIPSL Library

This paper introduces an innovative method for characterizing, implementing, and validating both three-phase and single-phase induction motor models, accompanied by a variable speed drive model. The primary goal is to investigate interactions between the electrical power grid and other dynamic domains (e.g., thermofluidic) that impact motor/load drive behavior. Our approach involves establishing a mechanical interface based on a physically meaningful equation linking motor torque/speed to the electrical model in the phasor domain. This allows seamless integration of diverse domain subsystems into a unified multi-domain model using Modelica v4.0.0 and the OpenIPSL library v3.0.1, overcoming co-simulation limitations. The proposed model, which requires only one Modelica-compliant tool for simulation, introduces additional dynamics through the mechanical interface, enabling explicit simulation of load disturbances based on constitutive physics. This deepens our understanding of dynamic interactions between the electrical power domain and other subsystems connected through the motor. We detail the modeled components using mathematical equations and textual descriptions, emphasizing the Modelica modeling approach. Simulation examples validate the implementation, demonstrating the multi-domain modeling capabilities of the newly developed components.

Computer Science↗
DOE OSTI2023

Difference in load predictions obtained with effective turbulence vs. a dynamic wake meandering modeling approach

Abstract. According to the international standard for wind turbine design, the effects of wind turbine wakes on structural loads can be considered in two ways: (1) by augmenting the ambient turbulence levels with the effective turbulence model (EFF) and then calculating the resulting loads and (2) by performing dynamic wake meandering (DWM) simulations, which compute wake effects and loads for all turbines on a farm at once. There is no definitive answer in scientific literature as to the consequences of choosing one model over the other, but the two approaches are unarguably very different. The work presented here expounds on these differences and investigates to what extent they affect the simulated structural loads. We consider an idealized 4×4 rectangular array of National Renewable Energy Laboratory 5 MW wind turbines with a spacing of 5 by 8 rotor diameters and three wind speed scenarios at high ambient turbulence. Load simulations are performed in OpenFAST with EFF and in FAST.Farm with the DWM implementation. We compare ambient turbulence, wind farm turbulence, and loads between both approaches. When omnidirectional results are compared, EFF wind farm turbulence intensity is consistently higher by 0.2 % (above-rated wind speed) to 2.7 % (below-rated wind speed). However, for certain wind directions, the EFF turbulence can be lower than FAST.Farm by almost 9 %. Wind speeds within the farm were found to differ by up to 3 m s−1 due to the lack of wake deficits in the EFF approach, leading to longer tails toward low values in the FAST.Farm mean load distributions. Consistent with the turbulence results, the median EFF load standard deviations are also consistently higher, by a maximum of 20 % and 17 % for blade-root out-of-plane and tower-base fore-aft moments, respectively.

17 WIND ENERGY↗
DOE OSTI2023

Uncertainty Error Modeling for Non-Linear State Estimation With Unsynchronized SCADA and µPMU Measurements

Distribution systems of the future smart grid require enhancements to the reliability of distribution system state estimation (DSSE) in the face of low measurement redundancy, unsynchronized measurements, and dynamic load profiles. Micro phasor measurement units (µPMUs) facilitate co-synchronized measurements with high granularity, albeit at an often prohibitively expensive installation cost. Supervisory control and data acquisition (SCADA) measurements can supplement µPMU data, although they are received at a slower sampling rate. Further complicating matters is the uncertainty associated with load dynamics and unsynchronized measurements–not only are the SCADA and µPMU measurements not synchronized with each other, but the SCADA measurements themselves are received at different time intervals with respect to one another. This paper proposes a non-linear state estimation framework which models dynamic load uncertainty error by updating the variances of the unsynchronized measurements, leading to a time-varying system of weights in the weighted least squares state estimator. Case studies are performed on the 33-Bus Distribution System in MATPOWER, using Ornstein–Uhlenbeck stochastic processes to simulate dynamic load conditions.

Cooper, Austin↗
DOE OSTI2023

ORBIT simulations of fast ion power loads on the wall of the Divertor Tokamak Test

Abstract Neutral beam injection in tokamaks produces a population of fast ions, which interact with 3D magnetic fields in a variety of ways, often resulting in energetic particle losses in very short times to the wall. Careful design of neutral beams and active control of error fields helps to keep these losses to a minimum. Nevertheless, past experience in tokamaks in the 1980s and detailed simulations for future machines, such as International Thermonuclear Experimental Reactor, suggest that very localized fast ion losses (‘hot spots’) can be present, even if the overall losses are low. In this paper, we discuss this issue in the Divertor Tokamak Test (DTT) project, and we show that in the standard single-null full-power scenario of the DTT, fast ions produce two hot spots, corresponding approximately to the beam injection and exit toroidal angles: the former being mainly due to prompt losses/passing particles, while the latter is due to trapped ions. However, the maximum power load in these spots is of the order ∼ 63 kW m −2 , below the tolerance of plasma-facing components of the machine.

Physics↗
DOE OSTI2026

Comparison between Alcator C-Mod ICRF experiments and 3D full wave simulations

Reliable modeling of ion cyclotron range of frequencies (ICRF) antenna performance is essential for interpreting present experiments and guiding the design of future reactors. In this work, a 3D model of the Alcator C-Mod field-aligned antenna is implemented in the Petra-M finite-element framework [S. Shiraiwa et al 2023 Nucl. Fusion 63 026024] and benchmarked against experimental results. Four experimental cases are examined. First, the simulated rectified sheath potentials on the antenna limiters are compared with measurements from a power tapering experiment. Second, proof-of-principle far-field sheath simulations are performed. In scenarios with low single-pass absorption, simulations predict enhanced sheath potentials on a distant poloidal limiter in the far field of the antenna, consistent with experimental observations. Third, the simulated antenna loading during edge localized modes follows experimental trends and appears to be dominated by the density gradient at the pedestal. Finally, Petra-M predicts the unintended excitation of high-$k$ ∥ modes during monopole phasing operation, in agreement with experimental evidence of poor wave coupling and accessibility to the plasma core. Overall, reasonable agreement is found between the simulations and experiments. At the same time, areas of imperfect agreement are identified. These provide important guidance on the limits of the current state-of-the-art modeling, which should be kept in mind when using it as a predictive tool for future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Porous Flow Modeling of Axial Gas Redistribution in Fragmented LWR Fuel Rods using MOOSE

Understanding how gas axially redistributes within fragmented fuel pellets is crucial for predicting the behavior of Light Water Reactor (LWR) fuel rods, particularly during transient and accidental scenarios. The time scale of this phenomenon plays a fundamental role in determining the progression and hazard of a Loss Of Coolant Accident (LOCA), especially when high burn-up fuel in a severe state of fragmentation is involved. Here, this study presents a Computational Fluid Dynamics (CFD) model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) to predict the time-scale of plenum depressurization in Light-Water Reactor (LWR) fuel rods driven by axial gas transport through fragmented pellets. The model examines the effects of incorporating non-linearities in the friction term by comparing the results with experimental data. The experiment employed surrogate fuel rods containing pellets subjected to mechanical and/or thermal loadings to simulate various severity of cracking, and aimed at studying the influence of fuel conditions on axial gas redistribution. The results of this analysis indicate that under certain flow regime conditions - determined by the value of an equivalent Reynolds number - accounting for the non-linear friction term in Navier-Stokes equations guarantees better predictions for the time-scale of plenum depressurization. Also, the model enabled the simulation of the pressure decay by assigning distinct permeability values to each pellet instead of a single uniform value. Multiple simulations were run across all possible pellet position combinations, having each pellet assigned with values of permeability extracted from the experimental data. This allows to quantify the impact of the considering various non-uniform distributions of permeability on the dynamics of axial gas redistribution. The present work findings enhance the understanding of axial gas transport, and provide valuable insights for the integration of a model for predicting the axial gas redistribution during a LOCA scenario into the BISON fuel performance code.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Control methods for mitigating flow oscillations in a supercritical CO 2 recompression closed Brayton cycle

A dynamic model of a 10 MWe supercritical CO 2 (sCO 2 ) recompression closed Brayton power generation cycle is used to investigate control methods for mitigating oscillations in flow conditions and power demand load, especially as the load ramp rate is increased and the load setpoint is more closely tracked. The focus is on control of the main compressor inlet temperature (MCIT) since oscillations are prompted by the strong nonlinearities in sCO 2 fluid properties near the critical point and by changes in system inventory for efficient load management. In this study, the oscillatory behavior during load ramping is first shown and then two potential control solutions for substantially reducing oscillations are given: 1) the use of the main compressor inlet guide vanes (IGV) and 2) bypassing a portion of the sCO 2 around the cooler to maintain the MCIT. These control methods reduce the effect of cycle feedback from interactions between the MCIT and load. When implementing either control method, a load setpoint ramp of 7.5%/min was maintainable while achieving closer load setpoint tracking with limited oscillations. Finally, the impact of sensor noise on the control was examined since derivative action was used in the IGV and cooler sCO 2 bypass control. Finally, while some effect was evident, signal noise was not problematic to the control.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Effects of upstream rotor tilt on a downstream floating wind turbine

This work quantifies the relationship between the design tilt of an upstream rotor and the structural response of a spar floating offshore wind turbine located in its wake. Three wind speed scenarios are considered: below, at, and above rated operation. The inflow is generated with the Mann model and the wake and loads are simulated with FAST.Farm. As the upstream rotor tilt goes from 0° to 6°, we find that the mean wake is displaced upward by more than 0.2 rotor diameters (D) by a downstream distance of 7 D. The vertical velocities increase by up to 35 cm/s in the center of the wake. As a result, the downstream rotor is partially waked and experiences a rotated velocity vector. With a higher upstream rotor tilt, the velocities and moments on the downstream turbine increase their mean axial value and their lateral and vertical standard deviation. These changes affect the blade and tower loading and the floater motion primarily in the out-of-plane direction: the damage-equivalent loads for the tower pitch moment and blade-root moment increase by up to 10% because of higher variability at the first mode for the tower and at one blade passing frequency for the blade root. Lesser effects are observed for the roll moments and for floater sway and heave. When the joint effect of rotor tilt and platform motion is considered, the load response on the downstream system is amplified primarily for tower pitch and blade out-of-plane moment.

17 WIND ENERGY↗
DOE OSTI2023

Molecular Dynamics Simulation of Pore-Size Effects on Gas Adsorption Kinetics in Zeolites

Strong gas-mineral interactions or slow adsorption kinetics require a molecular-level understanding of both adsorption and diffusion for these interactions to be properly described in transport models. In this combined molecular simulation and experimental study, noble gas adsorption and mobility is investigated in two naturally abundant zeolites whose pores are similar in size (clinoptilolite) and greater than (mordenite) the gas diameters. Simulated adsorption isotherms obtained from grand canonical Monte Carlo simulations indicate that both zeolites can accommodate even the largest gas (Rn). However, gas mobility in clinoptilolite is significantly hindered at pore-limiting window sites, as seen from molecular dynamics simulations in both bulk and slab zeolite models. Further, experimental gas adsorption isotherms for clinoptilolite confirm the presence of a kinetic barrier to Xe uptake, resulting in the unusual property of reverse Kr/Xe selectivity. Finally, a kinetic model is used to fit the simulated gas loading profiles, allowing a comparison of trends in gas diffusivity in the zeolite pores.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Decarbonizing all-electric communities via carbon-responsive control of behind-the-meter resources

The progression of electrification in the building and transportation sectors brings new opportunities for energy decarbonization. With higher dependence on the grid power supply, the variation of the grid carbon emission intensity can be utilized to reduce the carbon emissions from the two sectors. Existing coordinated control methods for buildings with distributed energy resources (DERs) either consider electricity price or renewable energy generation as the input signal, or adopt optimization in the decision-making, which is difficult to implement in the real-world environment. This paper aims to propose and validate an easy-to-deploy rule-based carbon responsive control framework that facilitates coordination between all-electric buildings and electric vehicles (EVs). The signals of the grid carbon emission intensity and the local photovoltaics (PV) generation are used for shifting the controllable loads. Extensive simulations were conducted using a model of an all-electric mixed-use community in a cold climate to validate the control performance with metrics such as emissions, energy consumption, peak demand, and EV end-of-day state-of-charge (SOC). Our study identifies that 4.5% to 27.1% of annual emission reduction can be achieved with limited impact on energy costs, peak demand, and thermal comfort. Additionally, up to 32.7% of EV emission reduction can be obtained if the EV owners reduce the target SOC by less than 21.2%.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Tailoring ductile-phase toughened tungsten hierarchical microstructures for plasma-facing materials

Biological materials, such as bones, nacre, etc., have been found to show attractive and unique combinations of stiffness, strength and fracture toughness. Research has been conducted in various engineering areas to mimic the naturally hierarchical microstructures or nanostructures of these materials to produce materials with optimum stiffness, high strength and fracture toughness for their intended structural applications. In the same objective, this work applies a multiscale microstructural approach recently developed to investigate the deformation and fracture behavior of ductile phase toughened tungsten (W) materials such as tungsten-nickel/iron (W-Ni/Fe) composites that possess lamellar-like and hierarchical “brick-and-mortar” (BAM) microstructures. First, the approach is used to simulate tensile loading of W-Ni/Fe specimens cut out from hot-rolled W-Ni/Fe plates which possess lamellar-like microstructures. The finite element model of the gage section of the specimen consists of a W-Ni/Fe dual-phase microstructural domain in which the constitutive behaviors of W and Ni/Fe phases are described by an elastic-plastic damage model. The predicted material stress-strain response and crack pattern development as a function of loading are then compared to the corresponding experimental results to determine the constitutive model parameters for W and Ni/Fe. Subsequently, the model parameters are used to analyze W-Ni/Fe BAM microstructures that are artificially created to investigate the effects of microstructural features and morphologies on the composite stress-strain response, damage, and fracture patterns. The modeling shows that regular BAM microstructures that experience bridging mechanism combined with crack penetration across W-phase regions exhibit significantly higher strengths than the rather random lamellar-like microstructures. More significantly, however adjusting the brick’s length-to-height ratio, the BAM microstructure can be designed to allow a more distributed damage zone which leads to increased strength, ductility and fracture energy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exploring the parameter space of MagLIF implosions using similarity scaling. I. Theoretical framework

Magneto-inertial fusion concepts, such as the magnetized liner inertial fusion (MagLIF) platform, constitute an alternative path for achieving ignition and significant fusion yields in the laboratory. The space of experimental input parameters defining a MagLIF load is highly multi-dimensional, and the implosion itself is a complex event involving many physical processes. In the first paper of this series, we develop a simplified analytical model that identifies the main physical processes at play during a MagLIF implosion. Using non-dimensional analysis, we determine the most important dimensionless parameters characterizing MagLIF implosions and provide estimates of such parameters using typical fielded or experimentally observed quantities for MagLIF. Here, we then show that MagLIF loads can be “incompletely” similarity scaled, meaning that the experimental input parameters of MagLIF can be varied such that many (but not all) of the dimensionless quantities are conserved. Based on similarity-scaling arguments, we can explore the parameter space of MagLIF loads and estimate the performance of the scaled loads. Then, in the follow-up papers of this series, we test the similarity-scaling theory for MagLIF loads against simulations for two different scaling “vectors,” which include current scaling and rise-time scaling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Commissioning plan for LCLS II 2 K cold box

Commissioning of the LCLS II cryogenic plant has been planned with respect to electrical and mechanical systems checkout, equipment cleanup, cool-down, and individual performance evaluation. This plant consists of warm and cold storage, helium compressors, 4.5 K and 2 K cold boxes, and all auxiliaries typical for helium refrigeration processes. The 2 K cold box, consisting of five cold compressors in series to produce the flow and pressures necessary for 2K operation, will be the last equipment installed and, along with an additional 10,000 L helium dewar, will allow the entire plant to be commissioned as a single entity. The pre-commissioning stages, commissioning plan for 2 K cold box without connecting to the LCLS II LINAC, additional test equipment needed to simulate the loads, all process studies and associated results are described herein.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mitigation of magneto-Rayleigh-Taylor instability growth in a triple-nozzle, neutron-producing gas-puff Z pinch

The gas-puff Z-pinch is a well-known source of x-rays and/or neutrons, but it is highly susceptible to the magneto-Rayleigh-Taylor instability (MRTI). Approaches to MRTI mitigation include density profile tailoring, in which nozzles are added or modified to alter the acceleration trajectory, and axial pre-magnetization, in which perturbations are smoothed out via magnetic field line tension. In this work, we present two-dimensional magnetohydrodynamic simulations of loads driven by an 850 kA,160 ns driver that suggest these mitigation strategies can be additive. The initial axial magnetic field, B z0 , to stabilize a 2.5-cm-radius Ne gas liner imploding onto an on-axis deuterium target can be reduced from 0.7 T to 0.3 T by adding a second liner with a radius of 1.25 cm. Because MRTI mitigation tends to increasingly lower yield with higher B z0 , the use of a lower field is advantageous. Here, we predict a reduction in yield penalty from > 100× with the single liner to < 10× with a double liner. A premagnetized, triple nozzle gas puff could therefore be an attractive source for intense neutrons or other fusion applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗