Search

In this work, we will give proper estimates for the discrete convolution complementary (DCC) kernels, which leads to the asymptotically compatible fractional Gr\"onwall inequality. The consequence can be applied in the analysis of the stability and pointwise-in-time error of difference-type schemes on a non-uniform mesh. The pointwise error is explicitly bound when a non-uniform time grid is given by a specific scale function e.g. graded mesh, can be given directly. Numerical experiments towards the conclusion of this work validate the error analysis., Comment: 22 pages, 4 figures

This paper explores the application of machine learning methods for classifying astronomical sources using photometric data, including normal and emission line galaxies (ELGs; starforming, starburst, AGN, broad line), quasars, and stars. We utilized samples from Sloan Digital Sky Survey (SDSS) Data Release 17 (DR17) and the ALLWISE catalog, which contain spectroscopically labeled sources from SDSS. Our methodology comprises two parts. First, we conducted experiments, including three-class, four-class, and seven-class classifications, employing the Random Forest (RF) algorithm. This phase aimed to achieve optimal performance with balanced datasets. In the second part, we trained various machine learning methods, such as $k$-nearest neighbors (KNN), RF, XGBoost (XGB), voting, and artificial neural network (ANN), using all available data based on promising results from the first phase. Our results highlight the effectiveness of combining optical and infrared features, yielding the best performance across all classifiers. Specifically, in the three-class experiment, RF and XGB algorithms achieved identical average F1 scores of 98.93 per~cent on both balanced and unbalanced datasets. In the seven-class experiment, our average F1 score was 73.57 per~cent. Using the XGB method in the four-class experiment, we achieved F1 scores of 87.9 per~cent for normal galaxies (NGs), 81.5 per~cent for ELGs, 99.1 per~cent for stars, and 98.5 per~cent for quasars (QSOs). Unlike classical methods based on time-consuming spectroscopy, our experiments demonstrate the feasibility of using automated algorithms on carefully classified photometric data. With more data and ample training samples, detailed photometric classification becomes possible, aiding in the selection of follow-up observation candidates., Comment: 13 pages, 10 figures, accepted for publication in MNRAS

Previous numerical simulations have shown that strong winds can be produced in the hot accretion flows around black holes. Most of those studies focus only on the region close to the central black hole, therefore it is unclear whether the wind production stops at large radii around Bondi radius. Bu et al. 2016 studied the hot accretion flow around the Bondi radius in the presence of nuclear star gravity. They find that when the nuclear stars gravity is important/comparable to the black hole gravity, winds can not be produced around the Bondi radius. However, for some galaxies, the nuclear stars gravity around Bondi radius may not be strong. In this case, whether winds can be produced around Bondi radius is not clear. We study the hot accretion flow around Bondi radius with and without thermal conduction by performing hydrodynamical simulations. We use the virtual particles trajectory method to study whether winds exist based on the simulation data. Our numerical results show that in the absence of nuclear stars gravity, winds can be produced around Bondi radius, which causes the mass inflow rate decreasing inwards. We confirm the results of Yuan et al. which indicates this is due to the mass loss of gas via wind rather convectional motions., Comment: 15 pages, 8 figures, accepted for publication in ApJ

We explore the effects of anisotropic thermal conduction, anisotropic pressure, and magnetic field strength on the hot accretion flows around black holes by solving the axisymmetric, steady-state magnetohydrodynamic equations. The anisotropic pressure is known as a mechanism for transporting angular momentum in weakly collisional plasmas in hot accretion flows with extremely low mass accretion rates. However, anisotropic pressure does not extensively impact the transport of the angular momentum, it leads to shrinkage of the wind region. Our results show that the strength of the magnetic field can help the Poynting energy flux overcomes the kinetic energy flux. This result may be applicable to understand the hot accretion flow in the Galactic Center Sgr A* and M87 galaxy., Comment: 13 pages, 6 figures, accepted for publication in ApJ

We investigate the dependency of the inflow-wind structure of the hot accretion flow on the kinematic viscosity coefficient. In this regard, we propose a model for the kinematic viscosity coefficient to mimic the behavior of the magnetorotational instability and would be maximal at the rotation axis. Then, we compare our model with two other prescriptions from numerical simulations of the accretion flow. We solve two-dimensional hydrodynamic equations of hot accretion flows in the presence of the thermal conduction. The self-similar approach is also adopted in the radial direction. We calculate the properties of the inflow and the wind such as velocity, density, angular momentum for three models of the kinematic viscosity prescription. On inspection, we find that in the model we suggested wind is less efficient than that in two other models to extract the angular momentum outward where the self-similar solutions are applied. The solutions obtained in this paper might be applicable to the hydrodynamical numerical simulations of the hot accretion flow., Comment: 13 pages, 5 figures, accepted for publication in ApJ

We present a residual-based a posteriori error estimator for the hybrid high-order (HHO) method for the Stokes model problem. Both the proposed HHO method and error estimator are valid in two and three dimensions and support arbitrary approximation orders on fairly general meshes. The upper bound and lower bound of the error estimator are proved, in which proof, the key ingredient is a novel stabilizer employed in the discrete scheme. By using the given estimator, adaptive algorithm of HHO method is designed to solve model problem. Finally, the expected theoretical results are numerically demonstrated on a variety of meshes for model problem.

We solve the two-dimensional hydrodynamic equations of hot accretion flow in the presence of the thermal conduction. The flow is assumed to be in steady-state and axisymmetric, and self-similar approximation is adopted in the radial direction. In this hydrodynamic study, we consider the viscous stress tensor to mimic the effects of the magnetorotational instability for driving angular momentum. We impose the physical boundary conditions at both the rotation axis and the equatorial plane and obtain the solutions in the full $ r-\theta $ space. We have found that thermal conduction is indispensable term for investigating the inflow-wind structure of the hot accretion flows with very low mass accretion rates. One of the most interesting results here is that the disc is convectively stable in hot accretion mode and in the presence of the thermal conduction. Furthermore, the properties of wind and also its driving mechanisms are studied. Our analytical results are consistent with previous numerical simulations of hot accretion flow., Comment: 20 pages, 10 figures, accepted for publication in ApJ

We present the two-dimensional inflow-outflow solutions of radiation hydrodynamic (RHD) equations of supercritical accretion flows. Compared with prior studies, we include all components of the viscous stress tensor. We assume steady state flow and use self-similar solutions in the radial direction to solve the equations in $ r-\theta $ domain of the spherical coordinates. The set of differential equations have been integrated from the rotation axis to the equatorial plane. We find that the self-similarity assumption requires that the radial profile of density is described by $ \rho(r) \propto r^{-0.5} $. Correspondingly, the radial profile of the mass inflow rate decreases with decreasing radii as $ \dot{M}_\mathrm{in} \propto r $. Inflow-outflow structure has been found in our solution. In the region $ \theta > 65^{\circ} $ there exist inflow while above that flow moves outward and outflow could launch. The driving forces of the outflow are analyzed and found that the radiation force is dominant and push the gas particles outwards with poloidal velocity $ \sim 0.25 c $. The properties of outflow are also studied. The results show that the mass flux weighted angular momentum of the inflow is lower than that of outflow, thus the angular momentum of the flow can be transported by the outflow. We also analyze the convective stability of the supercritical disk and find that in the absence of the magnetic field, the flow is convectively unstable. Our analytical results are fully consistent with the previous numerical simulations of the supercritical accretion flow., Comment: 14 pages, 5 figures, accepted for publication in ApJ

In this paper A Geometrical Design Method for Blade's surface $\Im$ in the impeller is provided here $\Im$ is a solution to a coupling system consisting of the well-known Navier-Stokes equations and a four order elliptic boundary value problem . The coupling system is used to describe the relations between solutions of Navier-Stokes equations and the geometry of the domain occupied by fluids, and also provides new theory and methods for optimal geometric design of the boundary of domain mentioned above. This coupling system is the Eular-Lagrange equations of the optimal control problem which is describing a new principle of the geometric design for the blade's surface of an impeller. The control variable is the surface of the blade and the state equations are Navier-Stokes equations with mixed boundary conditions in the channel between two blades. The objective functional depending on the geometry shape of blade's surface describes the dissipation energy of the flow and the power of the impeller. First we prove the existence of a solution of the optimal control problem. Then we use a special coordinate system of the Navier-Stokes equations to derive the objective functional which depends on the surface $\Theta$ explicitly. We also show the weakly continuity of the solution of the Navier-Stokes equations with respect to the geometry shape of the blade's surface., Comment: 27 pages