Search NASA⌕ Search

SEARCH · Search NASA

Results for “functions and mappings”

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.

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

At least 253 records · Page 14

Technology Candidates for Air-to-Air and Air-to-Ground Data Exchange

Technology Candidates for Air-to-Air and Air-to-Ground Data Exchange is a two-year research effort to visualize the U. S. aviation industry at a point 50 years in the future, and to define potential communication solutions to meet those future data exchange needs. The research team, led by XCELAR, was tasked with identifying future National Airspace System (NAS) scenarios, determining requirements and functions (including gaps), investigating technical and business issues for air, ground, & air-to-ground interactions, and reporting on the results. The project was conducted under technical direction from NASA and in collaboration with XCELAR's partner, National Institute of Aerospace, and NASA technical representatives. Parallel efforts were initiated to define the information exchange functional needs of the future NAS, and specific communication link technologies to potentially serve those needs. Those efforts converged with the mapping of each identified future NAS function to potential enabling communication solutions; those solutions were then compared with, and ranked relative to, each other on a technical basis in a structured analysis process. The technical solutions emerging from that process were then assessed from a business case perspective to determine their viability from a real-world adoption and deployment standpoint. The results of that analysis produced a proposed set of future solutions and most promising candidate technologies. Gap analyses were conducted at two points in the process, the first examining technical factors, and the second as part of the business case analysis. In each case, no gaps or unmet needs were identified in applying the solutions evaluated to the requirements identified. The future communication solutions identified in the research comprise both specific link technologies and two enabling technologies that apply to most or all specific links. As a result, the research resulted in a new analysis approach, viewing the underlying architecture of ground-air and air-air communications as a whole, rather than as simple "link to function" paired solutions. For the business case analysis, a number of "reference architectures" were developed for both the future technologies and the current systems, based on three typical configurations of current aircraft. Current and future costs were assigned, and various comparisons made between the current and future architectures. In general, it was assumed that if a future architecture offers lower cost than the current typical architecture, while delivering equivalent or better performance, it is likely that the future solution will gain industry acceptance. Conversely, future architectures presenting higher costs than their current counterparts must present a compelling benefit case in other areas or risk a lack of industry acceptance. The business case analysis consistently indicated lower costs for the proposed future architectures, and in most cases, significantly so. The proposed future solutions were found to offer significantly greater functionality, flexibility, and growth potential over time, at lower cost, than current systems. This was true for overall, fleet-wide equipage for domestic and oceanic air carriers, as well as for single, General Aviation (GA) aircraft. The overall research results indicate that all identified requirements can be met by the proposed solutions with significant capacity for future growth. Results also illustrate that the majority of the future communication needs can be met using currently allocated aviation RF spectrum, if used in more effective ways than it is today. A combination of such optimized aviation-specific links and commercial communication systems meets all identified needs for the 50-year future and beyond, with the caveat that a new, overall function will be needed to manage all information exchange, individual links, security, cost, and other factors. This function was labeled "Delivery Manager" (DM) within this research. DM employs a distributed client/server architecture, for both airborne and ground communications architectures. Final research results included identifying the most promising candidate technologies for the future system, conclusions and recommendations, and identifying areas where further research should be considered.

Communications↗

The LuNaMaps Project: Advancing Capabilities for Developing and Validating Digital Elevation Models of Rocky Surfaces from Orbital Data

Both navigation and surface science can benefit from the ability to generate high resolution and accurate maps of the surface of the Moon and other solar system bodies. The primary way these maps are generated is through the use of orbital imagery and ranging data. Traditionally, the process of using orbital imagery and ranging data is tedious and labor-intensive. Additionally, once maps have been built, there has generally been limited effort in developing standards by which to verify the accuracy and quality of the generated maps. The Lunar Navigation Maps (LuNaMaps) project is a NASA Game Changing Development (GCD) project which over the last 4 years has aimed to address these issues both for the Moon and for other rocky solar system bodies. This has been accomplished through development of new and existing capabilities including: a suite of methods and tools to combine all sources of orbital imagery; a benchmark data set as well as basic requirements for high-fidelity simulations of precision landing functions; tools to synthetically enhance map products with lander-scale features for use in the development and testing of hazard detection systems; methods and tools to evaluate the accuracy of developed digital elevation maps (DEMs) and their quality for use in terrain relative navigation scenarios; and tools to realistically render image and lidar data. In this work, we provide an overview of the capabilities developed through LuNaMaps, demonstrating its use for processing existing lunar data, and describing how it can be applied to other use cases. We additionally provide preliminary results showing the application of the developed tools and processes to the generation of elevation maps of the Lunar Surface Proving Grounds (LSPG) lunar analog at Astrobotic’s Mojave testing facility using “orbital imagery” captured by a drone. In this terrestrial demonstration, we have the benefit of being able to compare the results to a ground truth model of the LSPG. We finally describe plans to use the newly created maps in a terrestrial terrain relative navigation demonstration over the LSPG in early 2025.

optical navigation↗

Observations of the CO bulge on Venus and implications for mesospheric winds

Observations of CO at 2.6 mm (115.27 GHz) were made with the Owens Valley Radio Observatory (OVRO) millimeter interferometer in 1986 and 1988, yielding high-quality disk-resolved spectra which were inverted to determine the CO mixing ratio profile from distinct regions on the disk, allowing us to map the distribution of CO in the upper mesosphere of Venus both horizontally and vertically. The 1986 observations were of the morning terminator and were particularly useful in searching for a suspected CO maximum ('bulge') on the nightside. The resulting CO mixing ratio profiles were mapped for various altitudes as functions of latitude and local time, and we report that we have resolved the previously inferred CO bulge. The bulge increases in magnitude from a small day-night variation at 90 km to an extensive nightside peak at 100 km, the upper limit of our observations. The peak bulge-to-dayside ratio approached 20-30 at 100 km in 1986 and may have been as large as 50-100, assuming late-afternoon CO abundances found in 1988 were similar to those in 1986. Three-dimensional mapping shows that in the upper mesosphere the bulge was displaced from local midnight toward the morning equator, centered at 3:30 AM local time. Using the qualitative model of mesospheric circulation on Venus proposed by Clancy and Muhleman (1985, 1991), we explain this shift in terms of strong retrograde zonal winds throughout the mesosphere, matching the directly detected mesospheric circulation (Shah, K., D. O. Muhleman, and G. L. Berge 1991) observed with the same dataset in 1988.

Gurwell, Mark A.↗

Transition operators in acoustic-wave diffraction theory. I - General theory. II - Short-wavelength behavior, dominant singularities of Zk0 and Zk0 exp -1

A formal theory of the scattering of time-harmonic acoustic scalar waves from impenetrable, immobile obstacles is established. The time-independent formal scattering theory of nonrelativistic quantum mechanics, in particular the theory of the complete Green's function and the transition (T) operator, provides the model. The quantum-mechanical approach is modified to allow the treatment of acoustic-wave scattering with imposed boundary conditions of impedance type on the surface (delta-Omega) of an impenetrable obstacle. With k0 as the free-space wavenumber of the signal, a simplified expression is obtained for the k0-dependent T operator for a general case of homogeneous impedance boundary conditions for the acoustic wave on delta-Omega. All the nonelementary operators entering the expression for the T operator are formally simple rational algebraic functions of a certain invertible linear radiation impedance operator which maps any sufficiently well-behaved complex-valued function on delta-Omega into another such function on delta-Omega. In the subsequent study, the short-wavelength and the long-wavelength behavior of the radiation impedance operator and its inverse (the 'radiation admittance' operator) as two-point kernels on a smooth delta-Omega are studied for pairs of points that are close together.

Hahne, G. E.↗

Topographically Induced Thermal Effects on Lunar Hydrogen Distributions: Correlated Observations from the LRO LEND and LOLA Instruments

The question of whether water exists on the Moon's surface has long been an enigma to Lunar researchers. Largely, this was due to the thermally extreme lunar surface environment that would seem to preclude any long term maintenance, manufacture, transport or accumulation of hydrogen (H) volatiles over most of the lunar surface. As a result, for many years the cold permanent shadow regions (PSR) in the bottoms of craters near the lunar poles appeared to provide the basic conditions at least for maintenance of lunar hydrogen. Importantly, recent discoveries indicate that there is some hydrogen at the poles. However, the picture of the lunar hydrogen budget may be more complex than the PSR hypothesis has suggested. This evidence comes from observations by the Lunar Exploration Neutron Detector (LEND) onboard the Lunar Reconnaissance Orbiter (LRO) that inclici1te 1) some H concentrations lie outside PSR and 2) though a few of the larger PSR's have high hydrogen, PSR does not appear to be an independent factor influencing the large-scale suppression of polar epithermals observed by LEND and the Lunar Prospector Neutron Spectrometer. In this research we investigate the possibility that the thermal contrast between pole-facing and equator facing-slopes is a factor influencing the surface distributions of lunar H. We perform this bulk correlated observation and study by developing a thermal proxy from slope data of the Lunar Orbiting Laser Altimeter (LOLA) digital elevation model (DEM) which is registered with the collimated LEND epithermal map. From the LOLA transforms we impose a thermal functional decomposition and systematic statistical analysis of the LEND epithermal map. Our hypothesis testing suggests in most high latitude bands studied> +/- 45 deg: Epithermal rates in pole-facing slopes are significantly lower than epithermal rates in equivalent equator-facing slopes. As a control study, we find that there is no statistically significant difference between equivalent east and west facing slopes. This finding suggests topographic modulation of insolation is a factor influencing the lunar H budget. Importantly, this result is consistent with observations in terrestrial, Martian research.

McClanahan, T. P.↗

Modified Displacement Transfer Functions for Deformed Shape Predictions of Slender Curved Structures with Varying Curvatives

To eliminate the need to use finite-element modeling for structure shape predictions, a new method was invented. This method is to use the Displacement Transfer Functions to transform the measured surface strains into deflections for mapping out overall structural deformed shapes. The Displacement Transfer Functions are expressed in terms of rectilinearly distributed surface strains, and contain no material properties. This report is to apply the patented method to the shape predictions of non-symmetrically loaded slender curved structures with different curvatures up to a full circle. Because the measured surface strains are not available, finite-element analysis had to be used to analytically generate the surface strains. Previously formulated straight-beam Displacement Transfer Functions were modified by introducing the curvature-effect correction terms. Through single-point or dual-point collocations with finite-elementgenerated deflection curves, functional forms of the curvature-effect correction terms were empirically established. The resulting modified Displacement Transfer Functions can then provide quite accurate shape predictions. Also, the uniform straight-beam Displacement Transfer Function was applied to the shape predictions of a section-cut of a generic capsule (GC) outer curved sandwich wall. The resulting GC shape predictions are quite accurate in partial regions where the radius of curvature does not change sharply.

displacement transfer functions↗

Postural tachycardia syndrome: time frequency mapping

Orthostatic tachycardia is common but its specificity remains uncertain. Our preliminary work suggested that using autonomic function testing in conjunction with time-frequency mapping (TFM), it might be possible to characterize a subset of the postural tachycardia syndrome (POTS), that is due to a restricted autonomic neuropathy. We describe 20 patients (17 women and 3 men, aged 14-43 years) with florid POTS and 20 controls (14 women and 6 men, aged 20-41 years). Autonomic failure was quantified by its distribution (cardiovagal, adrenergic and sudomotor) and severity, a symptom profile was generated, and spectral indices, based on modified Wigner distribution during rest and head-up tilt (80 degrees) were evaluated. During tilt-up POTS patients differed from controls by an excessive heart rate (> 130 bpm) (P < 0.001), and higher diastolic pressure (P < 0.01). During rest, cardiovagal oscillations (at respiratory frequencies [RF]) and slow rhythms at nonrespiratory frequencies (NONRF) (from 0.01 to 0.07 Hz) in R-R intervals (RRI) (P < 0.01) were reduced. Both RF and NONRF rhythms in RRI were further blunted with tilt-up (P < 0.001). Slow adrenergic vasomotor rhythms in blood pressure (BP) (approximately 0.07 Hz) surged with tilt-up and returned to normal levels afterwards. The index of sympatho-vagal balance (NONRF-Systolic BP (SBP)/RF-RRI) was dramatically increased in POTS (P < 0.001). Distal postganglionic sudomotor failure was observed, and impairment of the BP responses to the Valsalva maneuver (phase II) suggested peripheral adrenergic dysfunction. Persistent orthostatic dizziness, tiredness, gastrointestinal symptoms and palpitations were common in POTS patients. It is possible to identify a subset of POTS patients who have a length-dependent autonomic neuropathy, affecting the peripheral adrenergic and cardiovagal fibers, with relative preservation of cardiac adrenergic fibers.

NASA Discipline Number 14-10↗

Wavespace-Based Coherent Deconvolution

Array deconvolution is commonly used in aeroacoustic analysis to remove the influence of a microphone array's point spread function from a conventional beamforming map. Unfortunately, the majority of deconvolution algorithms assume that the acoustic sources in a measurement are incoherent, which can be problematic for some aeroacoustic phenomena with coherent, spatially-distributed characteristics. While several algorithms have been proposed to handle coherent sources, some are computationally intractable for many problems while others require restrictive assumptions about the source field. Newer generalized inverse techniques hold promise, but are still under investigation for general use. An alternate coherent deconvolution method is proposed based on a wavespace transformation of the array data. Wavespace analysis offers advantages over curved-wave array processing, such as providing an explicit shift-invariance in the convolution of the array sampling function with the acoustic wave field. However, usage of the wavespace transformation assumes the acoustic wave field is accurately approximated as a superposition of plane wave fields, regardless of true wavefront curvature. The wavespace technique leverages Fourier transforms to quickly evaluate a shift-invariant convolution. The method is derived for and applied to ideal incoherent and coherent plane wave fields to demonstrate its ability to determine magnitude and relative phase of multiple coherent sources. Multi-scale processing is explored as a means of accelerating solution convergence. A case with a spherical wave front is evaluated. Finally, a trailing edge noise experiment case is considered. Results show the method successfully deconvolves incoherent, partially-coherent, and coherent plane wave fields to a degree necessary for quantitative evaluation. Curved wave front cases warrant further investigation. A potential extension to nearfield beamforming is proposed.

Bahr, Christopher J.↗

LET spectra of fragmented 600-MeV/A iron beam on aluminum and polyethylene: A preliminary analysis

An iron beam experiment recently conducted at the Lawrence Berkeley Laboratory Bevalac by Benton et al. provides an opportunity for verifying the new Green's function computer code (GRNTRN) and assessing the related nuclear database. The iron beam with 600 MeV/A at extraction traversed a series of beam transport elements, lead foil, and several triggering devices before impacting the target. Of these, only the 2.24 g/sq cm lead foil and target are considered in the transport analysis with an assumed (inferred from calibration) beam energy of 557 MeV/A. A thin layer of CR-39 plastic nuclear track detector (PNTD) was placed in front of the target to monitor the incident flux and a stack of four PNTD's placed behind to measure the linear energy transfer (LET) distribution of the transported beam. Test data are analyzed for three separate targets: 2 g/sq cm aluminum, 5-cm polyethylene, and 8-cm polyethylene. The two-layer GRNTRN results were mapped into the detector response function for comparison with the measured LET spectra. Reasonable agreement is obtained. Future research and analysis can be improved by using a more accurate isotope set or including other important media which significantly alter the beam. The assumption that the fragmentation cross sections are too small for aluminum is inconclusive because the production of fragments by 2 g/sq cm aluminum can be more significantly affected by the fragment contribution originating from the media in front of the target than that currently considered in the analysis.

Shinn, Judy L.↗

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗

Modeling Ring-Opening Polymerization with Machine-Learned Committors

Many important industrial and biological processes depend on ring-opening polymerization (ROP), including large-scale production polymers, such as nylon 6, and a variety of biodegradable polymers and high-performance resins. Large-scale simulations (>200K atoms) of ROP were performed using a deep neural network-based method that captured the effect of the local chemical environment of the reactive sites using training data from quantum chemical methods. The utility of this method was demonstrated by modeling the cationic ROP of bis-benzoxazine, a promising resin for use in advanced high-temperature application composites (>200°C). Density functional theory calculations were used to map out the various ring opening and polymerization reaction paths that occur during the thermal processing of bis-benzoxazine. The probabilities of these reactions, as a function of local atomic configurations, were used to construct committor functions using a neural network approach. The learned committor functions were then used to determine when a particular reaction occurred within the REACTER framework for reactive molecular dynamics simulations. The experimentally observed exotherm that occurs during bis-benzoxazine ROP was captured in an effective manner in the simulations by imposing a local heating of the system that corresponds to the enthalpy of reaction for each reaction. The dependence of the final polymer morphology on the ROP reaction temperature was analyzed at the molecular scale. The results presented in this work indicate that machine-learned committor functions are a promising approach for incorporating high-fidelity reactivity criteria that capture the effect of local chemical environments into large-scale polymerization simulations.

machine learning, polymer simulations, molecular d↗

Neural Development Under Conditions of Spaceflight

One of the key tasks the developing brain must learn is how to navigate within the environment. This skill depends on the brain's ability to establish memories of places and things in the environment so that it can form cognitive maps. Earth's gravity defines the plane of orientation of the spatial environment in which animals navigate, and cognitive maps are based on this plane of orientation. Given that experience during early development plays a key role in the development of other aspects of brain function, experience in a gravitational environment is likely to be essential for the proper organization of brain regions mediating learning and memory of spatial information. Since the hippocampus is the brain region responsible for cognitive mapping abilities, this study evaluated the development of hippocampal structure and function in rats that spent part of their early development in microgravity. Litters of male and female Sprague-Dawley rats were launched into space aboard the Space Shuttle Columbia on either postnatal day eight (P8) or 14 (P14) and remained in space for 16 days. Upon return to Earth, the rats were tested for their ability to remember spatial information and navigate using a variety of tests (the Morris water maze, a modified radial arm maze, and an open field apparatus). These rats were then tested physiologically to determine whether they exhibited normal synaptic plasticity in the hippocampus. In a separate group of rats (flight and controls), the hippocampus was analyzed using anatomical, molecular biological, and biochemical techniques immediately postlanding. There were remarkably few differences between the flight groups and their Earth-bound controls in either the navigation and spatial memory tasks or activity-induced synaptic plasticity. Microscopic and immunocytochemical analyses of the brain also did not reveal differences between flight animals and ground-based controls. These data suggest that, within the developmental window studied, microgravity has minimal long-term impact on cognitive mapping function and cellular substrates important for this function. Any differences due to development in microgravity were transient and returned to normal soon after return to Earth.

Kosik, Kenneth S.↗

Modern Era Retrospective Restrospective-Analysis for Research and Applications (MERRA) Data and Services at the GES DISC

The Modern Era Retrospective-analysis for Research and Applications (MERRA) dataset is a NASA satellite era, 30 year (1979 - present), reanalysis using the Goddard Earth Observing System Data Assimilation System, Version 5 (GEOS-5). The project, run out of NASA's Global Modeling and Assimilation Office at Goddard Space Flight Center, provides the science and application communities with a state-of-the-art global analysis with emphasis on improved estimates of the hydrological cycle over a broad range of weather and climate time scales. MERRA products are generated as a long-term synthesis that places the NASA EOS suite of observations in a climate context. The MERRA analysis is performed at a horizontal resolution of 2/3 longitude x 1/2 latitude (540x361 global gridpoints) with observational analyses every 6 hours. The MERRA output data will include 3 dimensional state fields for every 6 hourly analysis cycle on 42 pressure levels (or 72 terrain following model coordinate levels) from the surface through the stratosphere. Several data products are specifically designed to support chemistry and stratosphere transport modeling. The 2 dimensional surface and atmospheric diagnostics (numbering 259) are being stored on the native grid at 1 hourly intervals. These include radiation and vertical integrals of the atmosphere for water and energy budget studies and also surface diagnostics where the diurnal cycle is important. The one hourly surface and near surface data product will also facilitate research on the integrated analysis of Earth system observations in the land, ocean and cryosphere. The MERRA products are archived and distributed by the Goddard Earth Sciences Data and Information Services Center (GES DISC) through its Modeling DISC Web (MDISC) portal. Multiple data access methods and services are available for MERRA data through MDISC: (1) Mirador offers a quick, comprehensive search of MERRA and all GES DISC archived data holdings, allowing searches on keywords, location names or latitude/longitude box, and date/time, with responses within a few seconds. (2) Giovanni is a GES DISC developed Web application that provides data visualization and analysis online. Giovanni features popular visualizations such as latitude-longitude maps, animations, cross sections, profiles, time series, etc. and some basic statistical analysis functions such as scatter plots and correlation coefficient maps. Users are able to download results in several different formats, including Google Earth. (3) On-the-fly parameter subsetting of data within a spatial/temporal window is provided through a simple select and click Web page. (4) MERRA data are also available via OPeNDAP, GrADS Data Server (GDS) and can be converted to netCDF on the fly.

Berrick, Stephen W.↗

Photogrammetry and altimetry. Part A: Apollo 16 laser altimeter

The laser altimeter measures precise altitudes of the command and service module above the lunar surface and can function either with the metric (mapping) camera or independently. In the camera mode, the laser altimeter ranges at each exposure time, which varies between 20 and 28 sec (i.e., 30 to 43 km on the lunar surface). In the independent mode, the laser altimeter ranges every 20 sec. These altitude data and the spacecraft attitudes that are derived from simultaneous stellar photography are used to constrain the photogrammetric reduction of the lunar surface photographs when cartographic products are generated. In addition, the altimeter measurements alone provide broad-scale topographic relief around the entire circumference of the moon. These data are useful in investigating the selenodetic figure of the moon and may provide information regarding gravitational anomalies on the lunar far side.

Wollenhaupt, W. R.↗

Underground coal mining section data

A set of tables which display the allocation of time for ten personnel and eight pieces of underground coal mining equipment to ten function categories is provided. Data from 125 full shift time studies contained in the KETRON database was utilized as the primary source data. The KETRON activity and delay codes were mapped onto JPL equipment, personnel and function categories. Computer processing was then performed to aggregate the shift level data and generate the matrices. Additional, documented time study data were analyzed and used to supplement the KETRON databased. The source data including the number of shifts are described. Specific parameters of the mines from which there data were extracted are presented. The result of the data processing including the required JPL matrices is presented. A brief comparison with a time study analysis of continuous mining systems is presented. The procedures used for processing the source data are described.

Gabrill, C. P.↗

Barriers and dispersal surfaces in minimum-time interception

Minimum time interception of a target moving in a horizontal plane is analyzed as a one-player differential game. Dispersal points and points on the barrier are located for a class of pursuit evasion and interception problems. These points are determined by constructing cross sections of the isochrones and hence obtaining the barrier, dispersal, and control level surfaces. The game solution maps the controls as a function of the state within the capture region.

Rajan, N.↗

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

The paper describes a computer code for calculating the flow dynamics of a constant-density flow in the second-stage trumpet shaped nozzle section of a two-stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite-difference computer code has been developed for the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for the efficient solution of the equations. The code produces as output the flow field map of the nondimensional stream function, axial and swirl velocity. It was found that the best location for seed injection to obtain a uniform distribution at the combustor exit is in the central location for seed injected at the entrance to the second stage combustor.

Gupta, A. K.↗