Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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

X-ray transparent microfluidic chips for high-throughput screening and optimization of in meso membrane protein crystallization

Elucidating and clarifying the function of membrane proteins ultimately requires atomic resolution structures as determined most commonly by X-ray crystallography. Many high impact membrane protein structures have resulted from advanced techniques such as in meso crystallization that present technical difficulties for the set-up and scale-out of high-throughput crystallization experiments. In prior work, we designed a novel, low-throughput X-ray transparent microfluidic device that automated the mixing of protein and lipid by diffusion for in meso crystallization trials. Here, in this work, we report X-ray transparent microfluidic devices for high-throughput crystallization screening and optimization that overcome the limitations of scale and demonstrate their application to the crystallization of several membrane proteins. Two complementary chips are presented: (1) a high-throughput screening chip to test 192 crystallization conditions in parallel using as little as 8 nl of membrane protein per well and (2) a crystallization optimization chip to rapidly optimize preliminary crystallization hits through fine-gradient re-screening. We screened three membrane proteins for new in meso crystallization conditions, identifying several preliminary hits that we tested for X-ray diffraction quality. Further, we identified and optimized the crystallization condition for a photosynthetic reaction center mutant and solved its structure to a resolution of 3.5 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: synthetic experiments and Hussar land data set application

The problem of inverting for multiple physical parameters in the subsurface using seismic full-waveform inversion (FWI) is complicated by interparameter trade-off arising from inherent ambiguities between different physical parameters. Parameter resolution is often characterized using scattering radiation patterns, but these neglect some important aspects of interparameter trade-off. More general analysis and mitigation of interparameter trade-off in isotropic-elastic FWI is possible through judiciously chosen multiparameter Hessian matrix-vector products. We show that products of multiparameter Hessian off-diagonal blocks with model perturbation vectors, referred to as interparameter contamination kernels , are central to the approach. We then apply the multiparameter Hessian to various vectors designed to provide information regarding the strengths and characteristics of interparameter contamination, both locally and within the whole volume. With numerical experiments, we observe that S-wave velocity perturbations introduce strong contaminations into density and phase-reversed contaminations into P-wave velocity, but themselves experience only limited contaminations from other parameters. Based on these findings, we introduce a novel strategy to mitigate the influence of interparameter trade-off with approximate contamination kernels. Furthermore, we recommend that the local spatial and interparameter trade-off of the inverted models be quantified using extended multiparameter point spread functions (EMPSFs) obtained with pre-conditioned conjugate-gradient algorithm. Compared to traditional point spread functions, the EMPSFs appear to provide more accurate measurements for resolution analysis, by de-blurring the estimations, scaling magnitudes and mitigating interparameter contamination. Approximate eigenvalue volumes constructed with stochastic probing approach are proposed to evaluate the resolution of the inverted models within the whole model. With a synthetic Marmousi model example and a land seismic field data set from Hussar, Alberta, Canada, we confirm that the new inversion strategy suppresses the interparameter contamination effectively and provides more reliable density estimations in isotropic-elastic FWI as compared to standard simultaneous inversion approach.

58 GEOSCIENCES↗

Application of the weighted-density approximation to the accurate description of electron-positron correlation effects in materials

Here, we discuss positron-annihilation lifetimes for a set of illustrative bulk materials within the framework of the weighted-density approximation (WDA). The WDA can correctly describe electron-positron correlations in strongly inhomogeneous systems, such as surfaces, where the applicability of (semi-)local approximations is limited. We analyze the WDA in detail and show that the electrons which cannot screen external charges efficiently, such as the core electrons, cannot be treated accurately via the pair correlation of the homogeneous electron gas. We discuss how this problem can be addressed by reducing the screening in the homogeneous electron gas by adding terms depending on the gradient of the electron density. Further improvements are obtained when core electrons are treated within the LDA and the valence electron using the WDA. Finally, we discuss a semiempirical WDA-based approach in which a sum rule is imposed to reproduce the experimental lifetimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quasiparticle self-consistent GW band structures and phase transitions of LiAlO 2 in tetrahedrally and octahedrally coordinated structures

In this work, a first-principles computational study is presented of various phases of LiAlO 2 . The relative total energies and equations of state of the α, β, and γ phases are determined after structural relaxation of each phase. The β and γ tetrahedral phases are found to be very close in energy and lattice volume with the γ phase having the lowest energy. The octahedral α phase is a high-pressure phase and the transition pressure from the γ and β phases to α is determined to be about 1 GPa. The electronic band structures of each phase at their own equilibrium volume are determined using the quasiparticle self-consistent (QS) GW method as well as using the 0.8 Σ approach in which the QS GW self-energy is reduced by a factor of 0.8 to correct for the underscreening of W in QS GW. The effective masses of the band edges and the nature of the band gaps are presented. The lowest energy γ phase is found to have a pseudodirect gap of 7.69 eV. The gap is direct at Γ but corresponds to a dipole forbidden transition. The imaginary part of the dielectric function and the absorption coefficient are calculated in the long-wavelength limit and the random phase approximation, without local field or electron-hole interaction effects for each phase, and their anisotropies are discussed. Si doping on the Al site is investigated as a possible n-type dopant in γ-LiAlO 2 using a 128-atom supercell corresponding to 3.125% Si on the Al sublattice in the generalized gradient approximation and a smaller 16-atom cell with 25% Si in the QS GW approximation. The Si is found to significantly perturb the conduction band and lower the gap but a clearly separated deep donor defect level is not found. However, the donor binding energy is still expected to be relatively deep, on the order of a few tenths eV in the hydrogenic effective mass approximation.

36 MATERIALS SCIENCE↗

Combining synchrotron X-ray diffraction, mechanistic modeling and machine learning for in situ subsurface temperature quantification during laser melting

Laser melting, such as that encountered during additive manufacturing, produces extreme gradients of temperature in both space and time, which in turn influence microstructural development in the material. Qualification and model validation of the process itself and the resulting material necessitate the ability to characterize these temperature fields. However, well established means to directly probe the material temperature below the surface of an alloy while it is being processed are limited. To address this gap in characterization capabilities, a novel means is presented to extract subsurface temperature-distribution metrics, with uncertainty, from in situ synchrotron X-ray diffraction measurements to provide quantitative temperature evolution data during laser melting. Temperature-distribution metrics are determined using Gaussian process regression supervised machine-learning surrogate models trained with a combination of mechanistic modeling (heat transfer and fluid flow) and X-ray diffraction simulation. The trained surrogate model uncertainties are found to range from 5 to 15% depending on the metric and current temperature. The surrogate models are then applied to experimental data to extract temperature metrics from an Inconel 625 nickel superalloy wall specimen during laser melting. The maximum temperatures of the solid phase in the diffraction volume through melting and cooling are found to reach the solidus temperature as expected, with the mean and minimum temperatures found to be several hundred degrees less. The extracted temperature metrics near melting are determined to be more accurate because of the lower relative levels of mechanical elastic strains. However, uncertainties for temperature metrics during cooling are increased due to the effects of thermomechanical stress.

36 MATERIALS SCIENCE↗

Strategy to Develop a Control Scheme for Core Thermal Performance Optimization

This report outlines a three-year research plan for optimizing reactor core thermal performance through use of a Digital Twin (DT) model and a set of neutron detectors to correspondingly adapt the reactor control strategy. It identifies reactor features and achievable neutron measurements that factor into this optimization task. It considers spatial effects that are important and how they can be managed by a real-time algorithm that controls reactivity actuators. The report presents a set of tasks, along with corresponding methods, for accomplishing this objective. The final task in this set is to validate the combined optimization procedure in the Purdue University’s PUR-1 reactor, for which we have an agreement of understanding. Additionally, we describe an alternative approach for overcoming some of the limitations inherent in the above approach, which includes the computational costs of high-fidelity simulations and achievable integration of neutron measurements into the DT. The proposed algorithms lower the technology readiness level (TRL) of the overall approach as some additional analysis and the development of a new neutron detector concept are necessary. In particular, a sensor measuring the gradient of the neutron flux is required, and the prototype is expected to be validated in the PUR-1 reactor. The potential benefits and the opportunity to significantly push forward the state-of-the-art make this approach worthy of further exploration in the upcoming years.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimal detumbling of a large manned spacecraft using an internal moving mass

The use of a movable mass control system to stabilize an arbitrarily tumbling asymmetric vehicle about the maximum inertia axis is considered. A first-order gradient optimization technique is used to minimize angular velocity components along the intermediate and minimum inertia axes. This method permits a wide range of initial guesses for mass position history. Motion of the control mass is along a linear track fixed in the vehicle. The control variable is taken as mass acceleration with respect to body coordinates. Motion is limited to defined quantities and a penalty function is used to insure a given range of positions. Numerical solutions of the optimization equations verify that minimum time detumbling is achieved with the largest permissible movable mass, length of linear track, and positions of the mass on the two coordinates perpendicular to the linear motion.

Kunciw, B. G.↗

Theory of viscous transonic flow over airfoils at high Reynolds number

This paper considers viscous flows with unseparated turbulent boundary layers over two-dimensional airfoils at transonic speeds. Conventional theoretical methods are based on boundary layer formulations which do not account for the effect of the curved wake and static pressure variations across the boundary layer in the trailing edge region. In this investigation an extended viscous theory is developed that accounts for both effects. The theory is based on a rational analysis of the strong turbulent interaction at airfoil trailing edges. The method of matched asymptotic expansions is employed to develop formal series solutions of the full Reynolds equations in the limit of Reynolds numbers tending to infinity. Procedures are developed for combining the local trailing edge solution with numerical methods for solving the full potential flow and boundary layer equations. Theoretical results indicate that conventional boundary layer methods account for only about 50% of the viscous effect on lift, the remaining contribution arising from wake curvature and normal pressure gradient effects.

Melnik, R. E.↗

NDE of advanced ceramics

Radiographic, ultrasonic, and scanning laser acoustic microscopy (SLAM) techniques were used to characterize silicon nitride and silicon carbide modulus-of-rupture test specimens in various stages of fabrication. Conventional and microfocus X-ray techniques were found capable of detecting minute high-density inclusions in as-received powders, green compacts, and fully densified specimens. Significant density gradients in sintered bars were observed by radiography, ultrasonic velocity, and SLAM. Ultrasonic attenuation was found sensitive to microstructural variations due to grain and void morphology and distribution. SLAM was capable also of detecting voids, inclusions, and cracks in finished test bars. Consideration is given to the potential for applying thermoacoustic microscopy techniques to green and densified ceramics. Some limitations and the detection probability statistics of the aforementioned nondestructive evaluation (NDE) processes are also discussed.

Klima, S. J.↗

Time- and position-varying infrared scene simulation

Rapidly changing infrared scenes can be generated as a function of position using thermally produced irradiance. This can be done most conveniently by using a variable dwell-time, variable diameter, and scanning laser beam to deposit a predetermined amount of laser energy, as a function of position, on the target surface. An infrared target simulator has been designed and tested to perform within stringent temperature simulation requirements: heating rate of over 100 K/sec is achieved, and temperature gradients of 200 K/mm can be maintained. Problems with heat transfer between adjacent elements at different temperatures are solved by modeling the thermal layer on the target surface as a low-pass filter, which cuts off all the frequencies above the thermal cutoff frequency. When the system approach is taken to the infrared simulation design, the optical system limited performance results in the best realizable performance. In this case, the thermal layer on the target surface is designed so that the thermal cutoff frequency is above the optical cutoff frequency.

Scholl, Marija S.↗

Improved accuracy of radar WPMM estimated rainfall upon application of objective classification criteria

Application of the window probability matching method to radar and rain gauge data that have been objectively classified into different rain types resulted in distinctly different Z(sub e)-R relationships for the various classifications. These classification parameters, in addition to the range from the radar, are (a) the horizontal radial reflectivity gradients (dB/km); (b) the cloud depth, as scaled by the effective efficiency; (c) the brightband fraction within the radar field window; and (d) the height of the freezing level. Combining physical parameters to identify the type of precipitation and statistical relations most appropriate to the precipitation types results in considerable improvement of both point and areal rainfall measurements. A limiting factor in the assessment of the improved accuracy is the inherent variance between the true rain intensity at the radar measured volume and the rain intensity at the mouth of the rain guage. Therefore, a very dense rain gauge network is required to validate most of the suggested realized improvement. A rather small sample size is required to achieve a stable Z(sub e)-R relationship (standard deviation of 15% of R for a given Z(sub e)) -- about 200 mm of rainfall accumulated in all guages combined for each classification.

Rosenfeld, Daniel↗

High-Temperature, Thin-Film Strain Gages Improved

Conventional resistance strain gage technology uses "bonded" strain gages. These foil or wire gages are bonded onto the surface of the test article with glue, ceramic cements, or flame-sprayed ceramics. These bonding agents can, in some instances, limit both the degree of strain transmission from the test structure to the gage and the maximum working temperature of the gage. Also, the bulky, bonded gage normally disrupts aerodynamic gas flow on the surface of the test structure because of its intrusive character. To respond to the urgent needs in aeronautic and aerospace research where stress and temperature gradients are high, aerodynamic effects need to be minimized, and higher operational temperatures are required, the NASA Lewis Research Center developed a thin film strain gage. This gage, a vacuum-deposited thin film formed directly on the surface of a test structure, operates at much higher temperatures than commercially available gages do and with minimal disruption of the aerodynamic flow. The gage uses an alloy, palladium-13 wt % chromium (hereafter, PdCr), which was developed by United Technologies Research Center under a NASA contract. PdCr is structurally stable and oxidation resistant up to at least 1100 C (2000 F); its temperature-induced resistance change is linear, repeatable, and not sensitive to the rates of heating and cooling. An early strain gage, which was made of 25-micrometer-diameter PdCr wire and demonstrated to be useable to 800 C, won an R&D 100 award in 1991. By further improving the purity of the material and by developing gage fabrication techniques that use sputter-deposition, photolithography patterning, and chemical etching, we have made an 8- to 10-m PdCr thin-film strain gage that can measure dynamic and static strain to at least 1100 C. For static strain measurements, a 5-m-thick Pt element serves as a temperature compensator to further minimize the temperature effect of the gage. These thin-film gages provide the advantage of minimally intrusive surface strain measurements and give highly repeatable readings with low drift at temperatures from ambient to 1100 C. This is a 300 C advance in operating temperature over the PdCr wire gage and a 500 C advance over commercially available gages made of other materials.

Source record↗

Monitoring Atmospheric CO2 From Space: Challenge & Approach

Atmospheric CO2 is the key radiative forcing for the Earth's climate and may contribute a major part of the Earth's warming during the past 150 years. Advanced knowledge on the CO2 distributions and changes can lead considerable model improvements in predictions of the Earth's future climate. Large uncertainties in the predictions have been found for decades owing to limited CO2 observations. To obtain precise measurements of atmospheric CO2, certain challenges have to be overcome. For an example, global annual means of the CO2 are rather stable, but, have a very small increasing trend that is significant for multi-decadal long-term climate. At short time scales (a second to a few hours), regional and subcontinental gradients in the CO2 concentration are very small and only in an order of a few parts per million (ppm) compared to the mean atmospheric CO2 concentration of about 400 ppm, which requires atmospheric CO2 space monitoring systems with extremely high accuracy and precision (about 0.5 ppm or 0.125%) in spatiotemporal scales around 75 km and 10-s. It also requires a decadal-scale system stability. Furthermore, rapid changes in high latitude environments such as melting ice, snow and frozen soil, persistent thin cirrus clouds in Amazon and other tropical areas, and harsh weather conditions over Southern Ocean all increase difficulties in satellite atmospheric CO2 observations. Space lidar approaches using Integrated Path Differential Absorption (IPDA) technique are considered to be capable of obtaining precise CO2 measurements and, thus, have been proposed by various studies including the 2007 Decadal Survey (DS) of the U.S. National Research Council. This study considers to use the Intensity-Modulated Continuous-Wave (IM-CW) lidar to monitor global atmospheric CO2 distribution and variability from space. Development and demonstration of space lidar for atmospheric CO2 measurements have been made through joint adventure of NASA Langley Research Center and Exelis, Inc. As prototype space IPDA lidars, airborne laser absorption lidar systems operating in 1.57 CO2 absorption band have been developed and tested through lab, ground-based range, and flight campaigns. Very encouraging results have been obtained. The signal-to-noise ratio (SNR) for clear sky IPDA measurements of CO2 differential absorption optical depth (DAOD) for a 10-s integration over vegetated areas with about 10 km range was found to be as high as 1300, resulting in an error 0.077% or equivalent CO2 mixing ratio (XCO2) column precision of ~0.3 ppm. Precise range measurements using the IM-CW lidar approach were also achieved, and the uncertainties have been shown to be at sub meter level. Based on the airborne lidar development, space lidar and atmospheric CO2 observations are simulated. It shows that with the IM-CW approach, accurate atmospheric CO2 measurements can be achieved from space, and a space mission such as that proposed by the DS will meet science goals of atmospheric CO2 monitoring.

Lin, Bing↗

HFEM3D

Software solves the three-dimensional Poisson equation div(k(grad(u)) = f, by the finite element method for the case when material properties, k, are distributed over hierarchy of edges, facets and tetrahedra in the finite element mesh. Method is described in Weiss, CJ, Finite element analysis for model parameters distributed on a hierarchy of geometric simplices, Geophysics, v82, E155-167, doi:10.1190/GEO2017-0058.1 (2017). A standard finite element method for solving Poisson’s equation is augmented by including in the 3D stiffness matrix additional 2D and 1D stiffness matrices representing the contributions from material properties associated with mesh faces and edges, respectively. The resulting linear system is solved iteratively using the conjugate gradient method with Jacobi preconditioning. To minimize computer storage for program execution, the linear solver computes matrix-vector contractions element-by-element over the mesh, without explicit storage of the global stiffness matrix. Program output vtk compliant for visualization and rendering by 3rd party software. Program uses dynamic memory allocation and as such there are no hard limits on problem size outside of those imposed by the operating system and configuration on which the software is run. Dimension, N, of the finite element solution vector is constrained by the the addressable space in 32-vs-64 bit operating systems. Total storage requirements for the problem. Total working space required for the program is approximately 13*N double precision words.

Weiss, Chester J↗

A Novel Technique for Distributed Measuring of Strain in the Vadose Zone

Characterizing water content and pressure changes in the vadose zone is important to understanding a variety of geologic processes, ranging from permeability, to evapotranspiration, and aquifer recharge. Changes in water content, or pressure, cause strain in the solid porous medium of the vadose zone, and it is possible that measuring those strains could be used a characterization tool. The Coherence-length-gated Microwave Photonics Interferometry (CMPI) technique measures strain at high resolution along many intervals defined by pairs of reflectors distributed along an optical fiber. This technique has recently been developed at Clemson University, and it has the spatial, and temporal resolution to characterize the strains that are expected to occur with hydrologic changes in the vadose zone. However, the technique has never been used to measure strain in porous media, so its capabilities remain uncertain. The objective of this thesis is to evaluate the ability of using CMPI to measure strain changes in the vadose zone. The research approach consists of conducting laboratory tests using a column filled with sand that was subjected to changes in water content and pressure. An optical fiber with CMPI reflectors was deployed in a high surface area ribbon and used to measure strain along the axis of the column. The column was made from 8-inch, Schedule 40 PVC pipe (75 cm tall, 20 cm inner diameter) and filled with medium-grained sand (K = 2.5x10-6 m/s, porosity = 0.28, Coefficient of Uniformity = 2, van Genuchten (n = 4.34 and  = 0.00039 m-1), Young’s Modulus = 18 – 47 MPa). Five pressure and four temperature sensors attached to probes inserted into the wall. The optical fiber includes five reflectors spaced 10 cm apart on the inside and outside of the column along the vertical axis. Strain measurements are made between pairs of reflectors. The reflectors are created in 250-micron-diameter acrylate-coated single mode Corning SMF28e+ optical fiber using a femtosecond laser. The optical fiber was laminated between two pieces of polyester film creating a large surface area to transfer strain from the porous media to the optical fiber. The experiments were conducted by filling the column from the bottom or infiltrating water from the top. The water was allowed to equilibrate to room temperature prior to each test in order to limit thermoelastic strain. Five injection tests with a rate of 250 ml/min and three infiltration tests at varying rates were conducted, and the results show patterns of strain and pressure are generally similar. Hydrologic conditions define three zones based on the pressure magnitude and distribution. 1.) Ambient Zone where the pressure heads are quasi-static and the pressure gradient is roughly unity (head gradient of zero). This is the uppermost zone and is characterized by negative pressures. 2.) Transition Zone where the pressures increase from ambient to zero, are changing relatively rapidly and the pressure gradient is relatively steep (pressure head gradients of 2). The Transition zone is 10 to 15 cm thick. 3.) Positive Pressure Zone where the pressure is positive, the rate of change is slower than in the transition zone and the gradient is flatter (pressure head gradient 1.1 to 1.2). This is the lowest zone in the column. Injection of water causes the pressure to increase and the three zones to move upward. The Transition zone moves at a rate of approximately 0.0002 m/s +/- .00005, according to analyses of pressure profiles. This is the same as the average velocity of the water calculated as volumetric flux/effective porosity =5.4x10-5 m/s/0.28. Strain signals in the range of 10s  were observed with a noise level of generally less than 0.1  (signal to noise ratio of greater than 100) during injection and drainage. The spatial and temporal distributions of strain caused by injection depend on the location of the moving pressure zones. The locations of the different pressure zones were determined from pressure profiles at different times and these data were transferred to strain time series. This showed that strain in the Ambient pressure zone is either unchanged or slightly compressive, whereas strain in Transition zone is tensile and roughly proportional to the pressure change, and strain ranges from tensile to compressive in the Positive Pressure zone. In general, the strain is variable at the top of the Positive Pressure zone, but it appears to be consistently tensile lower in the zone. The spatial distribution of strain causes three distinct stages in the strain times series measured between pairs of CMPI reflectors. The strain is unchanged or slightly compressive during Stage 1, it increases (positive strain is tensile) approximately linearly at a rate of 0.02 to 0.03 μ/s during the Stage 2, and then flattens out in Stage 3. The strain increases again during Stage 4 of the time series. The pressure also changes in stages. It is unchanged during Stage 1 and then increases relatively rapidly at a rate of between 4 and 10 Pa/s during Stage 2 and continues to increase during Stages 3 and 4, but at a rate slightly slower than during Stage 2 (from 2.5 to 3.5 Pa/s). The rates of pressure and strain change are consistent with basic analyses. The rate of pressurization is similar to the calculated velocity of flow in the saturated zone, and this is similar to the measured velocity of the strain increase. The ratio of the pressurization rate and strain rate during Stage 2 (2.5 (Pa/s) /0.03 (μ/s)) is 80 MPa, which is approximately equal to the upper range expected considering the uniaxial compression and the Young’s Modulus of the sand. An interesting effect occurs when the upper surface of the sand becomes saturated. Significant compression (several 10s of ) occurs as the pressure and saturation increase at the upper surface. Compression occurs throughout the column, but the effect is greatest at the top of the column. The rate of compression is fastest slightly before ponding occurs, but it slows markedly and nearly stops when water starts to accumulate at the surface (ponding). This effect reverses when the surface of the sand is drained, resulting in tension throughout the column. The compression caused by this effect can be as large or larger than the tensile strain that accumulated during filling. This effect was unexpected because increasing pressure is normally associated with tensile strain. Nevertheless, this effect was observed consistently in all tests where the pressure changes at the upper surface, including tests where water was injected from below or infiltrated from above. This effect behaves as if the pore pressure at the upper surface of the sand exerts a normal force on the boundary (increasing pore pressure exerts a downward compression on the boundary). The results of these laboratory experiments indicate that CMPI can measure strain caused by fluid pressure changes in the vadose zone with a signal to noise ratio of 100 or more. Injection and drainage cause a strain signal that is complex, but repeatable. The magnitudes of the strain signal are consistent with magnitudes that are expected based on poroelastic calculations using independently measured properties of the sand. These results indicate that the CMPI technique with an optical fiber laminated in a polyester ribbon generates data that represent the strain distribution during hydrologic processes of imbibition and drainage in variably saturated sand. This suggests that distributed strain measurements using CMPI could be a viable approach for evaluating processes in the vadose zone, laying the groundwork for future field implementation.

47 OTHER INSTRUMENTATION↗

Modulated heat pulse propagation and partial transport barriers in chaotic magnetic fields

Direct numerical simulations of the time dependent parallel heat transport equation modeling heat pulses driven by power modulation in three-dimensional chaotic magnetic fields are presented. The numerical method is based on the Fourier formulation of a Lagrangian-Green's function method that provides an accurate and efficient technique for the solution of the parallel heat transport equation in the presence of harmonic power modulation. The numerical results presented provide conclusive evidence that even in the absence of magnetic flux surfaces, chaotic magnetic field configurations with intermediate levels of stochasticity exhibit transport barriers to modulated heat pulse propagation. In particular, high-order islands and remnants of destroyed flux surfaces (Cantori) act as partial barriers that slow down or even stop the propagation of heat waves at places where the magnetic field connection length exhibits a strong gradient. Results on modulated heat pulse propagation in fully stochastic fields and across magnetic islands are also presented. In qualitative agreement with recent experiments in large helical device and DIII-D, it is shown that the elliptic (O) and hyperbolic (X) points of magnetic islands have a direct impact on the spatio-temporal dependence of the amplitude of modulated heat pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mining Air for Fuels and Fine Chemicals

The Intergovernmental Panel on Climate Change’s 2021 report underscores that Direct Air Capture (DAC) is an essential component to limiting global warming to 2 °C, yet separating CO 2 from air and its subsequent concentration for conversion into carbon neutral fuels/chemicals or permanent storage remains energy intensive and costly. Further, most DAC technologies rely on repetitive capture and release cycles that reduces the amount of time the DAC materials are able to collect CO 2 and increases wear and tear. Moisture swing (MS) sorption within anion exchange materials utilizes differences in relative humidity to capture and concentrate CO 2 up to 500-fold with negligible energy inputs and low-cost sorbents. Molecular modeling predicts that membranes with a dry and wet side would actively transport CO 2 from its dry side to its wet side against a counter flow of H 2 O evaporating on the dry side, however, existing MS materials are too brittle to be processed into membranes to test this hypothesis. This project brought together experts in DAC materials and design, techno-economic analysis and project management at Arizona State University (ASU), polymer processing at University of Texas at Austin (UTA) and modeling gas transport and exchange at Norther Arizona University (NAU) to develop hollow fiber membranes that use energy from water evaporation to continuously pump CO 2 from air against a concentration gradient to the membranes interior and requiring significantly less energy than current DAC technologies. While this project did not meet its ultimate objectives of demonstrating an active CO 2 pumping membrane, a number of key accomplishments were made in developing low-cost, flexible anion exchange membranes (AEM) at 3.5 m 2 scale, characterizing the ionic and molecular transport within MS materials, and developing several analytical and mathematical models of the MS and pumping process leading to new fundamental knowledge about key rate limiting steps of CO 2 flux in anion exchange membranes at low water activity. Technoeconomic models show that if CO 2 fluxes ≥ 25 µmol CO 2 m -2 s -1 (membrane surface area) are achieved that cost ≤ $\$$100/tonne CO 2 captured, purified and compressed to pipeline quality is possible. The knowledge gained in this research will guide future research in developing new materials toward overcoming these rate limiting steps that in turn could enable transformative and disruptive DAC technologies needed for capturing gigatons of CO 2 per year needed to limiting global warming to 2 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multistage coupling of independent laser-plasma accelerators

Laser-plasma accelerators (LPAs) are capable of accelerating charged particles to very high energies in very compact structures. In theory, therefore, they offer advantages over conventional, large-scale particle accelerators. However, the energy gain in a single-stage LPA can be limited by laser diffraction, dephasing, electron-beam loading and laser-energy depletion. The problem of laser diffraction can be addressed by using laser-pulse guiding and preformed plasma waveguides to maintain the required laser intensity over distances of many Rayleigh lengths; dephasing can be mitigated by longitudinal tailoring of the plasma density; and beam loading can be controlled by proper shaping of the electron beam. To increase the beam energy further, it is necessary to tackle the problem of the depletion of laser energy, by sequencing the accelerator into stages, each powered by a separate laser pulse. In this work, we present results from an experiment that demonstrates such staging. Two LPA stages were coupled over a short distance (as is needed to preserve the average acceleration gradient) by a plasma mirror. Stable electron beams from a first LPA were focused to a twenty-micrometre radius-by a discharge capillary-based active plasma lens-into a second LPA, such that the beams interacted with the wakefield excited by a separate laser. Staged acceleration by the wakefield of the second stage is detected via an energy gain of 100 megaelectronvolts for a subset of the electron beam. Changing the arrival time of the electron beam with respect to the second-stage laser pulse allowed us to reconstruct the temporal wakefield structure and to determine the plasma density. Our results indicate that the fundamental limitation to energy gain presented by laser depletion can be overcome by using staged acceleration, suggesting a way of reaching the electron energies required for collider applications.

43 PARTICLE ACCELERATORS↗