Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Building MCP-native hierarchical AI scientist ecosystems: a perspective on scaling multi-agent scientific discovery

Large language models (LLMs) are evolving from chatbots with limited tool-using capabilities to agentic AI systems that can perform deep research, assist in proposing hypotheses, help design experiments, automate data analysis, and draft scientific reports. However, there are currently two bottlenecks limiting LLMs' real-world impact on the broader scientific research community beyond academic demonstrations: lack of interoperability (repetitive manual tool-integration is required across scenarios) and the need for scalable coordination (unstructured communication and memory become brittle as the number of agents grows). In this Perspective, we argue that the next phase of agentic scientific discovery requires the development of an ecosystem of protocol-native agents and tools organized through hierarchies inspired by human society, beyond the current paradigm of a single monolithic “AI scientist”. We use Model Context Protocol (MCP) as a concrete example of an emerging interoperability layer for scientific tool and context exchange, and we propose three complementary pathways to increase the scaling capabilities of an MCP-native scientific ecosystem by addressing the composability issues: (1) MCP servers for high-value scientific tools maintained by domain experts, (2) automated transformation of existing code repositories into MCP services, and (3) autonomous invention and evolution of new agents and workflows. Finally, we provide a practical roadmap for scaling AI-driven scientific discovery by expanding tool supply and coordination in MCP-native scientific ecosystems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Lidar clear atmosphere multiple scattering dependence on receiver range

The degree of multiple scattering encountered by nadir-directed lidar in clear atmosphere conditions was assessed by the Monte Carlo method based on a model of the vertical distribution of aerosol scattering in the atmosphere. The lowest 3 km of the atmosphere were regarded as containing 90% of the aerosol optical thickness, with the unity normalized value for the forward peak of the aerosol phase function being 4.38 at a 694.3 nm laser wavelength. Results were obtained for two lidar receiver heights for three receiver field of view (FOV) angular halfwidths. An increase in receiver height was determined to cause a significant increase in the amount of angular scattering of the lidar signal. A factor of 20 change in receiver height produced an order of magnitude change in the single to multiple scattering ratio.

Spinhirne, J. D.↗
NASA NTRS1983

Turbulent flame propagation and combustion in spark ignition engines

Pressure measurements synchronized with high-speed motion-picture records of flame propagation have been made in a transparent-piston engine. The data show that the initial expansion speed of the flame front is close to that of a laminar flame. As the flame expands, its speed rapidly accelerates to a quasi-steady value comparable with that of the turbulent velocity fluctuations in the unburned gas. During the quasi-steady propagation phase, a significant fraction of the gas behind the visible front is unburned. Final burnout of the charge may be approximated by an exponential decay in time. The data have been analyzed in a model-independent way to obtain a set of empirical equations for calculating mass burning rates in spark-ignition engines. The burning equations contain three parameters: the laminar burning speed, a characteristic speed (uT), and a characteristic length (lT). The laminar burning speed is known from laboratory measurements. Tentative correlations relating uT and lT to engine geometry and operating variables have been derived from the engine data.

Beretta, G. P.↗
NASA NTRS1989

Experiment plans to study preignition processes of a pool fire in low gravity

Science requirements are specified to guide experimental studies of transient heat transfer and fluid flow in an enclosure containing a two-layer gas-and-liquid system heated unevenly from above. Specifications are provided for experiments in three separate settings: (1) a normal gravity laboratory, (2) the NASA-LeRC Drop towers, and (3) a space-based laboratory (e.g., Shuttle, Space Station). A rationale is developed for both minimum and desired requirement levels. The principal objective of the experimental effort is to validate a computational model of the enclosed liquid fuel pool during the preignition phase and to determine via measurement the role of gravity on the behavior of the system. Preliminary results of single-phase normal gravity experiments and simulations are also presented.

Schiller, David N.↗
NASA NTRS1990

Current-driven plasma acceleration versus current-driven energy dissipation. I - Wave stability theory

The dominant unstable electrostatic wave modes of an electromagnetically accelerated plasma are investigated. The study is the first part of a three-phase program aimed at characterizing the current-driven turbulent dissipation degrading the efficiency of Lorentz force plasma accelerators such as the MPD thruster. The analysis uses a kinetic theory that includes magnetic and thermal effects as well as those of an electron current transverse to the magnetic field and collisions, thus combining all the features of previous models. Analytical and numerical solutions allow a detailed description of threshold criteria, finite growth behavior, destabilization mechanisms and maximized-growth characteristics of the dominant unstable modes. The lower hybrid current-driven instability is implicated as dominant and was found to preserve its character in the collisional plasma regime.

Kelly, A. J.↗
NASA NTRS1994

Charged dust in planetary magnetospheres: Hamiltonian dynamics and numerical simulations for highly charged grains

We use a combination of analytical and numerical methods to investigate the dynamics of charged dust grains in planetary magnetospheres. Our emphasis is on obtaining results valid for particles that are not necessarily dominated either by gravitational or electromagnetic forces. A Hamiltonian formulation of the problem yields exact results, for all values of charge-to-mass ratio, when we introduce two constraints: particles remain in the equatorial plane and the magnetic field is taken as axially symmetric. In particular, we obtain locations of equilibrium points, the frequencies of stable periodic orbits, the topology of separatrices in phase space, and the rate of longitudinal drift. These results are significant for specific applications: motion in the nearly aligned dipolar field of Saturn, and the trajectories of arbitrarily charged particles in complex magnetic fields for limited periods of time after ejection from parent bodies. Since the model is restrictive, we also use numerical integrations of the full three-dimensional equations of motion and illustrate under what conditions the constrained problem yields reasonable results. We show that a large fraction of the intermediately charged and highly charged (gyrating) particles will always be lost to a planet's atmosphere within a few hundred hours, for motion through tilted-dipole magnetic fields. We find that grains must have a very high charge-to-mass ratio in order to be mirrored back to the ring plane. Thus, except perhaps at Saturn where the dipole tilt is very small, the likely inhabitants of the dusty ring systems are those particles that are either nearly Keplerian (weakly charged) grains or grains whose charges place them in the lower end of the intermediate charge zone. Fianlly, we demonstrate the effect of plasma drag on the orbits of gyrating particles to be a rapid decrease in gyroradius followed by a slow radial evolution of the guiding center.

Schaffer, L.↗
NASA NTRS2001

Resolution Enhancement of Spaceborne Radiometer Images

Our progress over the last year has been along several dimensions: 1. Exploration and understanding of Earth Observatory System (EOS) mission with available data from NASA. 2. Comprehensive review of state of the art techniques and uncovering of limitations to be investigated (e.g. computational, algorithmic ...). and 3. Preliminary development of resolution enhancement algorithms. With the advent of well-collaborated satellite microwave radiometers, it is now possible to obtain long time series of geophysical parameters that are important for studying the global hydrologic cycle and earth radiation budget. Over the world's ocean, these radiometers simultaneously measure profiles of air temperature and the three phases of atmospheric water (vapor, liquid, and ice). In addition, surface parameters such as the near surface wind speed, the sea surface temperature, and the sea ice type and concentration can be retrieved. The special sensor microwaves imager SSM/I has wide application in atmospheric remote sensing over the ocean and provide essential inputs to numerical weather-prediction models. SSM/I data has also been used for land and ice studies, including snow cover classification measurements of soil and plant moisture contents, atmospheric moisture over land, land surface temperature and mapping polar ice. The brightness temperature observed by SSM/I is function of the effective brightness temperature of the earth's surface and the emission scattering and attenuation of the atmosphere. Advanced Microwave Scanning Radiometer (AMSR) is a new instrument that will measure the earth radiation over the spectral range from 7 to 90 GHz. Over the world's ocean, it will be possible to retrieve the four important geographical parameters SST, wind speed, vertically integrated water vapor, vertically integrated cloud liquid water L.

Krim, Hamid↗
NASA NTRS2010

Evaluation of Anomaly Detection Capability for Ground-Based Pre-Launch Shuttle Operations

This chapter will provide a thorough end-to-end description of the process for evaluation of three different data-driven algorithms for anomaly detection to select the best candidate for deployment as part of a suite of IVHM (Integrated Vehicle Health Management) technologies. These algorithms were deemed to be sufficiently mature enough to be considered viable candidates for deployment in support of the maiden launch of Ares I-X, the successor to the Space Shuttle for NASA's Constellation program. Data-driven algorithms are just one of three different types being deployed. The other two types of algorithms being deployed include a "nile-based" expert system, and a "model-based" system. Within these two categories, the deployable candidates have already been selected based upon qualitative factors such as flight heritage. For the rule-based system, SHINE (Spacecraft High-speed Inference Engine) has been selected for deployment, which is a component of BEAM (Beacon-based Exception Analysis for Multimissions), a patented technology developed at NASA's JPL (Jet Propulsion Laboratory) and serves to aid in the management and identification of operational modes. For the "model-based" system, a commercially available package developed by QSI (Qualtech Systems, Inc.), TEAMS (Testability Engineering and Maintenance System) has been selected for deployment to aid in diagnosis. In the context of this particular deployment, distinctions among the use of the terms "data-driven," "rule-based," and "model-based," can be found in. Although there are three different categories of algorithms that have been selected for deployment, our main focus in this chapter will be on the evaluation of three candidates for data-driven anomaly detection. These algorithms will be evaluated upon their capability for robustly detecting incipient faults or failures in the ground-based phase of pre-launch space shuttle operations, rather than based oil heritage as performed in previous studies. Robust detection will allow for the achievement of pre-specified minimum false alarm and/or missed detection rates in the selection of alert thresholds. All algorithms will also be optimized with respect to an aggregation of these same criteria. Our study relies upon the use of Shuttle data to act as was a proxy for and in preparation for application to Ares I-X data, which uses a very similar hardware platform for the subsystems that are being targeted (TVC - Thrust Vector Control subsystem for the SRB (Solid Rocket Booster)).

Martin, Rodney Alexander↗
NASA NTRS2010

GLGM-3: A Degree-ISO Lunar Gravity Model from the Historical Tracking Data of NASA Moon Orbiters

In preparation for the radio science experiment of the Lunar Reconnaissance Orbiter (LRO) mission, we analyzed the available radio tracking data of previous NASA lunar orbiters. Our goal was to use these historical observations in combination with the new low-altitude data to be obtained by LRO. We performed Precision Orbit Determination on trajectory arcs from Lunar Orbiter 1 in 1966 to Lunar Prospector in 1998, using the GEODYN II program developed at NASA Goddard Space Flight Center. We then created a set of normal equations and solved for the coefficients of a spherical harmonics expansion of the lunar gravity potential up to degree and order 150. The GLGM-3 solution obtained with a global Kaula constraint (2.5 x 10(exp -4)/sq l) shows good agreement with model LP150Q from the Jet Propulsion Laboratory, especially over the nearside. The levels of data fit with both gravity models are very similar (Doppler RMS of approx.0.2 and approx. 1-2 mm/s in the nominal and extended phases, respectiVely). Orbit overlaps and uncertainties estimated from the covariance matrix also agree well. GLGM-3 shows better correlation with lunar topography and admittance over the nearside at high degrees of expansion (l > 100), particularly near the poles. We also present three companion solutions, obtained with the same data set but using alternate inversion strategies that modify the power law constraint and expectation of the individual spherical harmonics coefficients. We give a detailed discussion of the performance of this family of gravity field solutions in terms of observation fit, orbit quality, and geophysical consistency.

Mazarico, E.↗
NASA NTRS2002

Re-Os Isotopic Constraints on the Chemical Evolution and Differentiation of the Martian Mantle

The (187)Re-187Os isotopic systematics of SNC meteorites, thought to be from Mars, provide valuable information regarding the chemical processes that affected the Martian mantle, particularly with regard to the relative abundances of highly siderophile elements (HSE). Previously published data (Birck and Allegre 1994, Brandon et al. 2000), and new data obtained since these studies, indicate that the HSE and Os isotopic composition of the Martian mantle was primarily set in its earliest differentiation history. If so, then these meteorites provide key constraints on the processes that lead to variation in HSE observed in not only Mars, but also Earth, the Moon and other rocky bodies in the Solar System. Processes that likely have an effect on the HSE budgets of terrestrial mantles include core formation, magma ocean crystallization, development of juvenile crust, and the addition of a late veneer. Each of these processes will result in different HSE variation and the isotopic composition of mantle materials and mantle derived lavas. Two observations on the SNC data to present provide a framework for which to test the importance of each of these processes. First, the concentrations of Re and Os in SNC meteorites indicate that they are derived from a mantle that has similar concentrations to the Earth's mantle. Such an observation is consistent with a model where a chondritic late veneer replenished the Earth and Martian mantles subsequent to core formation on each planet. Alternative models to explain this observation do exist, but will require additional data to test the limitations of each. Second, Re-Os isotopic results from Brandon et al. (2000) and new data presented here, show that initial yos correlates with variations in the short-lived systems of (182)Hf- (182)W and (142)Sm-142Nd in the SNC meteorites (epsilon(sub W) and epsilon(sub 142Nd)). These systematics require an isolation of mantle reservoirs during the earliest differentiation history of Mars, and subsequent inefficient mixing between these reservoirs. These data show that models for the origin of isotopic variation for SNC meteorites require at least two long-lived mantle reservoirs, and possibly three. The range in the projected present day gamma(sub Os) of these reservoirs is from -5.4+/-2.6, to +4+/-1. The isotopic systematics of these reservoirs may be linked to development of cumulate crystal piles in a Martian magma ocean and variable amounts of late stage intercumulus melt. In this model, fractional crystallization of olivine and possibly other phases with slightly subchondritic Re/Os, from a solidifying magma ocean, resulted in a lower Re/Os ratio early cumulates, and a resultant low gamma(sub Os). Later cumulates or evolved melts crystallized with higher Re/Os ratios to produce the mantle reservoir(s) with consequent higher gamma(sub Os). Crystallization of the Martian magma ocean followed earliest core formation, as indicated by the correlation of epsilon(sub W) with epsilon(sub 142Nd) and initial gamma(sub Os).

Brandon, Alan D.↗
NASA NTRS2019

initMIP-Antarctica: an Ice Sheet Model Initialization Experiment of ISMIP6

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

ice sheet model↗
NASA NTRS2021

Polarimetric sensitivity of light-absorbing carbonaceous aerosols over ocean: A theoretical assessment

Visible-light-absorbing carbonaceous aerosols within the boundary layer affect the radiance and polarization states of the radiation at the top of the atmosphere. Remote sensing from suborbital and satellite-based platforms utilizes these radiance and polarization signals to retrieve the key properties of these aerosols. Recent retrieval algorithms have shown a progressive trend toward including multi-angular and multi-spectral polarimetric measurements to produce better retrieval accuracy in comparison to those using measurements based on a single viewing angle. Here, we perform a theoretical investigation of the top of atmosphere (TOA) radiance-related reflectance factor (bidirectional reflectance factor (BRF)) and the two types of polarimetry-related factors (polarized bidirectional reflectance factor (pBRF) and the degree of linear polarization (DoLP)) for different types of atmospheric light-absorbing carbonaceous aerosols as a function of particle size distribution. We selected three polarimetric bands corresponding to those utilized by NASA's Airborne Multiangle SpectroPolarimetric Imager (AirMSPI)—near-UV (470 nm), visible (660 nm), and near-infrared (865 nm)—for our simulations which were performed over ocean surface using the successive order of scattering (SOS) algorithm coupled to a Lorenz-Mie aerosol optics model. The analysis of particle phase matrix elements indicates a close relationship between the angular dependencies of DoLP and associated phase matrix components at the shortest polarimetric band (470 nm). Using Jacobian analysis, we find that the radiance- and polarimetry-related reflectance factors of weakly light-absorbing aerosols, such as brown carbon, are more sensitive to changes in particle size and imaginary refractive index in comparison with those of black carbon, which is strongly light-absorbing. Our results suggest that the DoLP data could be used by future retrieval algorithms for reliably estimating microphysical properties of absorbing carbonaceous aerosols with imaginary refractive index less than 0.4.

refractive index↗
DOE OSTI2025

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2022

Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions

With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.

58 GEOSCIENCES↗
DOE OSTI2023

Modified Virtual Oscillator Based Operation of Grid Forming Converters with Single Voltage Sensor

Virtual oscillator based control architecture for grid forming converter with inner current control loop for additional objective accomplishments is proposed in this paper. In general, virtual oscillator is based on Lienard type of oscillator equation with a cubic nonlinearity. Virtual oscillator alone cannot accommodate unbalanced or harmonic mitigation while operating converters in grid forming mode. Therefore, to accomplish better tracking performance, incorporating the provision for nonlinear and/or unbalanced loads, an inner current control loop is utilized which is based upon Lyapunov energy function type control architecture due to its advantages. To reduce the overall cost of implementation and ensure proper phase difference in the generated voltages, the implemented virtual oscillator operates with only one voltage feedback from the point of common coupling as presented in this paper. To verify the effectiveness of the proposed approach, the overall system has been modeled in MATTLAB/Simulink and PLECS domain and case studies showing the successful production of low harmonic load voltages under nonideal loading condition along with other important results are presented in this paper.

dual second order generalized integrator↗
DOE OSTI2020

Complex Disorder in Type-I Clathrates: Synthesis and Structural Characterization of A8GaxSn46−x (A = Rb, Cs; 6.9 < x < 7.5)

Exploratory studies in the systems Rb–Ga–Sn and Cs–Ga–Sn yielded the cubic type-I clathrates with refined compositions Rb8GaxSn46−x and Cs8GaxSn46−x (6.9 < x < 7.5). Nearly single-phase materials with good crystallinity were obtained from stoichiometric reactions of the elements. The structures were characterized by means of single-crystal X-ray diffraction methods. Both Rb8GaxSn46−x and Cs8GaxSn46−x represents cases, where a Group 13 element randomly substitutes a Group 14 element in the structure. The extent of Ga/Sn mixing is apparently governed by the drive of the system to achieve an optimal valence electron count, and hence, Rb8GaxSn46−x and Cs8GaxSn46−x (x ≈ 8) can be regarded as Zintl phases. This notion is supported by structure refinements on a multitude of single-crystal X-ray diffraction data, which also confirm that both types of cages in the cubic type-I structure are fully occupied by Rb and Cs atoms. The open-framework, comprised of 46 nodes per formula unit, adapts to the incorporation of nearly eight Ga atoms within the matrix of Sn, whereby small, short-range distortions result. The exact nature of these effects is still unclear, as so far, the structural variations could only be modeled as both positional and occupational disorder at one of three framework sites. Since vacancies in the structures of the binary type-I clathrates A8Sn46−x☐x (A = Rb, Cs; ☐ = missing Sn atom) are also known to cause local distortions, the latter were also synthesized with the same protocols used for the synthesis of A8GaxSn46−x and structurally re-analyzed. The results from the latter studies confirm that homogeneity issues abound, and that the final structures/compositions are an intricate function of the experimental conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modified Virtual Oscillator-Based Operation of Grid-Forming Converters with Single Voltage Sensor: Preprint

This paper proposes a virtual oscillator-based control architecture for a grid-forming converter with an inner current control loop. In general, the virtual oscillator is based on a Lienard type of oscillator equation with a cubic nonlinearity. A virtual oscillator alone cannot accommodate unbalanced or harmonic mitigation while operating converters in grid-forming mode. Therefore, to accomplish better tracking performance, incorporating the provision for nonlinear and/or unbalanced loads, an inner current control loop is utilized that is based on the Lyapunov energy function type control architecture. To reduce the overall cost of implementation and ensure proper phase difference in the generated voltages, the implemented virtual oscillator operates with only one voltage feedback from the point of common coupling. To verify the effectiveness of the proposed approach, the overall system has been modeled in MATTLAB/Simulink and PLECS domain. This paper also presents case studies showing the successful production of low harmonic load voltages under nonideal loading conditions, along with other important results.

dual second order generalized integrator↗