SEARCH · Search NASA

Results for “spherical maps”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2000

Testing Models of Resonant Compton Scattering in X-Ray Pulsars

Over the performance period covered by the grant, the principal investigator modified a Monte Carlo Compton scattering code to model the propagation of x-rays through the magnetosphere of accreting neutron stars. These modifications were made to enable the author to compare the observations of x-ray pulsars to theoretical models of the system. The original code was designed to study relativistic plasmas with one of two geometries: a plane parallel plasma with a differential relativistic bulk velocity, and a static spherically symmetric plasma.- This code did not treat gravitational bending or bulk motion in the magnetosphere of a neutron star. Under the grant, the author incorporated code to trace light paths in a Schwarzschild metric. The code was modified to keep track of the photon polarization during propagati on. The investigator also modified the code so that bulk motion in an axisymmetric system is treated properly. An approximate treatment for resonant Compton scattering was added to the code. Finally, code was added that creates model observables that can be compared to observations, such as projected x-ray emission maps and energy-dependent light curves. Comparison to observations is now commencing.

Brainerd, Jerome J.↗

NASA NTRS2004

Impact Processes in the Solar System

Our laboratory has previously conducted impact fracture and dynamic failure tests. Polanskey and Ahrens [1990] mapped the fractures from a series of laboratory craters (Fig. 1) and Ahrens and Rubin [ 1993] inferred that the usually further extending radial cracks resulted from tensional failure during the compression of the shock propagation. The radial spreading induced by the particle velocity field caused the stresses perpendicular to the shock front to become sufficiently large and tensile. This induces "radial fractures." The concentric fractures are attributed to the tensional failure occurring after the initial compressive phase. Upon radial propagation of the stress wave the negative tension behind the stress-wave front caused failure along the quasi-spherical concentric fractures. The near-surface and spall fractures are attributed to the fractures described by Melosh [1984]. These are activated by impact and can launch relatively unshocked samples of planetary surfaces to speeds exceeding escape velocity. In the case of Mars, some of these surface samples presumably become the SNC (Mars) meteorites.

Ahrens, Thomas J.↗

NASA NTRS1991

Mimas: Tectonic structure and geologic history

Mimas, the innermost of the major saturnian satellites, occupies an important place in comparative studies of icy satellites. It is the smallest icy satellite known to have a mostly spherical shape. Smaller icy objects like Hyperion and Puck are generally irregular in shape, while larger ones like Miranda and Enceladus are spherical. Thus Mimas is near the diameter where the combination of increasing surface gravity and internal heating begin to have a significant effect on global structure. The nature and extent of endogenic surface features provide important constraints on the interior structure and history of this transitional body. The major landforms on Mimas are impact craters. Mimas has one of the most heavily cratered surfaces in the solar system. The most prominent single feature on Mimas is Herschel, an unrelaxed complex crater 130 km in diameter. The only other recognized landforms on Mimas are tectonic grooves and lineaments. Groove locations were mapped by Schenk, but without analysis of groove structures or superposition relationships. Mimas' tectonic structures are remapped here in more detail than previously has been done, as part of a general study of tectonic features on icy satellites.

Croft, Steven K.↗

NASA NTRS2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be-+12.390.981.30km and-+2.0720.0660.067Me respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is=-+TlogK5.99100.060.05([])for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters(within the prior support)and can thus be computed from the posterior samples.

Millisecond pulsars↗

DOE OSTI2026

Modeling of injected current stream-induced 3D perturbations in local helicity injection plasmas

Solenoid-free tokamak startup techniques are essential for spherical tokamaks and offer a pathway to cost reduction and design simplification in fusion energy systems. Local helicity injection (LHI) is one such approach, employing compact edge current sources to drive open field line current that initiates and sustains tokamak plasmas. The recently commissioned Pegasus-III Experiment provides a platform for advancing this and other solenoid-free startup methods. This study investigates the effect of LHI on magnetic topology in Pegasus-III plasmas. A helical filament model represents the injected current, and the linear plasma response to its three-dimensional field is calculated with M3D-C1. Poincaré mapping reveals substantial flux surface degradation in all modeled cases. The onset of overlapping magnetic structures and large-scale surface deformation begins near Ψ N ≈ 0.37, indicating a broad region of perturbed topology extending toward the edge. In rotating plasmas, both single-fluid and two-fluid models exhibit partial screening of the n = 1 perturbation, with two-fluid calculations showing stronger suppression near the edge. In contrast, the absence of rotation leads to strong resonant field amplification in the single-fluid case, while the two-fluid case with zero electron rotation mitigates this amplification and preserves edge screening. Magnetic probe measurements indicate that modeling the stream with spatial spreading—representing distributed current and/or oscillatory motion—better reproduces measured magnetic power profiles than a rigid filament model. The results underscore the role of rotation and two-fluid physics in screening stream perturbations and point to plasma flow measurements and refined stream models as key steps toward improving predictive fidelity.

Schaefer, Carolyn E. [Univ. of Wisconsin, Madison,↗

NASA NTRS2003

Modified Polar-Format Software for Processing SAR Data

HMPF is a computer program that implements a modified polar-format algorithm for processing data from spaceborne synthetic-aperture radar (SAR) systems. Unlike prior polar-format processing algorithms, this algorithm is based on the assumption that the radar signal wavefronts are spherical rather than planar. The algorithm provides for resampling of SAR pulse data from slant range to radial distance from the center of a reference sphere that is nominally the local Earth surface. Then, invoking the projection-slice theorem, the resampled pulse data are Fourier-transformed over radial distance, arranged in the wavenumber domain according to the acquisition geometry, resampled to a Cartesian grid, and inverse-Fourier-transformed. The result of this process is the focused SAR image. HMPF, and perhaps other programs that implement variants of the algorithm, may give better accuracy than do prior algorithms for processing strip-map SAR data from high altitudes and may give better phase preservation relative to prior polar-format algorithms for processing spotlight-mode SAR data.

Chen, Curtis↗

DOE OSTI2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be ${12.39}_{-0.98}^{+1.30}$ km and ${2.072}_{-0.066}^{+0.067}$ M ⊙ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is ${\mathrm{log}}_{10}(T\,[{\rm{K}}])={5.99}_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Furthermore, our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1978

The energetics of molecular clouds. II - The S140 molecular cloud

Techniques of analysis which allow estimates of kinetic temperature, gas density, gas cooling rate, and dust cooling rate are applied to the S140 molecular cloud. Maps of this cloud in various molecular lines are examined which reveal the existence of a strong temperature and density peak near the H-alpha rim and the presence of a strong near-IR source very near that peak. This source is found to have an energy distribution which is very similar to that of the Becklin-Neugebauer object, including the silicate feature at 9.5 microns. The heating and cooling rates for the gas and dust in S140 are estimated, and the predicted cooling rate for dust is shown to be much greater than that for gas. The molecular-line observations are compared with a spherical large-velocity-gradient trapping calculation. An optimum cloud model is developed which indicates that: (1) (C-13)O, HCO(+), and H2CO are about 10 times less abundant in S140 than in dark clouds; (2) the HCN abundance is approximately 4 times the H2CO abundance in S140, compared with roughly equal values in dark clouds; and (3) the (C-13)O/C(O-18) ratio is consistent with terrestrial isotope ratios.

Blair, G. N.↗

NASA NTRS1995

Solar Thermal Concept Evaluation

Concentrated solar thermal energy can be utilized in a variety of high temperature applications for both terrestrial and space environments. In each application, knowledge of the collector and absorber's heat exchange interaction is required. To understand this coupled mechanism, various concentrator types and geometries, as well as, their relationship to the physical absorber mechanics were investigated. To conduct experimental tests various parts of a 5,000 watt, thermal concentrator, facility were made and evaluated. This was in anticipation at a larger NASA facility proposed for construction. Although much of the work centered on solar thermal propulsion for an upper stage (less than one pound thrust range), the information generated and the facility's capabilities are applicable to material processing, power generation and similar uses. The numerical calculations used to design the laboratory mirror and the procedure for evaluating other solar collectors are presented here. The mirror design is based on a hexagonal faceted system, which uses a spherical approximation to the parabolic surface. The work began with a few two dimensional estimates and continued with a full, three dimensional, numerical algorithm written in FORTRAN code. This was compared to a full geometry, ray trace program, BEAM 4, which optimizes the curvatures, based on purely optical considerations. Founded on numerical results, the characteristics of a faceted concentrator were construed. The numerical methodologies themselves were evaluated and categorized. As a result, the three-dimensional FORTRAN code was the method chosen to construct the mirrors, due to its overall accuracy and superior results to the ray trace program. This information is being used to fabricate and subsequently, laser map the actual mirror surfaces. Evaluation of concentrator mirrors, thermal applications and scaling the results of the 10 foot diameter mirror to a much larger concentrator, were studied. Evaluations, recommendations and pit falls regarding the structure, materials and facility design are presented.

Hawk, Clark W.↗

NASA NTRS2002

Image Dis-Integration for Improved Plasmasphere Visualization

Traditionally, study of the plasmasphere has involved terrestrial observation of local characteristics. Global modeling of the plasmasphere in such an observation regime made use of an ensemble of (sparse) local measurements. Recently, sensors aboard the IMAGE (Imager for Magnetopause-to-Aurora Global Exploration) satellite (in particular, the EUV (Extreme Ultra Violet) Imager) have created the potential for truly global study of the plasmasphere. IMAGE was launched in spring of 2000 in an orbit with apogee altitude 7.2 R E (Earth radii) and perigee altitude 1000 km. IMAGE's EUV sensor allows an external view of the distribution of cold plasma in the plasmasphere to be acquired. EUV is designed to image light emission at 30.4 nanometers, which is the emission wavelength of the He + ion in the presence of solar radiation. He + makes up approximately 15-20% of the plasma in the plasmasphere, thus imaging of He + enables determination of plasma distribution. The EUV instrument provides a 90° by 84° field of view which is imaged as an equally spaced 150x140 pixel array on a spherical imaging surface. The EUV produces an image approximately every 10 minutes when the sensor is operating. Since EUV images contain line-of-sight integrations of plasma distributions, they do not directly express equatorial plane density (which would enable comparison of observed plasma distributions with predictions from models). Furthermore, the plasma density at any point in three-space is not known. The goal of our work was development of a technique that can enable plasma density to be determined throughout three-space. Our approach to creation of a three-space representation of the plasma distribution involves disintegrating the EUV lines of sight to form a volumetric map of plasma densities.

Timothy S. Newman↗

NASA NTRS1995

(abstract) Venus Gravity Field

A global gravity field model of Venus to degree and order 75 (5772 spherical harmonic coefficients) has been estimated from Doppler radio tracking of the orbiting spacecraft Pioneer Venus Orbiter (1979-1992) and Magellan (1990-1994). After the successful aerobraking of Magellan, a near circular polar orbit was attained and relatively uniform gravity field resolution (approximately 200 km) was obtained with formal uncertainties of a few milligals. Detailed gravity for several highland features are displayed as gravity contours overlaying colored topography. The positive correlation of typography with gravity is very high being unlike that of the Earth, Moon, and Mars. The amplitudes are Earth-like, but have significantly different gravity-topography ratios for different features. Global gravity, geoid, and isostatic anomaly maps as well as the admittance function are displayed.

Venus gravity field model gravity-topography ratio↗

NASA NTRS2022

A Study of Photoionized Gas in Two H װ Regions of the N44 Complex in the LMC Using MUSE Observations

We use the optical integral field observations with Multi-Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope, together with CLOUDY photoionization models, to study ionization structure and physical conditions of two luminous H װ regions in the N44 star-forming complex of the Large Magellanic Cloud. The spectral maps of various emission lines reveal a stratified ionization geometry in N44 D1. The spatial distribution of [O ו] λ6300 emission in N44 D1 indicates a partially covered ionization front at the outer boundary of the H װ region. These observations reveal that N44 D1 is a blister H װ region. The [O ו] λ6300 emission in N44 C does not provide a well-defined ionization front at the boundary, while patches of [S װ] λ6717 and [O ו] λ6300 emission bars are found in the interior. The results of spatially resolved MUSE spectra are tested with the photoionization models for the first time in these H װ regions. A spherically symmetric ionization-bounded model with a partial covering factor, which is appropriate for a blister H װ region, can well reproduce the observed geometry and most of the diagnostic line ratios in N44 D1. Similarly, in N44 C we apply a low-density and optically thin model based on the observational signatures. Our modeling results show that the ionization structure and physical conditions of N44 D1 are mainly determined by the radiation from an O5 V star. However, local X-rays, possibly from supernovae or stellar wind, play a key role. In N44 C, the main contribution is from three ionizing stars.

H II regions↗

NASA NTRS1988

An analytic-geometric model of the effect of spherically distributed injection errors for Galileo and Ulysses spacecraft - The multi-stage problem

In previous work the problem of injecting the Galileo and Ulysses spacecraft from low earth orbit into their respective interplanetary trajectories has been discussed for the single stage (Centaur) vehicle. The central issue, in the event of spherically distributed injection errors, is what happens to the vehicle? The difficulties addressed in this paper involve the multi-stage problem since both Galileo and Ulysses will be utilizing the two-stage IUS system. Ulysses will also include a third stage: the PAM-S. The solution is expressed in terms of probabilities for total percentage of escape, orbit decay and reentry trajectories. Analytic solutions are found for Hill's Equations of Relative Motion (more recently called Clohessy-Wiltshire Equations) for multi-stage injections. These solutions are interpreted geometrically on the injection sphere. The analytic-geometric models compare well with numerical solutions, provide insight into the behavior of trajectories mapped on the injection sphere and simplify the numerical two-dimensional search for trajectory families.

Longuski, James M.↗

NASA NTRS2003

Perfection of Apoferritin Crystals: An Advanced X-Ray Imaging and Diffraction Study

Ferritin is a well-known iron-storage protein, and is a spherical shell that consists of 24 identical subunits packed in a 432 symmetry. The typically large protein size and its distinction from lysozyme as to chemical and physical characteristics make ferritin an attractive model protein for crystal growth and perfection investigation-as an alternative to the most widely studied lysozyme. In this contribution, the latest results obtained from coherence-based x-ray diffraction imaging and diffraction experiments will be presented on octahedral apoferritin (a demetalized form of ferritin) crystals grown from various growth conditions. Crystal specimens, which have the measured rocking-curve widths varying from a few arcseconds to several tens arcseconds (or more), are comparatively examined by intrinsically highly sensitive mapping of lattice perfection and defects. The richness of the observed defects and growth features offers insight into perfection and growth of protein crystals. Beautiful interference fringe patterns formed in diffraction images and fine oscillation structure of rocking curves observed will be discussed for understanding of physical origins and the underlying impact.

Hu, Z. W.↗

NASA NTRS1993

Imaging Spectrometers Using Concave Holographic Gratings

Imaging spectroscopy combines the spatial attributes of imaging with the compositionally diagnostic attributes of spectroscopy. For spacebased remote sensing applications, mass, size, power, data rate, and application constrain the scanning approach. For the first three approaches, substantial savings in mass and size of the spectrometer can be achieved in some cases with a concave holographic grating and careful placement of an order-sorting filter. A hologram etched on the single concave surface contains the equivalent of the collimating, dispersing, and camera optics of a conventional grating spectrometer and provides substantial wavelength dependent corrections for spherical aberrations and a flat focal field. These gratings can be blazed to improve efficiency when used over a small wavelength range or left unblazed for broadband uniform efficiency when used over a wavelength range of up to 2 orders. More than 1 order can be imaged along the dispersion axis by placing an appropriately designed step order-sorting filter in front of the one- or two-dimensional detector. This filter can be shaped for additional aberration corrections. The VIRIS imaging spectrometer based on the broadband design provides simultaneous imaging of the entrance slit from lambda = 0.9 to 2.6 microns (1.5 orders) onto a 128 x 128 HgCdTe detector (at 77 K). The VIRIS spectrometer was used for lunar mapping with the UH 24.in telescope at Mauna Kea Observatory. The design is adaptable for small, low mass, space based imaging spectrometers.

Gradie, J.↗

DOE OSTI2021

Various Low-Symmetry Phases in High- χ and Conformationally Asymmetric PDMS- b -PTFEA Copolymers

Since the discovery of the Frank–Kasper (FK) phase a decade ago, single-component block copolymer systems have offered unique opportunities to explore the structural and phase-transition properties of topologically close-packed FK phases composed of point-particle motifs. Here, we present the formation of complex low-symmetry phases and their thermodynamic stability using high-χ and high-ε linear diblock copolymers, where χ and ε are the Flory–Huggins interaction parameter and conformational asymmetry between the two blocks, respectively. In this work, a series of polydimethylsiloxane-b-poly(2,2,2-trifluoroethyl acrylate) (PDMS-b-PTFEA) copolymers are prepared by tuning the volume fraction ($f_{\text{PTFEA}}$) of a minor PTFEA block. The packing structures of PDMS-b-PTFEAs are mapped in the space of temperature versus $f_{\text{PTFEA}}$ using X-ray scattering measurements. In addition to the packing structures reported earlier, we identify another A15 phase at higher $f_{\text{PTFEA}}$ close to a hexagonally packed cylinder (HEX) phase and establish the phase sequence of HEX–A15–σ–C14–BCC with decreasing $f_{\text{PTFEA}}$. Theoretical calculation of relative free-energy densities confirms the experimental phase sequence of the composition-dependent sphere-packing structures. Interestingly, the domain sizes of three-dimensional (3D) sphere-packing structures are nearly temperature-invariant compared to that of a two-dimensional HEX structure, and we attribute this dimensional stability to the strong enthalpic effects in the formation of 3D packing structures and the geometrical characteristics of spherical domains. Our results demonstrate that the interaction parameter χ, together with the known role of conformational asymmetry ε, serves as the primary thermodynamic parameter for selecting and stabilizing various low-symmetry packing structures in sphere-forming linear diblock copolymer systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

Computer program for analysis of high speed, single row, angular contact, spherical roller bearing, SASHBEAN. Volume 1: User's guide

The computer program SASHBEAN (Sikorsky Aircraft Spherical Roller High Speed Bearing Analysis) analyzes and predicts the operating characteristics of a Single Row, Angular Contact, Spherical Roller Bearing (SRACSRB). The program runs on an IBM or IBM compatible personal computer, and for a given set of input data analyzes the bearing design for it's ring deflections (axial and radial), roller deflections, contact areas and stresses, induced axial thrust, rolling element and cage rotation speeds, lubrication parameters, fatigue lives, and amount of heat generated in the bearing. The dynamic loading of rollers due to centrifugal forces and gyroscopic moments, which becomes quite significant at high speeds, is fully considered in this analysis. For a known application and it's parameters, the program is also capable of performing steady-state and time-transient thermal analyses of the bearing system. The steady-state analysis capability allows the user to estimate the expected steady-state temperature map in and around the bearing under normal operating conditions. On the other hand, the transient analysis feature provides the user a means to simulate the 'lost lubricant' condition and predict a time-temperature history of various critical points in the system. The bearing's 'time-to-failure' estimate may also be made from this (transient) analysis by considering the bearing as failed when a certain temperature limit is reached in the bearing components. The program is fully interactive and allows the user to get started and access most of its features with a minimal of training. For the most part, the program is menu driven, and adequate help messages were provided to guide a new user through various menu options and data input screens. All input data, both for mechanical and thermal analyses, are read through graphical input screens, thereby eliminating any need of a separate text editor/word processor to edit/create data files. Provision is also available to select and view the contents of output files on the monitor screen if no paper printouts are required. A separate volume (Volume-2) of this documentation describes, in detail, the underlying mathematical formulations, assumptions, and solution algorithms of this program.

Aggarwal, Arun K.↗