Search NASA⌕ Search

SEARCH · Search NASA

Results for “Instability waves”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Boundary Layer Receptivity due to Roughness and Freestream Sound for Supersonic Flows over Axisymmetric Cones

The receptivity of supersonic boundary layers due to the interaction of a two-dimensional acoustic wave with a three-dimensional roughness is numerically investigated over a 5-degree sharp tipped straight cone at a free stream Mach number of 3.5 and a high Reynolds number of 10(exp 6)/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. The simulations are performed with different roughness heights and roughness locations. The steady computations with the roughness reveal that the perturbations induced by the roughness remain very close to the wall and do not penetrate to the outer part of the boundary layer. The effects of the roughness persist for about 10 boundary layer thicknesses downstream of each element. The unsteady simulation shows the generation of unstable three-dimensional first mode disturbances within a short distance from the roughness. However, the amplitudes of the instability waves are very small near the neutral point, on the order of 10(exp -4) compared to the free stream disturbances. The receptivity coefficients due to the roughness elements located very close to the tip are one order of magnitude smaller than those due to the roughness elements located closer to the neutral point.

Balakumar, P.↗
NASA NTRSDate not supplied

Large Eddy Simulations of Upper-Ocean Turbulence Across the Equatorial Pacific Cold Tongue

Turbulence in the equatorial Pacific cold tongue plays multiple roles in the sea-surface temperature budget. However, observations of turbulence are limited to a handful of field experiments lasting from a few days to one month at 0 N , 140 W, and the key turbulent fluxes can only be inferred. Although explicit simulations of turbulence resolve these fluxes, the required large eddy simulations (LES) are costly and thus also limited to a few weeks. A two-fold hypothesis motivates this contribution to the PUMP pre-field modelling experiments in support of TPOS: 1) that there exist useful relationships between the surface wind, wave, and buoyancy forcing, oceanic shear and stratification, and turbulent mixing, but 2) these relationships differ between regions of equatorial dynamics (f~0) and higher latitudes, such that the state space is insufficiently sampled to understand these relationships and represent them in ocean general circulation models (OGCM). To address the lack of observations, we simulated 35 continuous days of upper ocean turbulence using LES during October-November 1985 at sites in two regions (70 days in total) of the cold tongue: 0 N, 140 W and 3 N, 140 W. The model reproduces many features of the observed turbulence at 0,140W, including the “deep-cycle” observed by Gregg et al. (1985) Nature and the turbulent pulses every few hours during the night and in the early morning. Furthermore, the results at 3 N support the hypothesis advanced by Cherian et al. (2021) JPO that strong, but intermittent deep-cycle turbulence is observed to the north of the equator in the presence of strong shear associated with Tropical Instability Waves. Comparisons with the K-profile parameterization (KPP) scheme show that additional tuning of the KPP parameters can reduce differences between the turbulent fluxes produced by KPP and LES at both sites. The results will suggest modifications to OGCMs and already emphasize the need for a more extensive measurements of vertical profiles of shear, stratification and turbulent mixing to reveal the spatio-temporal structure of turbulence and how it relates to forcing, shear and stratification across the equatorial Pacific Cold tongue.

Large Eddy Simulations↗
NASA NTRS2019

Disentangling Planets from Photoelectric Instability in Gas-rich Optically Thin Dusty Disks

Structures in circumstellar disks such as gaps and rings are often attributed to planets. This connection has been difficult to show unequivocally, as other processes may also produce these features. In particular, a photoelectric instability (PEI) has been proposed, operating in gas-rich optically thin disks, that generates structures predicted by planet–disk interactions. We examine the question of how to disentangle the planetary effects on disk structure from the effects of the PEI. We use the Pencil Code to perform 2D global hydrodynamical models of the dynamics of gas and dust in a thin disk with and without planetary perturbers. Photoelectric heating is modeled with an equation of state where pressure is proportional to dust surface density. The drag force on grains and its backreaction on the gas are included. Analyzing the situation without PEI, we find that gas–dust interactions alter the shape of the planetary gap from the dust-free case when the local dust-to-gas ratio ε approaches unity. This result also applies to primordial disks, because dust drifting inward accumulates at the edge of the planetary gap, and any initial dust-to-gas ratio eventually achieves ε = 1 if the dust reservoir is sufficient. We find a result particular to high dust-to-gas ratio disks as well: as dust drifts inward, the dust front becomes a sharp transition, and the backreaction triggers the Rossby wave instability. When PEI is included, we find that it obscures structures induced by planets unless the planet's mass is sufficiently large to carve a noticeable gap. Specifically, the instability generates arcs and rings of regular spacing: a planet is discernible when it carves a dust gap wider than the wavelength of the PEI.

Areli Castrejon↗
NASA NTRS2003

The Langley Stability and Transition Analysis Code (LASTRAC) : LST, Linear and Nonlinear PSE for 2-D, Axisymmetric, and Infinite Swept Wing Boundary Layers

During the past two decades, our understanding of laminar-turbulent transition flow physics has advanced significantly owing to, in a large part, the NASA program support such as the National Aerospace Plane (NASP), High-speed Civil Transport (HSCT), and Advanced Subsonic Technology (AST). Experimental, theoretical, as well as computational efforts on various issues such as receptivity and linear and nonlinear evolution of instability waves take part in broadening our knowledge base for this intricate flow phenomenon. Despite all these advances, transition prediction remains a nontrivial task for engineers due to the lack of a widely available, robust, and efficient prediction tool. The design and development of the LASTRAC code is aimed at providing one such engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. LASTRAC was written from scratch based on the state-of-the-art numerical methods for stability analysis and modem software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory (LST) or linear parabolized stability equations (LPSE) method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. Coupled with the built-in receptivity model that is currently under development, the nonlinear PSE method offers a synergistic approach to predict transition onset for a given disturbance environment based on first principles. This paper describes the governing equations, numerical methods, code development, and case studies for the current release of LASTRAC. Practical applications of LASTRAC are demonstrated for linear stability calculations, N-factor transition correlation, non-linear breakdown simulations, and controls of stationary crossflow instability in supersonic swept wing boundary layers.

Chang, Chau-Lyan↗
DOE OSTI2026

Electronic Layer Decoupling Driven by Density-Wave Order in La 4⁢ Ni 3 ⁢O 10

We probe the density-wave transition of the trilayer nickelate La 4⁢ Ni 3 ⁢O 10 with polarization-resolved infrared spectroscopy. The low-energy electrodynamics is strongly anisotropic, with metallic in-plane and insulating out-of-plane character. In the ordered phase, the anisotropy grows more than an order of magnitude as the out-of-plane conductivity is sharply suppressed. We interpret this enhancement as an effective electronic decoupling of the Ni-O layers driven by a spin-density-wave-induced redistribution of Ni−𝑑 𝑧 2 occupation within the trilayers. This electronic response is accompanied by clearly shifting and splitting out-of-plane phonons, compatible with a density-wave instability of electronic origin.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A comprehensive and fair comparison of two neural operators (with practical extensions) based on $\mathrm{FAIR}$ data

Neural operators can learn nonlinear mappings between function spaces and offer a new simulation paradigm for real-time prediction of complex dynamics for realistic diverse applications as well as for system identification in science and engineering. Herein, we investigate the performance of two neural operators, which have shown promising results so far, and we develop new practical extensions that will make them more accurate and robust and importantly more suitable for industrial-complexity applications. The first neural operator, DeepONet, was published in 2019 (Lu et al., 2019), and its original architecture was based on the universal approximation theorem of Chen & Chen (1995). The second one, named Fourier Neural Operator or FNO, was published in 2020, and it is based on parameterizing the integral kernel in the Fourier space. DeepONet is represented by a summation of products of neural networks (NNs), corresponding to the branch NN for the input function and the trunk NN for the output function; both NNs are general architectures, e.g., the branch NN can be replaced with a CNN or a ResNet. According to Kovachki et al. (2021), FNO in its continuous form can be viewed conceptually as a DeepONet with a specific architecture of the branch NN and a trunk NN represented by a trigonometric basis. In order to compare FNO with DeepONet computationally for realistic setups, we develop several extensions of FNO that can deal with complex geometric domains as well as mappings where the input and output function spaces are of different dimensions. We also develop an extended DeepONet with special features that provide inductive bias and accelerate training, and we present a faster implementation of DeepONet with cost comparable to the computational cost of FNO, which is based on the Fast Fourier Transform. Here we consider 16 different benchmarks to demonstrate the relative performance of the two neural operators, including instability wave analysis in hypersonic boundary layers, prediction of the vorticity field of a flapping airfoil, porous media simulations in complex-geometry domains, etc. We follow the guiding principles of FAIR (Findability, Accessibility, Interoperability, and Reusability) for scientific data management and stewardship. The performance of DeepONet and FNO is comparable for relatively simple settings, but for complex geometries the performance of FNO deteriorates greatly. We also compare theoretically the two neural operators and obtain similar error estimates for DeepONet and FNO under the same regularity assumptions.

42 ENGINEERING↗
DOE OSTI2024

Why Do Oceanic Nonlinearities Contribute Only Weakly to Extreme El Niño Events?

Abstract Extreme El Niño events have outsized global impacts and control the El Niño Southern Oscillation (ENSO) warm/cold phases asymmetries. Yet, a consensus regarding the relative contributions of atmospheric and oceanic nonlinearities to their genesis remains elusive. Here, we isolate the contribution of oceanic nonlinearities by conducting paired experiments forced with opposite wind stress anomalies in an oceanic general circulation model, which realistically simulates extreme El Niño events and oceanic nonlinearities thought to contribute to ENSO skewness (Tropical Instability Waves (TIWs), Nonlinear Dynamical Heating (NDH)). Our findings indicate a weak contribution of oceanic nonlinearities to extreme El Niño events in the eastern Pacific, owing to compensatory effects between lateral (NDH and TIWs) and vertical processes. These results hold across different vertical mixing schemes and modifications of the upper‐ocean heat budget mixed layer criterion. Our study reinforces previous research underscoring the pivotal role of atmospheric nonlinearities in shaping extreme El Niño events.

Geology↗
DOE OSTI2021

Thermal confinement and transport in spherical tokamaks: a review

Here, we review the thermal plasma confinement and transport properties observed and predicted in low aspect ratio tokamaks, or spherical tokamaks (STs), which can depart significantly from those observed at higher aspect ratio. In particular, thermal energy confinement scalings show a strong, near linear dependence of energy confinement time on toroidal magnetic field, while the dependence on plasma current is more modest, the opposite of what is seen at higher aspect ratio. STs have revealed a very strong improvement in normalized confinement with decreasing collisionality, much stronger than at higher aspect ratio, which bodes well for an ST-based fusion pilot plant should this trend continue at an even lower collisionality than has already been accessed. These differences arise because of fundamental differences in transport in STs due to the more extreme toroidicity (i.e. reduced region of bad curvature), and to the relatively larger $E_r \times B$ shearing rates, both of which can suppress electrostatic drift wave instabilities at both ion and electron gyroradius scales. In addition, electromagnetic effects are much stronger in STs because they operate at high β T . Gyrokinetic (GK) studies, coupled with low- and high-k turbulence measurements, have shed light on the underlying physics controlling transport. At lower β T , both ion- and electron-scale electrostatic drift turbulence may be responsible for transport. At higher β T , microtearing, kinetic ballooning, and hybrid trapped electron/kinetic ballooning modes increasingly play a role, and they have a much stronger impact in the core of ST plasmas than at higher aspect ratio. Flow shear affects the balance between ion- and electron-scale modes. Non-linear GK simulations find regimes where the electron heat flux decreases with decreasing collisionality, consistent with the experimental global normalized confinement scaling. The ST is unique in that the relatively low toroidal magnetic field allows for localized measurements of electron-scale turbulence, and this coupled with turbulence measurements at ion-scales has facilitated detailed comparisons with GK simulations. These data have provided compelling evidence for the presence of ion temperature gradient and electron temperature gradient turbulence in some plasmas, and direct experimental support for the impact of experimental actuators like rotation shear, density gradient and magnetic shear on turbulence and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Stability analysis of Alfvén eigenmodes excited by ion cyclotron resonance heating on EAST

Alfvén wave instabilities driven by energetic particles are common in fusion devices. Understanding their behavior is essential for the good confinement of fusion plasma. For this purpose, a series of experiments are performed on EAST to investigate the excitation of toroidal Alfvén eigenmodes (TAEs). Experimentally, it was found that AEs with frequencies around 80, 134, 157 kHz are excited by ion cyclotron resonance heating (ICRH) hydrogen minority heating scenario. Moreover, the excitation of AEs is independent of the wall-coating materials, but determined with the generation and confinement of the fast ions. Statistical analysis suggests that the TAEs can only be excited when the ICRH power is larger than 2 MW and the H 98 factor is larger than 1.15. In line with the experiments, a set of simulations using TORIC, ASCOT, and NOVA-C is performed. The simulation results show that with the experimental density profile and safety factor, the measured TAE with frequency of 134 kHz is well reproduced by taking into account the plasma toroidal rotation frequency. However, the appearance of the 80 kHz mode cannot be captured in the simulations. The ICRH-generated fast hydrogen ions mainly have banana orbits and that their perpendicular velocity is much larger than that of the parallel ions. The stability analysis with NOVA-C suggests that the radial gradient of the fast H ion distribution is the only driving source of TAEs. For TAEs locate in the center of Alfvén continuum, the D ion Landau damping is the dominant damping term; For most other TAEs which have intersections with the continuum boundary, the continuum damping becomes the dominant damping mechanism.

Alfvén eigenmodes↗
DOE OSTI2022

Explaining the lack of power degradation of energy confinement in wide pedestal quiescent H-modes via transport modeling

Wide pedestal quiescent H (WPQH)-mode is an attractive scenario for future burning plasmas as they operate without ELMs. WPQH is characterized by formation of a wider and higher pedestal (than quiescent H-mode), and broadband fluctuations in the pedestal. Unlike conventional H-modes, where the energy confinement time reduces with increasing heating power, the WPQH plasmas reported in this paper do not show power degradation of the energy confinement. As the injected neutral beam power was increased, reduced core (ρ ≤ 0.45) transport calculated by transp, as well as increased core temperatures, pressure gradient and diamagnetic E × B shear rate were observed. The reduction in the heat transport and rapid decrease in the ion temperature gradient scale length suggest the formation of an ion internal transport barrier (ITB) that was accompanied by increased stored energy in the core. Quasilinear turbulent transport modeling using the trapped gyro Landau fluid (tglf) code was used to predict the ITB and its turbulence stability properties. By using profiles and equilibria produced by matching the transp transport fluxes with the tglf transport model within the tgyro transport solver, the energy confinement time captures the experimentally observed insensitivity to the increased PNBI. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by E × B shear, T i /T e ratio and Shafranov shift; the latter was found to have the strongest effect on the turbulence suppression at the highest heating level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Visualizing the out-of-plane electronic dispersions in an intercalated transition metal dichalcogenide

Layered transition metal dichalcogenides have a rich phase diagram and they feature two-dimensionality in numerous physical properties. Co 1/3 NbS 2 is one of the newest members of this family where Co atoms are intercalated into the van der Waals gaps between NbS 2 layers. Here, we study the three-dimensional electronic band structure of Co 1/3 NbS 2 using both surface and bulk sensitive angle-resolved photoemission spectroscopy. We show that the electronic bands do not fit into the rigid band shift picture after the Co intercalation. Instead, Co 1/3 NbS 2 displays a different orbital character near the Fermi level compared to the pristine NbS 2 compound and has a clear band dispersion in the k z direction despite its layered structure. Our photoemission study demonstrates the out-of-plane electronic correlations introduced by the Co intercalation, thus offering a different perspective on this compound. Finally, we propose how Fermi level tuning could lead to exotic phases such as spin density wave instability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Superconductivity at carrier density 10 17 cm – 3 in quasi-one-dimensional Li 0.9 Mo 6 O 17

Superconductivity is a fascinating phenomenon that involves an attractive interaction by which electrons are paired and can move without electrical resistance below a critical temperature T c . A paradigm for understanding superconductive pairing, the successful “BCS” theory, developed more than 65 years ago, is challenged by certain classes of materials for which parameters (e.g. the value of T c , the density of electrons, etc.) fall outside the range for which BCS theory is applicable. The study of such materials is important because new physics may be revealed and because superconductivity has many potential applications in sensing and electronics. Superconductivity at very low carrier density is a rare and fascinating phenomenon that is of considerable importance and interest in condensed matter physics. The conditions of high density of states and effective screening of Coulomb repulsion, both key ingredients of stable Cooper pairing in the BCS theory, may not be met at low carrier density. Thus, such systems may also harbor unconventional electron pairing (e.g. non-phononic) or pairing without phase coherence as a precursor to Bose-Einstein condensation. An intriguing route toward extreme low-density (and possibly unconventional) superconductivity may be found in quasi-one-dimensional metals for which density-wave instabilities compete for the ground state and can lead to gapping of substantial portions of the Fermi surface (FS), leaving a residual (reconstructed) FS comprising a very small carrier density. Our work highlights a novel material in this unusual category: Li 0.9 Mo 6 O 17 (T c =2 K) (also known as "lithium purple bronze" or "LiPB"). We present both electrical and thermoelectric magnetotransport measurements that reveal superconductivity to occur at carrier density 2 x10 17 cm -3 , among the lowest known for any superconductor, and four orders of magnitude lower than is expected from its chemical valence and band structure. With its exceptionally large ratio of superconducting to Fermi temperature, T c /TF ≳ 0.1, LiPB is properly placed at the top of a very small group of superconductors (including cuprates, heavy-fermion and Fe-based superconductors) for which this ratio falls in the BCS-BEC crossover regime and for which pairing is believed to be spin mediated.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding Dynamics and Thermodynamics of ENSO and Its Complexity Simulated by E3SM and Other Climate Models

Despite the seeming success of most state-of-the-art climate models in simulating the El Niño-Southern Oscillation (ENSO), there is strong evidence that models achieve realistic levels of ENSO activity do so often for wrong reasons. This is owing to an often occurred near cancelation of large errors in terms of contributions to ENSO growth rate from coupled dynamic and thermodynamic feedback processes. Climate models remain deficient in simulating the observed ENSO’s spatial and temporal complexity that involves interplays of coupled dynamic and thermodynamic feedbacks, interactions across multiple scales, nonlinear processes in the tropical atmosphere and ocean system, biases in mean sate and physical processes, and influences external to equatorial Pacific coupled ENSO dynamics. Our proposed research aims at advancing predictive and process-level understandings of ENSO simulated in E3SM and other climate models under current and future climate conditions with two main objectives: (i) better understanding the aforementioned broad range interactive processes and sources that control fundamental properties of ENSO in E3SM and CIMP6 outputs as well as in observational (reanalysis) data sets, using a hierarchical of coupled dynamical frameworks consisting of theoretical analysis, intermediate complexity modeling, and coupled dynamic diagnostics; (ii) to use this understanding to explore pathways towards improving E3SM’s capability of simulating ENSO and its complexity. More specifically, we will focus on four main thrusts of research: (1) ENSO’s dynamic and thermodynamic feedbacks; (2) the across-scale interactions of ENSO with annual cycle and MJO/WWB/TIW (Madden Julian Oscillation/Westerly Wind Burst/Tropical Instability Wave) activity; (3) key nonlinear processes of ENSO involving atmospheric convective thresholds, nonlinear ocean dynamic heating, and thermocline outcropping; and (4) the impacts of climate mean-state biases/changes and perturbed physical processes on simulated ENSO and its complexity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Krash-ESL-Final report

During the period of the current Grant, the effort of the UCSD researchers was focused on: i) The application of the gyrokinetic code COGENT code to study the kinetic drift wave instability and its impact on blob dynamics, and ii) The modeling of edge plasma transport in collisional and semi-collisional regimes, and iii) Impurity transport in edge plasma and the incorporation of plasma-impurity collision operator (hereafter unlike particle collision operator) into COGENT. Over the years four UCSD graduate students (J. Angus, W. Lee, Y. Zhang, and A. Knyazev), funded under this Grant, contributed to the COGENT development/utilization and defended their PhD theses, which are based, partly, on their research related to the COGENT. Three of them (J. Angus, W. Lee, and Y. Zhang) have received “The Best MAE Graduate Student of the Year” awards. Whereas W. Lee after defending his PhD thesis returned back to his home country (South Korea), J. Angus, and Y. Zhang are continuing their research in plasma physics in the LLNL and LANL, whereas A. Knyazev currently is a postdoctoral fellow at Columbia University working on the theoretical aspects of advanced stellarators. The results obtained in our studies were disseminated in many papers published in first-class refereed scientific journals and presented on prestigious scientific meetings including the Annual APS DPP Meeting, Sherwood Fusion Theory Conference, IAEA Fusion Energy Conferences, International Conferences on Plasma-Surface Interactions in Controlled Fusion Devices, International Workshops on Edge Plasma Theory in Fusion Devices, etc. In what follows we highlight the major accomplishments of the UCSD Team in the main areas of our research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Structural Principles in Liquids and Glasses: Bottom-Up or Top-Down

The conventional approach to elucidate the atomic structure of liquid and glass is to start with local structural units made of several atoms, and to use them as building blocks to form a global structure, the bottom-up approach. We propose to add an alternative top-down approach in which we start with a global high-temperature gas state and then apply interatomic potentials to all atoms at once. This causes collective density wave instability in all directions with the same wavelength. These two driving forces, local and global, are in competition and are mutually frustrated. The final structure is determined through the compromise of frustration between these two, which creates the medium-range-order. This even-handed approach on global and local potential energy landscapes explains the distinct natures of short-range order and medium-range order, and strong temperature dependence of various properties of liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2019

The Observability of Vortex-driven Spiral Arms in Protoplanetary Disks: Basic Spiral Properties

Some circumstellar disks are observed to show prominent spiral arms in infrared scattered light or (sub-)millimeter dust continuum. The spirals might be formed from self-gravity, shadows, or planet–disk interactions. Recently, it was hypothesized that massive vortices can drive spiral arms in protoplanetary disks in a way analogous to planets. In this paper, we study the basic properties of vortex-driven spirals by the Rossby Wave Instability in 2D hydrodynamics simulations. We study how the surface density contrast, the number, and the shape of vortex-driven spirals depend on the properties of the vortex. We also compare vortex-driven spirals with those induced by planets. The surface density contrast of vortex-driven spirals in our simulations are comparable to those driven by a sub-thermal mass planet, typically a few to a few tens of Earth masses. In addition, different from the latter, the former is not sensitive to the mass of the vortex. Vortex-driven spiral arms are not expected to be detectable in current scattered light observations, and the prominent spirals observed in scattered light in a few protoplanetary disks, such as SAO 206462 (HD 135344B), MWC 758, and LkHα 330, are unlikely to be induced by the candidate vortices in them.

79 ASTRONOMY AND ASTROPHYSICS↗