Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated current”

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,063 records · Page 59

Novel Hybrid Current Limiter for Grid-Forming Inverter Control During Unbalanced Faults

Grid-forming (GFM) inverter controls have illustrated many desirable features to enable the bulk integration of renewable resources into the future power grid; however, the performance of GFM inverters during unbalanced faults remains underexplored. This paper proposes a novel current-limiting method for GFM inverters to handle unbalanced fault conditions while providing voltage support to the main grid. The proposed current limiter combines concepts of dynamic virtual-impedance and current-reference saturation limiting, all built in the stationary reference frame, to achieve improved current-limiting performance under unbalanced load/fault conditions. System-wide full-order transient simulations with multiple GFM inverters demonstrate the potential of the method and benchmark the performance against state-of-the-art current limiters. Simulation results portray improved voltage balancing performance of the proposed method compared to both current-reference saturation and virtual-impedance limiting.

controller anti windup↗

Novel Hybrid Current Limiter for Grid-Forming Inverter Control During Unbalanced Faults: Preprint

Grid-forming (GFM) inverter controls have illustrated many desirable features to enable bulk integration of renewable sources into the future power grid. However, the performance of GFM inverters during unbalanced faults remains under-explored. In this paper, a novel current-limiting method is proposed for GFM inverters to handle unbalanced fault conditions while providing voltage support to the main grid. The proposed current limiter combines concepts of dynamic virtual-impedance and current-reference saturation limiting, all built-up in the stationary reference frame to achieve improved current-limiting performance under unbalanced load/fault conditions. System-wide full-order transient simulations with multiple GFM inverters demonstrate the potential of the method and benchmark the performance against state-of-the-art current limiters. Simulation results portray improved voltage balancing performance of the proposed method compared to both current-reference saturation and virtual-impedance limiting.

controller anti windup↗

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comparison of laser models

Progress on the comparison of laser models is reported. Equations for an oscillatory model and a nonoscillatory model for the simulation of iodine laser operation are solved. Reaction rate coefficients used in both models are listed. Currently there are four models for the simulation of iodine laser operation are solved. Reaction rate coefficients used in both models are listed. Currently there are four models for the simulation of an iodine laser. They are: (1) a time dependent model; (2) a quasi-steady state model; (3) a noncompressible model; and (4) a compressible flow laser model. Current research is being directed toward: (1) parameter studies using the compressible flow laser model; (2) development of a two-pass amplifier model; and (3) solving a system of equations describing operation of the high powered iodine MOPA (master oscillator power amplifier).

Heinbockel, John H.↗

Comparison of 1D and 2D particle-in-cell simulations for DC magnetron sputtering discharges

One-dimensional and two-dimensional particle-in-cell (PIC)/Monte Carlo collision simulations are performed for a balanced direct current magnetron sputtering (DCMS) discharge. The plasma parameters obtained by both simulations above the target racetrack share similar features and magnitudes under similar discharge intensities, i.e., similar local discharge current densities. The commonly observed high-energy tail of electron energy distribution in DCMS discharges is reproduced in both simulations. Additionally, the accuracy and applicability of using onedimensional simulations to qualitatively investigate the characteristics of balanced DCMS discharge are confirmed. In terms of capturing the key physical mechanisms in DCMS, one-dimensional PIC simulation is an efficient method, which could largely alleviate the computational expense and preserve the physical fidelity of modeling results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Plasma dynamics on current-carrying magnetic flux tubes

A 1D numerical simulation is used to investigate the evolution of a plasma in a current-carrying magnetic flux tube of variable cross section. A large potential difference, parallel to the magnetic field, is applied across the domain. The result is that density minimum tends to deepen, primarily in the cathode end, and the entire potential drop becomes concentrated across the region of density minimum. The evolution of the simulation shows some sensitivity to particle boundary conditions, but the simulations inevitably evolve into a final state with a nearly stationary double layer near the cathode end. The simulation results are at sufficient variance with observations that it appears unlikely that auroral electrons can be explained by a simple process of acceleration through a field-aligned potential drop.

Swift, Daniel W.↗

Damping Effects of Drogue Parachutes on Orion Crew Module Dynamics

Because simulations of the Orion Crew Module (CM) dynamics with drogue parachutes deployed were under-predicting the amount of damping seen in free-flight tests, an attach-point damping model was applied to the Orion system. A key hypothesis in this model is that the drogue parachutes' net load vector aligns with the CM drogue attachment point velocity vector. This assumption seems reasonable and has historically produced good results, but has never been experimentally verified. The wake of the CM influences the drogue parachutes, which makes performance predictions of the parachutes difficult. Many of these effects are not currently modeled in the simulations. A forced oscillation test of the CM with parachutes was conducted in the NASA LaRC 20-Ft Vertical Spin Tunnel (VST) to gather additional data to validate and refine the attach-point damping model. A second loads balance was added to the original Orion VST model to measure the drogue parachute loads independently of the CM. The objective of the test was to identify the contribution of the drogues to CM damping and provide additional information to quantify wake effects and the interactions between the CM and parachutes. The drogue parachute force vector was shown to be highly dependent on the CM wake characteristics. Based on these wind tunnel test data, the attach-point damping model was determined to be a sufficient approximation of the parachute dynamics in relationship to the CM dynamics for preliminary entry vehicle system design. More wake effects should be included to better model the system.

Aubuchon, Vanessa V.↗

Plume Impingement Software Module for Real-Time Proximity Operations

Successfully executing proximity operations in space, such as docking or in-orbit servicing, requires sophisticated spacecraft design that accounts for induced environments. As a chaser vehicle’s attitude control thrusters fire, they create rarefied plumes that can impact the target vehicle, with the potential to overload components, exceed thermal limits, and spin the target vehicle out of control. High-fidelity simulations of the thruster plume impingement environment require the direct simulation Monte Carlo (DSMC) method, but DSMC is too computationally expensive to simulate proximity operations that involve thousands of thruster firings. For this analysis to be tractable, engineering models of the plume flowfield and impingement events are used to simulate these trajectories [1]. Currently, on-orbit plume impingement environments are modeled through an inefficient open-loop analysis cycle where the vehicle’s flight controller and plume impingement teams iterate on the trajectories until they pass the target vehicle’s plume requirements. As complex on-orbit missions evolve and become more frequent, lengthy design cycles will become operational bottlenecks. To address this gap, this work develops an advanced plume impingement module capable of operating at real-time scale that can be integrated with existing mission planning tools and onboard flight systems. The plume module leverages state-of-the-art plume simulation techniques [2] to deliver fast, physics-based impingement predictions in a software architecture that can be tailored to diverse proximity operations scenarios. A prototype of this plume impingement module is built to demonstrate the feasibility of real-time performance. This prototype completes plume impingement calculations in microseconds per target geometry mesh point. The software serves as a foundational capability for plume-aware trajectory design, operational risk assessment, and future autonomous decision-making systems.

Plume Impingement↗

Recent Dynamic Measurements and Considerations for Aerodynamic Modeling of Fighter Airplane Configurations

As airplane designs have trended toward the expansion of flight envelopes into the high angle of attack and high angular rate regimes, concerns regarding modeling the complex unsteady aerodynamics for simulation have arisen. Most current modeling methods still rely on traditional body axis damping coefficients that are measured using techniques which were intended for relatively benign flight conditions. This paper presents recent wind tunnel results obtained during large-amplitude pitch, roll and yaw testing of several fighter airplane configurations. A review of the similitude requirements for applying sub-scale test results to full-scale conditions is presented. Data is then shown to be a strong function of Strouhal number - both the traditional damping terms, but also the associated static stability terms. Additionally, large effects of sideslip are seen in the damping parameter that should be included in simulation math models. Finally, an example of the inclusion of frequency effects on the data in a simulation is shown.

Brandon, Jay M.↗

The spectral irradiance of some solar simulators and its effect on cell measurements

Moderate resolution spectral irradiance measurements in the range 300 nm to 1100 nm have been made of eight radiant sources which are currently being used as solar simulators. Spectral irradiance data are presented in graphical form. To demonstrate the interplay of source spectral distribution and cell spectral response, measurements of short circuit current of five cells of differing response characteristics have been made with these sources. These results are presented in tabular and graphical form.

Seaman, C. H.↗

Direct simulation for the instability and breakup of laminar liquid jets

A direct numerical simulation method is described for predicting the deformation of laminar liquid jets. In the present nonlinear direct simulation, the convective term, which has been discarded in past linear analyses by Rayleigh and others, is included in the hydrodynamic equations. It is shown that only by maintaining full complexity of the nonlinear surface tension term accurate drop formation can be predicted. The continuity and momentum equations in the transient form are integrated on an adaptive grid, conforming the jet and surface wave shape. The equations, which are parabolic in time and elliptic in space, are solved by a TVD scheme with characteristic flux splitting. The results of the present work are discussed and compared with available measurements and other analyses. The comparison shows that among the predictions, the current 1-D direct simulation results agree best with the experimental data. Furthermore, the computer time requirements are much (an order of magnitude) smaller than those of previously reported multidimensional analyses.

Chuech, S. G.↗

Qutrit and qubit circuits for three-flavor collective neutrino oscillations

We explore the utility of qutrits and qubits for simulating the flavor dynamics of dense neutrino systems. The evolution of such systems impacts some important astrophysical processes, such as core-collapse supernovae and the nucleosynthesis of heavy nuclei. Many-body simulations require classical resources beyond current computing capabilities for physically relevant system sizes. Quantum computers are therefore a promising candidate to efficiently simulate the many-body dynamics of collective neutrino oscillations. Previous quantum simulation efforts have primarily focused on properties of the two-flavor approximation due to their direct mapping to qubits. Furthermore, we present new quantum circuits for simulating three-flavor neutrino systems on qutrit- and qubit-based platforms, and demonstrate their feasibility by simulating systems of two, four, and eight neutrinos on IBM and Quantinuum quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Synergy of two lower hybrid waves with different frequencies on EAST

Lower hybrid current drive (LHCD) is one of the major approaches maintaining long pulse discharges on EAST. There are two lower hybrid (LH) systems on EAST launching waves at 4.6 GHz and 2.45 GHz into the plasma, respectively, and the input power of the 4.6 GHz wave is typically dominant. Here, the synergy of the two waves and the modification of the power deposition of the 4.6 GHz wave by coupling the 2.45 GHz wave are investigated theoretically. According to the phase space analysis, the wave propagation domains of the two LH waves always overlap substantially with each other for typical LHCD experimental parameters on EAST, indicating that the coupling between them might be strong. Ray-tracing/Fokker-Planck simulations for a LH current drive experiment on EAST show that the change of the power deposition profile due to the coupling could be understood by the rise and fall of the tails of the parallel electron distribution functions on different flux surfaces. If either the injected wave power at 2.45 GHz is comparable with that at 4.6 GHz or the incident N ∥ of the 2.45 GHz wave reaches a high value (e.g., $5.7/\sqrt{T_{e,0}}$), the power deposition profile of the 4.6 GHz wave can be modified greatly due to the coupling with the 2.45 GHz wave. With the parameters compatible with the ability of the LH wave system on EAST, injection of the 2.45 GHz wave can still modify the power deposition profile of the 4.6 GHz wave significantly by combining those two approaches. Finally, it is shown how a 4.6 GHz wave which cannot be Landau damped when injected alone is absorbed in low density plasmas with a magnetic shear reversal configuration in the presence of a 2.45 GHz wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Roadmap for IES Modeling and Simulation Activities: Update 2023

Since inception of the DOE-NE Integrated Energy Systems (IES) program, several program and modeling visions and roadmaps have been published. In 2020, a last comprehensive roadmap for the DOE-NE Integrated Energy Systems program was published. In 2017, a separate modeling and simulation capability development plan was published that was extensively updated in 2022. The current report incrementally updates the 2022 modeling and simulation capability roadmap. This is not a new report, but individual sections have been updated to reflect recent accomplishments and new directions of the program. The role of modeling and simulation within IES is to support the demonstration of new coupled integrated technologies along every step of the technology maturation from the strategic analysis of preferred system architecture to preliminary design, to laboratory testing, up to full commercial testing and integration. To achieve this role, modeling and simulation must be able to assess the technical performance and the economic viability of potential IES. Also, modeling and simulation must provide support for experimental evaluations (i.e., support component design and real-time operations of experimental demonstration systems). Finally, modeling and simulation can help scale-up experimental systems to the final commercial systems. This report details the current state of the modeling and simulation efforts within IES for the above-mentioned areas, provides a gap analysis, and proposes next steps for a time horizon over the next years.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Staggered solution procedures for multibody dynamics simulation

The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.

Park, K. C.↗

Effect Laws of Different Factors on Levitation Characteristics of High-T{sub c} Superconducting Maglev System with Numerical Solutions

The interaction between the high-temperature superconducting (HTS) bulk and permanent magnet guideway (PMG) has been a topic among the existing HTS magnetic levitation systems. Up to now, the theories about the electromagnetic and force characteristics have been relatively consummated. Here, using the finite element model by COMSOL Multiphysics 5.3a, we investigated the impact of the intrinsic and external factors on the levitation characteristics of HTS maglev systems, the application laws are concluded. Factors including the critical current density and geometric dimensions of the bulk superconductor, as well as the moving velocity, are discussed. The effects of the multiple back and forth movements are also analyzed. By comparing calculation results of the distributions of the induced current and the levitation force of the bulk superconductor, we can clearly confirm the effect laws of various factors on the electromagnetic characteristics and levitation force. These observations support the physical parameters that are difficult to be tested by experiments. And they also have a significant contribution to the experimental condition selections and optimization of the high-temperature superconducting magnetic levitation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Charge carrier motion and effect of fixed oxide charge in a microstructured silicon radiation detector

Signal formation in a microstructured semiconductor neutron detector is more complex than in planar diode geometry. Three-dimensional microstructures are laterally smaller than the ionization cloud length, and the electric fields may be weak enough to exhibit plasma time effects. This work is the first detailed treatment of charge carrier motion in these complex semiconductor devices to replicate the time profile and signal magnitude. Simulations were performed using COMSOL Multiphysics to investigate various parameters that affect the propagation of the charge cloud. It was observed that the size of the simulated three-dimensional structure had an impact on the induced current pulse, indicating the importance of simulation geometry optimization to accurately simulate charge cloud expansion. COMSOL Multiphysics was used to replicate accurate charge creation profiles using energy deposition information imported from radiation transport codes. A detailed simulation methodology is presented to benchmark preamplifier event pulses along with complexities in modeling the charge carrier motion along the etched microstructured trenches with Si–SiO 2 boundary conditions, including fixed oxide charge and interface trapping.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Gas-Grain Simulation Facility: Fundamental studies of particle formation and interactions. Volume 1: Executive summary and overview

An overview of the Gas-Grain Simulation Facility (GGSF) project and its current status is provided. The proceedings of the Gas-Grain Simulation Facility Experiments Workshop are recorded. The goal of the workshop was to define experiments for the GGSF--a small particle microgravity research facility. The workshop addressed the opportunity for performing, in Earth orbit, a wide variety of experiments that involve single small particles (grains) or clouds of particles. The first volume includes the executive summary, overview, scientific justification, history, and planned development of the Facility.

Fogleman, Guy↗