Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “projected gradient method”

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 379 records · Page 21

IRNet: A General Purpose Deep Risidual Regression Framework for Materials Discovery

Materials discovery is crucial for making scientific advances in many domains. Collections of data from experiments and first-principle computations have spurred interest in applying machine learning methods to create predictive models capable of mapping from composition and crystal structures to materials properties. Generally, these are regression problems with the input being a 1D vector composed of numerical attributes representing the material composition and/or crystal structure. While neural networks consisting of fully connected layers have been applied to such problems, their performance often suffers from the vanishing gradient problem when network depth is increased. Hence, predictive modeling for such tasks has been mainly limited to traditional machine learning techniques such as Random Forest. In this paper, we study and propose design principles for building deep regression networks composed of fully connected layers with numerical vectors as input. We introduce a novel deep regression network with individual residual learning, IRNet, that places shortcut connections after each layer so that each layer learns the residual mapping between its output and input. We use the problem of learning properties of inorganic materials from numerical attributes derived from material composition and/or crystal structure to compare IRNet’s performance against that of other machine learning techniques. Using multiple datasets from the Open Quantum Materials Database (OQMD) and Materials Project for training and evaluation, we show that IRNet provides significantly better prediction performance than the state-of-the-art machine learning approaches currently used by domain scientists. We also show that IRNet’s use of individual residual learning leads to better convergence during the training phase than when shortcut connections are between multi-layer stacks while maintaining the same number of parameters.

deep learning↗

Comparisons Between Array Derived Dynamic Strain Rate (ADDS) and Fiber‐Optic Distributed Acoustic Sensing (DAS) Strain Rate

Abstract Distributed acoustic sensing (DAS) strain rate and particle velocity can be compared through approximate scaling with medium velocity. We instead performed a direct comparison between array derived dynamic strain (ADDS) rate and DAS strain rate for six frequency bands. The PoroTomo project at Brady's Hot Springs, Nevada, deployed a 240‐geophone 3C array co‐located with fiber‐optic DAS system and 8.7 km of buried cable. We selected subsets of the geophone array to create four smaller arrays and computed ADDS. The horizontal components of the ADDS were rotated into the direction of the fiber‐optic cable and then compared with the observed DAS strain rates. From three example regional earthquakes of local magnitudes 2.9, 4.1, and 4.3, the ADDS are found to be coherent with DAS for frequencies ≤1 Hz. For frequencies >1‐Hz, this correlation decays quickly. Small differences between linear and areal dynamic strains at 1‐Hz suggest poor signal‐to‐noise or localized strain that is perturbed by shallow heterogeneities compare to the average strain propagating across the geophone array. The implication is that around 1‐Hz, straight fiber DAS is measuring axial strain along the fiber and can provide good approximations to translational particle motions. However, above 1‐Hz, DAS becomes more sensitive to shallow velocity gradients that can be beneficial for geophysical imaging yet becomes a limitation for traditional seismic analysis methods depending on absolute amplitude and phase from translational particle motions.

58 GEOSCIENCES↗

Back to the future: using historical climate variation to project near-term shifts in habitat suitable for coast redwood

Studies that mimic the effect of climate change on terrestrial ecosystems often use climate projections from downscaled global climate models (GCMs). These simulations are generally too coarse to capture patterns of fine-scale climate variation, such as the sharp coastal energy and moisture gradients associated with wind-driven upwelling of cold water. Coastal upwelling may limit future increases in coastal temperatures, compromising GCMs' ability to provide realistic scenarios of future climate in these coastal ecosystems. Taking advantage of naturally occurring variability in the high-resolution historic climatic record, we developed multiple fine-scale scenarios of California climate that maintain coherent relationships between regional climate and coastal upwelling. We compared these scenarios against coarse resolution GCM projections at a regional scale to evaluate their temporal equivalency. We utilized these historically based scenarios to estimate potential suitable habitat for coast redwood (Sequoia sempervirens D. Don) under 'normal' combinations of temperature and precipitation, and under anomalous combinations representative of potential future climates. We found that a scenario of warmer temperature with historically normal precipitation is equivalent to climate projected by GCMs for California by 2020-2030 and that under these conditions, climatically suitable habitat for coast redwood significantly contracts at the southern end of its current range. Our results indicate that historical climate data provide a high-resolution alternative to downscaled GCM outputs for near-term ecological forecasts. This method may be particularly useful in other regions where local climate is strongly influenced by ocean-atmosphere dynamics that are not represented by coarse-scale GCMs.

59 BASIC BIOLOGICAL SCIENCES↗

International Space Station Radiation Shielding Model Development

The projected radiation levels within the International Space Station (ISS) have been criticized by the Aerospace Safety Advisory Panel in their report to the NASA Administrator. Methods for optimal reconfiguration and augmentation of the ISS shielding are now being developed. The initial steps are to develop reconfigurable and realistic radiation shield models of the ISS modules, develop computational procedures for the highly anisotropic radiation environment, and implement parametric and organizational optimization procedures. The targets of the redesign process are the crew quarters where the astronauts sleep and determining the effects of ISS shadow shielding of an astronaut in a spacesuit. The ISS model as developed will be reconfigurable to follow the ISS. Swapping internal equipment rack assemblies via location mapping tables will be one option for shield optimization. Lightweight shield augmentation materials will be optimally fit to crew quarter areas using parametric optimization procedures to minimize the augmentation shield mass. The optimization process is being integrated into the Intelligence Synthesis Environment s (ISE s) immersive simulation facility at the Langley Research Center and will rely on High Performance Computing and Communication (HPCC) for rapid evaluation of shield parameter gradients.

Qualls, G. D.↗

Angular momentum desaturation for Skylab using gravity gradient torques

An angular momentum desaturation method for momentum exchange devices of orbiting spacecraft is described. The specific application of the method is to the Skylab which contains three double-gimbaled control moment gyros for precise attitude control and maneuvering. It is assumed that the attitude reference is inertially fixed and that two of the vehicle principal moments of inertia are much larger than the third. Gravity gradient torques and resultant angular momentum accumulation are developed for small deviations from the reference. The assumed moment-of-inertia distribution allows desaturation about all axes with only two attitude angles each for the two axes with large moments of inertia. The necessary desaturation maneuvers can be decoupled for a special set of orbital coordinates. All maneuvers are made during the night portion of the orbit, and the percentage utilized for desaturation is selectable. Expressions for the attitude angle commands are developed assuming infinite vehicle rates. The effect of finite rates introduces an efficiency into the desaturation. Expressions for this efficiency are developed and means for compensation are treated. Arbitrary misalignments between geometric vehicle axes and principal moment-of-inertia axes are permissible. An angle bias about the sun line minimizes the angular momentum accumulation about the sun line projection into the orbital plane. Adaptive desaturation maneuver limiting consistent with the available maneuver momentum is included.

Kennel, H. F.↗

Closures for Two-fluid Five-moment Equations for Magnetized Plasmas (Final Scientific Report)

In the proposed work, we will develop a complete set of closures for two-fluid five-moment equations for magnetized plasmas. Specifically, we propose to obtain explicit formulas for electrons and ions to express heat flow, viscosity, frictional force, and collisional heating in terms of density, and temperature, flow velocity. For this purpose, we solve a set of general moment equations for the relevant closure moments. An infinite hierarchy of general moment equations is completely equivalent to the original Landau-Fokker-Planck kinetic equation. In the small gyro-period expansion, the parallel components of the general moment equations agree with the moments of the drift-kinetic equation (DKE). The advantages of solving the moment equations include (i) an exact treatment of the collision operator, (ii) closure and transport relations with no flux-surface average, and (iii) a unified form of closures for general collisionality. Practically we solve a truncated set of moment equations and hence must verify the convergence of the solutions by increasing number of moments. The parallel integral (nonlocal) closures obtained by solving a set of linearized parallel moment equations are useful only when the nonlinear coupling of moments to the temperature and magnetic field gradient terms is negligible. When the nonlinear coupling terms are substantial, we propose a Fourier transform method to solve the parallel moment equations. By Fourier-expanding the magnetic field, we can convert the differential equations of parallel moments to linear algebraic equations for Fourier components of the moments. The convergence of the solutions can be confirmed by increasing the number of Fourier modes. Our preliminary study in axisymmetric circular magnetic geometry with a large aspect ratio shows that (i) the solutions for Fourier modes n=0 and n=1 converge when modes up to n=2 are included in the expansion, and (ii) lower collisionality requires more moments for convergent closures, as in the case of the integral closures. The proposed work will start from this simple magnetic geometry to find general parallel closures and then be generalized to explore a non-circular magnetic field with finite aspect-ratio. With parallel moments known, the perpendicular component of moment closures can be obtained by inverting the generalized cross product operator in the general moment equations. These efforts will be performed in conjunction with several multi-institutional projects, NIMROD, the Center for Extended MHD Modeling (CEMM), and the Plasma Science and Innovation (PSI) Center.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗

Describing rainfall in northern Australia using multiple climate indices

Savanna landscapes are globally extensive and highly sensitive to climate change, yet the physical processes and climate phenomena which affect them remain poorly understood and therefore poorly represented in climate models. Both human populations and natural ecosystems are highly susceptible to precipitation variation in these regions due to the effects on water and food availability and atmosphere–biosphere energy fluxes. Here we quantify the relationship between climate phenomena and historical rainfall variability in Australian savannas and, in particular, how these relationships changed across a strong rainfall gradient, namely the North Australian Tropical Transect (NATT). Climate phenomena were described by 16 relevant climate indices and correlated against precipitation from 1900 to 2010 to determine the relative importance of each climate index on seasonal, annual and decadal timescales. Precipitation trends, climate index trends and wet season characteristics have also been investigated using linear statistical methods. In general, climate index–rainfall correlations were stronger in the north of the NATT where annual rainfall variability was lower and a high proportion of rainfall fell during the wet season. This is consistent with a decreased influence of the Indian–Australian monsoon from the north to the south. Seasonal variation was most strongly correlated with the Australian Monsoon Index, whereas yearly variability was related to a greater number of climate indices, predominately the Tasman Sea and Indonesian sea surface temperature indices (both of which experienced a linear increase over the duration of the study) and the El Niño–Southern Oscillation indices. These findings highlight the importance of understanding the climatic processes driving variability and, subsequently, the importance of understanding the relationships between rainfall and climatic phenomena in the Northern Territory in order to project future rainfall patterns in the region.

54 ENVIRONMENTAL SCIENCES↗

A Measurement of Atmospheric Circular Polarization with POLARBEAR

At millimeter wavelengths, the atmospheric emission is circularly polarized owing to the Zeeman splitting of molecular oxygen by the Earth's magnetic field. We report a measurement of the signal in the 150 GHz band using 3 yr of observational data with the POLARBEAR project. Nonidealities of a continuously rotating half-wave plate (HWP) partially convert circularly polarized light to linearly polarized light. While POLARBEAR detectors are sensitive to linear polarization, this effect makes them sensitive to circular polarization. Although this was not the intended use, we utilized this conversion to measure circular polarization. We reconstruct the azimuthal gradient of the circular polarization signal and measure its dependency from the scanning direction and the detector bandpass. We compare the signal with a simulation based on atmospheric emission theory, the detector bandpass, and the HWP leakage spectrum model. We find the ratio of the observed azimuthal slope to the simulated slope is 0.92 ± 0.01(stat) ± 0.07(sys). This ratio corresponds to a brightness temperature of 3.8 mK at the effective band center of 121.8 GHz and bandwidth of 3.5 GHz estimated from representative detector bandpass and the spectrum of Zeeman emission. This result validates our understanding of the instrument and reinforces the feasibility of measuring the circular polarization using the imperfection of continuously rotating HWP. Continuously rotating HWP is popular in ongoing and future cosmic microwave background experiments to modulate the polarized signal. This work shows a method for signal extraction and leakage subtraction that can help measure circular polarization in such experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Transverse Emittance Reduction with Tapered Foil

The idea of reducing transverse emittance with tapered energy-loss foil is proposed by J.M. Peterson in 1980s and recently by B. Carlsten. In this paper, we present the physical model of tapered energy-loss foil and analyze the emittance reduction using the concept of eigen emittance. The study shows that, to reduce transverse emittance, one should collimate at least 4% of particles which has either much low energy or large transverse divergence. The multiple coulomb scattering is not trivial, leading to a limited emittance reduction ratio. Small transverse emittances are of essential importance for the accelerator facilities generating free electron lasers, especially in hard X-ray region. The idea of reducing transverse emittance with tapered energy-loss foil is recently proposed by B. Carlsten [1], and can be traced back to J.M. Peterson's work in 1980s [2]. Peterson illustrated that a transverse energy gradient can be produced with a tapered energy-loss foil which in turn leads to transverse emittance reduction, and also analyzed the emittance growth from the associated multiple coulomb scattering. However, what Peterson proposed was rather a conceptual than a practical design. In this paper, we build a more complete physical model of the tapered foil based on Ref. [2], including the analysis of the transverse emittance reduction using the concept of eigen emittance and confirming the results by various numerical simulations. The eigen emittance equals to the projected emittance when there is no cross correlation in beam's second order moments matrix [3]. To calculate the eigen emittances, it requires only to know the beam distribution at the foil exit. Thus, the analysis of emittance reduction and the optics design of the subsequent beam line section can be separated. In addition, we can combine the effects of multiple coulomb scattering and transverse energy gradient together in the beam matrix and analyze their net effect. We find that,when applied to an electron linac or electron beam line, the energy spread increase and angular growth due to multiple scattering are not trivial; as a result, the transverse emittance can only be reduced with a limited ratio, e.g. down to about 65% the original value. The contents of this paper are arranged as follows. In Sec. II, we build the physical model of the tapered foil, derive the transverse eigen emittance and discuss the emittance reduction criteria. In Sec. III, we implement numerical simulations to verify the physical model; and in Sec. IV, we present numerical experiments and subsequent beam line to remove the transverse energy gradient to demonstrate the applicability of such method. Conclusions are given in the last section.

43 PARTICLE ACCELERATORS↗

Stellar populations with MEGARA: The inner regions of NGC 7025

This paper aims to determine the capabilities of the MEGARA spectrograph at the Gran Telescopio Canarias (GTC), which is an optical integral-field unit, for studying stellar populations. We also aim to exploit its combination of high spectral (R ~ 6000, 12 000 and 20 000) and spatial (0.62") resolution within its 12".5 × 11"​​.3 field of view. We do this by analysing the commissioning data of the nearby S0a galaxy NGC 7025. We establish a systematic method through which we can determine the properties of the stellar populations in the observations made with MEGARA, more specifically, within the MEGADES legacy project. For this paper in particular, we determine the properties of the stellar populations of NGC 7025. We used MEGARA observations of galaxy NGC 7025 that were taken during the commissioning phase of the instrument. We applied different approaches to estimate the properties of the stellar populations with the highest possible certainty. In addition to the specific study of NGC 7025 and in the context of the MEGADES survey, we have carried out a number of tests to determine the expected errors (including potential biases) in these star formation history (SFH) derivations as a function of these parameters, namely spectral setup, signal-to-noise ratio, σ, and the SFH itself. All the studies we conduct (both full spectral fitting and absorption line indices) of the stellar populations of NGC 7025 indicate that the stars that form its bulge have supersolar metallicity and considerably old ages (~10 Gyr) in general. Using three different combinations of MEGARA spectral setups, we determined that the bulge of NGC 7025 has smild negative mass-weighted age gradient. For the more detailed SFH, our results indicate that in addition to a rather constant star formation at early epochs, a peak in the formation history of the stars in the bulge is also found 3.5–4.5 Gyr ago. This partly explains the mass-weighted age gradients we measured. The scenario presented in NGC 7025 is that of an isolated galaxy under secular evolution that about 3.5–4.5 Gyr ago likely experienced a minor merger (mass ratio 1/10) that induced an increase in star formation and also perturbed the morphology of its outer disc. In addition to these specific results for NGC 7025, we report different lessons learned for the ongoing exploitation of the MEGADES survey with the GTC, such as the need to obtain combined observations in the LR-B + LR-V setups and a signal-to-noise ratio of at least 20 per Å.

79 ASTRONOMY AND ASTROPHYSICS↗

The reliability of X-ray constraints of intrinsic cluster shapes

Using the simulation of Katz & White (1993) we have tested the viability of X-ray analysis for constraining the intrinsic shapes of clusters of galaxies considering the effects of both substructure and steep temperature gradients. We restrict our analysis to the aggregate shapes of clusters on scales of r approximately 1-2 Mpc in order to reduce our sensitivity to subclustering in the core. For low redshifts (z less than or approximately = 0.25) the X-ray method accurately measures the true ellipticity of the three-dimensional cluster dark matter provided the inclination of the cluster is known to within approximately 30 deg; assuming the gas is isothermal adds only small errors to the derived shapes. At higher redshifts the X-ray method yields unreliable results since the gas does not trace the cluster gravitational potential. We proffer some necessary conditions for the reliability of X-ray methods characterized by both the amount of substructure in the X-ray surface brightness images and the shapes of the isophotes. We conclude that measurements of the aggregate shapes of clusters on scales r approximately 1-2 Mpc are insensitive to core substructure representing scales of a few hundred kpc. Therefore our results suggest that the X-ray measurements of aggregate cluster shapes by Fabricant, Rybicki, & Gorenstein (1984) and Buote & Canizares (1992) are valid provided that they do not suffer from serious projection effects. A substantial number of Abell clusters observed with the ROSAT PSPC will be amenable to X-ray shape analysis.

Buote, David A.↗

Numerical and Experimental Investigation of Turbulent Transport Control via Shaping of Radial Plasma Flow Profiles

Turbulence, and turbulence-driven transport are ubiquitous in magnetically confined plasmas, where there is an intimate relationship between turbulence, transport, instability driving mechanisms (such as gradients), plasma flows, and flow shear. Though many of the detailed physics of the interrelationship between turbulence, transport, drive mechanisms, and flow remain unclear, there have been many demonstrations that transport and/or turbulence can be suppressed or reduced via manipulations of plasma flow profiles. This is well known in magnetic fusion plasmas [e.g., high confinement mode (H-mode) and internal transport barriers (ITB’s)], and has also been demonstrated in laboratory plasmas. However, it may be that the levels of particle transport obtained in such cases [e.g. H-mode, ITB’s] are actually lower than is desirable for a practical fusion device. Ideally, one would be able to actively feedback control the turbulent transport, via manipulation of the flow profiles. The purpose of this research was to investigate the feasibility of using both advanced model-based control algorithms, as well as non-model-based algorithms, to control cross-field turbulence-driven particle transport through appropriate manipulation of radial plasma flow profiles. The University of New Mexico was responsible for the experimental portion of the project, while our collaborators at the University of Montana provided plasma transport modeling, and collaborators at Lehigh University developed and explored control methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Seismic Monitoring at the Farnsworth CO2-EOR Field Using Time-Lapse Elastic-Waveform Inversion of 3D-3C VSP Data

During the Development Phase of the U.S. Southwest Regional Partnership on Carbon Sequestration, supercritical CO2 was continuously injected into the deep oil-bearing Morrow B formation of the Farnsworth Unit in Texas for Enhanced Oil Recovery (EOR). The project injected approximately 94 kilotons of CO2 to study geologic carbon storage during CO2-EOR. A three-dimensional (3D) surface seismic dataset was acquired in 2013 to characterize the subsurface structures of the Farnsworth site. Following this data acquisition, the baseline and three time-lapse three-dimensional three-component (3D-3C) vertical seismic profiling (VSP) data were acquired at a narrower surface area surrounding the CO2 injection and oil/gas production wells between 2014 and 2017 for monitoring CO2 injection and migration. With these VSP datasets, we inverted for subsurface velocity models to quantitatively monitor the CO2 plume within the Morrow B formation. We first built 1D initial P-wave (Vp) and S-wave (Vs) velocity models by upscaling the sonic logs. We improved the deep region of the Vp and Vs models by incorporating the deep part of a migration velocity model derived from the 3D surface seismic data. We improved the shallow region of 3D Vp and Vs models using 3D traveltime tomography of first arrivals of VSP downgoing waves. We further improved the 3D baseline velocity models using elastic-waveform inversion (EWI) of the 3D baseline VSP upgoing data. Our advanced EWI method employs alternative tomographic and conventional gradients and total-variation-based regularization to ensure the high-fidelity updates of the 3D baseline Vp and Vs models. We then sequentially applied our 3D EWI method to the three time-lapse datasets to invert for spatiotemporal changes of Vp and Vs in the reservoir. Our inversion results reveal the volumetric changes of the time-lapse Vp and Vs models and show the evolution of the CO2 plume from the CO2 injection well to the oil/gas production wells.

42 ENGINEERING↗

Large Coded Aperture Mask for Spaceflight Hard X-ray Images

The 2.6 square meter coded aperture mask is a vital part of the Burst Alert Telescope on the Swift mission. A random, but known pattern of more than 50,000 lead tiles, each 5 mm square, was bonded to a large honeycomb panel which projects a shadow on the detector array during a gamma ray burst. A two-year development process was necessary to explore ideas, apply techniques, and finalize procedures to meet the strict requirements for the coded aperture mask. Challenges included finding a honeycomb substrate with minimal gamma ray attenuation, selecting an adhesive with adequate bond strength to hold the tiles in place but soft enough to allow the tiles to expand and contract without distorting the panel under large temperature gradients, and eliminating excess adhesive from all untiled areas. The largest challenge was to find an efficient way to bond the > 50,000 lead tiles to the panel with positional tolerances measured in microns. In order to generate the desired bondline, adhesive was applied and allowed to cure to each tile. The pre-cured tiles were located in a tool to maintain positional accuracy, wet adhesive was applied to the panel, and it was lowered to the tile surface with synchronized actuators. Using this procedure, the entire tile pattern was transferred to the large honeycomb panel in a single bond. The pressure for the bond was achieved by enclosing the entire system in a vacuum bag. Thermal vacuum and acoustic tests validated this approach. This paper discusses the methods, materials, and techniques used to fabricate this very large and unique coded aperture mask for the Swift mission.

Vigneau, Danielle N.↗

Use of Microgravity to Control the Microstructure of Eutectics

This grant began in June of 1996. Its long term goal is to be able to control the microstructure of directionally solidified eutectic alloys, through an improved understanding of the influence of convection. The primary objective of the present projects is to test hypotheses for the reported influence of microgravity on the microstructure of three fibrous eutectics (MnBi-Bi, InSb-NiSb, Al3Ni-Al). A secondary objective is to determine the influence of convection on the microstructure of other eutectic alloys. Two doctoral students and a masters student supported as a teaching assistant were recruited for this research. Techniques were developed for directional solidification of MnBi-Bi eutectics with periodic application of current pulses to produce an oscillatory freezing rate. Image analysis techniques were developed to obtain the variation in MnBi fiber spacing, which was found to be normally distributed. The mean and standard deviation of fiber spacing were obtained for several freezing conditions. Eighteen ampoules were prepared for use in the gradient freeze furnace QUELD developed at Queen's University for use in microgravity. Nine of these ampoules will be solidified soon at Queen's in a ground-based model. We hope to solidify the other nine in the QUELD that is mounted on the Canadian Microgravity Isolation Mount on MIR. Techniques are being developed for directional solidification of the Al-Si eutectic at different freezing rates, with and without application of accelerated crucible rotation to induce convection. For the first time, theoretical methods are being developed to analyze eutectic solidification with an oscillatory freezing rate. In a classical sharp-interface model, we found that an oscillatory freezing rate increases the deviation of the average interfacial composition from the eutectic, and increases the undercooling of the two phases by different amounts. This would be expected to change the volume fraction solidifying and the fiber spacing. Because of difficulties in tracking the freezing interfaces of the two solid phases, a phase-field model is also being developed. A paper demonstrating application of phase field methods to periodic structures has been submitted for publication.

Wilcox, William R.↗

THREE-DIMENSIONAL NUMERICAL SIMULATIONS OF FAST-TO-ALFVEN CONVERSION IN SUNSPOTS

The conversion of fast waves to the Alfven mode in a realistic sunspot atmosphere is studied through three-dimensional numerical simulations. An upward propagating fast acoustic wave is excited in the high-{beta} region of the model. The new wave modes generated at the conversion layer are analyzed from the projections of the velocity and magnetic field in their characteristic directions, and the computation of their wave energy and fluxes. The analysis reveals that the maximum efficiency of the conversion to the slow mode is obtained for inclinations of 25 Degree-Sign and low azimuths, while the Alfven wave conversions peak at high inclinations and azimuths between 50 Degree-Sign and 120 Degree-Sign . Downward propagating Alfven waves appear at the regions of the sunspot where the orientation of the magnetic field is in the direction opposite to the wave propagation, since at these locations the Alfven wave couples better with the downgoing fast magnetic wave which is reflected due to the gradients of the Alfven speed. The simulations show that the Alfven energy at the chromosphere is comparable to the acoustic energy of the slow mode, being even higher at high inclined magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Final Technical Report for the SciDAC Center for the Study of Plasma Microturbulence

This final technical report for the SciDAC Center for the Study of Plasma Microturbulence documents contributions from MIT to the project over the period of the award (February 15, 2011 to February 14, 2020). This project focused on the development of neoclassical and gyrokinetic simulation capabilities relevant for the description of particle, momentum and energy transport in tokamak fusion plasmas. Extensive validation against experimental results from Alcator C-Mod and DIII-D experiments was performed. Overall areas of emphasis included gyrokinetic simulations, experiments, and observation of density gradient driven Trapped Electron Mode (TEM) turbulence relevant to burning plasmas; theory and gyrokinetic simulations of the parametric variation of the nonlinear upshift of the TEM critical density gradient; identification of a new mechanism for the isotope effect due to density gradient driven TEM turbulence; new synthetic diagnostics based on full-wave simulations which enabled unprecedented agreement between gyrokinetic simulations and measured fluctuation spectra with and without strong electron heating; development of new continuum pedestal neoclassical and turbulence codes; formulation and implementation of the first gyrokinetic exact Fokker-Planck collision operator in gyrokinetic turbulence simulations; development of new spectral approaches in velocity space for highly efficient treatments of collisions; new model collision operators implemented in the discontinuous Galerkin gyrokinetic edge turbulence code GKEYLL; and a new model and spectral code to study the interaction of ion- and electron-scale turbulence, which will serve as an efficient test-bed for multi-rate and multi-scale methods. These fundamental advances contribute to the development of a validated, predictive capability for turbulent and collisional particle, momentum and energy transport in both the core and edge of magnetically confined fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗