Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameter Identifiability”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

On the dose calculation at the cellular level and its implications for the RBE of {sup 99m}Tc and {sup 123}I

Purpose: Based on the authors’ previous findings concerning the radiotoxicity of{sup 99m}Tc, the authors compared the cellular survival under the influence of this nuclide with that following exposure to the Auger electron emitter {sup 123}I. To evaluate the relative biological effectiveness (RBE) of both radionuclides, knowledge of the absorbed dose is essential. Thus, the authors present the dose calculations and discuss the results based on different models of the radionuclide distribution. Both different target volumes and the influence of the uptake kinetics were considered. Methods: Rat thyroid PC Cl3 cells in culture were incubated with either{sup 99m}Tc or {sup 123}I or were irradiated using 200 kV x-rays in the presence or absence of perchlorate. The clonogenic cell survival was measured via colony formation. In addition, the intracellular radionuclide uptake was quantified. Single-cell dose calculations were based on Monte Carlo simulations performed using Geant4. Results: Compared with external radiation using x-rays (D{sub 37} = 2.6 Gy), the radionuclides {sup 99m}Tc (D{sub 37} = 3.5 Gy), and {sup 123}I (D{sub 37} = 3.8 Gy) were less toxic in the presence of perchlorate. In the absence of perchlorate, the amount of activity a{sub 37} that was necessary to reduce the surviving fraction (SF) to 0.37 was 22.8 times lower for {sup 99m}Tc and 12.4 times lower for {sup 123}I because of the dose increase caused by intracellular radionuclide accumulation. When the cell nucleus was considered as the target for the dose calculation, the authors found a RBE of 2.18 for {sup 99m}Tc and RBE = 3.43 for {sup 123}I. Meanwhile, regarding the dose to the entire cell, RBE = 0.75 for {sup 99m}Tc and RBE = 1.87 for {sup 123}I. The dose to the entire cell was chosen as the dose criterion because of the intracellular radionuclide accumulation, which was found to occur solely in the cytoplasm. The calculated number of intracellular decays per cell was (975 ± 109) decays/MBq for {sup 99m}Tc and (221 ± 82) decays/MBq for {sup 123}I. Conclusions: The authors’ data indicate that extra-nuclear targets to Auger electrons exist, which is obvious from our dose calculations. When considering the dose to the cell nucleus, the authors found an enhanced RBE for{sup 99m}Tc and {sup 123}I relative to acute x-ray irradiation and pure extracellular irradiation with both radionuclides. Surprisingly, the authors did not find any radionuclide accumulation in the cell nucleus, indicating that there are additional radiosensitive targets besides the DNA. In addition, the authors demonstrated the necessity of cellular dose calculations in radiobiological experiments using unsealed radionuclides and identified the relevant parameters.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

TH-A-18C-01: Design Optimization of Segmented Scintillators for Megavoltage Cone- Beam CT

Purpose: Active matrix flat-panel imagers incorporating thick, segmented scintillators for megavoltage cone-beam CT (MV CBCT) imaging have demonstrated strong potential for facilitating soft-tissue visualization at low, clinically practical doses. In order to identify scintillator design parameters that optimize performance for this purpose, a modeling technique which includes both radiation and optical effects and which lends itself to computationally practical implementation has been developed and explored. Methods: A hybrid modeling technique, based on Monte Carlo event-by-event simulation of radiation transport and separate determination of optical effects, was devised as an alternative to computationally prohibitive event-by- event simulations of both radiation and optical transport. The technique was validated against empirical results from a previously reported 1.13 cm thick, 1.016 mm element-to-element pitch BGO scintillator prototype. Using this technique, the contrast-to-noise ratio (CNR) and spatial resolution performance of numerous scintillator designs, with thicknesses ranging from 0.5 to 6 cm and pitches ranging from 0.508 to 1.524 mm, were examined. Results: CNR and spatial resolution performance for the various scintillator designs demonstrate complex behavior as scintillator thickness and pitch are varied - exhibiting a clear trade-off between these two imaging metrics up to a thickness of ~3 cm. Based on these results, an optimization map highlighting those regions of design that provide a balance between these metrics was created. The map indicates that, for a given set of optical parameters, scintillator thickness and pitch can be judiciously chosen to maximize performance without resorting to thicker, more costly scintillators. Conclusion: Modeling radiation and optical effects in thick, segmented scintillators through use of a hybrid modeling technique provides a practical way to gain insight as to how to optimize the performance of such devices for radiotherapy imaging. Assisted by such modeling, the development of practical designs should greatly facilitate low-dose, soft tissue visualization through MV CBCT imaging. This project was supported in part by NIH grant R01 CA051397.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

A strongly robust type II Weyl fermion semimetal state in Ta 3 S 2

Weyl semimetals are of great interest because they provide the first realization of the Weyl fermion, exhibit exotic quantum anomalies, and host Fermi arc surface states. The separation between Weyl nodes of opposite chirality gives a measure of the robustness of the Weyl semimetal state. To exploit the novel phenomena that arise from Weyl fermions in applications, it is crucially important to find robust separated Weyl nodes. Here, we propose a methodology to design robust Weyl semimetals with well-separated Weyl nodes. Using this methodology as a guideline, we search among the material parameter space and identify by far the most robust and ideal Weyl semimetal candidate in the single-crystalline compound tantalum sulfide (Ta 3 S 2 ) with new and novel properties beyond TaAs. Crucially, our results show that Ta 3 S 2 has the largest k-space separation between Weyl nodes among known Weyl semimetal candidates, which is about twice larger than the measured value in TaAs and 20 times larger than the predicted value in WTe 2 . Moreover, all Weyl nodes in Ta 3 S 2 are of type II. Therefore, Ta 3 S 2 is a type II Weyl semimetal. Furthermore, we predict that increasing the lattice by <4% can annihilate all Weyl nodes, driving a novel topological metal-to-insulator transition from a Weyl semimetal state to a topological insulator state. The robust type II Weyl semimetal state and the topological metal-to-insulator transition in Ta 3 S 2 are potentially useful in device applications. Our methodology can be generally applied to search for new Weyl semimetals.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Synthesis, X-ray Diffraction Characterization, Mössbauer Spectroscopy, and Dielectric Properties of Solid Solutions in the PbFe{sub 2/3}W{sub 1/3}O{sub 3}–PbSc{sub 2/3}W{sub 1/3}O{sub 3} System

Ceramic Pb(Fe{sub 1–x}Sc{sub x}){sub 2/3}W{sub 1/3}O{sub 3} samples with 0 ≤ x ≤ 1 have been prepared and characterized by X-ray diffraction, Mössbauer spectroscopy, and dielectric and pyroelectric measurements. The stoichiometry ranges of the perovskite solid solutions in this system have been identified, their structural parameters have been determined, and their dielectric permittivity, dielectric loss tangent, resistivity, and thermally stimulated depolarization current have been measured as functions of temperature. The composition dependences of the dielectric properties for the solid solutions have been obtained. The solid solutions have been shown to exhibit ferroelectric relaxor properties, with a well-defined maximum in their permittivity in the range 180–250 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Synergy of Coupled Kink and Tearing Modes in Fast Ion Transport

One aspect of the interaction between fast ions and magnetohydrodynamic (MHD) instabilities is the fast ion transport. Coupled kink and tearing MHD instabilities have also been reported to cause fast ion transport. Recently, the ''kick" model has been developed to compute the evolution of the fast ion distribution from the neutral beam injection using instabilities as phase space resonance sources. The goal of this paper is to utilize the kick model to understand the physics of fast ion transport caused by the coupled kink and tearing modes. Soft X-ray diagnostics are used to identify the mode parameters in NSTX. The comparison of neutron rates measured and computed from time-dependent TRANSP simulation with the kick model shows the coupling of kink and tearing mode is important in determination of the fast ion transport. The numerical scan of the mode parameters shows that the relative phase of the kink and tearing modes and the overlapping of kink and tearing mode resonances in the phase space can affect the fast ion transport, suggesting that the synergy of the coupled modes may be causing the fast ion transpor

Energetic Particles↗
DOE OSTI2020

Plasmonic nanoarcs: a versatile platform with tunable localized surface plasmon resonances in octave intervals

The tunability of the longitudinal localized surface plasmon resonances (LSPRs) of metallic nanoarcs is demonstrated with key relationships identified between geometric parameters of the arcs and their resonances in the infrared. The wavelength of the LSPRs is tuned by the mid-arc length of the nanoarc. The ratio between the attenuation of the fundamental and second order LSPRs is governed by the nanoarc central angle. Beneficial for plasmonic enhancement of harmonic generation, these two resonances can be tuned independently to obtain octave intervals through the design of a non-uniform arc-width profile. Because the character of the fundamental LSPR mode in nanoarcs combines an electric and a magnetic dipole, plasmonic nanoarcs with tunable resonances can serve as versatile building blocks for chiroptical and nonlinear optical devices.

42 ENGINEERING↗
DOE OSTI2024

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

Big Data Mining for Assessing Calibration of Building Energy Models

Residential and commercial buildings in China, India, the United States (US), United Kingdom (UK), and Italy consume 39-45% of each nation's primary energy. Building energy models can be used to automatically optimize the return-on-investment for retrofits to improve a building’s energy efficiency. However, with an average of 3,000 building descriptors necessary to accurately simulate a single building, there is a market need to reduce the transaction cost for creating a more robust model for simulating every building in a city and accurately estimate savings prior to capital expenditures.We used two of the world’s fastest supercomputers, assembled unique datasets, and developed innovative algorithms for big data mining to assess different methods to create an accurate building energy model. The team has leveraged a total of eight high performance computing resources and 218 servers to analyze the best methods, metrics, and algorithms for creating an accurate building energy model beyond current industry standards necessary for private-sector financing. The project developed the world’s fastest buildings simulator, has completed over 8 million simulations totaling over 300TB, and mined this data with over 130,000 parallel artificial intelligence algorithms. This was used to quantify accuracy of AI-developed calibration algorithms with results that surpass industry standard guidelines and can identify individual building parameters to between 15% and 32% of their actual value.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

FAST Model Calibration and Validation of the OC5-DeepCwind Floating Offshore Wind System Against Wave Tank Test Data

During the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project, which focused on the validation of numerical methods through comparison against tank test data, the authors created a numerical FAST model of the 1:50-scale DeepCwind semisubmersible system that was tested at the Maritime Research Institute Netherlands ocean basin in 2013. The OC5 project revealed a general underprediction of loads and motions by the participating numerical models. This paper discusses several model calibration studies that were conducted to identify potential model parameter adjustments that help to improve the agreement between the numerical simulations and the experimental test data. These calibration studies cover wind-field-specific parameters (coherence, turbulence), and hydrodynamic and aerodynamic modeling approaches, as well as rotor model (blade-pitch and blade-mass imbalances) and tower model (structural tower damping coefficient) adjustments. These calibration studies were conducted based on relatively simple calibration load cases (wave only/wind only). The agreement between the final FAST model and experimental measurements is then assessed based on more complex combined wind and wave validation cases. The analysis presented in this paper does not claim to be an exhaustive parameter identification study but is aimed at describing the qualitative impact of different model parameters on the system response. This work should help to provide guidance for future systematic parameter identification and uncertainty quantification efforts.

17 WIND ENERGY↗
DOE OSTI2018

Updates on Collective Effects Estimations for the JLEIC

JLEIC is the high luminosity and high polarization electron-ion collider (EIC) currently under active design at Jefferson Lab. It aims at high luminosity 1033~1034 cm-2s-1for a wide range of ion species and center-of-mass energies. This luminosity performance relies sensibly on beam stability with high intensity electron and ion beam operation. The impedance budget analysis and the estimations of the single and multibunch instabilities are currently underway. In this paper, we present the update status of estimations for the longitudinal and transverse coherent instabilities, and identify areas or parameter regimes where special attentions for instability mitigations are required.

Li, Rui↗
DOE OSTI2018

ISAP-MATLAB Package for Sensitivity Analysis of High-Dimensional Stochastic Chemical Networks

Stochastic simulation and modeling play an important role to elucidate the fundamental mechanisms in complex biochemical networks. The parametric sensitivity analysis of reaction networks becomes a powerful mathematical and computational tool, yielding information regarding the robustness and the identifiability of model parameters. However, due to overwhelming computational cost, parametric sensitivity analysis is a extremely challenging problem for stochastic models with a high-dimensional parameter space and for which existing approaches are very slow. Here we present an information-theoretic sensitivity analysis in path-space (ISAP) MATLAB package that simulates stochastic processes with various algorithms and most importantly implements a gradient-free approach to quantify the parameter sensitivities of stochastic chemical reaction network dynamics using the pathwise Fisher information matrix (PFIM; Pantazis, Katsoulakis, and Vlachos 2013). The sparse, block-diagonal structure of the PFIM makes its computational complexity scale linearly with the number of model parameters. As a result of the gradientfree and the sparse nature of the PFIM, it is highly suitable for the sensitivity analysis of stochastic reaction networks with a very large number of model parameters, which are typical in the modeling and simulation of complex biochemical phenomena. Finally, the PFIM provides a fast sensitivity screening method (Arampatzis, Katsoulakis, and Pantazis 2015) which allows it to be combined with any existing sensitivity analysis software.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Compaction and Morphology of Lost Circulation Materials

Lost circulation material (LCM) selection is critical to effectively and efficiently treating wellbore fluid losses in geothermal drilling where costs of treatment can be as much as 30% of the total drilling cost. We conducted several uniaxial compaction experiments on 10 different materials and several material mixtures to identify critical mechanical parameters of each. Materials degraded at 200°C were also investigated to understand how elevated temperatures in geothermal wells would degrade their compaction behavior. Granular materials tended to have lower compressibility and higher compression resistance, while more elongated and softer materials had less mechanical stiffness. Mixing materials tended to moderate the mechanical behaviors while heating universally increased the compaction of materials. Microscopy showed that particle strength tended to correlate positively with roundness and circularity and negatively with elongation of a material. Convexity of the degraded and undegraded materials showed heating may have increased the convexity or roughness of the individual particles. In conclusion, we concluded that granular materials are likely to provide the best seals in wells but that a mixture of size distribution, mechanical rigidity, and elongation is more likely to form a better seal for geothermal wells.

annular pressure drilling↗
DOE OSTI2011

FINAL SCIENTIFIC/TECHNICAL REPORT

Dynalene Inc has developed and patented a fuel cell coolant with the help of DOE SBIR Phase I and Phase II funding (Project DE-FG02-04ER83884). However, this coolant could only be produced in lab scale (500 ml to 2 L) due to problems in the optimization and scale-up of a nanoparticle ingredient. This project optimized the nanoparticle production process in 10 L and 100 L reactors (which translates to about 5000 gallons of coolant), optimized the filtration process for the nanoparticles, and develop a high throughput production as well as quality control method for the final coolant formulation. Scale-up of nanoparticle synthesis (using emulsion polymerization) is an extremely challenging task. Dynalene researchers, in collaboration with a university partner, identified all the parameters affecting the size, charge density and coagulation characteristics of the nanoparticles and then optimized these parameters to achieve the goals and the objectives of this project. Nanoparticle synthesis was demonstrated to be reproducible in the 10 L and 100 L scales.

08 HYDROGEN↗
DOE OSTI2012

Sensitivity Analysis of FEAST-Metal Fuel Performance Code: Initial Results

This memo documents the completion of the LANL milestone, M3FT-12LA0202041, describing methodologies and initial results using FEAST-Metal. The FEAST-Metal code calculations for this work are being conducted at LANL in support of on-going activities related to sensitivity analysis of fuel performance codes. The objective is to identify important macroscopic parameters of interest to modeling and simulation of metallic fuel performance. This report summarizes our preliminary results for the sensitivity analysis using 6 calibration datasets for metallic fuel developed at ANL for EBR-II experiments. Sensitivity ranking methodology was deployed to narrow down the selected parameters for the current study. There are approximately 84 calibration parameters in the FEAST-Metal code, of which 32 were ultimately used in Phase II of this study. Preliminary results of this sensitivity analysis led to the following ranking of FEAST models for future calibration and improvements: fuel conductivity, fission gas transport/release, fuel creep, and precipitation kinetics. More validation data is needed to validate calibrated parameter distributions for future uncertainty quantification studies with FEAST-Metal. Results of this study also served to point out some code deficiencies and possible errors, and these are being investigated in order to determine root causes and to improve upon the existing code models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Radiolysis Process Model

Assessing the performance of spent (used) nuclear fuel in geological repository requires quantification of time-dependent phenomena that may influence its behavior on a time-scale up to millions of years. A high-level waste repository environment will be a dynamic redox system because of the time-dependent generation of radiolytic oxidants and reductants and the corrosion of Fe-bearing canister materials. One major difference between used fuel and natural analogues, including unirradiated UO 2 , is the intense radiolytic field. The radiation emitted by used fuel can produce radiolysis products in the presence of water vapor or a thin-film of water (including OH• and H• radicals, O 2 -, eaq, H 2 O 2 , H 2 , and O 2 ) that may increase the waste form degradation rate and change radionuclide behavior. H 2 O 2 is the dominant oxidant for spent nuclear fuel in an O 2 depleted water environment, the most sensitive parameters have been identified with respect to predictions of a radiolysis model under typical conditions. As compared with the full model with about 100 reactions it was found that only 30-40 of the reactions are required to determine [H 2 O 2 ] to one part in 10–5 and to preserve most of the predictions for major species. This allows a systematic approach for model simplification and offers guidance in designing experiments for validation.

Hydrogen Peroxide↗
DOE OSTI2014

Analysis report for WIPP colloid model constraints and performance assessment parameters

An analysis of the Waste Isolation Pilot Plant (WIPP) colloid model constraints and parameter values was performed. The focus of this work was primarily on intrinsic colloids, mineral fragment colloids, and humic substance colloids, with a lesser focus on microbial colloids. Comments by the US Environmental Protection Agency (EPA) concerning intrinsic Th(IV) colloids and Mg-Cl-OH mineral fragment colloids were addressed in detail, assumptions and data used to constrain colloid model calculations were evaluated, and inconsistencies between data and model parameter values were identified. This work resulted in a list of specific conclusions regarding model integrity, model conservatism, and opportunities for improvement related to each of the four colloid types included in the WIPP performance assessment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Experimental deformation in sandstone, carbonates and quartz aggregate

The first part of my thesis is mainly focused on the effect of grain size distribution on compaction localization in porous sandstone. To identify the microstructural parameters that influence compaction band formation, I conducted a systematic study of mechanical deformation, failure mode and microstructural evolution in Bleurswiller and Boise sandstones, of similar porosity (~25%) and mineralogy but different sorting. Discrete compaction bands were observed to develop over a wide range of pressure in the Bleurswiller sandstone that has a relatively uniform grain size distribution. In contrast, compaction localization was not observed in the poorly sorted Boise sandstone. My results demonstrate that grain size distribution exerts important influence on compaction band development, in agreement with recently published data from Valley of Fire and Buckskin Gulch, as well as numerical studies. The second part aimed to improve current knowledge on inelastic behavior, failure mode and brittle-ductile transition in another sedimentary rock, porous carbonates. A micritic Tavel (porosity of ~13%) and an allochemical Indiana (~18%) limestones were deformed under compaction in wet and dry conditions. At lower confining pressures, shear localization occurred in brittle faulting regime. Through transitional regime, the deformation switched to cataclastic flow regime at higher confining pressure. Specifically in the cataclastic regime, the (dry and wet) Tavel and dry Indiana failed by distributed cataclastic flow, while in contrast, wet Indiana failed as compaction localization. My results demonstrate that different failure modes and mechanical behaviors under different deformation regimes and water saturation are fundamental prior to any geophysical application in porous carbonates. The third part aimed to focus on investigating compaction on quartz aggregate starting at low (MPa) using X-ray diffraction. We report the diffraction peak evolution of quartz with increasing pressures. Through evaluating the unit cell lattice parameters and the volume of the quartz sample, macroscopic stress and strain were resolved. Moreover, we observed quartz peak broadened asymmetrically at low pressure, such extent is more prominent in axial than in radial direction. Our evaluation on peak [101] (highest intensity among peaks) demonstrated that full width at half maximum can be a good proxy for microscopic stress distribution. We observed deviations in the pressurevolume curves at P = ~0.4 GPa and speculated that it was the point of which onset of grain crushing and pore collapse occur in quartz. This is on the same order of which onset of grain crushing (commonly known as P*) is observed in sandstones in the rock mechanics literature. This demonstrated that there is potential in estimating grain crushing and pore collapse pressure with our technique.

58 GEOSCIENCES↗
DOE OSTI2017

Stability and Control of Burning Tokamak Plasmas with Resistive Walls: Final Technical Report

This project is focused on theoretical and computational development for quantitative prediction of the stability and control of the equilibrium state evolution in toroidal burning plasmas, including its interaction with the surrounding resistive wall. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this project have performed theoretical and computational analyses, guided by analytic modeling, to address this physics in realistic configurations. The overall goal has been to understand the key physics mechanisms that describe stable toroidal burning plasmas under active feedback control. Several relevant achievements have occurred during this project, leading to publications and invited conference presentations. In theoretical efforts, with the physics of the resonant layers, resistive wall, and toroidal momentum transport included, this study has extended from cylindrical resistive plasma - resistive wall models with feedback control to toroidal geometry with strong shaping to study mode coupling effects on the stability. These results have given insight into combined tearing and resistive wall mode behavior in simulations and experiment, while enabling a rapid exploration of plasma parameter space, to identify possible domains of interest for large plasma codes to investigate in more detail. Resonant field amplification and quasilinear torques in the presence of error fields and velocity shear have also been investigated. Here it was found, surprisingly, that the Maxwell torque on resonant layers in the plasma which exhibit finite real frequencies ωr in their response is significantly different from the conventional results based on tearing layers with pure real growth (or damping) rates. This observation suggests the possibility that the torque on the tearing layers can lock the plasma rotation to this finite phase velocity, which may lead to locking in which velocity shear is maintained. More broadly, the sources of all torques driving flows in magnetic confinement experiments is not fully understood, and this theoretical work may shed light on puzzling experimental results. It was also found that real frequencies occur over a wide range of plasma response regimes, and are indeed the norm and not the exception, often leading to profound effects on the locking torque. Also, the influence of trapped energetic ions orbiting over the resistive plasma mode structure, a critical effect in burning plasmas, was investigated through analytic modeling and analysis of simulations and experiment. This effort has shown that energetic ions can drive the development of disruptive instabilities, but also damp and stabilize the instabilities, depending on the details of the shear in the equilibrium magnetic field. This finding could be critical to maintaining stable operations in burning plasmas. In the most recent work, a series of simulations have been conducted to study the effect of differential flow and energetic ion effects on entry to the onset of a disruptive instability in the most realistic conditions possible, with preexisting nonlinearly saturated benign instabilities. Throughout this work, the linear and quasilinear theory of resonant layers with differential flow between them, their interaction with resistive wall and error fields, and energetic ions effects, have been used to understand realistic simulations of mode onset and the experimental discharges they represent. These studies will continue to answer remaining questions about the relation between theoretical results obtained in this project and observations of the onset and evolution of disruptive instabilities in experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗