Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependent behavior”

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,747 records · Page 97

Data-Driven Performance Optimization of Gamma Spectrometers With Many Channels

In gamma spectrometers with variable spectroscopic performance across many channels (e.g., many pixels or voxels), a tradeoff exists between including data from successively worse-performing readout channels and increasing efficiency. Brute-force calculation of the optimal set of included channels is exponentially infeasible as the number of channels grows, and approximate methods are required. In this work, we present a data-driven framework for attempting to find near-optimal sets of included detector channels. The framework leverages non-negative matrix factorization (NMF) to learn the behavior of gamma spectra across the detector and clusters similarly-performing detector channels together. Performance comparisons are then made between spectra with channel clusters removed, which is more feasible than brute force. The framework is general and can be applied to arbitrary, user-defined performance metrics depending on the application. We apply this framework to optimizing gamma spectra measured by H3D M400 CdZnTe (CZT) spectrometers, which exhibit variable performance across their crystal volumes. In particular, we show several examples optimizing various performance metrics for uranium and plutonium gamma spectra in non-destructive assay (NDA) for nuclear safeguards, and explore trends in performance versus parameters such as clustering algorithm type. We also compare the NMF + clustering pipeline to several non-machine-learning (ML) algorithms, including several greedy algorithms. Although, we find that the NMF + clustering pipeline tends to find the best-performing set of detector voxels, significantly improving over the unoptimized spectra, but that a greedy accumulation of spectra segmented by detector depth can, in some cases, give similar performance improvements in much less computation time.

Energy resolution↗

Estimation of high temperature metal-silicate partition coefficients

It has been known for some time that abundances of siderophile elements in the upper mantle of the Earth are far in excess of those expected from equilibrium between metal and silicate at low pressures and temperatures. Murthy (1991) has re-examined this excess of siderophile element problem by estimating liquid metal/liquid silicate partition coefficients reduces from their measured values at a lower temperature, implying that siderophile elements become much less siderophilic at high temperatures. Murthy then draws the important conclusion that metal/silicate equilibrium at high temperatures can account for the abundances of siderophile elements in the Earth's mantle. Of course, his conclusion is critically dependent on the small values of the partition coefficients he calculates. Because the numerical values of most experimentally-determined partition coefficients increase with increasing temperature at both constant oxygen fugacity and at constant redox buffer, we think it is important to try an alternative extrapolation for comparison. We have computed high temperature metal/silicate partition coefficients under a different set of assumptions and show that such long temperature extrapolations yield values which are critically dependent upon the presumed chemical behavior of the siderophile elements in the system.

Jones, John H.↗

Long-term changes in amphetamine-induced reinforcement and aversion in rats following exposure to 56Fe particle

Exposing rats to heavy particles produces alterations in the functioning of dopaminergic neurons and in the behaviors that depend upon the integrity of the dopaminergic system. Two of these dopamine-dependent behaviors include amphetamine-induced reinforcement, measure using the conditioned place preference procedure, and amphetamine-induced reinforcement, measured using the conditioned place preference procedure, and amphetamine-induced aversion, measured using the conditioned taste aversion. Previous research has shown that exposing rats to 1.0 Gy of 1GeV/n 56Fe particles produced a disruption of an amphetamine-induced taste aversion 3 days following exposure, but produced an apparent enhancement of the aversion 112 days following exposure. The present experiments were designed to provide a further evaluation of these results by examining taste aversion learning 154 days following exposure to 1.0 Gy 56Fe particles and to establish the convergent validity of the taste aversion results by looking at the effects of exposure on the establishment of an amphetamine-induced conditioned place preference 3, 7, and 16 weeks following irradiation. The taste aversion results failed to confirm the apparent enhancement of the amphetamine-induced CTA observed in the prior experiment. However, exposure to 56Fe particles prevented the acquisition of amphetamine-induced place preference at all three-time intervals. The results are interpreted as indicating that exposure to heavy particles can produce long-term changes in behavioral functioning. c2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

Non-NASA Center↗

Characterization and Detailed Analysis of Regression Behavior for HTPB Solid Fuels Containing High Aluminum Loadings

NASA Marshall Space Flight Center's Materials and Processes Department, with support from the Propulsion Systems Department, has renewed the development and maintenance of a hybrid test bed for exposing ablative thermal protection materials to an environment similar to that seen in solid rocket motors (SRM). The Solid Fuel Torch (SFT), operated during the Space Shuttle program, utilized gaseous oxygen for oxidizer and an aluminized hydroxyl-terminated polybutadiene (HTPB) fuel grain to expose a converging section of phenolic material to a 400 psi, 2-phase flow combustion environment. The configuration allows for up to a 2 foot long, 5 inch diameter fuel grain cartridge. Wanting to now test rubber insulation materials with a turn-back feature to mimic the geometry of an aft dome being impinged by alumina particles, the throat area has now been increased by several times to afford flow similarity. Combined with the desire to maintain a higher operating pressure, the oxidizer flow rate is being increased by a factor of 10. Out of these changes has arisen the need to characterize the fuel/oxidizer combination in a higher mass flux condition than has been previously tested at MSFC, and at which the literature has little to no reporting as well. For (especially) metalized fuels, hybrid references have pointed out possible dependence of fuel regression rate on a number of variables: mass flux, G - oxidizer only (G0), or - total mass flux (Gtot), Length, L, Pressure, P, and Diameter, D.

Kibbey, Timothy P.↗

Crew collaboration in space: a naturalistic decision-making perspective

Successful long-duration space missions will depend on the ability of crewmembers to respond promptly and effectively to unanticipated problems that arise under highly stressful conditions. Naturalistic decision making (NDM) exploits the knowledge and experience of decision makers in meaningful work domains, especially complex sociotechnical systems, including aviation and space. Decision making in these ambiguous, dynamic, high-risk environments is a complex task that involves defining the nature of the problem and crafting a response to achieve one's goals. Goal conflicts, time pressures, and uncertain outcomes may further complicate the process. This paper reviews theory and research pertaining to the NDM model and traces some of the implications for space crews and other groups that perform meaningful work in extreme environments. It concludes with specific recommendations for preparing exploration crews to use NDM effectively.

Review, Tutorial↗

Thermo-Oxidation Analysis of Structural Adhesives: An Experimental Study

Abstract Cross-linked polymers play an important role in structural applications as adhesives. In this article, we performed an experimental study to understand the effects of high temperature i.e. thermo-oxidative aging on the mechanical behavior of polymers. Subsequently, the constitutive behavior, i.e. stress-strain relation of these adhesives before and after the accelerated aging tests were studied on samples aged at different temperature and different time durations at zero humidity. The aging conditions and periods selected for the tests have been specifically designed to simulate the effects of aging under the working conditions that the materials will be subjected to in their standard applications. The focus of this study is on the investigation of the effects of aging on the polymeric adhesives used in automotive applications. Consequently, four different polymeric adhesives (one silicone-based, two polyurethane-based, and one acrylic type), which are being used extensively in these applications, have been chosen. The quasi-static tensile loading was used on these four adhesives and the changes in the constitutive behavior were measured. The failure tensile tests as well as the cyclic tensile tests were performed in our lab using a Test-Resources tensile testing machine. The temperatures that were used in this work are 60, 80 and 95°C for varying aging periods of 24hrs, 240hrs, and 720hrs. Then, all the data was collected and consequently analyzed, and the resulting behavior was compared to that of the unaged material, showing varying results depending on the material, temperature, and aging period. Different behaviors have been observed after aging. The polymers’ behavior was expanding from no effect to significant damage caused by aging conditions. This work provides the core understanding required for modeling the material behavior during aging and lays down the groundwork for further data-infused model-free approaches.

Alazhary, Sharif↗

Analysis of Gaspra lightcurves using Galileo shape and photometric models

Galileo-based models for the shape of 951 Gaspra and the global-average photometric behavior of its surface have been used to model a representative subset of the asteroid's telescopic lightcurves. Fitting the synthetic lightcurves to the observed timing of lightcurve extrema, and knowing the orientation of Gaspra's axes at the time of the Galileo flyby, leads to a sidereal rotation period for the asteroid of 7.042024 +/- 0.000020 hr, a slight change from the period reported by Magnusson et al. (1992). Initially, the shapes, amplitudes, and absolute photometry of the synthetic and observed lightcurves agree with each other to within 0.05-0.1 mag. Small modifications to the Gaspra shape model on sides of the asteroid poorly imaged by Galileo (changes of 700 m or less in the southern hemisphere at longitudes 90 deg-270 deg W) reduce the typical discrepancies to approximately 0.05 mag in lightcurve shape and less than 0.03 mag in absolute photometry. The result demonstrates that Earth-based lightcurves can be used to refine the shape of a spacecraft-imaged irregular object in areas that are poorly constrained by the spacecraft observations. The consistency and phase-angle dependence of the Galileo-based model for Gaspra photometry, supports the accuracy of the absolute calibration of the Galileo SSI camera, and confirms the Earth-based determination of the V-filter geometric albedo of the asteroid (0.22 +/- 0.03; Tholen et al., submitted for publication). Remaining discrepancies between the synthetic and observed lightcurves show no indication of systematic latitudinal variations in albedo and also cannot be explained entirely by isolated albedo spots. These discrepancies are most likely caused by (1) small, remaining, hard-to-constrain errors in the Gaspra shape model and/or (2) moderate variations in macroscopic roughness across the asteroid's surface, in particular making longitudes 130 deg to 300 deg W moderately rougher than the opposite hemisphere.

Simonelli, Damon P.↗

Examining the Pilot and Controller Performance Data When in a Free Flight with Weather Phenomenon

The present study investigated effects of weather related factors on the performance of pilots under free flight. A weather scenario was defined by a combination of precipitation factors (light rain, moderate rain, and heavy rain or snow), visibility (1,4,8 miles), wind conditions (light, medium, or heavy), cloud ceiling (800ft. below, 1800ft above, and 4000ft horizontal). The performance of the aircraft self-separation was evaluated in terms of detection accuracy and detection times for student- and commercial (expert) pilots. Overall, the results obtained from a behavioral analysis showed that in general, the ability to recognize intruder aircraft conflict incidents, followed by the ability to acquire the spatial location of the intruder aircraft relative to ownership aircraft were judged to be the major cognitive tasks as perceived by the participants during self-separation. Further, the participants rarely used cockpit display of traffic information (CDTI) during conflict management related to aircraft separation, but used CDTI highly during decision-making tasks. In all weather scenarios, there were remarkable differences between expert and student pilots in detection times. In summary, weather scenarios were observed to affect intruder aircraft detection performance accuracies. There was interaction effects between weather Scenario-1 and Scenario-2 for climbing task data generated by both expert- and student- pilots at high traffic density. Scenario-3 weather condition provided an opportunity for poor detection accuracy as well as detection time increase. This may be attributed to low visibility. The intruder aircraft detection times were not affected by the weather conditions during climbing and descending tasks. The decision of pilots to fly into certain weather condition was dependent in part on the warning distance to the location of the weather. When pilots were warned of the weather conditions, they were more likely to fly their aircraft into it, but mostly when the warning was not close to the weather location.

Nituen, Celestine A.↗

Enhanced Diffusive Transport in Fluctuating Porous Media

Mass transport within porous structures is a ubiquitous process in biological, geological, and technological systems. Despite the importance of these phenomena, there is no comprehensive theory that describes the complex and diverse transport behavior within porous environments. While the porous matrix itself is generally considered a static and passive participant, many porous environments are in fact dynamic, with fluctuating walls, pores that open and close, and dynamically changing cross-links. While diffusion has been measured in fluctuating structures, notably in model biological systems, it is rarely possible to isolate the effect of fluctuations because of the absence of control experiments involving an identical static counterpart, and it is generally impossible to observe the dynamics of the structure. In this paper, we present a direct comparison of the diffusion of nanoparticles of various sizes within a trackable, fluctuating porous matrix and a geometrically equivalent static matrix, in conditions spanning a range of regimes from obstructed to highly confined. The experimental system comprised a close-packed layer of colloidal spheres that were either immobilized to a planar surface or allowed to fluctuate locally, within the space defined by their nearest neighbors. Interestingly, the effective long-time diffusion coefficient was approximately 35–65% greater in the fluctuating porous matrix than in the static one (depending on the size of the nanoparticle probes), regardless of the geometric regime. This was explained by considering the enhancing effects of matrix fluctuations on the short-time diffusion coefficient and cooperative “gate-opening” motions of matrix particles and nanoparticle probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fuel reid vapor pressure level and ethanol content on stochastic preignition, effects at steady and unsteady engine operation

The present work investigates relations between fuel Reid vapor pressure (RVP) and biofuel (ethanol) content on stochastic preignition (SPI) at both sustained steady-state engine operation and following load transients. This work stems from in-field observations that automotive original equipment manufacturers have observed consistent seasonal increases in United States customer drivability complaints and warranty claims during September and October where SPI is suspected to be responsible. The seasonal timing of these events coincides with the United States seasonal fuel property changeover initiating on September 15 each year, where fuel RVP increases. To explore potential linkage between fuel RVP and SPI the present study employs engine SPI experiments coupled with laboratory spray measurements of fuels with RVPs of 8, 12, and 16 psi in both E10 (10% ethanol) and E25 (25% ethanol) fuels. Engine results are partitioned into fuel RVP and ethanol content effects on SPI in steady-state, sustained high-load engine operation and unsteady-state low- to high-load transitions, where off-engine spray vessel patternation and tip penetration results help to elucidate the observed fuel effects on SPI. A boosted direct-injected, spark-ignition engine was fueled with three market relevant E10 and E25 fuels with RVPs of 8, 12, and 16 to characterize the interplay between winter fuels and abnormal combustion behavior. The steady-state work shows that for high-load, steady-state engine operation, SPI is directly linked to fuel retention, which was found to be dependent on fuel distillation. The unsteady-state engine operation work shows that following low-to high-load transitions, SPI can occur from a memory of fuel property effects at low-load operation. Specifically, the fuel RVP effect on fuel spray collapse at low loads was found to correlate with SPI with a more than 95% confidence interval following low- to high-engine-load transitions. Results suggest that fuel-wall impingement at low-load operation could carry over into high-load transitions and generate SPI events following low- to high-load transitions.

09 BIOMASS FUELS↗

Paradigm Shift in Data Content and Informatics Infrastructure Required for Generalized Constitutive Modeling of Materials Behavior

Materials property information such as composition and thermophysical/mechanical properties abound in the literature. Oftentimes, however, the corresponding response curves from which these data are determined are missing or at the very least difficult to retrieve. Further, the paradigm for collecting materials property information has historically centered on (1) properties for materials comparison/selection purposes and (2) input requirements for conventional design/analysis methods. However, just as not all materials are alike or equal, neither are all constitutive models (and thus design/ analysis methods) equal; each model typically has its own specific and often unique required materials parameters, some directly measurable and others indirectly measurable. Therefore, the type and extent of materials information routinely collected is not always sufficient to meet the current, much less future, needs of the materials modeling community. Informatics has been defined as the science concerned with gathering, manipulating, storing, retrieving, and classifying recorded information. A key aspect of informatics is its focus on understanding problems and applying information technology as needed to address those problems. The primary objective of this article is to highlight the need for a paradigm shift in materials data collection, analysis, and dissemination so as to maximize the impact on both practitioners and researchers. Our hope is to identify and articulate what constitutes "sufficient" data content (i.e., quality and quantity) for developing, characterizing, and validating sophisticated nonlinear time- and history-dependent (hereditary) constitutive models. Likewise, the informatics infrastructure required for handling the potentially massive amounts of materials data will be discussed.

Arnold, S. M.↗

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗

Characterization of single-shot attosecond pulses with angular streaking photoelectron spectra

Most of the traditional attosecond pulse retrieval algorithms are based on a so-called attosecond streak camera technique, in which the momentum of the electron is shifted by an amount depending on the relative time delay between the attosecond pulse and the streaking infrared pulse. Thus, temporal information of the attosecond pulse is encoded in the amount of momentum shift in the streaked photoelectron momentum spectrogram S(p,τ), where p is the momentum of the electron along the polarization direction and τ is the time delay. An iterative algorithm is then employed to reconstruct the attosecond pulse from the streaking spectrogram. This method, however, cannot be applied to attosecond pulses generated from free-electron x-ray lasers where each single shot is different and stochastic in time. However, using a circularly polarized infrared laser as the streaking field, a two (or three)-dimensional angular streaking electron spectrum can be used to retrieve attosecond pulses for each shot, as well as the time delay with respect to the circularly polarized IR field. Here we show that a retrieval algorithm previously developed for the traditional streaking spectrogram can be modified to efficiently characterize single-shot attosecond pulses. The methods have been applied to retrieve 188 single shots from recent experiments. We analyze the statistical behavior of these 188 pulses in terms of pulse duration, bandwidth, pulse peak energy, and time delay with respect to the IR field. Furthermore, the retrieval algorithm is efficient and can be easily used to characterize a large number of shots in future experiments for attosecond pulses at free-electron x-ray laser facilities.

74 ATOMIC AND MOLECULAR PHYSICS↗

Computer Modeling of Thermal Convection in Melts to Explain Glass Formation in Low Gravity and on Earth

Experiments conducted up to this time on glass forming melts in the low gravity environment of space show that glasses prepared in low-g are more chemically homogeneous and more resistant to crystallization than the comparable glasses prepared at 1-g on Earth. This result is somewhat surprising and opposite to the accepted concept on glass formation for a melt. A hypothesis based on "shear thinning" of a melt, a decrease in viscosity with increasing shear stress, is proposed as an explanation for the observed low-gravity results. This paper describes detailed simulation procedures to test the role of thermal convection in introducing shear stress in glass forming melts, using a lithium disilcate melt as a model. The simulation system in its idealized version consists of a cylinder that is heated at one end and cooled at the other with gravity acting in a transverse direction to the thermal gradient. The side wall of the cylinder is assumed to be insulating. The governing equations of motion and energy are solved using variable properties for viscosity (Arrehenius and non-Arrehenius behaviors) and density (constant and temperature dependent). Other parametric variables in the calculations include gravity level and gravity vector orientation. The shear stress in the system are then computed as a function of gravity from the calculated values of maximum melt velocity, and its effect on melt viscosity (shear thinning) is predicted. Also included and discussed are the modeling efforts related to other potential convective processes in glass forming melts and their possible effects on melt viscosity.

Ray, Chandra S.↗

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

15 GEOTHERMAL ENERGY↗

Nature of polariton transport in a Fabry-Perot cavity

Fabry-Perot microcavities can strongly enhance interactions between light and molecules, leading to the formation of hybrid light-matter states known as polaritons. Polaritons possess much smaller effective masses and much larger group velocities when the molecules are resonant with cavity modes that have finite (nonzero) in-plane wave vectors, giving rise to the possibilities of long-range and ultrafast ballistic transport. In this paper, we present the results of numerical simulations of the ultrafast ballistic transport phenomenon in real space and time during and after initialization with a short, spatially localized pulse. Here, we address this problem with two approaches: the standard transfer-matrix method (TMM) for planar structures and a second simulation based on the numerical solution of the Maxwell-Bloch equations that can be used for general configurations (complex metasurfaces, for instance) and excitation modes. The agreement between the TMM and the full numerical calculation with the Maxwell-Bloch equations when applied to multilayer planar configurations provides proof of the validity of both approaches for the present analyses. Overall, we find that the transport of the molecular excitons inside the cavity synchronize with the evolution of the enhanced electromagnetic field inside the cavity. Moreover, the synchronized transport rate is in good agreement with the group velocities predicted from a calculated dispersion relation across a wide range of frequencies. Finally, we relate the group velocity to the Hopfield coefficient obtained from quantum modeling and suggest that the dependence of light-matter coupling on the in-plane wave vector can be an important but overlooked factor for understanding the transport behavior of polaritons. These simulations provide an intuitive tool for understanding the collective motion of light and excitons and helps us to better understand how experimental observations of polaritons should be interpreted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗