Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 1,189 records · Page 66

First-principles investigation of high capacity, rechargeable CF x cathode batteries based on graphdiyne and “holey” graphene carbon allotropes

Batteries composed of CF x cathodes have high theoretical specific capacities (>860 mA h g -1 ). Attempts at realizing such batteries coupled with Li anodes have failed to deliver on this promise, however, due to a discharge voltage plateau below the theoretical maximum lowering the realized energy density and difficulties with recharging the system. Here, in this study, we use first-principles calculations to investigate novel carbon allotropes for these battery systems: graphdiyne and “holey” graphene. We first identify stable flourination structures and calculate their band gaps. We demonstrate that the holes in these carbon allotropes can induce the formation of an amorphous LiF network within the carbon and that this formation may, in fact, be kinetically favored. For structures where amorphous LiF forms within the carbon, we predict it is easier to recharge and higher discharge voltages can be achieved. If the LiF forms outside the carbon product, however, it will be crystalline in form and lead to lower discharge voltages and more difficulty in recharging the systems. Finally, we simulate XPS spectra of representative cases, demonstrating an experimental pathway for determining the reaction pathway of these systems. Our work suggests CF x allotropes with holes in them as potential targets for high capacity, rechargeable cathodes for Li batteries, provided they lead to the formation of amorphous LiF within the C structure.

25 ENERGY STORAGE↗

Improving Resilience of Bus Bunching Holding Strategy through a Rolling Horizon Approach

Providing public transportation with quality service is critical to attracting more passengers to the system. However, high-demand routes are prone to the so-called bus bunching -- a tendency of buses to group as a consequence of variations in travel times and demands. Bus holding is applied to overcome this effect. In this study, we present a novel method for bus holding in which the control law is based only on the buses' position using a computationally efficient rolling horizon approach. The method uses similar inputs as linear control approaches while not increasing significantly the computational time. However, the method overcomes key a weakness of the linear control approach thanks to the explicit constraint handling that always ensures the control action effectiveness. Simulation experiments in a validation case and a model-specific for a bus rapid transit line in Curitiba, Brazil showed a reduced holding time and improved resilience, delivering more than 20% reduction in delay time accounting for the on-board and station delays.

33 ADVANCED PROPULSION SYSTEMS↗

Electron temperature gradient instability and transport analysis in NSTX and NSTX-U plasmas

Extensive linear and nonlinear simulations to study electron temperature gradient (ETG) stability and thermal transport in National Spherical Torus Experiment (NSTX) and NSTX-U plasmas were performed using the fully electromagnetic gyrokinetic code CGYRO. Linear simulations were performed to determine ETG thresholds in different discharges, showing that ETG modes in spherical tokamaks can present different scalings compared to conventional aspect-ratio tokamaks. Nonlinear gyrokinetic simulations were conducted for selected cases to calculate electron thermal transport and compare to experimental values. Results are also compared with those of ETG modes in the multi-mode model and the Trapped-Gyro-Landau-Fluid reduced model codes, to better understand their applicability in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of Surfactants on Noble Metal Liquid-Gas Interface Mass Transfer in Molten Salt Reactors

Noble metals, a group of insoluble fission products (FPs), can accumulate in molten salt reactors (MSRs) and influence system performance and safety through deposition. Upon generation in fission, noble metals are transported by circulating salt and may either deposit on structural surfaces or enter the circulating gas bubble interfaces. While sparging was primarily designed for xenon and krypton removal, the mass transfer of noble metals into bubble interfaces highlights a potential pathway for their transport and partial extraction. The properties of the liquid-gas interface therefore play an important role in noble metal deposition and removal behaviors. Observations from earlier studies suggest that noble metals accumulated at the bubble interface can act as surfactants, making the bubble interface partially immobile. This behavior reduces the mass transfer capacity of circulating bubbles and affects both the overall distribution of noble metals in the loop and their removal efficiency. To investigate the influence of surfactants on noble metal transport, deposition, and extraction behaviors in MSRs, a previously developed and validated species transport model incorporating noble metal interphase mass transfer was used to simulate these behaviors. Two designed cases were studied to predict the distribution of 132Te along the Molten Salt Reactor Experiment (MSRE) loop under different bubble interface conditions. The results show that the mass transfer coefficients at the liquid-gas interface, determined by bubble interface characteristics, significantly influence noble metal distribution within the reactor loop. These findings emphasize the importance of continuously introducing fresh circulating bubbles into the molten salt to improve the removal efficiency of insoluble FPs from the MSR primary loop.

MSR↗

Analysis of Dynamic Behavior of the Target and the Deuterium-Tritium Ice in Magnetic-Field Assisted Implosions

The use of strong magnetic fields to augment the output energy of inertial confinement fusion experiments at the National Ignition Facility is of high interest. It offers the potential of reducing electron thermal conduction and increasing hot-spot alpha heating with little to no change in hohlraum behavior. In these magnetically assisted ignition experiments, the ultimate goal is to add a B-field in the form of a pulse ranging from 25 to 60 T to a high-performing hohlraum implosion several microseconds before impingement of the laser beams. This requires eliminating metallic components in the target and replacing them with electrically nonconducting materials. However, the strong eddy currents generated by the rapidly increasing high B-field, which were calculated to be as high as 2000 K, can heat the hohlraum. Here in this paper, we examine the transient effects of this rapid temperature change on the behavior of the target as well as the fuel layer composed typically of deuterium and tritium. Using simulations and calculations for limiting case scenarios, we find that the effect of the heating is not restrictive toward the performance of the target or the quality of the deuterium-tritium ice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An integrated manifold learning approach for high-dimensional data feature extractions and its applications to online process monitoring of additive manufacturing

As an effective dimension reduction and feature extraction technique, manifold learning has been successfully applied to high-dimensional data analysis. With the rapid development of sensor technology, a large amount of high-dimensional data such as image streams can be easily available. Thus, a promising application of manifold learning is in the field of sensor signal analysis, particular for the applications of online process monitoring and control using high-dimensional data. The objective of this study is to develop a manifold learning-based feature extraction method for process monitoring of Additive Manufacturing (AM) using online sensor data. Due to the non-parametric nature of most existing manifold learning methods, their performance in terms of computational efficiency, as well as noise resistance has yet to be improved. To address this issue, this study proposes an integrated manifold learning approach termed multi-kernel metric learning embedded isometric feature mapping (MKML-ISOMAP) for dimension reduction and feature extraction of online high-dimensional sensor data such as images. Based on the extracted features with the utilization of supervised classification and regression methods, an online process monitoring methodology for AM is implemented to identify the actual process quality status. Finally, in the numerical simulation and real-world case studies, the proposed method demonstrates excellent performance in both prediction accuracy and computational efficiency.

36 MATERIALS SCIENCE↗

Neural network based fast prediction of β N limits in HL-2M

Artificial neural networks (NNs) are trained, based on the numerical database, to predict the no-wall and ideal-wall β N limits, due to onset of the n = 1 ( n is the toroidal mode number) ideal external kink instability, for the HL-2M tokamak. The database is constructed by toroidal computations utilizing both the equilibrium code CHEASE (Lütjens et al 1992 Comput. Phys. Commun. 69 287) and the stability code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681). The stability results show that (1) the plasma elongation generally enhances both β N limits, for either positive or negative triangularity plasmas; (2) the effect is more pronounced for positive triangularity plasmas; (3) the computed no-wall β N limit linearly scales with the plasma internal inductance, with the proportionality coefficient ranging between 1 and 5 for HL-2M; (4) the no-wall limit substantially decreases with increasing pressure peaking factor. Furthermore, both the NN model and the convolutional neural network (CNN) model are trained and tested, producing consistent results. The trained NNs predict both the no-wall and ideal-wall limits with as high as 95% accuracy, compared to those directly computed by the stability code. Additional test cases, produced by the Tokamak Simulation Code (Jardin et al 1993 Nucl. Fusion 33 371), also show reasonable performance of the trained NNs, with the relative error being within 10%. The constructed database provides effective references for the future HL-2M operations. The trained NNs can be used as a real-time monitor for disruption prevention in the HL-2M experiments, or serve as part of the integrated modeling tools for ideal kink stability analysis.

Physics↗

Evidence of electron microbunching in laser-driven modulated downramp injection and prospects for beam-driven implementation

Plasma accelerators can generate high-energy, high-brightness electron beams over centimeter-scale distances, offering novel pathways to compact x-ray free-electron lasers. Generating beams pre-bunched at the desired radiation wavelength would significantly enhance longitudinal coherence and reduce saturation length. Plasma density-modulated downramp injection offers an in-situ way to generate such beams with nanometer-scale bunching. Here we report the first experimental evidence of this mechanism in a laser-driven wakefield accelerator, showing that modulated density downramps generate modulated electron energy spectra absent in unmodulated cases. Particle-in-cell simulations reproduce these observations and reveal bunching factors of 0.05 at 0.4 μm, with a compression factor of approximately 7. Building on this demonstration, we propose a beam-driven implementation for FACET-II to generate multi-GeV beams pre-bunched at hundreds of nanometers wavelength, with sub-micrometer emittance, kiloampere peak current, and sub-percent slice energy spread. Two-stage magnetic compression enables tunable bunching from optical to extreme ultraviolet wavelengths while achieving peak currents exceeding 100 kA. Coherent transition radiation calculations confirm diagnostic feasibility. This approach offers a promising path towards compact, high-energy pre-bunched electron sources for advanced photon science applications.

electron microbunching↗

Progress in extending high poloidal beta scenarios on DIII-D towards a steady-state fusion reactor and impact of energetic particles

To prepare for steady-state operation of future fusion reactors (e.g. the International Thermonuclear Experimental Reactor and China Fusion Engineering Test Reactor (CFETR)), experiments on DIII-D have extended the high poloidal beta ( β P ) scenario to reactor-relevant edge safety factor q 95 ~ 6.0, while maintaining a large-radius internal transport barrier (ITB) using negative magnetic shear. Excellent energy confinement quality ( H 98y2 > 1.5) is sustained at high normalized beta ( β N ~ 3.5). This high-performance ITB state with Greenwald density fraction near 100% and q min ≥ 3 is achieved with toroidal plasma rotation V tor ~ 0 at ρ ≥ 0.6. This is a key result for reactors expected to have low V tor . At high β P (≥1.9), large Shafranov shift can stabilize turbulence leading to a high confinement state with a low pedestal and an ITB. At lower β P (<1.9), negative magnetic shear in the plasma core contributes to turbulence suppression and can compensate for reduced Shafranov shift to continue to access a large-radius ITB and excellent confinement with low V tor , consistent with the results of gyrofluid transport simulations. These high- β P cases are characterized by weak/no Alfvén eigenmodes (a.e.) and classical fast-ion transport. At high density, the fast-ion deceleration time decreases and Δ β fast is lower; these reduce a.e. drive. The reverse-shear Alfvén eigenmodes are weaker or stable because the negative magnetic shear region is located at higher radius, away from the peaked fast-ion profile. Resistive wall modes can be a limitation at simultaneous high β N , low internal inductance, and low rotation. Analysis suggests that additional off-axis external current drive could provide a more stable path at reduced q 95 . Based on a DIII-D high- β P plasma with large-radius ITB, two scenarios are proposed for CFETR Q = 5 steady-state operation with ~1 GW fusion power: a lower-$l_i$($l_i$ ~ 0.66) and a higher-$l_i$($l_i$ ~ 0.75) case. Using a Landau closure model, multiple energetic particle (EP) effects on the a.e. stability are analyzed modifying the growth rate of the a.e.s triggered by the neutral-beam-injection EPs and alpha particles, although the stabilizing/destabilizing effect is weak for the cases analyzed. The stabilizing effects of the combined EP species β , energy, and density profile in CFETR need further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Search for interacting galaxy clusters from SDSS DR-17 employing optimized friends-of-friends algorithm and multimessenger tracers

ABSTRACT In the theoretical framework of hierarchical structure formation, galaxy clusters evolve through continuous accretion and mergers of substructures. Cosmological simulations have revealed the best picture of the universe as a 3D filamentary network of dark-matter distribution called the cosmic web. Galaxy clusters are found to form at the nodes of this network and are the regions of high merging activity. Such mergers being highly energetic, contain a wealth of information about the dynamical evolution of structures in the Universe. Observational validation of this scenario needs a colossal effort to identify numerous events from all-sky surveys. Therefore, such efforts are sparse in literature and tend to focus on individual systems. In this work, we present an improved search algorithm for identifying interacting galaxy clusters and have successfully produced a comprehensive list of systems from SDSS DR-17. By proposing a set of physically motivated criteria, we classified these interacting clusters into two broad classes, ‘merging’ and ‘pre-merging/postmerging’ systems. Interestingly, as predicted by simulations, we found that most cases show cluster interaction along the prominent cosmic filaments of galaxy distribution (i.e. the proxy for dark matter filaments), with the most violent ones at their nodes. Moreover, we traced the imprint of interactions through multiband signatures, such as diffuse cluster emissions in radio or X-rays. Although we could not find direct evidence of diffuse emission from connecting filaments and ridges; our catalogue of interacting clusters will ease locating such faintest emissions as data from sensitive telescopes such as eROSITA or SKA, becomes accessible.

Oak, Tejas↗

Operation and Control of Electric Vehicle Charger with Enhanced Dynamic Performance Under Non-Ideal Grid Voltage Condition

This paper presents a three-phase electric vehicle charger connected to the grid, featuring multiple boost converters on the DC side, specifically designed to ensure smooth, oscillation-free power transfer during unsymmetrical voltage sags. Precise control mechanisms are implemented on the boost converter side to regulate both voltage and current on the electric vehicle side, thereby maintaining optimal charging conditions. The control architecture for both the grid-connected and boost converter components is based on the Lyapunov energy function, which is employed to achieve superior dynamic performance and stability. The system's robustness and reliability are demonstrated through its ability to maintain stable operation and efficient power transfer despite fluctuations in grid conditions. Furthermore, the implementation of Lyapunov-based control ensures rapid response and minimal energy loss, enhancing the overall efficiency of the system. To validate the effectiveness of this approach, a comprehensive model of the system was developed and tested using MATLAB/Simulink, with detailed computer simulations conducted across various significant case studies.

DC-DC boost converter↗

A 0.2–2 GHz Time-Interleaved Multistage Switched-Capacitor Delay Element Achieving 2.55–448.6 ns Programmable Delay Range and 330 ns/mm 2 Area Efficiency

Simulation of radar returns, full-duplex systems, and signal repeaters require hundreds of ns of programmable broadband radio frequency (RF) delay in the signal path to simulate large distances in the case of radar returns, for signal cancellation in full-duplex, and for isolation from reflections in signal repeaters. However, programmable broadband RF delay has been limited to ones of ns due to challenges in miniaturization with low loss and low power consumption. In this work, we present a 0.2–2 GHz digitally programmable RF delay element based on a time-interleaved multistage switched-capacitor (TIMS-SC) approach. The proposed approach enables hundreds of ns of broadband RF delay by employing sample time expansion in multiple stages of switched-capacitor storage elements. Further, the delay element was implemented in a 45 nm SOI CMOS process and achieves a 2.55–448.6 ns programmable delay range with < 0.12% delay variation across 1.8 GHz of bandwidth at maximum delay, 2.42 ns programmable delay steps, and 330 ns/mm 2 area efficiency. Through the proposed approach, the device shows minimal delay change across a -40 °C to 85 °C temperature range and < 0.25 dB gain variation across delay settings. The device achieves 26 dB gain, 7.4 dB noise figure, and consumes 74 mW from a 1 V supply with an active area of 1.36 mm 2 .

42 ENGINEERING↗

Proactive Frequency Stability Scheme via Bayesian Filters and Synchrophasors

Underfrequency (UF) load shedding schemes are traditionally implemented in two ways: One approach is based on manual load shedding, with system operators requesting loads to be shed ahead of anticipated stressful operating conditions. Manual load shedding is usually done through phone calls. The second method is automatic load shedding via underfrequency relays. Using static static settings, these schemes can be designed to operate in stages and drop previously identified loads. The main limitation of traditional load shedding schemes is that they are reactive and leave little room for optimized corrective actions. This work presents a proactive and automatic underfrequency load shedding solution for power systems. Measurements are captured via phasor measurement units (PMUs) at relatively low sampling rates of 30 Hz. These measurements are then processed by particle filters who predict the future state of the system's frequency. Based on these predictions excess load is determined and shed. Comparative case studies are performed in simulated environments. Easy-to-implement models, without hard-to-derive parameters, highlight potential aspects for real-life implementation.

Paramo, Gian↗

Modeling the Strategic Behavior of an Active Distribution Network in the ISO Markets

With increasing integration of distributed energy resources (DERs), active distribution networks (ADNs) can actively participate in the electricity markets by dispatching their DERs, which can change the existing electricity market paradigm. It is essential to investigate the strategic behaviors of ADNs and their DER dispatch when they participate in the wholesale market as price-makers. This paper proposes a bi-level optimization model to study the strategic behavior of an ADN in both energy and reserve markets. The optimal scheduling of DERs in the ADN is modeled as the upper level problem and the joint energy and reserve market-clearing of the ISO is modeled as the lower-level problem. The two-level optimization models exchange bidding information and energy/reserve prices with each other. The proposed bi-Ievel optimization problem is converted to a mathematical programming with equilibrium constraints (MPEC) by using Karush-Kuhn Tucker (KKT) conditions and strong duality theory. Further, the MPEC problem is reformulated as a computationally-solvable mixed integer second order cone programming (MISOCP) model. The simulation results on an illustrative case demonstrate the impact of the strategic bidding of the ADN on the day-ahead energy and reserve market prices.

active distribution network↗

Benchmark High-Fidelity EMT Models for Power Grid with PV Plants

In recent times electromagnetic transient (EMT) modeling tools have been identified as one of the most important requirements in replicating, analyzing, and investigating the dynamics of the power grid with photovoltaic (PV) plants. However, there are no benchmark models for power grid with PVs to investigate emerging challenges with higher penetration of PVs (like trips and momentary cessations during faults from a region far away). To this end, in this paper, synthetic benchmark high-fidelity EMT dynamic models of power grid with large-scale PV plants are presented. The models are developed in PSCAD and PSCAD/Fortran. Simulation results for different use cases (events) and scenarios are presented.

Marthi, Phani Ratna Vanamali↗

Ensemble Learning Based Convex Approximation of Three-Phase Power Flow

Though the convex optimization has been widely used in power systems, it still cannot guarantee to yield a tight (accurate) solution to some problems. To mitigate this issue, this paper proposes an ensemble learning based convex approximation for alternating current (AC) power flow equations that differs from the existing convex relaxations. The proposed approach is based on three-phase quadratic power flow equations in rectangular coordinates. To develop this data-driven convex approximation of power flows, the polynomial regression (PR) is first deployed as a basic learner to fit convex relationships between the independent and dependent variables. Then, ensemble learning algorithms such as gradient boosting (GB) and bagging are introduced to combine learners to boost model performance. Based on the learned convex approximation of power flow, optimal power flow (OPF) is formulated as a convex quadratic programming problem. The simulation results on IEEE standard cases of both balanced and unbalanced systems show that, in the context of solving OPF, the proposed data-driven convex approximation outperforms the conventional semi-definite programming (SDP) relaxation in both accuracy and computational efficiency, especially in the cases that the conventional SDP relaxation fails

Convex approximation↗

Software for Transformative Remedial Action Scheme Tool (TRAST)

The transformative remedial action scheme tool (TRAST) can be applied to improve and validate the power system remedial action scheme (RAS), and further improve the performance of power system operation and control. This tool provides a full suite of advanced functionalities, which are given as follows: 1. Advanced statistical data analysis; 2. OPF-based automated power flow case generation; 3. Customized dynamic simulation in HPC/cloud platform; 4. Machine learning based RAS coefficient prediction; 5. A reliable RAS validation strategy in multiple commercial platforms.

Fan, Xiaoyuan↗

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗