Engineering topics

Lai, Dong

Publications and source records attributed to Lai, Dong.

13 records

DOE OSTI2024

Dynamical instability in multi-orbiter systems with gas friction

ABSTRACT Closely packed multiplanet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in active galactic nuclei discs around supermassive black holes. We use N-body integrations to evaluate how the frictional damping of orbital eccentricity affects the growth of dynamical instability for a wide range of K (the difference in the planetary semimajor axes in units of the mutual Hill radius) and (unequal) planet masses. We find that, in general, the stable region (large K) and unstable region (small K) are separated by a “grey zone”, where the (in)stability is not guaranteed. We report the numerical values of the critical spacing for stability Kcrit and the “grey zone” range in different systems, and provide fitting formulae for arbitrary frictional forcing strength. We show that the stability of a system depends on the damping time-scale τ relative to the zero-friction instability growth time-scale tinst: two-planet systems are stable if tinst ≳ τ; three-planet systems require tinst ≳ 10τ−100τ. When K is sufficiently small, tinst can be less than the synodic period between the planets, which makes frictional stabilization unlikely to occur. As K increases, tinst tends to grow exponentially, but can also fluctuate by a few orders of magnitude. We also devise a linear map to analyse the dynamical instability of the “planet + test mass” system, and find qualitative agreement with N-body simulations.

Li, Jiaru (ORCID:0000000155507421)↗

NASA NTRS2013

XIPE the X-Ray Imaging Polarimetry Explorer

X-ray polarimetry, sometimes alone, and sometimes coupled to spectral and temporal variability measurements and to imaging, allows a wealth of physical phenomena in astrophysics to be studied. X-ray polarimetry investigates the acceleration process, for example, including those typical of magnetic reconnection in solar flares, but also emission in the strong magnetic fields of neutron stars and white dwarfs. It detects scattering in asymmetric structures such as accretion disks and columns, and in the so-called molecular torus and ionization cones. In addition, it allows fundamental physics in regimes of gravity and of magnetic field intensity not accessible to experiments on the Earth to be probed. Finally, models that describe fundamental interactions (e.g. quantum gravity and the extension of the Standard Model) can be tested. We describe in this paper the X-ray Imaging Polarimetry Explorer (XIPE), proposed in June 2012 to the first ESA call for a small mission with a launch in 2017. The proposal was, unfortunately, not selected. To be compliant with this schedule, we designed the payload mostly with existing items. The XIPE proposal takes advantage of the completed phase A of POLARIX for an ASI small mission program that was cancelled, but is different in many aspects: the detectors, the presence of a solar flare polarimeter and photometer and the use of a light platform derived by a mass production for a cluster of satellites. XIPE is composed of two out of the three existing JET-X telescopes with two Gas Pixel Detectors (GPD) filled with a He-DME mixture at their focus. Two additional GPDs filled with a 3-bar Ar-DME mixture always face the Sun to detect polarization from solar flares. The Minimum Detectable Polarization of a 1 mCrab source reaches 14 in the 210 keV band in 105 s for pointed observations, and 0.6 for an X10 class solar flare in the 1535 keV energy band. The imaging capability is 24 arcsec Half Energy Width (HEW) in a Field of View of 14.7 arcmin 14.7 arcmin. The spectral resolution is 20 at 6 keV and the time resolution is 8 s. The imaging capabilities of the JET-X optics and of the GPD have been demonstrated by a recent calibration campaign at PANTER X-ray test facility of the Max-Planck-Institut fr extraterrestrische Physik (MPE, Germany). XIPE takes advantage of a low-earth equatorial orbit with Malindi as down-link station and of a Mission Operation Center (MOC) at INPE (Brazil). The data policy is organized with a Core Program that comprises three months of Science Verification Phase and 25 of net observing time in the following 2 years. A competitive Guest Observer program covers the remaining 75 of the net observing time.

Astronomy X-ray Polarimetry↗

NASA NTRS1995

Hydrodynamics of coalescing binary neutron stars: Ellipsoidal treatment

We employ an approximate treatment of dissipative hydrodynamics in three dimensions to study the coalescence of binary neutron stars driven by the emission of gravitational waves. The stars are modeled as compressible ellipsoids obeying a polytropic equation of state; all internal fluid velocities are assumed to be linear functions of the coordinates. The hydrodynamics equations then reduce to a set of coupled ordinary differential equations for the evolution of the principal axes of the ellipsoids, the internal velocity parameters, and the binary orbital parameters. Gravitational radiation reaction and viscous dissipation are both incorporated. We set up exact initial binary equilibrium configurations and follow the transition from the quasi-static, secular decay of the orbit at large separation to the rapid dynamical evolution of the configurations just prior to contact. A hydrodynamical instability resulting from tidal interactions significantly accelerates the coalescence at small separation, leading to appreciable radial infall velocity and tidal lag angles near contact. This behavior is reflected in the gravitational waveforms and may be observable by gravitational wave detectors under construction. In cases where the neutron stars have spins which are not aligned with the orbital angular momentum, the spin-induced quadrupole moment can lead to precession of the orbital plane and therefore modulation of the gravitational wave amplitude even at large orbital radius. However, the amplitude of the modulation is small for typical neutron star binaries with spins much smaller than the orbital angular momentum.

Lai, Dong↗

NASA NTRS1995

Gravitational radiation from rapidly rotating nascent neutron stars

We study the secular evolution and gravitational wave signature of a newly formed, rapidly rotating neutron star. The neutron star may arise from core collapse in a massive star or from the accretion-induced collapse of a white dwarf. After a brief dynamical phase, the nascent neutron star settles into an axisymmetric, secularly unstable state. Gravitational radiation drives the star to a nonaxisymmetric, stationary equilibrium configuration via the bar-mode instability. The emitted quasi-periodic gravitational waves have a unique signature: the wave frequency sweeps downward from a few hundred Hertz to zero, while the wave amplitude increase from zero to a maximum and then decays back to zero. Such a wave signal could detected by broadband gravitational wave interferometers currently being constructed. We also characterize two other types of gravitational wave signals that could arise in principle from a rapidly rotating, secularly unstable neutron star: a high-frequency (f greater than or approximately = 1000 Hz) wave which increases the pattern-speed of the star, and a wave that actually increases the angular momentum of the star.

Lai, Dong↗

NASA NTRS1995

Spin-Orbit Interaction in Neutron Star/Main Sequence Binaries and Implications for Pulsar Timing

The spin-induced quadrupole moment of a rapidly rotating star changes the orbital dynamics in a binary system, giving rise to advance (or regression) of periastron and precession of the orbital plane. We show that these effects are important in the recently discovered radio pulsar/mail sequence star binary system PSR J0045-7319, and reliably account for the observed peculiar timing residuals. Precise measurements of the apsidal motion and orbital plane precession can yield valuable information on the internal structure and rotation of the star. The detection of orbital precession implies that the spin of the companion star is not aligned with the orbital angular momentum, and suggests that the supernova gave the pulsar a kick out of the original orbital plane. Excitations of g-mode oscillations near periastron (the dynamical tide) can induce measurable changes in the orbital period and eccentricity at each passage for this system. We also discuss the spin-orbit coupling effects for the accreting X-ray pulsars and the other know radio pulsar/main sequence binary, PSR B1259-63.

neutron stars↗

NASA NTRS1994

Hydrodynamics of rotating stars and close binary interactions: Compressible ellipsoid models

We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids, and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest, such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.

Lai, Dong↗

NASA NTRS1994

Equilibrium, stability, and orbital evolution of close binary systems

We present a new analytic study of the equilibrium and stability properties of close binary systems containing polytropic components. Our method is based on the use of ellipsoidal trial functions in an energy variational principle. We consider both synchronized and nonsynchronized systems, constructing the compressible generalizations of the classical Darwin and Darwin-Riemann configurations. Our method can be applied to a wide variety of binary models where the stellar masses, radii, spins, entropies, and polytropic indices are all allowed to vary over wide ranges and independently for each component. We find that both secular and dynamical instabilities can develop before a Roche limit or contact is reached along a sequence of models with decreasing binary separation. High incompressibility always makes a given binary system more susceptible to these instabilities, but the dependence on the mass ratio is more complicated. As simple applications, we construct models of double degenerate systems and of low-mass main-sequence star binaries. We also discuss the orbital evoltuion of close binary systems under the combined influence of fluid viscosity and secular angular momentum losses from processes like gravitational radiation. We show that the existence of global fluid instabilities can have a profound effect on the terminal evolution of coalescing binaries. The validity of our analytic solutions is examined by means of detailed comparisons with the results of recent numerical fluid calculations in three dimensions.

Lai, Dong↗

NASA NTRS1994

Hydrodynamic instability and coalescence of binary neutron stars

We study the importance of hydrodynamic effects on the evolution of coalescing binary neutron stars. Using an approximate energy functional constructed from equilibrium solutions for polytropic binary configuration, we incorporate hydrodynamic effects into the calculation of the orbital decay driven by gravitational wave emission. In particular, we follow the transition between the quasi-static, secular decay of the orbit at large separation and the rapid dynamical evolution of confirgurations approaching contact. We show that a purely Newtonian hydrodynamic instability can significantly accelerate the coalescence at small separation. Such an instability occurs in all close binary configurations containing sufficiently incompressible stars. Calculations are performed for various neutron star masses, radii, and spins. The influence of the stiffness of the equation of state is also explored by varying the effective polytopic index. Typically, we find that the radial infall velocity just prior to contact is about 10% of the tangential orbital velocity. Once the stability limit is reached, the final evolution only takes another orbit. Post-Newtonian effects can move the stability limit to a larger binary separation, and may induce an even larger radial velocity. We also consider the possibiltiy of mass transfer form one neutron star to the other. We show that stable mass transfer is unlikely except when the mass of one of the components is very small (M is less than or approximately 0.4 solar mass) and the viscosity is high enough to maintain corotation. Otherwise, either the two stars come into contact or the dynamical instability sets in before a Roche limit can be reached.

Lai, Dong↗

NASA NTRS1993

Ellipsoidal figures of equilibrium - Compressible models

The results of Chandrasekhar (1969) are generalized to polytropes, using a formalism based on ellipsoidal energy variational principle to construct approximate stellar equilibrium solutions and study their stability. After reviewing the energy variational method and describing the approach, several equivalent stability conditions are established and secular vs. dynamical instabilities are discussed. Then, the equilibrium structure equations are derived for isolated, rotating polytropes, and axisymmetric configurations (compressible Maclaurin spheroids) are considered. Particular attention is given to triaxial configurations, either in a state of uniform rotation (generalizing the classical Jacobi ellipsoids) or with internal fluid motions of uniform vorticity (the compressible analogues of Riemann-S ellipsoids) and to the stability of these single star configurations. The compressible generalizations of the Roche and Roche-Riemann problems for a polytrope in orbit about a point-mass companion are solved, and the generalized Darwin problem for two identical polytropes in a binary is considered.

Lai, Dong↗

NASA NTRS1993

Collisions and close encounters between massive main-sequence stars

Smooth-particle hydrodynamics is used to investigate collisions and close encounters between two massive main-sequence stars. Approximate scaling relations and fitting formulas are given for the amount of mass loss and for the critical impact parameters for capture or merging. The more massive stars, which have smaller ratios of specific binding energy to the square of escape velocity, are more easily disrupted in collisions. As the relative velocity increases, the orbital energy of the colliding stars can only be dissipated by shock heating, and the critical capture radius decreases much faster than predicted by linear theory. The average fractional mass loss per collision in a cluster is found not to depend sensitively on the stellar velocity dispersion.

Lai, Dong↗

NASA NTRS1993

Hydrodynamic instability and coalescence of close binary systems

The meaning of hydrodynamic instability of close binary systems is clarified, and the resulting large (but finite) value of the radial infall velocity near the end of the coalescence is estimated. Emphasis is placed on the case where the orbital decay is driven by gravitational radiation. A simple calculation of the orbital decay is presented for a simple generic model of a close binary.

Lai, Dong↗

NASA NTRS1992

Hydrogen molecules and chains in a superstrong magnetic field

The electronic structures of hydrogen polymolecules H(n) (n = 2,3,4,...) is studied in a superstrong magnetic field (B greater than about 10 exp 12 G) typically found on the surface of a neutron star. Simple analytical scaling relations for several limiting cases (e.g., large n, high B field) are derived. The binding energies of H(n) molecules are numerically calculated for various magnetic-field strengths. For a given magnetic-field strength, the binding energy per atom in the H(n) molecules is found to approach a constant value as n increases. For typical field strengths of interest, energy saturation is essentially achieved once n exceeds 3 to 4. Also considered is the structure of negative H ions in a high magnetic field. For B about 10 exp 12 G, the dissociation energy of an atom in a hydrogen chain and the ionization potential of H(-) are smaller than the ionization potential of neutral atomic hydrogen.

Lai, Dong↗

NASA NTRS1991

Cold equation of state in a strong magnetic field - Effects of inverse beta-decay

The influence of a high magnetic field (B is greater than 10 exp 12 G) on the degenerate matter equation of state appropriate to a neutron star is studied. The regime dominated by relativistic electrons up to the neutron drip density is highlighted. The equilibrium matter composition and equation of state, allowing for inverse beta-decay. Two different equilibrium models are determined: an ideal neutron-proton-electron (npe) gas and the more realistic model of Baym, Pethick, and Sutherland (1971) consisting of a Coulomb lattice of heavy nuclei embedded in an electron gas. For a sufficiently high field strength, the magnetic field has an appreciable effect, changing the adiabatic index of the matter and the nuclear transition densities. The influence of a strong field on some simple nonequilibrium processes, including beta-decay and inverse beta-decay (electron capture) is also considered. The effects produced by the magnetic field are mainly due to the changes in the transverse electron quantum orbits and the allowed electron phase space induced by the field.

Lai, Dong↗