Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transition path”

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 883 records · Page 49

A rheological model for loose sands with insights from DEM

A rheological model for loose granular media is developed to capture both solid-like and fluid-like responses during shearing. The proposed model is built by following the mathematical structure of an extended Kelvin–Voigt model, where an elastic spring and plastic slider act in parallel to a viscous damper. This arrangement requires the partition of the total stress into rate-independent and rate-dependent stress components. To model the solid-like behavior, a simple frictional plasticity model is adopted without modifications, thus contributing to the rate-independent stress. Instead, the fluid-like or rate-dependent stress is further decomposed into deviatoric and volumetric parts, by proposing a new formulation based on a combination of the μ(I) relation, originally developed under pressure-controlled shear, with a pressure-shear rate relation derived under volume-controlled shear. The proposed formulation allows the model to capture both the increase in the friction coefficient and the enhanced dilation at high shear rates. High-fidelity simulation data, obtained from discrete element method and multiscale modelling, are used to evaluate the performance of the proposed constitutive model. The model provides accurate results under both drained and undrained simple shear paths across a wide range of shear rates. Furthermore, it successfully reproduces at much lower computational cost the flowslide mobility computed through multiscale simulations, which is primarily regulated by the shear rate dependence of the material properties during the dynamic runout stage.

Elasticity↗

Reducing On-Board Computer Propagation Errors Due to Omitted Geopotential Terms by Judicious Selection of Uploaded State Vector

Several future, and some current missions, use an on-board computer (OBC) force model that is very limited. The OBC geopotential force model typically includes only the J(2), J(3), J(4), C(2,2) and S(2,2) terms to model non-spherical Earth gravitational effects. The Tropical Rainfall Measuring Mission (TRMM), Wide-field Infrared Explorer (WIRE), Transition Region and Coronal Explorer (TRACE), Submillimeter Wave Astronomy Satellite (SWAS), and X-ray Timing Explorer (XTE) all plan to use this geopotential force model on-board. The Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX) is already flying this geopotential force model. Past analysis has shown that one of the leading sources of error in the OBC propagated ephemeris is the omission of the higher order geopotential terms. However, these same analyses have shown a wide range of accuracies for the OBC ephemerides. Analysis was performed using EUVE state vectors that showed the EUVE four day OBC propagated ephemerides varied in accuracy from 200 m. to 45 km. depending on the initial vector used to start the propagation. The vectors used in the study were from a single EUVE orbit at one minute intervals in the ephemeris. Since each vector propagated practically the same path as the others, the differences seen had to be due to differences in the inital state vector only. An algorithm was developed that will optimize the epoch of the uploaded state vector. Proper selection can reduce the previous errors of anywhere from 200 m. to 45 km. to generally less than one km. over four days of propagation. This would enable flight projects to minimize state vector uploads to the spacecraft. Additionally, this method is superior to other methods in that no additional orbit estimates need be done. The definitive ephemeris generated on the ground can be used as long as the proper epoch is chosen. This algorithm can be easily coded in software that would pick the epoch within a specified time range that would minimize the OBC propagation error. This techniques should greatly improve the accuracy of the OBC propagation on-board future spacecraft such as TRMM, WIRE, SWAS, and XTE without increasing complexity in the ground processing.

Greatorex, Scott↗

Charged Bosons Made of Fermions in Bilayer Structures with Strong Metallic Screening

Two-dimensional monolayer structures of transition metal dichalogenides (TMDs) have been shown to allow many higher-order excitonic bound states, including trions (charged excitons), biexcitons (excitonic molecules), and charged biexcitons. Here, we report here experimental evidence and the theoretical basis for a new bound excitonic complex, consisting of two free carriers bound to an exciton in a bilayer structure. Our experimental measurements on structures made using two different materials show a new spectral line at the predicted energy with two different TMD materials (MoSe 2 and WSe 2 ) with both n- and p-doping, if and only if all the required theoretical conditions for this complex are fulfilled, in particular, only in the presence of a parallel metal layer that significantly screens the repulsive interaction between the like-charge carriers. Because these four-carrier bound states are charged bosons, they could eventually be the basis for a new path to superconductivity without Cooper pairing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optical vortex manipulation for topological quantum computation

Topological quantum computation based on Majorana bound states may enable new paths to fault-tolerant quantum computing. Several recent experiments have suggested that the vortex cores of topological superconductors, such as iron-based superconductors, may host Majorana bound states at zero energy. However, quantum computation with these zero-energy vortex bound states requires precise and fast manipulation of individual vortices, which is difficult to do in a scalable manner. To address this issue, in this study we propose a control scheme based on local heating via, for example, scanning optical microscopy to braid vortex-bound Majorana zero modes in a two-dimensional topological superconductor. First, we derive the conditions required for transporting a single vortex between two defects in the superconducting material by trapping it with a hot spot generated by local optical heating. Equipped with critical conditions for the vortex motion, we then establish the ideal material properties for vortex braiding and describe how transition errors resulting from finite speed and/or temperature can be minimized. Our work paves the way toward optical or microscopic control of zero-energy vortex bound states in two-dimensional topological superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Digital Twin Model for Advanced Manufacture of a Rotating Detonation Engine Injector

A digital twin material model (DTMM) of an additive manufacturing (AM) process was created to advance the state of the art in rotating detonation engine (RDE) injector design. Current RDE injectors are designed with large pressure drops, enabling a stable and repeatable combustion process. However, this comes at the cost of system efficiency. For the technology to transition to commercial fossil-based power generation, it is important to develop injectors with reduced flow losses. Low-loss injectors are difficult to design and manufacture with conventional manufacturing techniques. AM enables new design options, but the AM manufacturing process must be thoroughly understood to result in a robust design. A DTMM provides the necessary insight by defining the cause-effect relationships between process parameters, microstructure features, and properties. Therefore, a DTMM to support the design and manufacturing process was developed and applied to the design of a new additively manufactured low-loss injector. The injector combustion behavior was characterized through hot-fire tests, and mechanical performance was compared to the DTMM predictions. The two project goals were the successful development of the DTMM and the demonstration of an improved RDE injector design. The RDE injector design and DTMM developments occurred on parallel but dependent paths. The injector was designed to reduce pressure drop by increasing the cross-sectional flow area ratio between the injector air passages and the combustor annulus. This resulted in less structural material, raising the concern that thin members would be susceptible to high-cycle fatigue (HCF) under the periodic loading inherent to an RDE. It was most important for the DTMM to predict behavior in these features; therefore, the injector design concept guided the material thicknesses used in fatigue tests. The DTMM development started by manufacturing a series of coupons over the range of possible AM process variations. A design-of-experiment approach was used to select which process variable combinations gave the most efficient coverage relevant to the injector design space. The microstructure in each of these coupons was characterized, and then computational methods were used to create a numerical model of the correlation between process variables and microstructure. Next, a set of HCF samples were tested to calibrate existing models that map microstructure to HCF performance. Together, these two links formed the DTMM that calculates HCF behavior from AM process variables. Two injector prototypes were additively manufactured. The first injector design strategy aggressively pursued low-loss performance by substantially increasing the oxidizer flow area. The combination of manufacturing lead times and the fatigue testing schedule meant that the DTMM was not available when building this first prototype. Therefore, its process parameters were chosen based on a manual review of the available coupon data. This prototype was built successfully and evaluated in 58 combustion tests. Sustained detonation was achieved with remarkably reduced pressure loss, and some tests even displayed pressure loss characteristics similar to conventional gas turbine combustors. This achieved the project goal of improving RDE injector design. The second injector was manufactured according to the optimized parameters predicted by the DTMM. The flow area modifications of this injector were less aggressive than the first injector since demonstrating low pressure loss was not an objective of the second hot-fire test series. Rather, the test objective was to cause high cycle fatigue failure in the part due to periodic loading from the rotating detonation wave. The observed number of cycles to failure was to be compared to the number predicted by the DTMM and thereby assess the utility of the DTMM in component design. However, the required level of vibration was not obtained during combustion. Therefore, high cycle fatigue was not experienced in the hot-fire tests of the second injector. Fatigue data was obtained by further testing the second injector in a conventional HCF test apparatus. The injector demonstrated HCF strength above the DTMM prediction. In fact, it did not fail and testing was only discontinued due to reaching the end of the period of performance. This points to some success in the project’s primary goal of successfully developing and applying the DTMM to a component design. Implementing the DTMM recommendations for optimal processing parameters led to a part with acceptable properties. The DTMM was also shown to be an efficient correlator of data and to provide insight into the relationship between process settings, microstructure, and property performance. However, the failure of the DTMM prediction to match the experimental result of the injector fatigue test also points to the need to include significantly more data in the model development. In this project, coupons made with identical processing parameters exhibited drastically different properties from each other and from the injector part, which clearly influences the accuracy of a model that predicts performance based on parameters. Uncertainties in the build process must be quantified to develop more robust models. A denser and broader matrix of coupon process and geometry variations, several repeated builds of every point, more in-situ build process measurements, and direct observation of tensile and HCF sample microstructure (as opposed to separate microstructure specimens) are recommendations to improve future AM modeling efforts.

20 FOSSIL-FUELED POWER PLANTS↗

Assessment of Cloud-Based Applications Enabling a Scalable Risk-Informed Predictive Maintenance Strategy Across the Nuclear Fleet

The current light water reactor fleet uses time-based or failure-based maintenance strategies to achieve high-capacity factors. But to make nuclear more competitive in the energy market, these reactors could utilize emerging technologies in terms of artificial intelligence (AI) and cloud computing to enable a cost-effective, predictive maintenance strategy. This report examines the feasibility of cloud computing for the nuclear industry’s needs in terms of the cloud’s computing capabilities, feasibility, and regulatory concerns. The technical viability of cloud computing was analyzed using one year worth of data from a boiling water reactor’s safety relief valve. Models were hosted on a local desktop, Idaho National Laboratory’s high-performance computer, and Microsoft Azure. Data was loaded, processed, and two types of models were trained in an A/B fashion. Based on the speed at which these actions were completed, it was used to determined that cloud computing has adequate computing resources. Additionally, the computing power can scale with the demanded load. To enable cloud computing in the existing fleet, additional sensors, networks, and other requirements must be implemented to ensure a smooth transition from current maintenance strategies. However, there is a benefit as the plant no longer needs manage their own servers, software, cybersecurity, and IT support staff. Many of these features can be offloaded on to the cloud provider. A comprehensive analysis was completed that showed the current annual cost of operating is more expensive than using cloud computing resources. Lastly, the regulatory framework does not explicitly address AI or autonomous control. Currently, the NRC and other regulatory bodies are evaluating providing guidance to address gaps rather than new regulations to address the use of AI and ML. But since many of the AI applications are focused on non-safety related applications, such as balance-of-plant components, they will likely have little or no regulatory restrictions or necessary approvals. Demonstrating how AI can improve maintenance and operation of these non-safety related systems seems like the likely path forward for implementing AI and cloud computing resources inside nuclear power plants (NPPs).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Environmental fatigue of an Al-Li-Cu alloy. Part 2: Microscopic hydrogen cracking processes

Based on a fractographic analysis of fatigue crack propagation (FCP) in Al-Li-Cu alloy 2090 stressed in a variety of inert and embrittling environments, microscopic crack paths are identified and correlated with intrinsic da/dN-delta K kinetics. FCP rates in 2090 are accelerated by hydrogen producing environments (pure water vapor, moist air, and aqueous NaCl), as defined in Part 1. For these cases, subgrain boundary fatigue cracking (SGC) dominates for delta K values where the crack tip process zone, a significant fraction of the cyclic plastic zone, is sufficiently large to envelop 5 micron subgrains in the unrecrystallized microstructure. SGC may be due to strong hydrogen trapping at T1 precipitates concentrated at sub-boundaries. At low delta K, the plastic zone diameter is smaller than the subgrain size and FCP progresses along (100) planes due to either local lattice decohesion or aluminum-lithium hydride cracking. For inert environments (vacuum, helium, and oxygen), or at high delta K where the hydrogen effect on da/dN is small, FCP is along (111) slip planes; this mode does not transition with increasing delta K and plastic zone size. The SGC and (100) crystallographic cracking modes, and the governing influence of the crack tip process zone volume (delta K), support hydrogen embrittlement rather than a surface film rupture and anodic dissolution mechanism for environmental FCP. Multi-sloped log da/dN-log delta K behavior is produced by changes in process zone hydrogen-microstructure interactions, and not by purely micromechanical-microstructure interactions, in contradiction to microstructural distance-based fatigue models.

Robert S Piascik↗

Towards controlled synthesis of 2D crystals by chemical vapor deposition (CVD)

We report the emergence of two-dimensional (2D) materials has captured the imagination of researchers since graphene was first exfoliated from graphite in 2004. Their exotic properties give rise to many exciting potential applications in advanced electronic, optoelectronic, energy and biomedical technologies. Scalable growth of high quality 2D materials is crucial for their adoption in technological applications the same way the arrival of high quality silicon single crystals was to the semiconductor industry. A huge amount of effort has been devoted to grow large-area, highly crystalline 2D crystals such as graphene and transition metal dichalcogenides (TMDs) through various methods. While CVD growth of wafer-scale monolayer graphene and TMDs has been demonstrated, considerable challenges still remain. In this perspective, we advocate for the focus on the crystal growth morphology as an underpinning for understanding, diagnosing and controlling the CVD process and environment for 2D material growth. Like snowflakes in nature, 2D crystals exhibit a rich variety of morphologies under different growth conditions. The mapping of crystal shapes in the growth parameter space “encodes” a wealth of information, the deciphering of which will lead to better understanding of the fundamental growth mechanism and materials properties. To this end, we envision a collective effort by the 2D materials community to establish the correlation between crystal shapes and the intrinsic thermodynamic and kinetic parameters for CVD reactions through integrated crystal growth experiment, database development and machine learning assisted predictive modeling, which will pave a robust path towards controlled synthesis of 2D materials and heterostructures.

36 MATERIALS SCIENCE↗

Phonon spectrum in the spin-Peierls phase of CuGeO 3

CuGeO 3 has long been studied as a prototypical example of the spin-Peierls transition in a 𝑆 = 1/2 Heisenberg chain. Despite intensive investigation of this quasi-one-dimensional material, systematic measurements and calculations of the phonon excitations in the dimerized phase have not to date been possible, leaving certain aspects of the spin-Peierls phenomenon unresolved. We perform state-of-the-art density functional theory (DFT) calculations to compute the electronic structure and phonon dynamics in the low-temperature dimerized phase. We also perform high-resolution neutron spectroscopy to measure the full phonon spectrum over multiple Brillouin zones. We find excellent agreement between our numerical and experimental results that extend to all measurement temperatures. Notable features of our phonon spectra include a number of steeply dispersive modes, nonmonotonic dispersion features, and specific phonon anticrossings, which we relate to the mode eigenvectors. By calculating the magnetic interactions within DFT and studying the effects of different phonon modes on the superexchange paths, we discuss the possibility of observing spin-phonon hybridization effects in experiments performed both in and out of equilibrium.

density functional theory↗

Zero-point motion of liquid and solid hydrogen

We present an inelastic neutron scattering study of liquid and solid hydrogen carried out using the wide Angular Range Chopper Spectrometer at Oak Ridge National Laboratory. From the observed dynamic structure factor, we obtained empirical estimates of the molecular mean-squared displacement and average translational kinetic energy. Here we find that the former quantity increases with temperature, indicating that a combination of thermal and quantum effects is important near the liquid-solid phase transition, contrary to previous measurements. We also find that the kinetic energy drops dramatically on melting of the crystals, a consequence of the large increase in molar volume together with the Heisenberg indeterminacy principle. Our results are compared with quantum Monte Carlo simulations based on different model potentials. In general, there is good agreement between our findings and theoretical predictions based on the Silvera-Goldman and Buck potentials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Equation-of-state, sound speed, and reshock of shock-compressed fluid carbon dioxide

Mechanical equation-of-state data of initially liquid and solid CO2 shock-compressed to terapascal conditions are reported. Diamond-sapphire anvil cells were used to vary the initial density and state of CO2 samples that were then further compressed with laser-driven shock waves, resulting in a data set from which precise derivative quantities, including Grüneisen parameter and sound speed, are determined. Reshock states are measured to 800 GPa and map the same pressure-density conditions as the single shock using different thermodynamic paths. The compressibility data reported here do not support current density-functional-theory calculations, but are better represented by tabular equation-of-state models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Moving-Base Simulation Evaluation of Control/Display Integration Issues for ASTOVL Aircraft

A moving-base simulation has been conducted on the Vertical Motion Simulator at Ames Research Center using a model of an advanced, short takeoff and vertical landing (STOVL) lift fan fighter aircraft. This experiment expanded on investigations during previous simulations with this STOVL configuration with the objective of evaluating (1) control law modifications over the low speed flight envelope, (2) integration of the throttle inceptor with flight control laws that provide direct thrust command for conventional flight, vertical and short takeoff, and flightpath or vertical velocity command for transition, hover, and vertical landing, (3) control mode blending for pitch, roll, yaw, and flightpath control during transition from wing-borne to jet-borne flight, and (4) effects of conformal versus nonconformal presentation of flightpath and pursuit guidance symbology on the out-the-window display for low speed STOVL operations. Assessments were made for takeoff, transition, hover, and landing, including precision hover and landing aboard an LPH-type amphibious assault ship in the presence of winds and rough seas. Results yielded Level 1 pilot ratings for the flightpath and vertical velocity command modes for a range of land-based and shipboard operation and were consistent with previous experience with earlier control laws and displays for this STOVL concept. Control mode blending was performed over speed ranges in accord with the pilot's tasks and with the change of the basic aircraft's characteristics between wing-borne and hover flight. Blending of yaw control from heading command in hover to sideslip command in wing-borne flight performed over a broad speed range helped reduce yaw transients during acceleration through the low speed regime. Although the pilots appreciated conformality of flightpath and guidance symbols with the external scene during the approach, increased sensitivity of the symbols for lateral path tracking elevated the pilots' control activity in the presence of turbulence. The pilots preferred the choice of scaling that was originally established during the display development and in-flight evaluations.

Franklin, James A.↗

Temperature-induced densification in compressed basaltic glass revealed by in-situ ultrasonic measurements

Abstract Acoustic velocities of a model basalt glass (64 mol% CaMgSi2O6 + 36 mol% CaAl2Si2O8) were measured along different pressure-temperature (P-T) paths. One set of experiments involved isothermal compression-decompression cycles, performed at temperatures of 300, 641, 823, and 1006 K and pressures up to 12.2 GPa. The other set of experiments involved constant-load heating-cooling cycles at temperatures up to 823 K and pressures up to 7.5 GPa. Both sets of experiments were performed in a multi-anvil apparatus using a synchrotron-based ultrasonic technique. Our results show that the glass compressed isothermally at 300 K (cold-compression) displays anomalously decreasing compressional (VP) and shear (VS) wave velocities with increasing pressure until ~8 GPa. Beyond 8 GPa, both VP and VS start to increase sharply with pressure and irreversible densification of the glass occurred, producing large hysteresis loops of velocities upon decompression. However, for the glass compressed isothermally at increasingly higher temperatures (hot-compression), the velocity minima gradually shift to lower pressures. At temperature close to the glass transition temperature Tg, the velocity minima disappear completely, displaying a monotonic increase of velocities during compression and higher VP and VS during decompression. In addition, constant-load heating-cooling experiments show that velocities generally decrease slightly with increasing temperature, but start to increase once heated above a threshold temperature (~650 K). During cooling the velocities increase almost linearly with decreasing temperature, resulting in higher velocities (~1.5–2.5% higher) when returned to 300 K. This implies that a temperature-induced densification may have occurred in the glass at high pressures. Raman spectra on recovered samples show that the hot-compressed and high-P heated glasses contain distinctly densified and depolymerized structural signatures compared to the initial glass and the cold-compressed glass below the velocity transition pressure PT (~8 GPa). Such densification may be attributed to the breaking of bridging oxygen bonds and compaction in the intermediate-range structure. Our results demonstrate that temperature can facilitate glass densification at high pressures and point out the importance of P-T history in understanding the elastic properties of silicate glasses. Comparison with melt velocity suggests that hot-compressed glasses may better resemble the pressure dependence of velocity of silicate melts than cold-compressed glasses, but still show significantly higher velocities than melts. If the abnormal acoustic behaviors of cold-compressed glasses were used to constrain melt fractions in the mantle low-velocity regions, the melt fractions needed to explain a given velocity reduction would be significantly underestimated at high pressures.

Geochemistry & Geophysics↗

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gelation and Re-entrance in Mixtures of Soft Colloids and Linear Polymers of Equal Size

Liquid mixtures composed of colloidal particles and much smaller non-adsorbing linear homopolymers can undergo a gelation transition due to polymer-mediated depletion forces. We now show that the addition of linear polymers to suspensions of soft colloids having the same hydrodynamic size yields a liquid-togel-to-re-entrant liquid transition. In particular, the dynamic state diagram of 1,4-polybutadiene star-linear polymer mixtures was determined with the help of linear viscoelastic and small-angle Xray scattering experiments. While keeping the star polymers below their nominal overlap concentration, a gel was formed upon increasing the linear polymer content. Further addition of linear chains yielded a re-entrant liquid. This unexpected behavior was rationalized by the interplay of three possible phenomena: (i) depletion interactions, driven by the size disparity between the stars and the polymer length scale which is the mesh size of its entanglement network; (ii) colloidal deswelling due to the increased osmotic pressure exerted onto the stars; and (iii) a concomitant progressive suppression of the depletion efficiency on increasing the polymer concentration due to reduced mesh size, hence a smaller range of attraction. Our results unveil an exciting new way to tailor the flow of soft colloids and highlight a largely unexplored path to engineer soft colloidal mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The advanced tokamak path to a compact net electric fusion pilot plant

Abstract Physics-based simulations project a compact net electric fusion pilot plant with a nuclear testing mission is possible at modest scale based on the advanced tokamak concept, and identify key parameters for its optimization. These utilize a new integrated 1.5D core-edge approach for whole device modeling to predict performance by self-consistently applying transport, pedestal and current drive models to converge fully non-inductive stationary solutions, predicting profiles and energy confinement for a given density. This physics-based approach leads to new insights and understanding of reactor optimization. In particular, the levering role of high plasma density is identified, which raises fusion performance and self-driven ‘bootstrap currents’, to reduce current drive demands and enable high pressure with net electricity at a compact scale. Solutions at 6–7 T, ∼4 m radius and 200 MW net electricity are identified with margins and trade-offs possible between parameters. Current drive comes from neutral beam and ultra-high harmonic (helicon) fast wave, though other advanced approaches are not ruled out. The resulting low recirculating power in a double null configuration leads to a divertor heat flux challenge that is comparable to ITER, though reactor solutions may require more dissipation. Strong H-mode access (x2 margin over L–H transition scalings) and ITER-like heat fluxes are maintained with ∼20%–60% core radiation, though effects on confinement need further analysis. Neutron wall loadings appear tolerable. The approach would benefit from high temperature superconductors, as higher fields would increase performance margins while potential for demountability may facilitate nuclear testing. However, solutions are possible with conventional superconductors. An advanced load sharing and reactive bucking approach in the device centerpost region provides improved mechanical stress handling. The prospect of an affordable test device which could close the loop on net-electric production and conduct essential nuclear materials and breeding research is compelling, motivating research to validate the techniques and models employed here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Boundary layer receptivity and control

Receptivity processes initiate natural instabilities in a boundary layer. The instabilities grow and eventually break down to turbulence. Consequently, receptivity questions are a critical element of the analysis of the transition process. Success in modeling the physics of receptivity processes thus has a direct bearing on technological issues of drag reduction. The means by which transitional flows can be controlled is also a major concern: questions of control are tied inevitably to those of receptivity. Adjoint systems provide a highly effective mathematical method for approaching many of the questions associated with both receptivity and control. The long term objective is to develop adjoint methods to handle increasingly complex receptivity questions, and to find systematic procedures for deducing effective control strategies. The most elementary receptivity problem is that in which a parallel boundary layer is forced by time-harmonic sources of various types. The characteristics of the response to such forcing form the building blocks for more complex receptivity mechanisms. The first objective of this year's research effort was to investigate how a parallel Blasius boundary layer responds to general direct forcing. Acoustic disturbances in the freestream can be scattered by flow non-uniformities to produce Tollmien-Schlichting waves. For example, scattering by surface roughness is known to provide an efficient receptivity path. The present effort is directed towards finding a solution by a simple adjoint analysis, because adjoint methods can be extended to more complex problems. In practice, flows are non-parallel and often three-dimensional. Compressibility may also be significant in some cases. Recent developments in the use of Parabolized Stability Equations (PSE) offer a promising possibility. By formulating and solving a set of adjoint parabolized equations, a method for mapping the efficiency with which external forcing excites the three-dimensional motions of a non-parallel boundary layer was developed. The method makes use of the same computationally efficient formulation that makes the PSE currently so appealing. In the area of flow control, adjoint systems offer a powerful insight into the effect of control forces. One of the simplest control strategies for boundary layers involves the application of localized mean wall suction.

Hill, D. C.↗

Toward nonthermal control of excited quantum materials: framework and investigations by ultrafast electron scattering and imaging

Quantum material systems upon applying ultrashort laser pulses provide a rich platform to access excited material phases and their transformations that are not entirely like their equilibrium counterparts. The addressability and potential controls of metastable or long-trapped out-of-equilibrium phases have motivated interests both for the purposes of understanding the nonequilibrium physics and advancing the quantum technologies. Thus far, the dynamical spectroscopic probes eminently focus on microscopic electronic and phonon responses. For characterizing the long-range dynamics, such as order parameter fields and fluctuation effects, the ultrafast scattering probes offer direct sensitivity. Bridging the connections between the microscopic dynamics and macroscopic responses is central toward establishing the nonequilibrium physics behind the light-induced phases. Here, we present a path toward such understanding by cross-examining the structure factors associated with different dynamical states obtained from ultrafast electron scattering, imaging, and modeling. We give the basic theoretical framework on describing the non-equilibrium scattering problems and briefly describe how such framework relates to the out-of-equilibrium phenomena. We give effective models outlining the emergences of nonthermal critical points, hidden phases, and non-equilibrium relaxational responses from vacuum-suspended rare-earth tritellurides, tantalum disulfides thin films, and vanadium dioxide nanocrystalline materials upon light excitations.

79 ASTRONOMY AND ASTROPHYSICS↗