Researcher profile

Hsiang-Hsu Wang

Hsiang-Hsu Wang contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2020arXiv

Self-gravitational Force Calculation of Second Order Accuracy Using Multigrid Method on Nested Grids

We present a simple and effective multigrid-based Poisson solver of second-order accuracy in both gravitational potential and forces in terms of the one, two and infinity norms. The method is especially suitable for numerical simulations using nested mesh refinement. The Poisson equation is solved from coarse to fine levels using a one-way interface scheme. We introduce anti-symmetrically linear interpolation for evaluating the boundary conditions across the multigrid hierarchy. The spurious forces commonly observed at the interfaces between refinement levels are effectively suppressed. We validate the method using two- and three-dimensional density-force pairs that are sufficiently smooth for probing the order of accuracy.

preprint2016arXiv

Self-Gravitational Force Calculation of Infinitesimally Thin Gaseous Disks on Nested Grids

We extend the work of Yen et al. (2012) and develop 2nd order formulae to accommodate a nested grid discretization for the direct self-gravitational force calculation for infinitesimally thin gaseous disks. This approach uses a two-dimensional kernel derived for infinitesimally thin disks and is free of artificial boundary conditions. The self-gravitational force calculation is presented in generalized convolution forms for a nested grid configuration. A numerical technique derived from a fast Fourier transform is employed to reduce the computational complexity to be nearly linear. By comparing with analytic potential-density pairs associated with the generalized Maclaurin disks, the extended approach is verified to be of second order accuracy using numerical simulations. The proposed method is accurate, computationally fast and has the potential to be applied to the studies of planetary migration and the gaseous morphology of disk galaxies.

preprint2015arXiv

Self-Gravitational Force Calculation of Second Order Accuracy for Infinitesimally Thin Gaseous Disks in Polar Coordinates

Investigating the evolution of disk galaxies and the dynamics of proto-stellar disks can involve the use of both a hydrodynamical and a Poisson solver. These systems are usually approximated as infinitesimally thin disks using two- dimensional Cartesian or polar coordinates. In Cartesian coordinates, the calcu- lations of the hydrodynamics and self-gravitational forces are relatively straight- forward for attaining second order accuracy. However, in polar coordinates, a second order calculation of self-gravitational forces is required for matching the second order accuracy of hydrodynamical schemes. We present a direct algorithm for calculating self-gravitational forces with second order accuracy without artifi- cial boundary conditions. The Poisson integral in polar coordinates is expressed in a convolution form and the corresponding numerical complexity is nearly lin- ear using a fast Fourier transform. Examples with analytic solutions are used to verify that the truncated error of this algorithm is of second order. The kernel integral around the singularity is applied to modify the particle method. The use of a softening length is avoided and the accuracy of the particle method is significantly improved.

preprint2014arXiv

Formation of Isothermal Disks around Protoplanets. I. Introductory Three-Dimensional Global Simulations for Sub-Neptune-Mass Protoplanets

The regular satellites found around Neptune ($\approx 17~M_{\Earth}$) and Uranus ($\approx 14.5~M_{\Earth}$) suggest that past gaseous circumplanetary disks may have co-existed with solids around sub-Neptune-mass protoplanets ($< 17~M_{\Earth}$). These disks have been shown to be cool, optically thin, quiescent, with low surface density and low viscosity. Numerical studies of the formation are difficult and technically challenging. As an introductory attempt, three-dimensional global simulations are performed to explore the formation of circumplanetary disks around sub-Neptune-mass protoplanets embedded within an isothermal protoplanetary disk at the inviscid limit of the fluid in the absence of self-gravity. Under such conditions, a sub-Neptune-mass protoplanet can reasonably have a rotationally supported circumplanetary disk. The size of the circumplanetary disk is found to be roughly one-tenth of the corresponding Hill radius, which is consistent with the orbital radii of irregular satellites found for Uranus. The protoplanetary gas accretes onto the circumplanetary disk vertically from high altitude and returns to the protoplanetary disk again near the midplane. This implies an open system in which the circumplanetary disk constantly exchanges angular momentum and material with its surrounding prenatal protoplanetary gas.

preprint2014arXiv

Spiral Density Waves in M81. I. Stellar Spiral Density Waves

Aside from the grand-design stellar spirals appearing in the disk of M81, a pair of stellar spiral arms situated well inside the bright bulge of M81 has been recently discovered by Kendall et al. (2008). The seemingly unrelated pairs of spirals pose a challenge to the theory of spiral density waves. To address this problem, we have constructed a three component model for M81, including the contributions from a stellar disk, a bulge, and a dark matter halo subject to observational constraints. Given this basic state for M81, a modal approach is applied to search for the discrete unstable spiral modes that may provide an understanding for the existence of both spiral arms. It is found that the apparently separated inner and outer spirals can be interpreted as a single trailing spiral mode. In particular, these spirals share the same pattern speed 25.5 km s$^{-1}$ kpc$^{-1}$ with a corotation radius of 9.03 kpc. In addition to the good agreement between the calculated and the observed spiral pattern, the variation of the spiral amplitude can also be naturally reproduced.

preprint2014arXiv

Spiral Density Waves in M81. II. Hydrodynamic Simulations for Gas Response to Stellar Spiral Density Waves

Gas response to the underlying stellar spirals is explored for M81 using unmagnetized hydrodynamic simulations. Constrained within the uncertainty of observations, 18 simulations are carried out to study the effects of selfgravity and to cover the parameter space comprising three different sound speeds and three different arm strengths. The results are confronted with those data observed at wavelengths of 8 $μ$m and 21 cm. In the outer disk, the ring-like structure observed in 8 $μ$m image is consistent with the response of cold neutral medium with an effective sound speed 7 km s$^{-1}$, while for the inner disk, the presence of spiral shocks can be understood as a result of 4:1 resonances associated with the warm neutral medium with an effective sound speed 19 km s$^{-1}$. Simulations with single effective sound speed alone cannot simultaneously explain the structures in the outer and inner disks. This justifies the coexistence of cold and warm neutral media in M81. The anomalously high streaming motions observed in the northeast arm and the outward shifted turning points in the iso-velocity contours seen along the southwest arm are interpreted as signatures of interactions with companion galaxies. The level of simulated streaming motions narrows down the uncertainty of observed arm strength toward larger amplitudes.

preprint2013arXiv

Hydrodynamical Simulations of the Barred Spiral Galaxy NGC 1097

NGC 1097 is a nearby barred spiral galaxy believed to be interacting with the elliptical galaxy NGC 1097A located to its northwest. It hosts a Seyfert 1 nucleus surrounded by a circumnuclear starburst ring. Two straight dust lanes connected to the ring extend almost continuously out to the bar. The other ends of the dust lanes attach to two main spiral arms. To provide a physical understanding of its structural and kinematical properties, two-dimensional hydrodynamical simulations have been carried out. Numerical calculations reveal that many features of the gas morphology and kinematics can be reproduced provided that the gas flow is governed by a gravitational potential associated with a slowly rotating strong bar. By including the self-gravity of the gas disk in our calculation, we have found the starburst ring to be gravitationally unstable which is consistent with the observation in \citet{hsieh11}. Our simulations show that the gas inflow rate is 0.17 M$_\sun$ yr$^{-1}$ into the region within the starburst ring even after its formation, leading to the coexistence of both a nuclear ring and a circumnuclear disk.

preprint2010arXiv

Equilibrium Initialization and Stability of Three-Dimensional Gas Disks

We present a new systematic way of setting up galactic gas disks based on the assumption of detailed hydrodynamic equilibrium. To do this, we need to specify the density distribution and the velocity field which supports the disk. We first show that the required circular velocity has no dependence on the height above or below the midplane so long as the gas pressure is a function of density only. The assumption of disks being very thin enables us to decouple the vertical structure from the radial direction. Based on that, the equation of hydrostatic equilibrium together with the reduced Poisson equation leads to two sets of second-order non-linear differential equation, which are easily integrated to set-up a stable disk. We call one approach `density method&#39; and the other one `potential method&#39;. Gas disks in detailed balance are especially suitable for investigating the onset of the gravitational instability. We revisit the question of global, axisymmetric instability using fully three-dimensional disk simulations. The impact of disk thickness on the disk instability and the formation of spontaneously induced spirals is studied systematically with or without the presence of the stellar potential. In our models, the numerical results show that the threshold value for disk instability is shifted from unity to 0.69 for self-gravitating thick disks and to 0.75 for combined stellar and gas thick disks. The simulations also show that self-induced spirals occur in the correct regions and with the right numbers as predicted by the analytic theory.