Researcher profile

Chien-Chang Yen

Chien-Chang Yen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

6 published item(s)

preprint2022arXiv

A Comprehensive Analysis of the Gravitational Wave Events with the Hilbert-Huang Transform: From Compact Binary Coalescence to Supernova

We analyze the gravitational wave signals with a model-independent time-frequency analysis, which is improved from the Hilbert-Huang transform (HHT) and optimized for characterizing the frequency variability on the time-frequency map. Except for the regular HHT algorithm, i.e., obtaining intrinsic mode functions with ensemble empirical mode decomposition and yielding the instantaneous frequencies, we propose an alternative algorithm that operates the ensemble mean on the time-frequency map. We systematically analyze the known gravitational wave events of the compact binary coalescence observed in LIGO O1 and O2, and in the simulated gravitational wave signals from core-collapse supernovae (CCSNe) with our method. The time-frequency maps of the binary black hole coalescence cases show much better details compared to those wavelet spectra. Moreover, the oscillation in the instantaneous frequency caused by mode-mixing could be reduced with our algorithm. For the CCSNe data, the oscillation from the proto-neutron star and the radiation from the standing accretion shock instability can be precisely determined with the HHT in great detail. More importantly, the initial stage of different modes of oscillations can be clearly separated. These results provide new hints for further establishment of the detecting algorithm, and new probes to investigate the underlying physical mechanisms.

preprint2022arXiv

Revealing the drag instability in one-fluid nonideal MHD simulations of a 1D isothermal C-shock

C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in non-ideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the postshock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed.

preprint2020arXiv

Direct Calculation of Self-Gravitational Force for Infinitesimally Thin Gaseous Disks Using Adaptive Mesh Refinement

Yen et al. (2012) advanced a direct approach for the calculation of self-gravitational force to second order accuracy based on uniform grid discretization. This method improves the accuracy of N-body calculation by using exact integration of kernel functions and employing the Fast Fourier Transform (FFT) to reduce complexity of computation to nearly linear. This direct approach is free of artificial boundary conditions, however, the applicability is limited by the uniform discretization of grids. We report here an advancement in the direct method with the implementation of adaptive mesh refinement (AMR) and maintaining second-order accuracy, which breaks the barrier set by uniform grid discretization. The adoption of graphic process units (GPUs) can significantly speed up the computation and make application of this method possible to astrophysical systems of gaseous disk galaxies and protoplanetary disks.

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.

preprint2019arXiv

Modeling the Thermal Bulge of A Hot Jupiter with the Two-Stream Approximation

We revisit the problem of thermal bulge of asynchronous hot Jupiters, using HD 209458 b as a fiducial study. We improve upon previous works by using a double-gray atmosphere model and interior structure from MESA as the background state, and then solve for the thermal bulge in response to the semidiurnal component of stellar insolation. The atmosphere model is based on the radiative transfer with Eddington's two-stream approximation. Two opacity cases are considered: the first introduces a greenhouse effect and the second exhibits a strong temperature inversion. We find that for the predominant thermal bulges excited by g-modes of lower orders, our results are qualitatively similar to the adiabatic results from Arras and Socrates (2010). It arises because the perturbed heating due to self-absorption of thermal emissions can be significant (i.e., greenhouse effect) against Newtonian damping, thereby leading to almost undamped thermal bulges. We also find that the contribution to the thermal bulge from the evanescent waves in the convective zone is not negligible, implying that the thermal bulge is not merely confined in the atmosphere and radiative envelope. Assuming the torque balance between the thermal and gravitational bulges, we estimate the tidal quality factor of the planet for gravitational tides to match the observed radius. The limitations of our model are also briefly discussed.

preprint2015arXiv

Global Transonic Solutions of Planetary Atmospheres in Hydrodynamic Region $\mbox{--}$ Hydrodynamic Escape Problem due to Gravity and Heat

The hydrodynamic escape problem (HEP), which is characterized by a free boundary value problem of Euler equation with gravity and heat, is crucial for investigating the evolution of planetary atmospheres. In this paper, the global existence of transonic solutions to the HEP is established using the generalized Glimm method. The new version of Riemann and boundary-Riemann solvers, are provided as building blocks of the generalized Glimm method by inventing the contraction matrices for the homogeneous Riemann (or boundary-Riemann) solutions. The extended Glimm-Goodman wave interaction estimates are investigated for obtaining a stable scheme and positive gas velocity, which matches the physical observation. The limit of approximation solutions serves as an entropy solution of bounded variations. Moreover, the range of the hydrodynamical region is also obtained.