Pith. sign in

REVIEW 4 major objections 5 minor 83 references

Novel methodology to obtain transonic solutions for dissipative flows around compact objects

T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper claims that an implicit-explicit (ImEx) integration scheme overcomes the stiffness in the angular-momentum equation and regenerates global transonic accretion and wind solutions around compact objects.

desk verdict A useful ImEx fix for a known viscous-transonic integration bottleneck, but the printed equations have a factor-of-two inconsistency and the method lacks a convergence test. read the letter →

arxiv 2505.24839 v1 pith:WD5ABKRO submitted 2025-05-30 astro-ph.HE hep-thphysics.comp-phphysics.plasm-ph

Shilpa Sarkar This is my paper
classification astro-ph.HEhep-thphysics.comp-phphysics.plasm-ph
keywords transonicaccretionwindsolutionsimplicit-explicit(ImEx)schemesviscousflowsstiffdifferentialequationsblackholemultiplesonicpointsshocks
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper claims to solve a long-standing numerical obstruction in astrophysical fluid dynamics: with a realistic viscous stress proportional to $d\Omega/dr$, the angular-momentum equation becomes stiff, and explicit integration from the sonic point fails in the inward direction. The author proposes a partitioned implicit-explicit (ImEx) scheme—implicit trapezoidal update for the angular momentum $\lambda$, explicit adaptive Runge-Kutta for velocity and temperature—and reports that this scheme successfully regenerates the supersonic accretion branch and, for the first time, obtains global subsonic wind branches. If correct, the method yields global transonic solutions of every topology, including multiple sonic points and shocks, for both accretion and winds, and would unify the treatment of inflow and outflow around black holes. The approach is demonstrated across a broad parameter space with a variable-adiabatic-index equation of state and realistic cooling.

What carries the argument

The load-bearing object is the differential equation for specific angular momentum, $d\lambda/dr = 2r\Omega + r^2\,d\Omega/dr$, with $d\Omega/dr$ given by Eq. (13) and containing the stiff gradient term proportional to $\Omega_K(\lambda-\lambda_0)/(\alpha a_s^2 r^2)$. This equation is the only part treated implicitly. The ImEx partition pairs an implicit trapezoidal step (Butcher tableau) for $\lambda$ with an explicit Heun-Euler adaptive step for $v$ and $\Theta$, chosen because the two tableaus share the same nodes. The sonic-point location uses the extended IRM-SP-RV relaxation method, and at the critical point the two possible slopes of the transonic solution (accretion and wind) come from L'Hôpital's rule on the $0/0$ form of $dv/dr$.

What would settle it

Directly rederive Eq. (13) from Eqs. (6), (9), and (12) symbolically: if the result differs by a factor of 1/2 from the printed form, the stiff equation being solved is not the model's equation. Alternatively, for a parameter set where the flow is not stiff, compare the ImEx output with a converged high-accuracy explicit or fully implicit reference solution; a discrepancy beyond the integration tolerance would falsify the claim that the scheme reproduces the true global solutions.

Watch

Extended reading notes

Core claim

The central claim is that the integration failure reported in earlier work on viscous transonic flows is purely a stiffness problem in the coupled ordinary differential equations, not an actual absence of physically admissible solutions. The paper asserts that treating only the $d\lambda/dr$ equation implicitly—using an implicit trapezoidal method—while integrating $dv/dr$ and $d\Theta/dr$ explicitly with an adaptive RK method, regenerates the branches that explicit methods lose. With the sonic point located by the extended IRM-SP-RV procedure, the same scheme produces all topologies of global transonic solutions, the multiple-sonic-point regime, and shocks, for accretion and winds alike.

Load-bearing premise

The load-bearing premise is that Eq. (13) is the exact azimuthal momentum balance, although a direct derivation from the stated model appears to introduce a factor of 1/2 that the paper does not discuss.

Editorial extensions

If this is right

  • If correct, global wind solutions become obtainable for dissipative flows with realistic viscosity, removing the previous restriction to accretion-only solutions.
  • The same method can produce multiple-sonic-point and shock topologies for both accretion and winds, extending the older single-sonic-point treatments.
  • A unified parameter-space study of inflow and outflow becomes feasible, including bremsstrahlung, synchrotron, and inverse-Compton cooling.
  • The long-reported 'inward integration fails' issue is resolved, so solutions no longer depend on integrating only in the outward direction.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If validated, the partitioned ImEx strategy could generalize to other astrophysical stiff systems, such as two-temperature plasmas or flows with thermal conduction, where similar gradient-coupling stiff terms appear.
  • The paper leaves a formal stability and order analysis of the partitioned ImEx scheme implicit; an error study comparing against a fully implicit or high-order reference solver would be the natural next test.
  • The derivation of Eq. (13) from Eqs. (6), (9), and (12) appears to yield an extra factor of 1/2; if this check stands, the reported solutions may correspond to a slightly different viscosity normalization than the one stated.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper reports a numerical methodology for obtaining global transonic accretion and wind solutions around compact objects when the viscous stress is proportional to dΩ/dr. The author combines the IRM-SP iterative relaxation method (extended to this viscous prescription as IRM-SP-RV) with an implicit-explicit (ImEx) integration scheme in which the angular momentum equation is treated implicitly and the velocity and temperature equations explicitly. The central claim is that this ImEx approach cures the failure of inward integration from the sonic point noted by Becker & Le (2003), thereby allowing both accretion and wind solutions, including multiple sonic points and shocks, to be obtained for the first time for this viscosity prescription. The paper includes a broad parameter-space survey with realistic cooling processes.

Significance. If the claims hold, the paper addresses a long-standing technical barrier in transonic accretion/wind modeling and provides a unified framework for both inflows and outflows under a realistic viscous stress. The author deserves credit for explicitly identifying the BL03 failure, extending the IRM-SP scheme to the dΩ/dr viscosity, and including detailed cooling processes. The sonic-point location procedure is a legitimate iterative shooting method, not an empirical fit. However, the central numerical claim is not yet established: the printed equations contain a factor-of-two inconsistency, and the ImEx scheme lacks convergence, stiffness, and stability analyses. The paper's strength is its qualitative demonstration over a wide parameter space; its weakness is the absence of quantitative numerical verification and the unresolved consistency of the governing ODE.

major comments (4)
  1. [Section II.B.2, Eqs. (9) and (13)] As printed, Eq. (9) is not the integral of Eq. (8) with Eq. (6): because Σ v r = Mdot/π from Eq. (6), the integrated azimuthal momentum equation is Mdot(λ − λ0) = −π r² t_rφ, not the printed Mdot(λ − λ0) = −2π r² t_rφ. Substituting the realistic stress, Eq. (12), into the printed Eq. (9) gives dΩ/dr = −Γ v Ω_K (λ − λ0) / (2 α a_s² r²), which differs from Eq. (13) by a factor of 2. The printed Eq. (13) is actually the correct result if Eq. (9) is corrected to have π. Since Eq. (13) is the central ODE whose inward integration is the entire subject of the paper, this inconsistency is load-bearing: the reported solutions may be computed from a different viscous torque than the stated model. Please correct Eq. (9) or Eq. (13), and re-check the consistency of Eq. (11) and the generalized Bernoulli constant, Eq. (19), under the corrected definition.
  2. [Section III.B.2] No convergence test is presented for the ImEx scheme. The claim that the scheme "successfully as well as efficiently regenerates" the transonic solutions is supported only by visual overlap in Figs. 2 and 5. The manuscript should report a step-size or tolerance study showing that the inward branches converge to the outward branches, or to a high-accuracy reference, with decreasing step size, and should present error norms for representative accretion and wind solutions.
  3. [Section III.B.2] The identification of Eq. (14) as "stiff" is asserted, not demonstrated. The failure of explicit methods shown in Fig. 2 could also arise from the critical-point singularity or from accuracy loss near the sonic point. A stiffness ratio, a controlled step-size study, or a demonstration that the explicit method's error grows as h→0 is needed to support the central numerical diagnosis.
  4. [Section III.B.2 and Table I] The partitioned ImEx update is not analyzed. The implicit trapezoidal rule is applied to the λ equation while velocity and temperature are advanced explicitly, but those equations depend on λ; nothing establishes the stability, order, or error balance of the coupled composite scheme. In particular, the Butcher tableaus define an explicit and an implicit method of different orders, and the composite integrator's effective order is never stated. Please provide a local truncation-error or linear-stability analysis, or at least an empirical order verification.
minor comments (5)
  1. [Section IV.B.2 and Fig. 7] The flow parameters listed in the text for Fig. 7 (E = 1.001, λ0 = 1.5, α = 0.01) do not match the figure caption, which states λ0 = 1.6 and α = 0.03. Please make the text and figure consistent.
  2. [Section III.A, step 2] The expression "E≃E" in the description of the inner boundary condition is unclear; it should be written as E ≃ Ẽ or another notation indicating the algebraic Bernoulli constant without dissipation.
  3. [Fig. 9 caption] The caption "variation of λ0 with E" is ambiguous; the figure sweeps both λ0 and E, so the caption should state "variation with λ0 and E".
  4. [Table I] The Butcher tableau for the implicit method is written in an unusual layout; in particular, the lower-order b* row is difficult to distinguish from the a-coefficients. Adding explicit row labels (c, a, b, b*) would improve clarity.
  5. [General] The paper does not include a data or code availability statement; please add one, as reproducibility would strengthen the numerical claims.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the ImEx scheme is validated by independent directional consistency checks, not by fitting to the quantity it predicts.

full rationale

The central claim is that an ImEx integration scheme cures the failure of inward integration for flows with t_rphi proportional to dOmega/dr. This claim is checked against the paper's own equations and against independently obtained branches: the ImEx solutions retrace the supersonic accretion branch obtained separately by outward integration from the inner boundary (Figs. 1c, 2, 5), and the accretion/wind branches are also compared with the traditional explicit CK integration where that integration succeeds. This is a self-contained numerical consistency test, not a fit of a parameter to the quantity being predicted. The IRM-SP and IRM-SHOCK schemes are cited from Paper 1 and reused, but they are auxiliary locating procedures whose output is a sonic-point boundary condition, not the claimed result; the claimed result is the ability to integrate from that boundary, and that is demonstrated by the direction-independence of the solutions. The paper contains self-citations, but none is load-bearing in the sense of importing a forbidden uniqueness theorem or hiding an ansatz. The apparent factor-of-two discrepancy between Eq. (9) and the integration of Eq. (8) is a possible correctness issue, not a circularity, because the equations are stated rather than being defined in terms of the results they are used to produce.

Assumptions & free parameters 0 free parameters · 9 assumptions · 0 invented entities

The central claim rests on the Becker-Le viscous flow model and the CR equation of state rather than on new physics. No empirical data are fitted; E, λ0, α, β, Mdot, MBH, and ξ are physical inputs. The most paper-specific assumptions are the inner-boundary simplification and the unproved partitioned ImEx stability, which are the main fragility.

assumptions (9)
  • domain assumption Paczyński-Wiita pseudo-potential Φ = -1/(r-1) approximates Schwarzschild gravity for the transonic region.
    Introduced in Section II to avoid full general relativity while retaining strong-gravity effects; all solutions inherit this approximation.
  • domain assumption Steady, axisymmetric, height-integrated thin-disk flow with vertical hydrostatic equilibrium, giving Eq. (7) for H.
    Section II.B; the entire one-dimensional ODE system is built on this geometry.
  • domain assumption Single-temperature, charge-neutral electron-proton fluid closed by the variable adiabatic index CR equation of state, Eq. (1).
    Section II.A; thermodynamics of all solutions depends on this closure.
  • domain assumption Viscous stress vanishes at the horizon, giving the integration constant λ0 and Eq. (9).
    Section II.B, taken from Becker and Le 2003; it determines the angular momentum equation and is not derived in this paper.
  • domain assumption Inner boundary at r_in = 1.001 is radiatively inefficient so E reduces to its inviscid form and λ_in ≈ λ0.
    Section III.A step 2; this is the boundary condition on which IRM-SP-RV is built.
  • domain assumption Synchrotron, bremsstrahlung, Comptonization, and Compton heating terms from Sarkar et al. 2020 and 2022 are correct for stochastic fields specified by β.
    Section II.B; cooling and heating enter Eq. (15) and affect the thermodynamics.
  • domain assumption Maximum-entropy principle selects the physical solution when multiple sonic points exist.
    Section II.C; used to choose among solutions and to justify shocks.
  • standard math Rankine-Hugoniot jump conditions, mass, energy, angular momentum, and radial momentum flux, apply to the height-integrated one-dimensional flow.
    Section II.D; used to locate shocks.
  • ad hoc to paper The partitioned ImEx update remains stable and accurate when the implicit stage uses only the λ equation and explicit stages provide v and Θ.
    Section III.B.2; no stability or order analysis is given for the coupled partitioned scheme.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Novel methodology to obtain transonic solutions for dissipative flows around compact objects." pith.science (2026). https://pith.science/paper/WD5ABKRO

@misc{pith2026250524839,
  author       = {Pith},
  title        = {Pith review of: Novel methodology to obtain transonic solutions for dissipative flows around compact objects},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WD5ABKRO}},
  note         = {Machine review of arXiv:2505.24839}
}
read the original abstract

A novel methodology to obtain global transonic solutions around compact objects is reported here. A unified methodology to obtain accretion as well as wind solutions around these objects has been presented. Flows around compact objects are dissipative, and the conservation equations are therefore stiff. In such conditions, obtaining of sonic point(s) and hence, the transonic solution is not trivial. The conserved equations of motion fail to integrate in the presence of realistic viscosity, thereby making it difficult to obtain a global solution. This inhibits one from getting an actual picture of an astrophysical flow. The current work addresses this long-standing issue of obtaining solutions for both accretion and wind. The methodology developed utilises the inner boundary conditions and takes recourse to implicit-explicit (ImEx) integration schemes, to obtain general global transonic accretion and wind solutions. This is the first time such an attempt has been made. Current work considers the different cooling processes like bremsstrahlung, synchrotron and their inverse-Comptonizations, which are found to affect the thermodynamics of the flow. This methodology could successfully generate all topologies of global solutions, multiple sonic point regime, as well as shocks. A broad parameter space study has been done in this work. In an upcoming part II of the paper, a detailed discussion on the spectra and luminosity of the accretion and wind solutions has been presented.

Figures

Figures reproduced from arXiv: 2505.24839 by the authors.

Figure 1
Figure 1. FIG. 1: Finding of SPs (left panel) and global transonic [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Failure of traditional explicit integration [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Typical global transonic solution for accretion (solid curves) and wind (dashed curves). The direction of [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Finding multiple sonic points using [PITH_FULL_IMAGE:figures/full_fig_p011_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: Typical global transonic solutions for accretion [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6: Flow variables in accreting systems harbouring shock. Couloured star represents [PITH_FULL_IMAGE:figures/full_fig_p012_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7: Typical global transonic solutions for accretion [PITH_FULL_IMAGE:figures/full_fig_p013_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8: Flow variables associated with the shocked [PITH_FULL_IMAGE:figures/full_fig_p013_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9: Typical global transonic solutions of accretion (solid curves) and wind (dashed curves) for variation of [PITH_FULL_IMAGE:figures/full_fig_p014_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10: Typical global transonic solution for accretion (solid curves) and wind (dotted curves). The direction of [PITH_FULL_IMAGE:figures/full_fig_p015_10.png]

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

83 extracted references · 72 canonical work pages

  1. [1]

    This is famously used by many authors, because it reduces Eqs

    First is the general,t rϕ ∝αpprescription, whereαis the Shakura and Sunyaev viscos- ity parameter [2]. This is famously used by many authors, because it reduces Eqs. 8–9 to an algebraic form [19, 54, Paper 1]. The exact expression is given below: trϕ,αp =−αW,(10) where,Wis the height-integrated gas pressure W= 2Hp. Then, Eq. 9 reduces to: λ=λ 0 + 2πr2αW ˙...

  2. [2]

    The second prescription is more realistic and involves the gradient of angular velocity [25, 27, 34, 39]. It is given as, trϕ,Real =ηr dΩ dr ,(12) where,Ωis the angular velocity,ηis the dynamical viscosity parameter and is given byη= Σν= 2ρH[αa 2 s /(ΓΩK)], whereν is the kinematic viscosity parameter, Ω K = 1/[2r(r−1) 2] is the Keplerian angular velocity ...

    1984

  3. [3]

    LOCATE SP(s): Revisit the IRM-SP scheme, where a short explanation on the method to ob- tain the SP in the presence of realistic viscosity prescription is presented

  4. [4]

    This part is subdivided into two further subsections

    OBTAIN TSs: Discuss the methodology to obtain the global TSs for accretion and winds. This part is subdivided into two further subsections. •One section discusses the issues associated with using the traditional methods of integra- tion, to obtain global TS. •Second part discusses the scheme which can remove this issue and allow one to obtain global TSs f...

    work page

  5. [5]

    Here,βis the inverse of plasma beta parameter and is used to find the magnitude of stochastic magnetic field required for computing the synchrotron emission

    FLOW PARAMETERS: SupplyE,λ 0,α,β, ˙M, MBH andξ. Here,βis the inverse of plasma beta parameter and is used to find the magnitude of stochastic magnetic field required for computing the synchrotron emission. This step defines or char- acterises the system

    work page

  6. [6]

    BOUNDARY V ALUES, finding appropriate values ofr in,v in, Θ in andλ in: Selectr in →r g = 1.001 to be the boundary and focus on the accretion branch. The reason for selecting this inner bound- ary is that asymptotically close to the central ob- ject, gravity is strong, and infall timescales are much shorter than any other timescales. Hence, matter does no...

    work page

  7. [7]

    14),dv/dr(Eq

    IRM-SP-R V TECHNIQUE: Usingv in, Θ in,λ in (from step 2) and the flow parameters supplied (in step 1) atr in, as boundary conditions (BCs), an integration technique is employed (preferably an explicit adaptive RK method) to solve the set of coupled differential equationsi.e.,the EoM:dλ/dr (Eq. 14),dv/dr(Eq. 17) anddΘ/dr(Eq. 15). It is to note that, in the...

    work page

  8. [8]

    For the present case, in Fig

    CHECKING OF MSP: To check whether there are any other SPs present (apart from the one obtained above) for the same set of flow parameters, the Θ in used to obtain SP in step 3 is changed by a large factor, and the same IRM-SP-R V technique is fol- lowed. For the present case, in Fig. 1, only one SP is present. However, if the parameters are changed, then ...

    work page

Show all 83 references
  1. [9]

    1 (a1, b1, c1) are overplotted on the global TS plotted in panels a2, b2 and c2, respectively

    The issue of integration The thick-dashed black curves in Fig. 1 (a1, b1, c1) are overplotted on the global TS plotted in panels a2, b2 and c2, respectively. It is seen that the supersonic portion of the accretion solutions, in panels a2 and b2 (represented by thin, solid, col...

  2. [10]

    Implicit-Explicit (ImEx) scheme Since thedλ/drequation is stiff, an implicit scheme is experimented upon and included in the picture of in- tegration methodology [48, 51, 61], unlike the sole use of traditional explicit schemes. It is seen and reported in this work, that an Im...

  3. [11]

    The solution passing throughr ci is non-global (NG), andr co is global (G)

    Accretion flows with shocks As have been mentioned before in detail, accretion flows passing through bothr ci andr co are connected to the inner boundary. The solution passing throughr ci is non-global (NG), andr co is global (G). Since both these G and NG solutions are connec...

  4. [12]

    While the solution passing throughr ci is always connected to the inner boundary the solution passing throughr co is not necessarily so

    Wind flows with shocks Unlike accretion flows, MSP for winds are difficult to obtain. While the solution passing throughr ci is always connected to the inner boundary the solution passing throughr co is not necessarily so. In the section before, it was seen that if the wind so...

  5. [13]

    Otherwise, the wind solution would pass through eitherr ci orr co whichever is global and has a higher entropy

    If a wind flow harbours a shock, then only it is pos- sible to locater co and the global TS would pass through both the SPs via a shock transition. Otherwise, the wind solution would pass through eitherr ci orr co whichever is global and has a higher entropy. The IRM-SHOCK tec...

  6. [14]

    Use IRM-SP-R V to findr ci for a given set of flow parameters:E,λ 0,α,β, ˙M,M BH andξ

  7. [15]

    Obtain the global transonic wind solution passing throughr ci using ImEx scheme

  8. [16]

    Using these values, we assume a shock jump at that location, and the EoM are inte- grated further outwards

    At each point of the supersonic branch of wind solution, the shock conditions are utilised to find post-shock values. Using these values, we assume a shock jump at that location, and the EoM are inte- grated further outwards. The post-shock solution thus obtained need not be transonic

  9. [17]

    In this wayr co is located

    Step 3 is iteratedi.e.,relaxation method is ap- plied by iterating onruntil the post-shock solution passes through a SP. In this wayr co is located. It could be possible that the post-shock solution does not pass through an SP, even after multiple iterations ofr. This would su...

    2003

  10. [18]

    Bondi, On spherically symmetrical accretion, Monthly Notices of the Royal Astronomical Society112, 195 (1952)

    H. Bondi, On spherically symmetrical accretion, Monthly Notices of the Royal Astronomical Society112, 195 (1952)

    1952

  11. [19]

    N. I. Shakura and R. A. Sunyaev, Reprint of 1973A&A....24..337S. Black holes in binary systems. Ob- servational appearance., A&A500, 33 (1973)

    1973

  12. [20]

    J. E. Pringle and M. J. Rees, Accretion Disc Models for Compact X-Ray Sources, A&A21, 1 (1972)

    1972

  13. [21]

    A. V. Koldoba, R. V. E. Lovelace, G. V. Ustyugova, and M. M. Romanova, Funnel Flows from Disks to Magnetized Stars, AJ123, 2019 (2002), arXiv:astro- ph/0104058 [astro-ph]

    2002

  14. [22]

    S. K. Chakrabarti, Grand Unification of Solutions of Accretion and Winds around Black Holes and Neu- tron Stars, Astrophys. J.464, 664 (1996), arXiv:astro- ph/9606145 [astro-ph]

    1996

  15. [23]

    Mukhopadhyay, Hydrodynamical study of advective accretion flow around neutron stars, Int

    B. Mukhopadhyay, Hydrodynamical study of advective accretion flow around neutron stars, Int. J. Mod. Phys. D11, 1305 (2002), arXiv:astro-ph/0203438

    2002 arXiv

  16. [24]

    D.-F. Bu, E. Qiao, and X.-H. Yang, Hot Accretion Flow around Neutron Stars, Astrophys. J.875, 147 (2019), arXiv:1903.10708 [astro-ph.HE]

    2019 arXiv

  17. [25]

    R. D. Blandford and R. L. Znajek, Electromagnetic ex- traction of energy from Kerr black holes., Monthly No- tices of the Royal Astronomical Society179, 433 (1977)

    1977

  18. [26]

    R. D. Blandford and D. G. Payne, Hydromagnetic flows from accretion disks and the production of radio jets., Monthly Notices of the Royal Astronomical Society199, 883 (1982)

    1982

  19. [27]

    Meliani, C

    Z. Meliani, C. Sauty, K. Tsinganos, and N. Vlahakis, Rel- ativistic Parker winds with variable effective polytropic index, A&A425, 773 (2004), arXiv:astro-ph/0407100 [astro-ph]

    2004 arXiv

  20. [28]

    S. K. Chakrabarti, Accretion processes on a black hole., PhR266, 229 (1996), arXiv:astro-ph/9605015 [astro-ph]

    1996 arXiv

  21. [29]

    E. N. Parker, Dynamics of the Interplanetary Gas and Magnetic Fields., Astrophys. J.128, 664 (1958)

    1958

  22. [30]

    E. J. Weber and J. Davis, Leverett, The Angular Momen- tum of the Solar Wind, Astrophys. J.148, 217 (1967)

    1967

  23. [31]

    T. E. Holzer and W. I. Axford, The Theory of Stellar Winds and Related Flows, Annual Review of Astronomy and Astrophysics8, 31 (1970)

    1970

  24. [32]

    G. R. Blumenthal and W. G. Mathews, Spherical winds and accretion in general relativity., Astrophys. J.203, 714 (1976)

    1976

  25. [33]

    E. P. T. Liang and K. A. Thompson, Transonic disk ac- cretion onto black holes., Astrophys. J.240, 271 (1980)

    1980

  26. [34]

    R. A. Flammang, Stationary spherical accretion into black holes II. Theory of optically thick accretion., Monthly Notices of the Royal Astronomical Society199, 833 (1982)

    1982

  27. [35]

    Fukue, Transonic disk accretion revisited, PASJ39, 309 (1987)

    J. Fukue, Transonic disk accretion revisited, PASJ39, 309 (1987)

    1987

  28. [36]

    Sarkar and I

    S. Sarkar and I. Chattopadhyay, Viscous dissipative two- temperature accretion flows around black holes, Journal of Astrophysics and Astronomy43, 34 (2022)

    2022

  29. [37]

    Sarkar and I

    S. Sarkar and I. M. Kulikov, Iterative Relaxation Method to Obtain Global Transonic Flows around Compact Objects (accepted in LOBACHEVSKII JOURNAL OF MATHEMATICS), arXiv e-prints , arXiv:2505.07397 (2025), arXiv:2505.07397 [astro-ph.HE]

    2025 arXiv

  30. [38]

    S. K. Chakrabarti,Theory of Transonic Astrophysical Flows(1990)

    1990

  31. [39]

    S. K. Chakrabarti and D. Molteni, Viscosity prescrip- tions in accretion discs with shock waves, Monthly Notices of the Royal Astronomical Society272, 80 (1995), https://academic.oup.com/mnras/article- pdf/272/1/80/2897406/mnras272-0080.pdf

    1995

  32. [40]

    Matsumoto, S

    R. Matsumoto, S. Kato, J. Fukue, and A. T. Okazaki, Viscous transonic flow around the inner edge of geomet- rically thin accretion disks, PASJ36, 71 (1984)

    1984

  33. [41]

    K. E. Nakamura, R. Matsumoto, M. Kusunose, and S. Kato, Global Structures of Advection-Dominated Two-Temperature Accretion Disks, PASJ48, 761 (1996)

    1996

  34. [42]

    Debnath, I

    S. Debnath, I. Chattopadhyay, and R. Joshi, Oscillating shocks in the transonic viscous, variable gamma accretion flows around black holes, Monthly Notices of the Royal Astronomical Society528(2024)

    2024

  35. [43]

    Papaloizou and E

    J. Papaloizou and E. Szuszkiewicz, A Comparison of One-Dimensional and Two-Dimensional Models of Tran- sonic Accretion Discs around Collapsed Objects, Monthly Notices of the Royal Astronomical Society268, 29 (1994), arXiv:astro-ph/9401025 [astro-ph]

    1994 arXiv

  36. [44]

    Takahashi, Advection-dominated accretion flows with causal viscosity, A&A461, 393 (2007), arXiv:0704.3554 [astro-ph]

    R. Takahashi, Advection-dominated accretion flows with causal viscosity, A&A461, 393 (2007), arXiv:0704.3554 [astro-ph]

    2007 arXiv

  37. [45]

    Gu and J.-F

    W.-M. Gu and J.-F. Lu, Causal viscosity in global dy- namics of advection-dominated accretion flows around black holes, Chinese Physics Letters17, 151 (2000)

    2000

  38. [46]

    C. F. Gammie and R. Popham, Advection-dominated Ac- cretion Flows in the Kerr Metric. I. Basic Equations, Astrophys. J.498, 313 (1998), arXiv:astro-ph/9705117 [astro-ph]

    1998 arXiv

  39. [47]

    Popham and R

    R. Popham and R. Narayan, Supersonic Infall and Causality in Accretion Disk Boundary Layers, Astrophys. J.394, 255 (1992)

    1992

  40. [48]

    Lu, W.-M

    J.-F. Lu, W.-M. Gu, and F. Yuan, Global dynamics of advection-dominated accretion revisited, The Astrophys- ical Journal523, 340 (1999)

    1999

  41. [49]

    Chen and R

    X. Chen and R. E. Taam, The Structure and Stability of Transonic Accretion Disks Surrounding Black Holes, Astrophys. J.412, 254 (1993)

    1993

  42. [50]

    Narayan, S

    R. Narayan, S. Kato, and F. Honma, Global Struc- ture and Dynamics of Advection-dominated Accretion Flows around Black Holes, Astrophys. J.476, 49 (1997), arXiv:astro-ph/9607019 [astro-ph]

    1997 arXiv

  43. [51]

    P. A. Becker and T. B. Le, Inner Boundary Conditions for Advection-dominated Accretion onto Black Holes, Astro- phys. J.588, 408 (2003), arXiv:astro-ph/0301151 [astro- ph]

    2003 arXiv

  44. [52]

    Palit, A

    I. Palit, A. Janiuk, and P. Sukova, Effects of adia- batic index on the sonic surface and time variability of low angular momentum accretion flows, Monthly No- tices of the Royal Astronomical Society487, 755 (2019), arXiv:1905.02289 [astro-ph.HE]

    2019 arXiv

  45. [53]

    S. K. Chakrabarti, Spectral properties of accretion disks around black holes. ii. sub-keplerian flows with and with- out shocks, The Astrophysical Journal484, 313 (1997)

    1997

  46. [54]

    Chattopadhyay and D

    I. Chattopadhyay and D. Ryu, Effects of Fluid Composi- tion on Spherical Flows Around Black Holes, Astrophys. J.694, 492 (2009), arXiv:0812.2607 [astro-ph]

    2009 arXiv

  47. [55]

    Chandrasekhar,An Introduction to the Study of Stel- 18 lar Structure(Dover Publications, New York, 1957) orig- inally published by University of Chicago Press, 1939

    S. Chandrasekhar,An Introduction to the Study of Stel- 18 lar Structure(Dover Publications, New York, 1957) orig- inally published by University of Chicago Press, 1939

    1957

  48. [56]

    Kumar and I

    R. Kumar and I. Chattopadhyay, Dissipative advective accretion disc solutions with variable adiabatic index around black holes, Monthly Notices of the Royal As- tronomical Society443, 3444 (2014), arXiv:1407.2130 [astro-ph.HE]

    2014 arXiv

  49. [57]

    Sarkar, I

    S. Sarkar, I. Chattopadhyay, and P. Laurent, Two- temperature solutions and emergent spectra from rela- tivistic accretion discs around black holes, A&A642, A209 (2020), arXiv:2007.00919 [astro-ph.HE]

    2020 arXiv

  50. [58]

    Sarkar, K

    S. Sarkar, K. Singh, I. Chattopadhyay, and P. Lau- rent, Two-temperature accretion flows around strongly magnetized stars and their spectral analysis, Monthly Notices of the Royal Astronomical Society522, 3735 (2023), https://academic.oup.com/mnras/article- pdf/522/3/3735/5024...

    2023

  51. [59]

    Mitra, S

    S. Mitra, S. M. Ghoreyshi, A. Mosallanezhad, S. Ab- bassi, and S. Das, Global transonic solution of hot accretion flow with thermal conduction, Monthly Notices of the Royal Astronomical Society523, 4431 (2023), https://academic.oup.com/mnras/article- pdf/523/3/4431/50660094/st...

    2023

  52. [60]

    A. A. Esin, Heating and Cooling of Hot Accretion Flows by Nonlocal Radiation, Astrophys. J.482, 400 (1997), arXiv:astro-ph/9701039 [astro-ph]

    1997 arXiv

  53. [61]

    Martin-Vaquero, A 17th-order radau iia method for package radau

    J. Martin-Vaquero, A 17th-order radau iia method for package radau. applications in mechanical systems, Com- puters and Mathematics with Applications59, 2464 (2010)

    2010

  54. [62]

    Martin-Vaquero and A

    J. Martin-Vaquero and A. Kleefeld, Eserk5: A fifth-order extrapolated stabilized explicit runge–kutta method, Journal of Computational and Applied Mathematics 356, 22 (2019)

    2019

  55. [63]

    Maurya, V

    P. Maurya, V. Yadav, B. Mahato, G. Naveen, M. Ra- jpoot, and Y. Bhumkar, New optimized implicit-explicit runge-kutta methods with applications to the hyperbolic conservation laws, Journal of Computational Physics 446, 110650 (2021)

    2021

  56. [64]

    Hairer and G

    E. Hairer and G. Wanner, Stiff differential equations solved by radau methods, Journal of Computational and Applied Mathematics111, 93 (1999)

    1999

  57. [65]

    Pareschi and G

    L. Pareschi and G. Russo, Implicit-explicit runge-kutta schemes for stiff systems of differential equations (2001) pp. 269–288

    2001

  58. [66]

    Butcher, Numerical methods for ordinary differential equations in the 20th century (2001)

    J. Butcher, Numerical methods for ordinary differential equations in the 20th century (2001)

    2001

  59. [67]

    Kennedy and M

    C. Kennedy and M. Carpenter, Diagonally implicit runge-kutta methods for ordinary differential equations. a review (2016)

    2016

  60. [68]

    Takei and Y

    Y. Takei and Y. Iwata, Numerical scheme based on the implicit runge-kutta method and spectral method for cal- culating nonlinear hyperbolic evolution equations, Ax- ioms11, 28 (2022)

    2022

  61. [69]

    Paczy´ nsky and P

    B. Paczy´ nsky and P. J. Wiita, Thick Accretion Disks and Supercritical Luminosities, A&A88, 23 (1980)

    1980

  62. [70]

    Sarkar and I

    S. Sarkar and I. Chattopadhyay, General relativistic two-temperature accretion solutions for spherical flows around black holes, International Journal of Modern Physics D28, 1950037 (2019), arXiv:1811.05947 [astro- ph.HE]

    2019 arXiv

  63. [71]

    Sadowski, Slim disks around kerr black holes revisited, The Astrophysical Journal Supplement Series183, 171 (2009)

    A. Sadowski, Slim disks around kerr black holes revisited, The Astrophysical Journal Supplement Series183, 171 (2009)

    2009

  64. [72]

    Kumar, C

    R. Kumar, C. B. Singh, I. Chattopadhyay, and S. K. Chakrabarti, Effect of the flow composition on outflow rates from accretion discs around black holes, Monthly Notices of the Royal Astronomical Society436, 2864 (2013), arXiv:1310.0144 [astro-ph.HE]

    2013 arXiv

  65. [73]

    Maity, M

    S. Maity, M. A. Shaikh, P. Tarafdar, and T. K. Das, Carter-penrose diagrams for emergent spacetime in ax- isymmetrically accreting black hole systems, Physical Re- view D106, 10.1103/physrevd.106.044062 (2022)

    2022 doi

  66. [74]

    K. Pal, S. Ghose, S. Sarkar, and T. K. Das, Effect of spin on the dynamics of multi-component trans- relativistic accretion flows around Kerr black holes, arXiv e-prints , arXiv:2503.04321 (2025), arXiv:2503.04321 [astro-ph.HE]

    2025 arXiv

  67. [75]

    L. D. Landau and E. M. Lifshitz,Fluid mechanics(1959)

    1959

  68. [76]

    J. R. Cash and A. H. Karp, A variable order runge-kutta method for initial value problems with rapidly varying right-hand sides, ACM Trans. Math. Softw.16, 201–222 (1990)

    1990

  69. [77]

    W. H. Press, B. P. Flannery, and S. A. Teukolsky,Nu- merical recipes. The art of scientific computing(1986)

    1986

  70. [78]

    Naldi and L

    G. Naldi and L. Pareschi, Numerical schemes for hyper- bolic systems of conservation laws with stiff diffusive re- laxation, SIAM Journal on Numerical Analysis (2000)

    2000

  71. [79]

    Hairer and G

    E. Hairer and G. Wanner,Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems, Vol. 14 (1996)

    1996

  72. [80]

    K. E. Atkinson,An Introduction to Numerical Analysis, 2nd ed. (John Wiley & Sons, New York, 1989)

    1989

  73. [81]

    S¨ uli and D

    E. S¨ uli and D. F. Mayers,An Introduction to Numerical Analysis(Cambridge University Press, 2003)

    2003

  74. [82]

    U. M. Ascher and L. R. Petzold,Computer Methods for Ordinary Differential Equations and Differential- Algebraic Equations(Society for Industrial and Applied Mathematics, Philadelphia, 1998)

    1998

  75. [83]

    Butcher,Numerical Methods for Ordinary Differential Equations(John Wiley & Sons, Chichester, UK, 2003)

    J. Butcher,Numerical Methods for Ordinary Differential Equations(John Wiley & Sons, Chichester, UK, 2003)

    2003

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.