The energy-dependent radiation transfer in a homogeneous finite sphere

Radiation transfer problem in the slowing-down region is solved for a spherical geometry. Modelled kernels are used to represent the slowing-down kernel. Calculations are performed for the radiation flux at the boundary of the sphere for two cases: (i) in the presence of internal source, (ii) if the radiation is incident on the sphere. Numerical results are obtained by using the Padé approximant technique.     

Energy-dependent transport problem with generalized boundary conditions

The slowing-down Boltzmann equation for generalized boundary conditions is considered and transformed to one-speed equation in Laplace space. Exact relations between energy reflection and transmission coefficients for a problem with diffuse reflecting boundary conditions and the albedos for the problem with isotropic boundary conditions are obtained. The Galerkin method is used to calculate the energy reflection coefficient for a finite slab for different thicknesses at different mass ratiosA, target to projectile mass, at different synthetic-scattering kernels. The results for partial heat fluxes for isotropic and anisotropic-scattering dispersive medium are given. The results obtained for isotropic boundary conditions are compared well with the exact results.     

Braking index diagnostics of pulsars. I. alignment,Counteralignment and slowing-down noise

We show that the strong correlation observed between the braking indices (n) and the slowing-down ages (Τ) of pulsars is inconsistent with counteralignment between their rotation and magnetic axes, but that the data on pulsars with positive braking indices is consistent with alignment. Alternatively, slowing-down noise can quantitatively account for the data on all pulsars except the Crab and the Vela, and so for the apparent|n| ∼ Τ² correlation observed for the older pulsars.     

Maximum entropy for half-space energy albedo

The maximum entropy method is used to find the energy albedo for half-space, i.e., the classical albedo problem but extended to include energy dependence. The slowing-down kernel is described by elastic scattering and various model representations for dealing with anisotropic scattering. The energy albedo for different values of secondary particles per collisionC is calculated and compared with exact and variational values given by Williams and Padé's approximant given by Machaliet al. (1981). Also the energy albedo for different mass-ratios is found and compared with Williams and Padé's results.     

Particle transfer in decelerating regions

Equations connecting the transmission and reflection functions of a finite medium to the reflection function of a semi-infinite one are used to get the albedos in the slowing-down region. The transport equation is solved by the modified Eddington method using Fermi's backward-forward scattering model, modulated to allow for different orders of backward-forward and isotropic scattering. Numerical results for the energy albedos are obtained and compared.On leave of absence from Atomic Energy Center, Inchass, Egypt.     

A fourier transform solution of radiation transfer in a spherical-inhomogeneous medium

The radiative heat flux at the boundary of a sphere containing an internal energy source and subject to general boundary conditions (problem 1) is obtained in terms of the albedo of the corresponding source-free problem with isotropic boundary condition (problem 2). The solution of problem 2 is performed on the basis of the integral Fourier transforms method. Numerical results for the partial heat flux and emissivity for a given internal energy source and inhomogeneous medium, isotropic scattering are obtained.     

Energy albedo for multiregion slabs

An asymptotic solution for the equation of radiative transfer in an inhomogeneous medium was obtained on the basis of the corresponding solutions for homogeneous sub-layers in the slowing down region. Function relations between the reflection and transmission coefficients for the whole slab and those of the sublayers are given. The invariant embedding concepts are used to get the reflection and transmission coefficients for the sub-layers. We assumed different models for the slowing-down kernels. Laplace transform was used to transform the Boltzmann equation to one velocity approximation with re-scaled mean-free path and single-scattering albedo. Numerical results are given for energy albedo as a function of the mass number of the host medium.     

Radiation transfer through atmospheric aerosol media with anisotropic scattering

Radiation transfer in atmospheric aerosol media with general boundary conditions has been studied for anisotropic scattering. The considered aerosol medium assumed to have specular and diffused reflecting boundary surfaces and in the presence of internal source. The radiation transfer scattering parameters as single scattering albedo, asymmetry factor, scattering, absorption, extinction efficiencies and anisotropic scattering coefficient have been calculated using the Mie theory. The problem with general boundary conditions is solved in terms of the solution of source-free problem with simply boundary conditions. Pomraning-Eddington approximation is used to solve the source-free problem. For the sake of comparison, a weight function is introduced and used in two special forms. The calculated partial heat fluxes with the two methods are compared and showed good agreement. Some of our results are found in a good agreement with published data.     

Surface heat current in an inhomogeneous slab with reflecting walls

The problem of radiative transfer in an inhomogeneous finite medium with boundary surfaces which reflects both diffusely and specularly is connected with a source-free specular boundary condition radiation tramsfer. Numerical results are obtained for the partial heat flux using the bi-variational technique.     

A new approach to impulsive rendezvous near circular orbit

A new approach is presented for the problem of planar optimal impulsive rendezvous of a spacecraft in an inertial frame near a circular orbit in a Newtonian gravitational field. The total characteristic velocity to be minimized is replaced by a related characteristic-value function and this related optimization problem can be solved in closed form. The solution of this problem is shown to approach the solution of the original problem in the limit as the boundary conditions approach those of a circular orbit. Using a form of primer-vector theory the problem is formulated in a way that leads to relatively easy calculation of the optimal velocity increments. A certain vector that can easily be calculated from the boundary conditions determines the number of impulses required for solution of the optimization problem and also is useful in the computation of these velocity increments. Necessary and sufficient conditions for boundary conditions to require exactly three nonsingular non-degenerate impulses for solution of the related optimal rendezvous problem, and a means of calculating these velocity increments are presented. A simple example of a three-impulse rendezvous problem is solved and the resulting trajectory is depicted. Optimal non-degenerate nonsingular two-impulse rendezvous for the related problem is found to consist of four categories of solutions depending on the four ways the primer vector locus intersects the unit circle. Necessary and sufficient conditions for each category of solutions are presented. The region of the boundary values that admit each category of solutions of the related problem are found, and in each case a closed-form solution of the optimal velocity increments is presented. Similar results are presented for the simpler optimal rendezvous that require only one-impulse. For brevity degenerate and singular solutions are not discussed in detail, but should be presented in a following study. Although this approach is thought to provide simpler computations than existing methods, its main contribution may be in establishing a new approach to the more general problem.     

A Shotgun Model for Gamma-Ray Bursts

We propose that gamma-ray bursts (GRBs) are produced by a shower of heavy blobs running into circumstellar material at highly relativistic speeds. The gamma-ray emission is produced in the shocks these bullets drive into the surrounding medium. The short-term variability seen in GRBs is set by the slowing-down time of the bullets, while the overall duration of the burst is set by the lifetime of the central engine. A requirement of this model is that the ambient medium be dense, consistent with a strong stellar wind. The efficiency of the burst can be relatively high.     

Non-axisymmetric magnetostatic equilibrium

This paper considers the structural properties of a sunspot-like magnetic flux tube which lacks perfect axisymmetry. The flux tube is taken to be in static equilibrium with an atmosphere in a uniform gravity. Assuming the departure from axisymmetry to be slight, the equations for the first order non-axisymmetric part of the equilibrium are derived in cylindrical coordinates. These first order equations reduce to a linear second order hyperbolic partial differential equation in the r-z plane. Whereas Cauchy type boundary conditions are appropriate for hyperbolic equations, physical considerations dictate the specification of boundary conditions on a closed surve for our problem of interest. The construction of solutions to this boundary value problem is illustrated with three analytically soluble cases, where the zero-order axisymmetric equilibria are chosen to have magnetic field geometry of different complexity. A physical discussion of the results is given.     

The heat flux calculation for the laser-fusion problem

The problem of heat flux at the critical surfaces and the surfaces of a pellet of deuterium and tritium (conduction zone) heated by laser have been considered. Ion-electron collisions are only allowed for: i.e., the linear transport equation is used to describe the problem with boundary conditions consists of isotropic and diffuse boundary conditions. The bi-variational technique has been used to calculate the electron density and temperature across the conduction zone as well as the heat flux. Numerical results are given and compared with those of Rouse and Williams (1981) results.     

Symplectic adaptive algorithm for solving nonlinear two-point boundary value problems in Astrodynamics

In this paper, from a Hamiltonian point of view, the nonlinear optimal control problems are transformed into nonlinear two-point boundary value problems, and a symplectic adaptive algorithm based on the dual variational principle is proposed for solving the nonlinear two-point boundary value problem. The state and the costate variables within a time interval are approximated by using the Lagrange polynomial and the costate variables at two ends of the time interval are taken as independent variables. Then, based on the dual variational principle, the nonlinear two-point boundary value problems are replaced by a system of nonlinear equations which can preserve the symplectic structure of the nonlinear optimal control problem. Furthermore, the computational efficiency of the proposed symplectic algorithm is improved by using the adaptive multi-level iteration idea. The performance of the proposed algorithm is tested by the problems of Astrodynamics, such as the optimal orbital rendezvous problem and the optimal orbit transfer between halo orbits.     

Anisotropic radiation transfer in a sphere with internal source

In this work we connect the problem of anisotropic radiation transfer in a sphere with diffuse reflectivity and containing an energy source with the problem of source free amisotropic radiation transfer with isotropic boundary condition. Exact equation for radiation heat flux for the first problem is obtained in terms of the source, the flux and the albedo of the second problem. Modeled kernel is used to represent the anisotropy of the phase function. Numerical results are given.     

On some problems concerning the calculation of the form of the magnetopause and the electric field on the magnetopause in the framework of the continuum medium model

In order to understand the reason of the existence of the electric field in the magnetosphere, and for the theoretical evaluation of its value, it is necessary to find the solution of the problem of determination of the magnetosphere boundary form in the frameworks of the continuum medium model which takes into account part of the magnetospheric plasma movement in supporting the magnetospheric boundary equilibrium. A number of problems for finding the distribution of the pressure, the density, the magnetic field and the electric field on the particular tangential discontinuity is considered in the case when the form of discontinuity is set (the direct problem) and a number of problems for finding the form of the discontinuity and the distribution of the above-mentioned physical quantities on the discontinuity is considered when the law of the change of the external pressure along the boundary is set (for example, with the help of the approximate Newton equation). The problem which is considered here, which deals with the calculation of the boundary form and with the calculation of the distribution of the corresponding physical quantities on the discontinuity of the 1st kind for the compressible fluid with the magnetic field with field lines which are perpendicular to the plane of the flow in question, concerns the last sort of problems. The comparison of the results of the calculation with the data in the equatorial cross-section of the magnetosphere demonstrates that the calculated form of the boundary, the value of the velocity of the return flow and the value of the electric field on the magnetopause, agree satisfactorily with the observational data.     

An estimate of the contribution of high nergy cosmic-ray protons to the absorbed dose inventory of a cometary nucleus

We discuss the contribution of cosmic-ray protons at all energies above 1 MeV to the absorbed doses of the surface layers of a comet. Since there exists no calculation which takes into account proton energy losses by means of losses to electrons and nuclear collisions (in a cascade process), and losses due to the low energy end-products of the cascade, we have made a rough estimate of all of these contributions. An analytical formula is proposed that allows a rapid estimate of ranges and the dose absorbed at any depth. We give dose-depth curves for two extreme values of the energy at which nuclear collisions begin to dominate the slowing-down process, and for an intermediate value we display the dose-depth curve down to 20 m from the surface. The relevance of these findings to dosimetry in comets and some alterations of cometary material are considered. The need for improving the analytical expression proposed is stressed.     

An intermediate matching technique for solving two point boundary value problems using the perturbation method

The perturbation method, a numerical method for solving two point boundary value problems (TPBVP), is modified to attempt to improve inherent instability and sensitivity problems associated with the method. The desired solution to the TPBVP is divided into two time intervals. The differential equations required to define a solution to the two point boundary value problem are integrated independently over these shorter segments rather than consecutively over the entire trajectory. The independent integration of the differential equations over approximately half of the trajectory instead of the entire trajectory substantially decreases sensitivity and stability properties associated with the numerical integration. The equations for both time segments can be integrated simultaneously. By this procedure, a system of twice the dimension of the original problem is integrated for a period of time equal to half of the time interval for the original problem. To show the effectiveness of the method, two impulse trajectories which minimize the total velocity increment required to transfer a spacecraft from an Earth orbit into a lunar orbit are calculated.     

Boltzmann-Fokker-Planck model for the electron distribution function in the Earth's ionosphere

Using Boltzmann-Fokker-Planck methods and the diffusion approximation, we derive coupled non-linear equations for the first two angular moments of the electron distribution function in the Earth's ionosphere. The theory includes a phenomenological treatment of photionization of the neutral species by an externally produced photon flux; electron-ion recombination; electronneutral particle attachment; elastic, excitation, deexcitation, and ionizing electron-neutral particle collisions; and elastic electron-electron and electron-ion collisions. At high and low energies, we obtain approximate analytic solutions for the steady-state electron distribution function. Under certain conditions we also obtain the standard continuous slowing-down formulae for the steady-state electron distribution function in the diffusion approximation.     

Secular Love numbers of the major planets

Secular Love numbers of the major planets have been determined by solving the first (Dirichlet's) boundary — value problem for the level ellipsoid as the boundary surface. It has been demonstrated that the rotational distortions are responsible for the actual figure of the bodies above. The 6th degree Stokes zonal parameter for Uranus and for Neptune have been predicted.