SEARCH · Engineering Papers

Results for “invariant”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,477 records · Page 82

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Ordinary and exotic mesons in the extended Linear Sigma Model

The extended Linear Sigma Model (eLSM) is a hadronic model based on the global symmetries of QCD and the corresponding explicit, anomalous, and spontaneous breaking patterns. In its basic three-flavor form, its mesonic part contains the dilaton/glueball as well as the nonets of pseudoscalar, scalar, vector, and axial–vector mesons, thus chiral symmetry is linearly realized. In the chiral limit and neglecting the chiral anomaly, only one term – within the dilaton potential – breaks dilatation invariance, and all terms are chirally symmetric. Spontaneous symmetry breaking is implemented by a generalization of the Mexican-hat potential, with explicit symmetry breaking responsible for its tilting. The overall mesonic phenomenology up to ~2 GeV is in agreement with the PDG compilation of masses and partial and total decay widths. The eLSM was enlarged in a straightforward way to include other conventional quark–antiquark nonets (pseudovector and orbitally excited vector mesons, tensor and axial-tensor mesons, radially excited (pseudo)scalar mesons, etc.), as well as two nonets of hybrid mesons, the lightest one with exotic quantum numbers J $\mathscr{PC}$ = 1 -+ not allowed for $\bar{q}q$ objects, such as the resonance π 1 (1600) and the recently discovered η 1 (1855) . In doing so, different types of chiral multiplets are introduced: heterochiral and homochiral multiplets, which differ in the way they transform under chiral transformations. Moreover, besides the scalar glueball that is present from the beginning as dilaton, other glueballs, the tensor, the pseudoscalar and the vector glueballs were coupled to the eLSM: the scalar resonance ƒ 0 (1710) turns out to be mostly gluonic, the tensor glueball couples strongly to vector mesons, and the pseudoscalar glueball couples sizably to $ππη'$ and can be assigned to X (2370) or X (2600). In all cases above, masses and decays can be analyzed allowing for a better understanding of both conventional and non-conventional mesons: whenever data are available, a comparison is performed and, when this is not the case, predictions of decay widths and decay ratios are outlined. The eLSM contains chiral partners on an equal footing and is therefore well suited for studies of chiral symmetry restoration at nonzero temperature and densities: this is done by coupling it to the Polyakov loop. Finally, the QCD phase diagram and the location of the critical endpoint were investigated within this framework.

Chiral symmetry↗

DOE OSTI2020

Observation of structure in the J/ψ-pair mass spectrum

Using proton-proton collision data at centre-of-mass energies of $\sqrt{s}$ = 7, 8 and 13 TeV recorded by the LHCb experiment at the Large Hadron Collider, corresponding to an integrated luminosity of 9 fb -1 , the invariant mass spectrum of J/ψ pairs is studied. A narrow structure around 6.9 GeV/c 2 matching the lineshape of a resonance and a broad structure just above twice the J/ψ mass are observed. The deviation of the data from nonresonant J/ψ-pair production is above five standard deviations in the mass region between 6.2 and 7.4 GeV/c 2 , covering predicted masses of states composed of four charm quarks. The mass and natural width of the narrow X(6900) structure are measured assuming a Breit-Wigner lineshape.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Evidence of a $J/ψΛ$ structure and observaton of excited $Ξ_{b}^{-}$ → $J/ψΛΚ^{-}$ decay

First evidence of a structure in the $J/ψΛ$ invariant mass distribution is obtained from an amplitude analysis of $Ξ_{b}^{-}$ → $J/ψΛΚ^{-}$ decays. The observed structure is consistent with being due to a charmonium pentaquark with strangeness with a significance of 3.1σ including systematic uncertainties and look elsewhere effect. Its mass and width are determined to be 4458.8 ± $2.9_{-1.1}^{+4.7}$ MeV and 17.3 ± $6.5_{-5.7}^{+8.0}$ MeV respectively, where the quoted uncertainties are statistical and systematic. The structure is also consistent with being due to two resonances. In addition, the narrow excited $Ξ^{-}$ states, $Ξ(1690)^{-}$ and $Ξ(1820)^{-}$, are seen for the first time in a $Ξ_{b}^{-}$ decay, and their masses and widths are measured with improved precision. The analysis is performed using $pp$ collision data corresponding to a total integrated luminosity of 9 fb -1 , collected with the LHCb experiment at centre-of-mass energies of 7, 8 and 13 TeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Variability and associated uncertainty in image analysis for soiling characterization in solar energy systems

The accumulation of soiling on photovoltaic modules and on the mirrors of concentrating solar power systems causes non-negligible energy losses with economic consequences. These challenges can be mitigated, or even prevented, through appropriate actions if the magnitude of soiling is known. Particle counting analysis is a common procedure to characterize soiling, as it can be easily performed on micrographs of glass coupons or solar devices that have been exposed to the environment. Particle counting does not, however, yield invariant results across institutions. The particle size distribution analysis is affected by the operator of the image analysis software and the methodology utilized. The results of a round-robin study are presented in this work to explore and elucidate the uncertainty related to particle counting and its effect on the characterization of the soiling of glass surfaces used in solar energy conversion systems. An international group of soiling experts analyzed the same 8 micrographs using the same open-source ImageJ software package. The variation in the particle analyses results were investigated to identify specimen characteristics with the lowest coefficient of variation (CV) and the least uncertainty among the various operators. The mean particle diameter showed the lowest CV among the investigated characteristics, whereas the number of particles exhibited the largest CV. Additional parameters, such as the fractional area coverage by particles and parameters related to the distribution's shape yielded intermediate CV values. These results can provide insights on the magnitude inter-lab variability and uncertainty for optical and microscope-based soiling monitoring and characterization.

14 SOLAR ENERGY↗

DOE OSTI2021

Model-based and Model-free Designs for an Extended Continuous-time LQR with Exogenous Inputs

We present an extended linear quadratic regulator (LQR) design for continuous-time linear time-invariant (LTI) systems in the presence of exogenous inputs. We first propose a model-based solution with cost minimization guarantees for states and inputs using dynamic programming (DP). The control law consists of a combination of the optimal state feedback and an additional optimal term dependent on the exogenous inputs. The control gains for the two components are obtained by solving a set of matrix differential equations. We provide these solutions for both finite horizons and steady-state cases. In the second part of the paper, we formulate a reinforcement learning (RL) based algorithm which does not need any model information except the input matrix, and can compute an approximate steady-state LQR gain using measurements of the states, the control inputs, and the exogenous inputs. Both model-based and data-driven optimal control algorithms are tested with a numerical example under different exogenous inputs showcasing the effectiveness of the designs.

Mukherjee, Sayak↗

DOE OSTI2026

Sliding friction and wear behavior of nuclear graphite in high temperature inert environment: Influence of contact load, speed and temperature

Repeated dynamic interactions of graphitic components in pebble-bed gas-cooled nuclear reactors can cause abrasive wear-induced pebble surface damage, generate hazardous fine graphite debris, and alter fuel circulation dynamics due to changes in friction behavior. Comprehensive tribological characterization of nuclear graphitic materials in conditions relevant to reactor operation is needed to assess reactor long-term safety and performance. This work reports sliding friction and wear behavior of self-mated nuclear graphite ET-10 at various elevated temperatures (650 °C and 750 °C), sliding speeds (1 and 10 mm/s) and contact loads (20 and 40 N) in a controlled argon environment. The results revealed nonmonotonic frictional behavior with a higher running-in coefficient of friction (COF) followed by a lower steady-state COF, as a result of transition from two-body abrasion to three-body abrasion along with formation of a tribofilm. A key finding of this work is the sensitivity of the running-in COF to experimental conditions; maximum running-in values were lower at either elevated temperature (0.52–0.54) or reduced sliding speed (0.51–0.54). Conversely, the steady-state COF remained invariant at approximately 0.3 across all tested parameters. Transmission electron microscopy revealed a 0.5–2.0 μm thick nanocrystalline tribofilm that was thought to be formed by the compaction of the graphitic wear debris on the contact surface during the sliding process. The nanocrystalline nature of the tribofilm was further confirmed by Raman spectroscopy. As a result, the combination of tribological testing and morphological characterization provided a mechanistic understanding of the frictional behavior of nuclear graphite upon sliding.

Friction↗

DOE OSTI2023

Comparative study of two Saccharomyces cerevisiae strains with kinetic models at genome-scale

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by 13 C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

On the effectiveness of neural operators at zero-shot weather downscaling

Machine-learning (ML) methods have shown great potential for weather downscaling. These data-driven approaches provide a more efficient alternative for producing high-resolution weather datasets and forecasts compared to physics-based numerical simulations. Neural operators, which learn solution operators for a family of partial differential equations, have shown great success in scientific ML applications involving physics-driven datasets. Neural operators are grid-resolution-invariant and are often evaluated on higher grid resolutions than they are trained on, i.e., zero-shot super-resolution. Given their promising zero-shot super-resolution performance on dynamical systems emulation, we present a critical investigation of their zero-shot weather downscaling capabilities, which is when models are tasked with producing high-resolution outputs using higher upsampling factors than are seen during training. To this end, we create two realistic downscaling experiments with challenging upsampling factors (e.g., 8x and 15x) across data from different simulations: the European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) and the Wind Integration National Dataset Toolkit. While neural operator-based downscaling models perform better than interpolation and a simple convolutional baseline, we show the surprising performance of an approach that combines a powerful transformer-based model with parameter-free interpolation at zero-shot weather downscaling. We find that this Swin-Transformer-based approach mostly outperforms models with neural operator layers in terms of average error metrics, whereas an Enhanced Super-Resolution Generative Adversarial Network-based approach is better than most models in terms of capturing the physics of the ground truth data. We suggest their use in future work as strong baselines.

17 WIND ENERGY↗

DOE OSTI2023

Towards critical and supercritical electromagnetic fields

Abstract The availability of ever stronger, laser-generated electromagnetic fields underpins continuing progress in the study and application of nonlinear phenomena in basic physical systems, ranging from molecules and atoms to relativistic plasmas and quantum electrodynamics. This raises the question: how far will we be able to go with future lasers? One exciting prospect is the attainment of field strengths approaching the Schwinger critical field ${E}_{\mathrm{cr}}$ in the laboratory frame, such that the field invariant ${E}^2-{c}^2{B}^2>{E}_{\mathrm{cr}}^2$ is reached. The feasibility of doing so has been questioned, on the basis that cascade generation of dense electron–positron plasma would inevitably lead to absorption or screening of the incident light. Here we discuss the potential for future lasers to overcome such obstacles, by combining the concept of multiple colliding laser pulses with that of frequency upshifting via a tailored laser–plasma interaction. This compresses the electromagnetic field energy into a region of nanometre size and attosecond duration, which increases the field magnitude at fixed power but also suppresses pair cascades. Our results indicate that laser facilities with peak power of tens of PW could be capable of reaching ${E}_{\mathrm{cr}}$ . Such a scenario opens up prospects for the experimental investigation of phenomena previously considered to occur only in the most extreme environments in the universe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

The effect of initial conditions on mixing transition of the Richtmyer–Meshkov instability

We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

A scale-wise analysis of intermittent momentum transport in dense canopy flows

We investigate the intermittent dynamics of momentum transport and its underlying time scales in the near-wall region of the neutrally stratified atmospheric boundary layer in the presence of a vegetation canopy. This is achieved through an empirical analysis of the persistence time scales (periods between successive zero-crossings) of momentum flux events, and their connection to the ejection–sweep cycle. Using high-frequency measurements from the GoAmazon campaign, spanning multiple heights within and above a dense canopy, the analysis suggests that, when the persistence time scales ( $t_p$ ) of momentum flux events from four different quadrants are separately normalized by $\varGamma _{w}$ (integral time scale of the vertical velocity), their distributions $P(t_p/\varGamma _{w})$ remain height-invariant. This result points to a persistent memory imposed by canopy-induced coherent structures, and to their role as an efficient momentum-transporting mechanism between the canopy airspace and the region immediately above. Moreover, $P(t_p/\varGamma _{w})$ exhibits a power-law scaling at times $t_{p}<\varGamma _{w}$ , with an exponential tail appearing for $t_{p} \geq \varGamma _{w}$ . By separating the flux events based on $t_p$ , we discover that around 80 % of the momentum is transported through the long-lived events ( $t_{p} \geq \varGamma _{w}$ ) at heights immediately above the canopy, while the short-lived ones ( $t_{p} < \varGamma _{w}$ ) only contribute marginally ( $\approx 20\,\%$ ). To explain the role of instantaneous flux amplitudes in momentum transport, we compare the measurements with newly developed surrogate data and establish that the range of time scales involved with amplitude variations in the fluxes tends to increase as one transitions from within to above the canopy.

Mechanics↗

DOE OSTI2022

Analysis of scale-dependent kinetic and potential energy in sheared, stably stratified turbulence

Budgets of turbulent kinetic energy (TKE) and turbulent potential energy (TPE) at different scales $\ell$ in sheared, stably stratified turbulence are analysed using a filtering approach. Competing effects in the flow are considered, along with the physical mechanisms governing the energy fluxes between scales, and the budgets are used to analyse data from direct numerical simulation at buoyancy Reynolds number $Re_b=O(100)$ . The mean TKE exceeds the TPE by an order of magnitude at the large scales, with the difference reducing as $\ell$ is decreased. At larger scales, buoyancy is never observed to be positive, with buoyancy always converting TKE to TPE. As $\ell$ is decreased, the probability of locally convecting regions increases, though it remains small at scales down to the Ozmidov scale. The TKE and TPE fluxes between scales are both downscale on average, and their instantaneous values are correlated positively, but not strongly so, and this occurs due to the different physical mechanisms that govern these fluxes. Moreover, the contributions to these fluxes arising from the sub-grid fields are shown to be significant, in addition to the filtered scale contributions associated with the processes of strain self-amplification, vortex stretching and density gradient amplification. Probability density functions (PDFs) of the $Q,R$ invariants of the filtered velocity gradient are considered and show that as $\ell$ increases, the sheared-drop shape of the PDF becomes less pronounced and the PDF becomes more symmetric about $R=0$ .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Onset of instabilities in rotating flows by direct numerical simulation

A rotating disk is the canonical experiment for measuring surface reaction rates in geochemical and electrochemical systems. Using the similarity solution for laminar flow around an infinite disk, the mass transfer coefficient can be simply related to the intrinsic reaction rate at the surface. However, measurements of mass transfer rates use a finite-size disk within a larger container of solution; here the flow is no longer strictly laminar, but there must always be some recirculation. Our interest was initially in the assumption of a uniform radial concentration field, how this breaks down near the perimeter of the disk, and what effect that might have on the measured mass transfer rates. However, our numerical simulations suggest that the flow around a finite-size disk becomes time dependent at Reynolds number ($\textit{Re}$) below 1000, which is much smaller than the typical values in mass-transfer measurements ($\textit{Re}$~10 4 ). In this work, we observe the formation of coherent structures in the flow, which suggest a non-uniform mass transfer at the disk surface. The rotating-disk flow follows a similar sequence of instabilities to the Taylor–Couette flow: a centrifugal instability leading an axisymmetric, time-invariant flow, followed by a Hopf bifurcation to a time-periodic flow. To minimise the possibility that our results are a numerical artefact, we have also simulated the instability in the stationary boundary layer of a rotor–stator flow, comparing with self-similar solutions at low $\textit{Re}$ and with spectral methods near the critical Reynolds number.

42 ENGINEERING↗

DOE OSTI2024

Reynolds-number scaling of wall-pressure–velocity correlations in wall-bounded turbulence

Wall-pressure fluctuations are a practically robust input for real-time control systems aimed at modifying wall-bounded turbulence. The scaling behaviour of the wall-pressure–velocity coupling requires investigation to properly design a controller with such input data so that it can actuate upon the desired turbulent structures. A comprehensive database from direct numerical simulations (DNS) of turbulent channel flow is used for this purpose, spanning a Reynolds-number range$Re_\tau \approx 550\unicode{x2013}5200$. Spectral analysis reveals that the streamwise velocity is most strongly coupled to the linear term of the wall pressure, at a Reynolds-number invariant distance-from-the-wall scaling of$\lambda _x/y \approx 14$(and$\lambda _x/y \approx 8$for the wall-normal velocity). When extending the analysis to both homogeneous directions in$x$and$y$, the peak coherence is centred at$\lambda _x/\lambda _z \approx 2$and$\lambda _x/\lambda _z \approx 1$for$p_w$and$u$, and$p_w$and$v$, respectively. A stronger coherence is retrieved when the quadratic term of the wall pressure is concerned, but there is only little evidence for a wall-attached-eddy type of scaling. An experimental dataset comprising simultaneous measurements of wall pressure and velocity complements the DNS-based findings at one value of$Re_\tau \approx 2$k, with ample evidence that the DNS-inferred correlations can be replicated with experimental pressure data subject to significant levels of (acoustic) facility noise. It is furthermore shown that velocity-state estimations can be achieved with good accuracy by including both the linear and quadratic terms of the wall pressure. An accuracy of up to 72 % in the binary state of the streamwise velocity fluctuations in the logarithmic region is achieved; this corresponds to a correlation coefficient of$\approx$0.6. This thus demonstrates that wall-pressure sensing for velocity-state estimation – e.g. for use in real-time control of wall-bounded turbulence – has merit in terms of its realization at a range of Reynolds numbers.

Mechanics↗

DOE OSTI2025

Asymptotic coefficients of the attached-eddy model derived from an adiabatic atmosphere

The attached-eddy model (AEM) predicts that the mean streamwise velocity and streamwise velocity variance profiles follow a logarithmic shape, while the vertical velocity variance remains invariant with height in the overlap region of high Reynolds number wall-bounded turbulent flows. Moreover, the AEM coefficients are presumed to attain asymptotically constant values at very high Reynolds numbers. Here, the AEM predictions are examined using sonic anemometer measurements in the near-neutral atmospheric surface layer, with a focus on the logarithmic behaviour of the streamwise velocity variance. Utilizing an extensive 210-day dataset collected from a 62 m meteorological tower located in the Eastern Snake River Plain, Idaho, USA, the inertial sublayer is first identified by analysing the measured momentum flux and mean velocity profiles. The logarithmic behaviour of the streamwise velocity variance and the associated ‘−1’ scaling of the streamwise velocity energy spectra are then investigated. The findings indicate that the Townsend–Perry coefficient (A 1 ) is influenced by mild non-stationarity that manifests itself as a Reynolds number dependence. After excluding non-stationary runs, and requiring the bulk Reynolds number defined using the atmospheric boundary layer height to be larger than 4 × 10 7 , the inferred A 1 converges to values ranging between 1 and 1.25, consistent with laboratory experiments. Furthermore, nine benchmark cases selected through a restrictive quality control reveal a close relation between the ‘−1’ scaling in the streamwise velocity energy spectrum and the logarithmic behaviour of streamwise velocity variance. However, additional data are required to determine whether the plateau value of the pre-multiplied streamwise velocity energy spectrum is identical to A 1 .

atmospheric flows↗

DOE OSTI2026

Modulation induced by very-large-scale motions on the inclination angle of wall-attached eddies: an atmospheric surface layer experiment

The forward leaning inclination angle, 𝛾, of coherent turbulent structures is a well-known feature of wall-bounded turbulent flows. Although invariant across friction Reynolds numbers within the range 𝑅𝑒 𝜏 =10 3 −10 6 , 𝛾 can vary significantly across turbulent scales within a high-Reynolds-number flow. Very-large-scale motions (VLSMs) are known to induce significant changes in the instantaneous shear profile, which is a conditioning event that could trigger variability in the inclination angle of smaller coherent turbulent structures. Although this aspect has been extensively studied via numerical and laboratory experiments, few studies have explored this feature for a very-high-Reynolds-number atmospheric flow. In this work, the inclination angle of turbulent structures within the atmospheric surface layer at a very high Reynolds number (𝑅𝑒 𝜏 =7.9 ×10 5 ) is investigated by deploying a scanning Doppler light detection and ranging and a super large particle image velocimetry (SLPIV) apparatus. The inclination angle of wall-attached eddies is inferred either from the two-point correlation of streamwise velocity (𝛾 =41.1∘) or with a scale-dependent approach through the spectral linear stochastic estimator (SLSE). The SLSE (and, thus, the scale-dependent inclination angle) is conditionally evaluated based on the high- and low-momentum events induced by VLSMs, both in the streamwise (𝑢′𝑉𝐿𝑆𝑀) and in the vertical (𝑤′ 𝑉𝐿𝑆𝑀) velocity components. As a result, lower inclination angles (𝛾 =30° −50°) are found for 𝑢$^{'}_{VLSM}$ >0 (𝑤$^{'}_{VLSM}$ <0), while higher values (50° −85°) are ascribed to 𝑢$^{'}_{VLSM}$ <0 (𝑤$^{'}_{VLSM}$ >0). This result emphasises the primary role that VLSMs play in shaping the wall-attached eddy geometry, which, in turn, is crucial to determine the Reynolds stress balance within the wall-attached eddy range.

42 ENGINEERING↗

DOE OSTI2022

Discovery of peculiar radio morphologies with ASKAP using unsupervised machine learning

Abstract We present a set of peculiar radio sources detected using an unsupervised machine learning method. We use data from the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to train a self-organizing map (SOM). The radio maps from three ASKAP surveys, Evolutionary Map of Universe pilot survey (EMU-PS), Deep Investigation of Neutral Gas Origins pilot survey (DINGO), and Survey With ASKAP of GAMA-09 + X-ray (SWAG-X), are used to search for the rarest or unknown radio morphologies. We use an extension of the SOM algorithm that implements rotation and flipping invariance on astronomical sources. The SOM is trained using the images of all ‘complex’ radio sources in the EMU-PS which we define as all sources catalogued as ‘multi-component’. The trained SOM is then used to estimate a similarity score for complex sources in all surveys. We select 0.5% of the sources that are most complex according to the similarity metric and visually examine them to find the rarest radio morphologies. Among these, we find two new odd radio circle (ORC) candidates and five other peculiar morphologies. We discuss multiwavelength properties and the optical/infrared counterparts of selected peculiar sources. In addition, we present examples of conventional radio morphologies including: diffuse emission from galaxy clusters, and resolved, bent-tailed, and FR-I and FR-II type radio galaxies. We discuss the overdense environment that may be the reason behind the circular shape of ORC candidates.

Astronomy & Astrophysics↗

DOE OSTI2020

Multi-scale dynamics of magnetic flux tubes and inverse magnetic energy transfer

Here, we report on an analytical and numerical study of the dynamics of a three-dimensional array of identical magnetic flux tubes in the reduced-magnetohydrodynamic description of the plasma. We propose that the long-time evolution of this system is dictated by flux-tube mergers, and that such mergers are dynamically constrained by the conservation of the pertinent (ideal) invariants, viz. the magnetic potential and axial fluxes of each tube. We also propose that in the direction perpendicular to the merging plane, flux tubes evolve in a critically balanced fashion. These notions allow us to construct an analytical model for how quantities such as the magnetic energy and the energy-containing scale evolve as functions of time. Of particular importance is the conclusion that, like its two-dimensional counterpart, this system exhibits an inverse transfer of magnetic energy that terminates only at the system scale. We perform direct numerical simulations that confirm these predictions and reveal other interesting aspects of the evolution of the system. We find, for example, that the early time evolution is characterized by a sharp decay of the initial magnetic energy, which we attribute to the ubiquitous formation of current sheets. We also show that a quantitatively similar inverse transfer of magnetic energy is observed when the initial condition is a random, small-scale magnetic seed field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗