Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary effects”

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 703 records · Page 39

Numerical Evaluation of Effective Thermal Conductivity of PCM with Metal Foam Incorporating Buoyancy Effects for Thermal Energy Storage

The thermal energy storage (TES) system has the capability to efficiently preserve thermal energy directly derived from the energy source, minimizing any conversion losses. Especially latent heat storage offers distinct advantages, including a substantial increase in energy storage density and minimization of temperature fluctuations within the plants. However, the phase change material (PCM) employed in latent heat storage has low thermal conductivity. Consequently, various studies are being conducted to enhance heat transfer. One approach to enhance heat transfer involves utilizing metal foam to maximize the heat transfer area. However, modeling metal foam with its intricate structure is a challenging task in numerical analysis. For this reason, ongoing research focuses on simplifying the modeling of metal foam. Nevertheless, fully encompassing all the characteristics of actual metal foam proves to be a challenging task for the simplified analytical model. The objective of this paper is to interpret the simple lattice metal foam analysis model from the perspective of behavior induced by buoyancy. When comparing the analysis results of solid PCM and liquid PCM with the same thermal conductivity under changes in porosity and gravity direction, we conducted an analysis to discern the trends in effective thermal conductivity that are overestimated due to convection. In the analysis, a constant heat flux of 10 kW and a constant surface boundary condition of 350 K were applied, and a sensitivity study regarding the mesh was conducted. The results indicate that, from the perspective of gravity in the simple lattice model, the solid analysis yields an effective thermal conductivity 29-47% higher compared to the liquid analysis. Additionally, as porosity increases, there is an observed increase of 24-33% in effective thermal conductivity.

25 ENERGY STORAGE↗

Cement stress and microstructure evolution during curing in semi-rigid high-pressure environments

The process of cement slurry curing into a solid has a direct effect on the initial state of stress in a composite concrete and steel structure. Failure models for critical infrastructure and wellbores could benefit from measurements of this initial stress because it greatly influences structural performance. Here, we present measurements of the initial stress state (total stress, effective stress, and pore pressure) of multiple cement formulations cured under a high-pressure of 40 MPa within a semi-rigid steel pipe. Our results show that cements cured in this restrained environment exhibit dense microstructures with low permeability (nano-Darcy), develop an anisotropic stress state greater than zero and less than the curing pressure, and undergo a staged curing process where hydration and mechanical behaviors evolve from a slurry into a pre-stressed porous medium. Both drained and undrained scenarios were investigated, revealing a feedback from the semi-rigid boundary towards improved cement performance.

36 MATERIALS SCIENCE↗

High-temperature deformation mechanisms in a BCC+B2 refractory complex concentrated alloy

The current study reports the analysis of the deformation mechanisms at 600 °C in a two-phase, BCC+B2, refractory complex concentrated alloy (RCCA) Al 0.5 NbTa 0.8 Ti 1.5 V 0.2 Zr. At this temperature, the alloy microstructure is unstable and dynamic coarsening of B2 precipitates is evidenced during the mechanical testing. After true plastic strain of 0.030 at strain rate of 10 -4 s -1 , the deformation becomes highly localized in wavy bands reflecting the profusion of cross-slip. Scanning transmission electron microscopy (STEM) observations highlight the presence of paired a/2$\langle$111$\rangle$ dislocations that shear the B2 precipitates in a cooperative process. Additionally, some chemical segregation effect is observed along the narrow dislocation bands likely induced to decrease the antiphase boundary (APB) energy of the system.

ordering↗

Combined Cr and S poisoning behaviors of La 1-x Sr x MnO 3±δ and La 1-x Sr x Co 1-y FeyO 3-δ cathodes in solid oxide fuel cells

Although the individual effects of airborne Cr and S contaminants on SOFC cathode performance degradation have been extensively studied, the combined effects of Cr and S contaminants remain largely unexplored. Under the real SOFC operating condition where the Cr and S species coexist, their effects may compete, affecting the poisoning behavior. Here, our investigation reveals that the combined Cr and S poisoning behavior of LSCF remains different from those of individual Cr and S poisonings, while the combined poisoning mechanism of LSM is equivalent to the sum of those of individual Cr and S effects. For LSCF electrode, gaseous Cr species are deposited mainly at the LSCF/GDC interface by electrochemical reduction, rather than forming SrCrO 4 on LSCF surfaces (as in the case of Cr-only poisoning), indicating no reaction between Cr vapors and SrO on the LSCF surface. Thermodynamic analysis demonstrates that the SrO on LSCF surface absorbs SO 2 (g) and thereby loses the Cr-gettering effect, allowing Cr vapors to flow through LSCF and to reach the LSCF/GDC interface where the Cr deposition occurs. Unlike LSCF, LSM electrode shows cumulative effects of Cr and S, as Cr accumulation occurs at the triple-phase boundary and S absorption takes place at localized Sr-rich regions.

36 MATERIALS SCIENCE↗

Synthesis and morphological characterization of linear and miktoarm star poly(solketal methacrylate)- block -polystyrene copolymers

The synthesis, molecular characterization, and morphological evaluation of AB n, (n = 2,3) miktoarm star block copolymers consisting of poly(glycerol monomethacrylate) (PGMA) and polystyrene (PS) with varying molecular weights and compositions is described. The system is known to demonstrate a remarkably high Flory-Huggins interaction parameter. The corresponding linear diblock copolymer analogues, were synthesized as well, and comparisons with regard to feature dimensions and morphologies are provided. Well-ordered nanostructures of various morphologies were formed with domain spacing as low as 7.2 nm. Different morphologies were attained by some of the topological isomers indicating that in miktoarm star block copolymers the phase boundaries were strongly shifted. Additionally, a triblock ABA analogue was studied to investigate the effect of triblock copolymer conformations. Noteworthy is that for copolymers with different macromolecular architecture leading to similar morphology; different domain spacings were obtained. The synthesis of all the samples was carried out by high-vacuum anionic polymerization techniques. Molecular characterization with Size Exclusion Chromatography (SEC) and Proton Nuclear Magnetic Resonance Spectroscopy ( 1 H NMR) confirmed well-defined copolymers obtained. The morphological characterization was accomplished by Small-Angle X-ray Scattering (SAXS). In conclusion, the observations from this study highlight the potential of incorporating macromolecular architecture in the self-assembly of strongly immiscible block copolymers to attain ultra-small nanofeatures with desired morphologies.

36 MATERIALS SCIENCE↗

Defect cluster and nonequilibrium gas bubble associated growth in irradiated UMo fuels – a cluster dynamics and phase field model

Irradiation examination shows that gas bubble swelling kinetics is much faster after irradiation-induced recrystallization than that prior recrystallization in UMo fuels. It implies that gas bubbles in coarse grains and small recrystallized grains have different growth behavior. In this work, a phase-field model of gas bubble evolution integrating microstructure dependent cluster dynamics has been developed, for the first time, to study the gas bubble swelling behavior in the recrystallization zone of UMo fuels. Generation, diffusion, reaction, sink, emission and clustering of vacancies and interstitials are described by the cluster dynamics model while a phase-field model is used to describe the evolution of non-equilibrium gas bubbles including nucleation and growth. With the coupled model, the effect of defect generation rate, clustering rate, interstitial emission and sink rates on grain boundaries on the gas bubble evolution are systematically simulated. A set of model parameters (defect generation rate, clustering rate, interstitial emission and sink rates) is determined by comparing measured and simulated gas bubble swelling kinetics. The results demonstrate that interstitial clustering is one of the important physical mechanisms which results in a fast gas bubble swelling kinetics in the recrystallization zone. The developed model can also be extended to study the associated growth of defect and second phase precipitates often observed in irradiated materials.

Hu, Shenyang↗

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗

Large-Eddy Simulations of Marine Boundary Layer Clouds Associated with Cold-Air Outbreaks during the ACTIVATE Campaign. Part II: Aerosol–Meteorology–Cloud Interaction

Aerosol effects on micro/macrophysical properties of marine stratocumulus clouds over the western North Atlantic Ocean (WNAO) are investigated using in situ measurements and large-eddy simulations (LES) for two cold-air outbreak (CAO) cases (28 February and 1 March 2020) during the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE). The LES is able to reproduce the vertical profiles of liquid water content (LWC), effective radius r eff and cloud droplet number concentration N c from fast cloud droplet probe (FCDP) in situ measurements for both cases. Furthermore, we show that aerosols affect cloud properties (N c , r eff , and LWC) via the prescribed bulk hygroscopicity of aerosols ($\bar{κ}$) and aerosol size distribution characteristics. N c , r eff , and liquid water path (LWP) are positively correlated to $\bar{κ}$ and aerosol number concentration (N a ) while cloud fractional cover (CFC) is insensitive to $\bar{κ}$ and aerosol size distributions for the two cases. The realistic changes to aerosol size distribution (number concentration, width, and the geometrical diameter) with the same meteorology state allow us to investigate aerosol effects on cloud properties without meteorological feedback. We also use the LES results to evaluate cloud properties from two reanalysis products, ERA5 and MERRA-2. Compared to LES, the ERA5 is able to capture the time evolution of LWP and total cloud coverage within the study domain during both CAO cases while MERRA-2 underestimates them.

54 ENVIRONMENTAL SCIENCES↗

Boundary-induced classical generalized Gibbs ensemble with angular momentum

We investigate how confinement geometry leads to the emergence of a Generalized Gibbs Ensemble (GGE) in classical systems. Unlike the standard Gibbs ensemble, the GGE includes additional conserved quantities, such as angular momentum, that arise from boundary-induced symmetries. Using analytical arguments based on the maximum entropy principle, we show that circular boundaries preserve angular momentum and drive the system toward a chiral, non-ergodic GGE that violates time-reversal symmetry. This ensemble differs fundamentally from the Gibbs case, producing near-boundary condensation and revealing how geometry alone can alter thermal equilibration. To quantify these effects, we introduce an order parameter measuring deviations from Gibbs behavior and demonstrate that conventional Monte Carlo methods must incorporate angular momentum conservation under such conditions. Our study highlights how geometric constraints shape non-equilibrium statistical ensembles and lead to subtle departures from the Bohr-van Leeuwen theorem. These predictions are validated through detailed simulations of confined classical hard-disk gases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Implications of Microstructure in Helium-Implanted Nanocrystalline Metals

Helium bubbles are known to form in nuclear reactor structural components when displacement damage occurs in conjunction with helium exposure and/or transmutation. If left unchecked, bubble production can cause swelling, blistering, and embrittlement, all of which substantially degrade materials and—moreover—diminish mechanical properties. On the mission to produce more robust materials, nanocrystalline (NC) metals show great potential and are postulated to exhibit superior radiation resistance due to their high defect and particle sink densities; however, much is still unknown about the mechanisms of defect evolution in these systems under extreme conditions. Here, the performances of NC nickel (Ni) and iron (Fe) are investigated under helium bombardment via transmission electron microscopy (TEM). Bubble density statistics are measured as a function of grain size in specimens implanted under similar conditions. While the overall trends revealed an increase in bubble density up to saturation in both samples, bubble density in Fe was over 300% greater than in Ni. To interrogate the kinetics of helium diffusion and trapping, a rate theory model is developed that substantiates that helium is more readily captured within grains in helium-vacancy complexes in NC Fe, whereas helium is more prone to traversing the grain matrices and migrating to GBs in NC Ni. Our results suggest that (1) grain boundaries can affect bubble swelling in grain matrices significantly and can have a dominant effect over crystal structure, and (2) an NC-Ni-based material can yield superior resistance to irradiation-induced bubble growth compared to an NC-Fe-based material and exhibits high potential for use in extreme environments where swelling due to He bubble formation is of significant concern.

36 MATERIALS SCIENCE↗

Radiation Damage Mitigation in FeCrAl Alloy at Sub-Recrystallization Temperatures

Traditional defect recovery methods rely on high-temperature annealing, often exceeding 750 °C for FeCrAl. In this study, we introduce electron wind force (EWF)-assisted annealing as an alternative approach to mitigate irradiation-induced defects at significantly lower temperatures. FeCrAl samples irradiated with 5 MeV Zr 2+ ions at a dose of 10 14 cm −2 were annealed using EWF at 250 °C for 60 s. We demonstrate a remarkable transformation in the irradiated microstructure, where significant increases in kernel average misorientation (KAM) and low-angle grain boundaries (LAGBs) typically indicate heightened defect density; the use of EWF annealing reversed these effects. X-ray diffraction (XRD) confirmed these findings, showing substantial reductions in full width at half maximum (FWHM) values and a realignment of peak positions toward their original states, indicative of stress and defect recovery. To compare the effectiveness of EWF, we also conducted traditional thermal annealing at 250 °C for 7 h, which proved less effective in defect recovery as evidenced by less pronounced improvements in XRD FWHM values.

FeCrAl alloys↗

Response of marine post-frontal clouds to Gulf Stream variability

Understanding how Gulf Stream variation influences cloud morphology is critical for evaluating cloud feedback in the western North Atlantic Ocean and beyond, where mesoscale air-sea interactions dominate. This study investigates the impact of altered mean sea surface temperature (SST) and SST gradients on post-frontal cloud characteristics during cold-air outbreaks, using the Weather Research and Forecasting (WRF) model. Three sensitivity experiments are conducted: a control simulation (default SST), Plus4 (uniform SST increase of 4 K), and Gradplus (SST gradient enhanced by 25 %, centered around mean SST). Results reveal distinctly different responses in boundary layer dynamics and cloud macro-physics. In Plus4, a warmer and moister boundary layer reduces total cloud cover but promotes larger cloud sizes and elongated cloud streets, with diminished liquid water and enhanced ice-phase hydrometeors. Conversely, Gradplus amplifies impacts in the upwind colder SST regions, yielding a drier, colder boundary layer, weaker energy transport, and higher liquid water path but reduced ice water content and cloud lines. Tracer analysis highlights that SST modifications alter airmass sources near cloud tops due to the entrainment of ambient air, with Plus4 amplifying boundary layer contributions to cloud-top regions. These findings underscore the spatially varying effects of SST gradients and mean SST on cloud organization and microphysics, emphasizing the need to resolve ocean-atmosphere coupling in global models to improve the prediction of marine cloud feedback under warming scenarios.

AIR-SEA INTERACTION↗

Investigating the physical mechanisms that modify wind plant blockage in stable boundary layers

Abstract. Wind plants slow down the approaching wind, a phenomenon known as blockage. Wind plant blockage undermines turbine performance for front-row turbines and potentially for turbines deeper into the array. We use large-eddy simulations to characterize blockage upstream of a finite-size wind plant in flat terrain for different atmospheric stability conditions and investigate the physical mechanisms modifying the flow upstream of the turbines. To examine the influence of atmospheric stability, we compare simulations of two stably stratified boundary layers using the Weather Research and Forecasting model in large-eddy simulation mode, representing wind turbines using the generalized actuator disk approach. For a wind plant, a faster cooling rate at the surface, which produces stronger stably stratified flow in the boundary layer, amplifies blockage. As a novelty, we investigate the physical mechanisms amplifying blockage by evaluating the different terms in the momentum conservation equation within the turbine rotor layer. The velocity deceleration upstream of a wind plant is caused by an adverse pressure gradient and momentum advection out of the turbine rotor layer. The cumulative deceleration of the flow upstream of the front-row turbines instigates vertical motions. The horizontal flow is diverted vertically, reducing momentum availability in the turbine rotor layer. Although the adverse pressure gradient upstream of the wind plant remains unchanged with atmospheric stability, vertical advection of horizontal momentum is amplified in the more strongly stable boundary layer, mainly by larger shear of the horizontal velocity, thus increasing the blockage effect.

17 WIND ENERGY↗

Expression of interfacial Seebeck coefficient through grain boundary engineering with multi-layer graphene nanoplatelets

Energy filtering has been a long-sought strategy to enhance a thermoelectric material's figure of merit zT through improving its power factor. Here we show a composite of multi-layer graphene nanoplatelets (GNP) and n-type Mg 3 Sb 2 leads to the expression of an energy filtering like effect demonstrated by an increase in the material's Seebeck coefficient and maximum power factor, without impact on the material's carrier concentration. We analyse these findings from the perspective of a heterogeneous material consisting of grain and grain boundary phases, instead of a more traditional and common analysis that assumes a homogeneously transporting medium. An important implication of this treatment is that it leads to the development of an interfacial Seebeck coefficient term, which can explain the observed increase in the material's Seebeck coefficient. The contribution of this interfacial Seebeck coefficient to the overall Seebeck coefficient is determined by the relative temperature drop across the grain boundary region compared to that of the bulk material. In Te doped Mg 3 Sb 2 we show the introduction of GNP increases the interfacial thermal resistance of grain boundaries, enhancing the contribution of the interfacial Seebeck coefficient arising from grain boundaries to the overall Seebeck coefficient. Without significant detriment to the electrical conductivity this effect results in a net increase in maximum power factor. This increased interfacial thermal resistance also leads to the synergistic reduction of the total thermal conductivity. As a result, we enhance zT of the Mg 3 Sb 2 to a peak value of 1.7 near 750 K. Considering the two-dimensional nature of the grain boundary interface, this grain boundary engineering strategy could be applied to a few thermoelectric systems utilizing various two-dimensional nanomaterials.

42 ENGINEERING↗

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗

The effect of remnant CdSe layers on the performance of CdSeTe/CdTe photovoltaic devices

Thin film CdTe-based photovoltaic devices have achieved high efficiency above 22 %. The recent improvement in efficiency is due to Se alloying in the CdTe absorbers to form a CdSeTe/CdTe structure. The subsequent band gap grading increases the short circuit current density. The Se can be introduced by depositing a precursor thin film of either CdSe or a CdSeTe alloy and then diffusing the Se into the CdTe during the high temperature cadmium chloride activation process. Using CdSe is preferred because it is easier to control the Se concentration. However, during fabrication of the CdSeTe/CdTe devices, the CdSe thickness needs to be precisely controlled to prevent the retention of a CdSe remnant layer after the activation treatment. Retention of a remnant CdSe layer causes a dramatic reduction in device efficiency. In this work, we show that the reduction in efficiency is caused by a number of factors. The remnant CdSe layer is n-type which moves the position of the p-n junction. Also, it is widely thought that the CdSe remnants are photo-inactive. In this work, we clarify that the individual CdSe grains are actually highly photo-active. However, the grain sizes in the CdSe remnant and the adjacent CdSeTe layer are very small resulting in a high grain boundary area. Although the grain boundaries are passivated with chlorine, cathodoluminescence imaging and electrical measurements show that this is only partially effective. Also, EQE measurements show that the remnant CdSe causes parasitic absorption. Overall, the remnant CdSe layer causes a reduction in short circuit current density and device efficiency. The thickness of the CdSe precursor layer and the cadmium chloride activation process conditions must be precisely optimised to ensure that all the CdSe is consumed and inter-diffused to form the CdSeTe alloy for highest efficiency devices.

CdSe remnant↗

Focusing of laser-accelerated ions with spherically curved targets

A theoretical model for laser-accelerated ion focusing from spherically curved targets using a Gaussian self-similar solution is presented. This model describes the evolution of the focal location and focal radius with final ion energy and the initial radial boundary of the ions, the latter being a function of both target geometry and the ion acceleration radius for a given ion energy. The theory is supported by particle-in-cell simulations of a variety of target shapes with varying radii of curvature and target opening angles, as well as variations in the injected electron beam radius and energy spectrum. The theory and simulations suggest that the focal location varies linearly with the radius of curvature, with the square root of the ratio of ion energy to effective electron temperature, and monotonically increases with the energy-dependent ratio of the initial ion radial boundary to the radius of curvature. Considering ponderomotive acceleration, this suggests that the focal length should scale inversely with intensity, $(I_Lλ_L^2)$ −1/4 , suggesting a 10× increase in intensity will reduce the focal distance by ̃1.8× for a given ion energy.

Physics - Plasma physics↗

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

ABSTRACT The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

Schoeffler, K. M. (ORCID:0000000276447381)↗