Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear Parameter Varying”

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 163 records · Page 9

Nonlinear modeling of the scaling law for the m/n = 3/2 error field penetration threshold

The scaling law for the m/n=3/2 error field (EF) penetration threshold is predicted numerically based on nonlinear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field b r at the plasma edge is scaled to the electron density n e , temperature T e , viscous time τ μ , toroidal field B t and the natural frequency ω in the form of b r /B t ∝n e αn T e αT τ μ αμ B t αB ω αω by scanning these parameters separately. Here, α n , α T , α μ , α B and α ω are the scaling coefficients on n e , T e , τ μ , B t and ω, respectively. Single-fluid modeling shows that the 3/2 EF threshold scales as b r /B t ∝n e 0.56 T e 0.6 τ μ -0.59 B t -1.15 ω, which is similar with the analytical scaling law in both the Rutherford and visco-resistive regimes. Yet, two-fluid modeling shows that the scaling law differs significantly in particular regarding the dependence on plasma rotation. In detail, the scaling coefficient α n on density decreases from 0.67 to 0.56 and α T on temperature decreases from 0.67 to 0.32, while α μ on viscous time is around -0.45 and α B on toroidal field decreases slightly from -1.15 to -1, when the ratio |ω E /ω *e | between plasma rotation frequency ω E and diamagnetic drift frequency ω *e varies from 0 to 10. Scans of the plasma rotation reveals that the penetration threshold linearly depends on the perpendicular rotation frequency (or natural frequency) ω ⊥e =ω E +ω *e , and there is a minimum in the required field amplitude when ω ⊥e 0. In addition, the enduring mystery of non-zero penetration threshold at zero plasma natural frequency in EF experiments is resolved by two-fluid simulations. We report that the very small island and smooth bifurcation in EF penetration near zero frequency is hard to detect in the experiment, leading to a finite penetration threshold within the capability of the experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗

Experimental Safe Extremum Seeking for Accelerators

Here, we demonstrate the recent designs of safe extremum seeking (Safe ES) on the 1-km-long charged particle accelerator at the Los Alamos Neutron Science Center (LANSCE). Safe ES is a modification of extremum seeking (ES) which, in addition to minimizing an analytically unknown cost, also employs a safety filter based on an analytically unknown control barrier function (CBF) safety metric. Tuning is necessitated by accelerators being large complex systems, with many drifting parameters due to thermal effects and degradation. At the same time, safe operation (the maintenance of state constraints) is crucial, as damage brings astronomical costs, both financially and in operation downtime. Our measured (but analytically unknown) safety metric is the beam current. We perform multivariable Safe ES on three accelerator applications, in which we adapt 4, 6, and 3 magnet strength parameters, respectively. Two of the three applications are for validated simulation models of beamlines at LANSCE: the first for the proton radiography (pRad) beamline of 800-MeV protons for spot size tuning; the second on a high-performance code, HPSim, for tuning the low-energy beam transport (LEBT) region that contains a beam of 750-keV protons. The third is an experimental tuning of the steering magnets in the LEBT at LANSCE.

43 PARTICLE ACCELERATORS↗

High-Speed OH-PLIF Diagnostics of Flame Flashback in Low Swirl Hydrogen-Enriched Flames

This paper reports flashback events observed in hydrogen-enriched flames stabilized in a low swirl burner (LSB) at atmospheric pressure and temperature conditions. The fundamentals of hydrogen-rich stable flames and spatiotemporal investigation of flashback phenomena were observed experimentally using a high-repetition-rate nanosecond (ns)-duration hydroxyl radical planar laser-induced fluorescence (OH-PLIF) diagnostic. Testing was conducted in an optically accessible pre-mixing section of a LSB with inlet pre-mixing velocities ranging from 5 to 10 m/s for methane and hydrogen (50 – 90% by mole) blends. Swirlers with two different turning angles, 26o and 33o, and three different perforated plate hole diameters, 1.08, 1.12 and 1.16 mm, with measured swirl numbers varying from 0.43 to 0.49 were used in this study at atmospheric temperature and pressure inlet conditions. The flashback propensity showed a dependence on the proximity of the lifted flame to the burner exit (termed as lift-off length, L) which was dependent on the pre-mixer velocity (V), hydrogen content (XH2) and equivalence ratio (ϕ) at constant temperature and pressure conditions. High-speed OH-PLIF images revealed lifted flames were first observed at low ϕ/ XH2 condition which then changed to a M-shaped flame attached to burner rim with the increase in ϕ/ XH2. Further increases in ϕ/XH2, depending upon other inlet parameters, triggered flame flashback into the pre-mixing section. Flashback occurred when burning began in the outer shear layer and the leading flame brush propagated into the nozzle. Such spatiotemporally resolved high-speed OH-PLIF imaging provides insights on different stages of flashback in a LSB – flashback initiation, transition of flame propagation from burner exit to the swirler and flashback to flame holding transition. Flame ϕ at flashback showed an expected linearly increasing trend with increasing V and decreasing XH2 and the conclusions drawn aligned well with detailed investigations. For identical inlet conditions, flashback propensity decreased with increasing perforated-plate hole diameter and increasing swirler vane angle.

Parajuli, Pradeep↗

The Davis–Chandrasekhar–Fermi method revisited

ABSTRACT Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical scales, because direct measurement of the field strength is observationally difficult. The Davis–Chandrasekhar–Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength from polarization data. It is based on the assumption that gas turbulent motion is the driving source of field distortion via linear Alfvén waves. In this work, using MHD simulations of star-forming clouds, we test the validity of the assumption underlying the DCF method by examining its accuracy in the real 3D space. Our results suggest that the DCF relation between turbulent kinetic energy and magnetic energy fluctuation should be treated as a statistical result instead of a local property. We then develop and investigate several modifications to the original DCF method using synthetic observations, and propose new recipes to improve the accuracy of DCF-derived magnetic field strength. We further note that the biggest uncertainty in the DCF analysis may come from the linewidth measurement instead of the polarization observation, especially since the line-of-sight gas velocity can be used to estimate the gas volume density, another critical parameter in the DCF method.

79 ASTRONOMY AND ASTROPHYSICS↗

Command of active and responsive elastomers by topological defects and patterns

The project resulted in the development of stimuli-responsive elastomer coatings formed by photopolymerized liquid crystal molecules. The molecular orientation of the liquid crystal elastomers is coupled to rubber-like elasticity. The project developed an approach to produce complex patterns of molecular orientation in elastomers by employing photopatterning technique based on plasmonic metamasks that convert unpolarized incoming light into a transmitted light beam with spatially-varying linear polarization. The polarization-modulated light beam aligns the film of azodye molecules, which serves as the substrate for liquid crystal elastomer coating. The alignment of the azodye molecules follows the polarization pattern with a phase shift of the alignment direction by 90 degrees. The azodye substrate imposes orientation of molecules in the adjacent liquid crystal layer which is preserved after the monomers are polymerized into a liquid crystal elastomer. By using differently aligned substrates one could produce various patterns of molecular orientation in the elastomer, which dictate the mechanical properties of the elastomers and their response to external stimuli such as temperature, humidity, and ultraviolet irradiation. In particular, the predesigned molecular orientation produces deterministically defined topography response of the liquid crystal coatings, in which their free surface changes from flat to either locally elevated or locally compressed, depending on whether the pre-inscribed molecular orientation is predesigned with a splay or bend deformation. The research established that the mechanism of this deterministic relation between the dynamic topography and pre-inscribed molecular orientation is rooted in the change of the tensorial order parameter that characterizes the degree of orientational order, and in the contraction/expansion of polymer molecules in response to the order parameter changes. In particular, lowering the order causes contraction of polymer globules. The work demonstrated that the dynamic topography of elastomer coatings with three-dimensional pattern of orientational order includes up and down motion of the material (along the normal to the coating) and also shifts in lateral directions. We demonstrate that the deterministic relationship between the orientational patterns and variations of coatings’ profile is caused by forces that are mathematically equivalent to active forces in the system of “pullers” or “pushers” in active matter, demonstrating universality of the out-of-equilibrium behavior. The project resulted in theoretical models that predictively describe the dynamic response of coatings. The topography of liquid crystal elastomers is sensitive to stimuli such as temperature, ultraviolet irradiation, and humidity. The applicability of the patterned approach to produce dynamic coatings was expanded from nematic elastomers to their smectic counterparts. The project combined a battery of imaging techniques to unveil the dynamic properties of the coatings, ranging from the state-of-the-art dynamic holographic microscopy to fluorescent confocal polarizing microscopy and polscope microscopy. The research uncovered potential applications of dynamic elastomer coatings, such as control placement of colloidal particles, engineering of living tissues with alignment patterns of cells that follow the orientational order of elastomer coatings. The approach developed in the project can be used to design programmable dynamic coatings with functionalities that mimic biological tissues such as skin.

36 MATERIALS SCIENCE↗

Driven resonant current amplification in self-organized plasma configurations with uniform λ and plasma pressure confinement

Driven resonant current amplification in uniform λ=$^{μoJ}_{B}$, self-organized plasma configurations with plasma pressure confinement is studied. This work extends previous analyses of force-free states (J ∥ B) to equilibria with plasma pressure gradients (∇P = J × B ≠ 0). An analytical treatment shows that in equilibria with uniform λ and a plasma pressure gradient that varies linearly with poloidal magnetic flux, driven resonances are either downshifted or upshifted from their force-free values. Additionally, a modification of off resonance current amplification is observed in pressure confining equilibria. Computational results add support to analytical predictions, and also include driven resonance behavior with equilibrium profiles not considered analytically. Furthermore, these results suggest not only persistence, but also modification of driven resonant current amplification in sustained spheromak configurations with plasma pressure confinement when compared to force-free states. By providing a basis for the existence of driven, high current amplification, pressure confining equilibria, this work provides insight into the accessibility of sustained spheromak states of most interest for potential fusion energy applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reciprocal swimming at intermediate Reynolds number

In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers Re≈1–1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, in this work, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate Re. In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.

42 ENGINEERING↗

Theoretical investigation of the triggering of neoclassical tearing modes by transient resonant magnetic perturbations in NSTX

The extended perturbed equilibrium code asymptotic matching code is used to simulate the triggering of n = 1 neoclassical tearing modes (NTMs) by a pulsed, rotating, n = 1, resonant magnetic perturbation (RMP) in two example NSTX discharges. Although the two discharges are significantly different from one another, the results of the two sets of simulations are quite similar. The critical n = 1 RMP pulse amplitude required to trigger an n = 1 NTM is minimized when the RMP pulse rotation frequency matches the linear natural frequency of an n = 1 tearing mode, resonant within the plasma, that is metastable to an NTM. However, if there is a frequency mismatch, then the seed magnetic island chain driven at the relevant resonant surface is forced to rotate with respect to the RMP, because the RMP pulse amplitude is nowhere near sufficient to lock the island chain to the RMP. This rotation causes the critical RMP pulse amplitude required to trigger an NTM to oscillate as the RMP pulse duration is varied. Furthermore, the critical amplitude is minimized when the RMP pulse duration is such that seed island chain executes a half-integer number of rotations with respect to the pulse. All of the minima have the same value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance Analysis and Limiting Parameters of Cross-flow Membrane-based Liquid-desiccant Air Dehumidifiers

We report that to dehumidify a humid air stream, existing air conditioning (AC) systems substantially overcool the outdoor humid air below its dew point, thereby significantly reducing energy efficiency. Directly capturing humidity, membrane-based liquid-desiccant dehumidification systems separate sensible and latent cooling (SSLC) loads and thus offer a promising pathway for a high-performance AC solution. Design of an energy-efficient SSLC-AC system, however, rests largely on detailed understating of the dehumidification process. While some studies have identified the dehumidification process mainly depends on membrane characteristics, other studies have argued that the process is limited by desiccant liquid or alternatively air thermo-hydraulic physics for typical humid climate conditions. The present study examines performance and physics of the membrane-based liquid-desiccant dehumidification process over a wide range of climate conditions through a novel 3D, two-phase, multi-species CFD model. Decoupling the thermodynamic and hydraulic effects, the study reveals that the dehumidification rate is a linear function of the water vapor pressure potential ($J=α ΔP$) summarizing the system's thermodynamic state. The slope of the curve (i.e., α) depends on hydraulic transport characteristics of the membrane pores, air stream, and desiccant solution. More importantly, it was found that the air dehumidification process is mainly limited by the air-side transport physics for thin liquid-desiccant films and commonly used porous superhydrophobic membranes. Additionally, results show that, depending on ambient/desiccant conditions and physical dehumidifier characteristics, energy effectiveness and dehumidification rate vary from 13 to 34% and from 0.13 to 1.4 g m -2 s -1 , respectively. Therefore, the present study allows to efficiently design future SSLC-based AC systems exhibiting high performance energy metrics.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Photo-z outlier self-calibration in weak lensing surveys

Calibrating photometric redshift errors in weak lensing surveys with external data is extremely challenging. In this work, we show that both Gaussian and outlier photo-z parameters can be self-calibrated from the data alone. This comes at no cost for the neutrino masses, curvature and dark energy equation of state w0, but with a 65% degradation when both w0 and wa are varied. We perform a realistic forecast for the Vera Rubin Observatory (VRO) Legacy Survey of Space and Time (LSST) 3× 2 analysis, combining cosmic shear, projected galaxy clustering and galaxy - galaxy lensing. We confirm the importance of marginalizing over photo-z outliers. We examine a subset of internal cross-correlations, dubbed "null correlations", which are usually ignored in 3× 2 analyses. Despite contributing only ~ 10% of the total signal-to-noise, these null correlations improve the constraints on photo-z parameters by up to an order of magnitude. Using the same galaxy sample as sources and lenses dramatically improves the photo-z uncertainties too. Together, these methods add robustness to any claim of detected new Physics, and reduce the statistical errors on cosmology by 15% and 10% respectively. Finally, including CMB lensing from an experiment like Simons Observatory or CMB-S4 improves the cosmological and photo-z posterior constraints by about 10%, and further improves the robustness to systematics. To give intuition on the Fisher forecasts, we examine in detail several toy models that explain the origin of the photo-z self-calibration. Our Fisher code LaSSI (Large-Scale Structure Information), which includes the effect of Gaussian and outlier photo-z, shear multiplicative bias, linear galaxy bias, and extensions to LCDM, is publicly available at \href{https://github.com/EmmanuelSchaan/LaSSI}{https://github.com/EmmanuelSchaan/LaSSI}.

79 ASTRONOMY AND ASTROPHYSICS↗

Collisionless adiabatic afterglow

In this work, we study, by numerical and analytical means, the evolution of a collisionless plasma initiated between absorbing walls. The ensuing flow is described by rarefaction waves that travel inward from the boundaries, interact, and eventually vanish after crossing through, leading up to the asymptotic stage of the decay. Particle simulations indicate that the kinetic evolution strongly resembles one found in isentropic gas dynamics. Namely, a very gradual density profile forms in the expanding central region where the rarefaction waves interact, with an accompanying linear velocity profile. Asymptotically, the density falls off as 1/t. The density and the flux at the boundary show little variation over the period when rarefaction waves still exist. Plasma potential, on the other hand, drops quite rapidly (on the underlying ion-acoustic timescale) to less than initial electron temperature Te when over 70% of the particles still remain in the system. This is due to electron kinetics being governed by conservation of adiabatic invariant in a slowly varying potential well. Analytical model of the velocity distribution is presented to explain the simulations. The results have implications for afterglow plasmas used in material processing and also for ion-extraction devices. One property of potential interest is good uniformity of the decaying plasma that occurs after approximately one ion-acoustic time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data-Informed Evaluation Framework for Integrated Energy Systems: Insights from Power, Process Heat, and Hydrogen Production Applications

The multi-criteria decision analysis (MCDA) framework provides a systematic evaluation of the diverse preferences and performance metrics associated with alternative solutions. This approach is advantageous over a single-criterion methodology, which are only valid under conditions that assume ceteris paribus or an "apples-to-apples" comparison. However, selecting suitable technologies for integrated energy systems (IES) can be likened to an "apples-to-oranges" comparison, given the heterogeneous factors at stake. These factors include economics and performance parameters, geological compatibility, and environmental impacts. Consequently, past research has often employed a mixture of qualitative and quantitative criteria tailored to the specific interests of each study. While the method proves effective in handling the intricate interplay of criteria, the resulting rankings and scores can vary from study to study. This inconsistency is introduced from the use of subjectively defined thresholds and weights. As a result, decision-makers frequently find it challenging to establish clear connections between specific criteria and the resulting scores, as the transformation of criteria into ordinal scores results in a substantial loss of information. To address this challenge, we introduce a data-informed IES evaluation framework that offers comprehensive, interpretable, and traceable evaluations backed by quantifiable rationale. First, we identified key IES evaluation criteria from a decade of literature, focusing on relevant IES applications in power, process heat, and hydrogen production. We leveraged state-of-the-art cost estimates from the Idaho National Laboratory (INL) and technical data from 78 reactor designs from the International Atomic Energy Agency (IAEA) and the Organization for Economic Co-operation and Development - Nuclear Energy Agency (OECD-NEA). Lastly, we established thresholds by analyzing the mean, variance, root mean square, and slope of values across alternatives, categorizing the preferences of decision-makers into distinct utility functions, such as linear, saturating, exponential, and stepwise. Our approach yielded two main outcomes: (1) it provided consistent assessments across different stakeholder groups and (2) it visualized uncertainties in the decision-making context via comprehensive sensitivity analysis. To demonstrate the impact of our framework, we conducted case studies on 6 reactor designs (AP1000, NuScale, BWRX-300, Xe-100, eVinci, iMSR) for the three applications. Our data-driven framework proved to be highly effective in addressing heterogenous uncertainties faced by varied decision-makers? preferences and IES applications, as well as cost and technical estimates of advanced reactors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Data-Informed Evaluation Framework for Integrated Energy Systems: Insights from Power, Process Heat, and Hydrogen Production Applications

The multi-criteria decision analysis (MCDA) framework provides a systematic evaluation of the diverse preferences and performance metrics associated with alternative solutions. This approach is advantageous over a single-criterion methodology, which are only valid under conditions that assume ceteris paribus or an "apples-to-apples" comparison. However, selecting suitable technologies for integrated energy systems (IES) can be likened to an "apples-to-oranges" comparison, given the heterogeneous factors at stake. These factors include economics and performance parameters, geological compatibility, and environmental impacts. Consequently, past research has often employed a mixture of qualitative and quantitative criteria tailored to the specific interests of each study. While the method proves effective in handling the intricate interplay of criteria, the resulting rankings and scores can vary from study to study. This inconsistency is introduced from the use of subjectively defined thresholds and weights. As a result, decision-makers frequently find it challenging to establish clear connections between specific criteria and the resulting scores, as the transformation of criteria into ordinal scores results in a substantial loss of information. To address this challenge, we introduce a data-informed IES evaluation framework that offers comprehensive, interpretable, and traceable evaluations backed by quantifiable rationale. First, we identified key IES evaluation criteria from a decade of literature, focusing on relevant IES applications in power, process heat, and hydrogen production. We leveraged state-of-the-art cost estimates from the Idaho National Laboratory (INL) and technical data from 78 reactor designs from the International Atomic Energy Agency (IAEA) and the Organization for Economic Co-operation and Development - Nuclear Energy Agency (OECD-NEA). Lastly, we established thresholds by analyzing the mean, variance, root mean square, and slope of values across alternatives, categorizing the preferences of decision-makers into distinct utility functions, such as linear, saturating, exponential, and stepwise. Our approach yielded two main outcomes: (1) it provided consistent assessments across different stakeholder groups and (2) it visualized uncertainties in the decision-making context via comprehensive sensitivity analysis. To demonstrate the impact of our framework, we conducted case studies on 6 reactor designs (AP1000, NuScale, BWRX-300, Xe-100, eVinci, iMSR) for the three applications. Our data-driven framework proved to be highly effective in addressing heterogenous uncertainties faced by varied decision-makers? preferences and IES applications, as well as cost and technical estimates of advanced reactors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Effects of internal swelling on residual elasticity of a quasi-brittle material through a composite sphere model

This work describes the development of a micromechanical-based constitutive model accounting for the effect of internal expansion on the residual elasticity of a Hashin composite material. This material is made of spherical inclusions that are subjected to gradual swelling within a quasi-brittle matrix. The main focus of this work is to describe and analyze the material mechanical response, with an additional focus on the internal swelling’s effect on the stress–strain response and residual elasticity. Microstructural features and parameters of major importance for the mechanical responses were identified. The innovative characteristics of the proposed approach are summarized as follows: (1) a full determination of the physics of a complete-damage problem throughout the whole process of inclusion swelling with upscaling techniques, which transfers the microcrack-related properties from the lower scale to upper scale; and (2) an evolution of the mechanical fields and the corresponding residual elasticity for various inclusion swelling levels. Concerning the matrix–inclusion composite, it was hypothesized that only the matrix was susceptible to cracking, with varied degrees of damage, whereas the inclusions behave elastically and the elastic modulus of the expanding inclusions remains constant. It is to emphasize that the gradual swelling of inclusions is modeled by an increasing strain in the current micromechanical-based constitutive model, and the microcracks are represented by a set of randomly oriented penny-shaped microcracks with identical radii (namely, closed cracks). The main contribution of the current research is to establish the exact mathematical solutions for the mechanical fields (stress, strain) caused by the swelling of the inclusions (mechanical loading) and derive the effective residual elasticity of a composite (structural response) subject to internal expansion and quasi-brittle damage. Based on the assumption of closed cracks that could be extended to open cracks in the upcoming work, the results proposed by the present paper help to understand the non-linear mechanical behavior of quasi-brittle materials subject to microcracking and provide a theoretical framework to be used as academical benchmark for numerical simulations.

36 MATERIALS SCIENCE↗

Temperature effects in plasma-based positron acceleration schemes using electron filaments

Preserving the quality of a positron beam in a plasma-based accelerator, where a wakefield suitable for positron transport and acceleration is generated by means of an electron filament, is challenging. This is due to the nature of the wakefields, characterized by focusing fields that vary nonlinearly in the transverse direction, and by accelerating fields that are non-uniform. These fields also change slice-by-slice along the beam. Maintaining a high beam quality is pivotal for application of positron beams in a plasma-based collider. In this paper, we show that an initial background plasma temperature can help mitigate the positron beam quality degradation in plasma-based accelerators that rely on electron filaments. We show that temperature effects broaden the electron filament and smooth radially both the non-linear transverse and the non-uniform longitudinal wakefields. Using warm plasmas opens up new possibilities to improve beam quality in several plasma-based positron acceleration concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of a linear surface wave plasma source for tin removal in an extreme ultraviolet source

Tin deposition mitigation employs hydrogen radicals and ions, formed in a hydrogen plasma, to interact with tin to form tin hydride (SnH 4 ) in the gaseous state, which is then pumped away. Surface wave plasma (SWP) technology developed at Illinois generates hydrogen radicals and ions, resulting in tin etch rates that are high enough to keep extreme ultraviolet (EUV) lithographic tools clean. An advantage of an SWP antenna is the ability to generate a high density of hydrogen radicals and hydrogen ions directly at the desired etching location. In situ etching of tin enables high availability EUV tools by maintaining high reflectivity of the multilayer mirror of the collector. Additionally, the SWP is characterized with low ion energies and low electron temperature, such that the multilayer mirror does not suffer any damage from sputtering or implantation of hydrogen ions during operation. Here, experiments elucidating the fundamental processes of tin removal are conducted by varying pressure, power, surface temperature, and gas flow rate in order to observe the etch rate behavior. Our results have shown that the presence of hydrogen ions increases etch rates because ion bombardment weakens Sn–Sn bonds, which, in turn, allows for a higher rate of chemical etching by the radicals. The ion bombardment reduces the number of radicals needed to etch a single tin atom to the range of 10 2 –10 3 . The linear SWP antenna yields plasma densities on the order of 10 17 to 10 18 m -3 and radical densities on the order of 10 18 to 10 19 m -3 , allowing for greater utilization of ion etch enhancement. Etch rates of up to 200 nm/min have been achieved. The surface temperature of the samples is an important factor in the etching process such that the decrease in the surface temperature increases the etch rates and decreases the hydrogen desorption rates. In addition, a kinetic etch model is developed to explain the behavior of etch rates as a function of surface temperature. Furthermore, results from experiments performed in an Illinois NXE:3100 chamber will be discussed.

42 ENGINEERING↗