Search NASA⌕ Search

SEARCH · Search NASA

Results for “geometrical sum”

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.

49 records · Page 3

Bounding solutions of geometrically nonlinear viscoelastic problems

Integral transform techniques, such as the Laplace transform, provide simple and direct methods for solving viscoelastic problems formulated within a context of linear material response and using linear measures for deformation. Application of the transform operator reduces the governing linear integro-differential equations to a set of algebraic relations between the transforms of the unknown functions, the viscoelastic operators, and the initial and boundary conditions. Inversion either directly or through the use of the appropriate convolution theorem, provides the time domain response once the unknown functions have been expressed in terms of sums, products or ratios of known transforms. When exact inversion is not possible approximate techniques may provide accurate results. The overall problem becomes substantially more complex when nonlinear effects must be included. Situations where a linear material constitutive law can still be productively employed but where the magnitude of the resulting time dependent deformations warrants the use of a nonlinear kinematic analysis are considered. The governing equations will be nonlinear integro-differential equations for this class of problems. Thus traditional as well as approximate techniques, such as cited above, cannot be employed since the transform of a nonlinear function is not explicitly expressible.

Stubstad, J. M.↗

Investigation of the external flow analysis for density measurements at high altitude

The results of analysis performed on the external flow around the shuttle orbiter nose regions at the Shuttle Upper Atmosphere Mass Spectrometer (SUMS) inlet orifice are presented. The purpose of the analysis is to quantitatively characterize the flow conditions to facilitate SUMS flight data reduction and subsequent determination of orbiter aerodynamic force coefficients in the hypersonic rarefied flow regime. Experimental determination of aerodynamic force coefficients requires accurate simultaneous measurement of forces (or acceleration) and dynamic pressure along with independent knowledge of density and velocity. The SUMS provides independent measurement of dynamic pressure; however, it does so indirectly and requires knowledge of the relationship between measured orifice conditions and the dynamic pressure which can only be determined on the basis of molecule or theory for a winged configuration. Monte Carlo direct simulation computer codes were developed for both the flow field solution at the orifice and for the internal orifice flow. These codes were used to study issues associated with geometric modeling of the orbiter nose geometry and the modeling of intermolecular collisions including rotational energy exchange and a preliminary analysis of vibrational excitation and dissociation effects. Data obtained from preliminary simulation runs are presented.

Bienkowski, G. K.↗

A uniform geometrical theory of diffraction for an edge in a perfectly conducting surface

A compact dyadic diffraction coefficient for electromagnetic waves obliquely incident on a curved edge formed by perfectly conducting curved or plane surfaces is obtained. This diffraction coefficient remains valid in the transition regions adjacent to shadow and reflection boundaries, where the diffraction coefficients of Keller's original theory fail. Our method is based on Keller's method of the canonical problem, which in this case is the perfectly conducting wedge illuminated by plane, cylindrical, conical, and spherical waves. When the proper ray-fixed coordinate system is introduced, the dyadic diffraction coefficient for the wedge is found to be the sum of only two dyads, and it is shown that this is also true for the dyadic diffraction coefficients of higher order edges. One dyad contains the acoustic soft diffraction coefficient; the other dyad contains the acoustic hard diffraction coefficient. The expressions for the acoustic wedge diffraction coefficients contain Fresnel integrals, which ensure that the total field is continuous at shadow and reflection boundaries.

Kouyoumjian, R. G.↗

Applicability of the single equivalent point dipole model to represent a spatially distributed bio-electrical source

Although the single equivalent point dipole model has been used to represent well-localised bio-electrical sources, in realistic situations the source is distributed. Consequently, position estimates of point dipoles determined by inverse algorithms suffer from systematic error due to the non-exact applicability of the inverse model. In realistic situations, this systematic error cannot be avoided, a limitation that is independent of the complexity of the torso model used. This study quantitatively investigates the intrinsic limitations in the assignment of a location to the equivalent dipole due to distributed electrical source. To simulate arrhythmic activity in the heart, a model of a wave of depolarisation spreading from a focal source over the surface of a spherical shell is used. The activity is represented by a sequence of concentric belt sources (obtained by slicing the shell with a sequence of parallel plane pairs), with constant dipole moment per unit length (circumferentially) directed parallel to the propagation direction. The distributed source is represented by N dipoles at equal arc lengths along the belt. The sum of the dipole potentials is calculated at predefined electrode locations. The inverse problem involves finding a single equivalent point dipole that best reproduces the electrode potentials due to the distributed source. The inverse problem is implemented by minimising the chi2 per degree of freedom. It is found that the trajectory traced by the equivalent dipole is sensitive to the location of the spherical shell relative to the fixed electrodes. It is shown that this trajectory does not coincide with the sequence of geometrical centres of the consecutive belt sources. For distributed sources within a bounded spherical medium, displaced from the sphere's centre by 40% of the sphere's radius, it is found that the error in the equivalent dipole location varies from 3 to 20% for sources with size between 5 and 50% of the sphere's radius. Finally, a method is devised to obtain the size of the distributed source during the cardiac cycle.

NASA Discipline Cardiopulmonary↗

Light aircraft lift, drag, and moment prediction: A review and analysis

The historical development of analytical methods for predicting the lift, drag, and pitching moment of complete light aircraft configurations in cruising flight is reviewed. Theoretical methods, based in part on techniques described in the literature and in part on original work, are developed. These methods form the basis for understanding the computer programs given to: (1) compute the lift, drag, and moment of conventional airfoils, (2) extend these two-dimensional characteristics to three dimensions for moderate-to-high aspect ratio unswept wings, (3) plot complete configurations, (4) convert the fuselage geometric data to the correct input format, (5) compute the fuselage lift and drag, (6) compute the lift and moment of symmetrical airfoils to M = 1.0 by a simplified semi-empirical procedure, and (7) compute, in closed form, the pressure distribution over a prolate spheroid at alpha = 0. Comparisons of the predictions with experiment indicate excellent lift and drag agreement for conventional airfoils and wings. Limited comparisons of body-alone drag characteristics yield reasonable agreement. Also included are discussions for interference effects and techniques for summing the results above to obtain predictions for complete configurations.

Smetana, F. O.↗

Calculation of Heavy Ion Inactivation and Mutation Rates in Radial Dose Model of Track Structure

In the track structure model, the inactivation cross section is found by summing an inactivation probability over all impact parameters from the ion to the sensitive sites within the cell nucleus. The inactivation probability is evaluated by using the dose response of the system to gamma rays and the radial dose of the ions and may be equal to unity at small impact parameters. We apply the track structure model to recent data with heavy ion beams irradiating biological samples of E. Coli, B. Subtilis spores, and Chinese hamster (V79) cells. Heavy ions have observed cross sections for inactivation that approach and sometimes exceed the geometric size of the cell nucleus. We show how the effects of inactivation may be taken into account in the evaluation of the mutation cross sections in the track structure model through correlation of sites for gene mutation and cell inactivation. The model is fit to available data for HPRT (hypoxanthine guanine phosphoribosyl transferase) mutations in V79 cells, and good agreement is found. Calculations show the high probability for mutation by relativistic ions due to the radial extension of ions track from delta rays. The effects of inactivation on mutation rates make it very unlikely that a single parameter such as LET (linear energy transfer) can be used to specify radiation quality for heavy ion bombardment.

Cucinotta, Francis A.↗

Computer prediction of large reflector antenna radiation properties

A FORTRAN program for calculating reflector antenna radiation patterns was rewritten and extended to include reflectors composed of a number of panels. These individual panels must be analytic surfaces. The theoretical foundation for the program is as follows: Geometrical optics techniques are used to trace rays from a feed antenna to the reflector surface and back to a mathematical plane just in front of the reflector. The resulting tangential electric field distribution, assumed to be the only source of forward radiation, is integrated numerically to calculate the radiation pattern for a desired set of angles. When the reflector is composed of more than one panel, each panel is treated as a separated antenna, the ray-tracing procedure and integration being repeated for each panel. The results of the individual aperture plane integrations are stored and summed to yield the relative electric field strength over the angles of interest. An example and several test cases are included to demonstrate the use of the program and verify the new method of computation.

Botula, A.↗

Using Correlated Photons to Suppress Background Noise

A proposed method of suppressing the effect of background noise in an optical communication system would exploit the transmission and reception of correlated photons at the receiver. The method would not afford any advantage in a system in which performance is limited by shot noise. However, if the performance of the system is limited by background noise (e.g., sunlight in the case of a free-space optical communication system or incoherently scattered in-band photons in the case of a fiber-optic communication system), then the proposed method could offer an advantage: the proposed method would make it possible to achieve a signal-to-noise ratio (S/N) significantly greater than that of an otherwise equivalent background- noise-limited optical communication system based on the classical transmission and reception of uncorrelated photons. The figure schematically depicts a classical optical-communication system and a system according to the proposed method. In the classical system, a modulated laser beam is transmitted along an optical path to a receiver, the optics of which include a narrow-band-pass filter that suppresses some of the background noise. A photodetector in the receiver detects the laser-beam and background photons, most or all of which are uncorrelated. In the proposed system, correlated photons would be generated at the transmitter by making a modulated laser beam pass through a nonlinear parametric down-conversion crystal. The sum of frequencies of the correlated photons in each pair would equal the frequency of the incident photon from which they were generated. As in the classical system, the correlated photons would travel along an optical path to a receiver, where they would be band-pass filtered and detected. Unlike in the classical system, the photodetector in the receiver in this system would be one that intrinsically favors the detection of pairs of correlated photons over the detection of uncorrelated photons. Even though there would be no way of knowing the precise location and time of creation of a given pair of correlated signal photons in the nonlinear down-conversion crystal, the fact that the photons are necessarily created at the same time and place makes it possible to utilize conventional geometrical imaging optics to reunite the photons in coincidence in the receiving photodetector. Because most or all of the signal photons would be correlated while most or all of the noise photons would be uncorrelated, the S/N would be correspondingly enhanced in the photodetector output. An additional advantage to be gained by use of a correlated-photon detector is that it could be capable of recovering the signal even in the presence of background light so bright that a classical uncorrelated-photon detector would be saturated. A blocked-impurity-band (BIB) photodetector that preferentially detects pairs of correlated photons over uncorrelated ones and that operates at a quantum efficiency of 88 percent is commercially available. This detector must be cooled to the temperature of liquid helium to obtain the desired low-noise performance. It is planned to use this detector in a proof-of-principle demonstration. In addition, it may be possible to develop GaN-based photodetectors that could offer the desired low-noise performance at room temperature.

Jackson, Deborah↗

System Finds Horizontal Location of Center of Gravity

An instrumentation system rapidly and repeatedly determines the horizontal location of the center of gravity of a laboratory vehicle that slides horizontally on three air bearings (see Figure 1). Typically, knowledge of the horizontal center-of-mass location of such a vehicle is needed in order to balance the vehicle properly for an experiment and/or to assess the dynamic behavior of the vehicle. The system includes a load cell above each air bearing, electronic circuits that generate digital readings of the weight on each load cell, and a computer equipped with software that processes the readings. The total weight and, hence, the mass of the vehicle are computed from the sum of the load-cell weight readings. Then the horizontal position of the center of gravity is calculated straightforwardly as the weighted sum of the known position vectors of the air bearings, the contribution of each bearing being proportional to the weight on that bearing. In the initial application for which this system was devised, the center- of-mass calculation is particularly simple because the air bearings are located at corners of an equilateral triangle. However, the system is not restricted to this simple geometry. The system acquires and processes weight readings at a rate of 800 Hz for each load cell. The total weight and the horizontal location of the center of gravity are updated at a rate of 800/3 approx. equals 267 Hz. In a typical application, a technician would use the center-of-mass output of this instrumentation system as a guide to the manual placement of small weights on the vehicle to shift the center of gravity to a desired horizontal position. Usually, the desired horizontal position is that of the geometric center. Alternatively, this instrumentation system could be used to provide position feedback for a control system that would cause weights to be shifted automatically (see Figure 2) in an effort to keep the center of gravity at the geometric center.

Johnston, Albert S.↗

Contact and pressure balance structures in two-fluid cosmic-ray hydrodynamics

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

Webb, G. M.↗

Model of Image Artifacts from Dust Particles

A mathematical model of image artifacts produced by dust particles on lenses has been derived. Machine-vision systems often have to work with camera lenses that become dusty during use. Dust particles on the front surface of a lens produce image artifacts that can potentially affect the performance of a machine-vision algorithm. The present model satisfies a need for a means of synthesizing dust image artifacts for testing machine-vision algorithms for robustness (or the lack thereof) in the presence of dust on lenses. A dust particle can absorb light or scatter light out of some pixels, thereby giving rise to a dark dust artifact. It can also scatter light into other pixels, thereby giving rise to a bright dust artifact. For the sake of simplicity, this model deals only with dark dust artifacts. The model effectively represents dark dust artifacts as an attenuation image consisting of an array of diffuse darkened spots centered at image locations corresponding to the locations of dust particles. The dust artifacts are computationally incorporated into a given test image by simply multiplying the brightness value of each pixel by a transmission factor that incorporates the factor of attenuation, by dust particles, of the light incident on that pixel. With respect to computation of the attenuation and transmission factors, the model is based on a first-order geometric (ray)-optics treatment of the shadows cast by dust particles on the image detector. In this model, the light collected by a pixel is deemed to be confined to a pair of cones defined by the location of the pixel s image in object space, the entrance pupil of the lens, and the location of the pixel in the image plane (see Figure 1). For simplicity, it is assumed that the size of a dust particle is somewhat less than the diameter, at the front surface of the lens, of any collection cone containing all or part of that dust particle. Under this assumption, the shape of any individual dust particle artifact is the shape (typically, circular) of the aperture, and the contribution of the particle to the attenuation factor for a given pixel is the fraction of the cross-sectional area of the collection cone occupied by the particle. Assuming that dust particles do not overlap, the net transmission factor for a given pixel is calculated as one minus the sum of attenuation factors contributed by all dust particles affecting that pixel. In a test, the model was used to synthesize attenuation images for random distributions of dust particles on the front surface of a lens at various relative aperture (F-number) settings. As shown in Figure 2, the attenuation images resembled dust artifacts in real test images recorded while the lens was aimed at a white target.

Willson, Reg↗

Magic-T Junction using Microstrip/Slotline Transitions

An improved broadband planar magic-T junction that incorporates microstrip/slotline transitions has been developed. In comparison with a prior broadband magic-T junction incorporating microstrip/slotline transitions, this junction offers superior broadband performance. In addition, because this junction is geometrically simpler and its performance is less affected by fabrication tolerances, the benefits of the improved design can be realized at lower fabrication cost. There are potential uses for junctions like this one in commercial microwave communication receivers, radar and polarimeter systems, and industrial microwave instrumentation. A magic-T junction is a four-port waveguide junction consisting of a combination of an H-type and an E-type junction. An E-type junction is so named because it includes a junction arm that extends from a main waveguide in the same direction as that of the electric (E) field in the waveguide. An H-type junction is so named because it includes a junction arm parallel to the magnetic (H) field in a main waveguide. A magic-T junction includes two input ports (here labeled 1 and 2, respectively) and two output ports (here labeled E and H, respectively). In an ideal case, (1) a magic-T junction is lossless, (2) the input signals add (that is, they combine in phase with each other) at port H, and (3) the input signals subtract (that is, they combine in opposite phase) at port E. The prior junction over which the present junction is an improvement affords in-phase-combining characterized by a broadband frequency response, and features a small slotline area to minimize in-band loss. However, with respect to isolation between ports 1 and 2 and return loss at port E, it exhibits narrowband frequency responses. In addition, its performance is sensitive to misalignment of microstrip and slotline components: this sensitivity is attributable to a limited number of quarter-wavelength (lambda/4) transmission-line sections for matching impedances among all four ports, and to strong parasitic couplings at the microstrip/slotline T junction, where four microstrip lines and a slotline are combined. The present improved broadband magic-T junction (see figure) includes a microstrip ring structure and two microstrip- to-slotline transitions. One of the microstrip/slotline transitions is a small T junction between the ring and a slotline; the other microstrip/slotline transition effects coupling between the slotline and port E. The smallness of the T junction and the use of minimum-size slotline terminations help to minimize radiation loss. An impedance-transformation network that includes multiple quarter-wavelength sections is used to increase the operating bandwidth and minimize the parasitic coupling around the microstrip/slotline T junction. As a result, the improved junction has greater bandwidth and lower phase imbalance at the sum and difference ports than did the prior junction.

U-yen, Kongpop↗

Combined Experimental and Modeling Study of the Interactions of Acid Gas with Common Spacecraft Surfaces for Fire Safety Applications

A fire in a spacecraft poses detrimental consequences and risks mission success in addition to crew safety. This is compounded during long-duration missions when the crew has limited options to recover from a fire. A common spacecraft fire concern is the smoldering of wire insulation, typically made from Polyvinyl chloride (PVC) or Polytetrafluoroethylene (PTFE). This creates acid gases such as Hydrogen Chloride (HCl), Hydrogen Fluoride (HF) and Hydrogen Cyanide (HCN). These poisonous gases are hazardous to the crew. They also interact with common surfaces within the spacecraft more than dominant combustion products such as CO2 and H2O. This makes them more difficult to track for potential fire detection techniques, or for postfire clean-up. It is imperative to be able to understand and predict the fate of these poisonous species in a microgravity environment in order to design a safe vehicle. HCl interacts with a number of materials inside a spacecraft. Primary among these materials is aluminum, which is abundantly used due to its strong and light weight nature. Aluminum has a natural oxide layer that protects it from corrosion but is typically treated to enhance this oxide layer. Among these treatments is a chromate conversion coating (CCC), which provides a thin enough protective oxide layer to still conduct electricity, and a traditional anodized material that has a thicker oxide layer that does not conduct electricity. Nomex is another common material found inside a spacecraft. It is a flame-resistant woven polymer that is related to nylon. This commercially available material is used for cargo storage bags and as a fire barrier. Physics-based models were developed to predict the uptake of HCl by these materials. The ultimate objective of these models is to predict the fate of HCl within the spacecraft so that sensors can be placed in meaningful locations in future missions based on the model predictions. To support these modeling efforts, experiments were performed in a cast acrylic test cell that measured the difference between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. Different uptake capacities were realized for each type of sample tested. A computational fluid dynamics model (CFD) model of the reactor was then constructed that used a one-step global reaction rate with calibratable reaction (or kinetic) constants. These constants were calibrated to match the HCl uptake on the CCC aluminum samples, and the same kinetic constants were then tested for the stock and anodized aluminum samples. Model predictions matched the experimental data for the stock aluminum, and to a much lesser extent, the anodized aluminum. The model was additionally validated at different flow rates, sample surface areas, and inlet concentrations, and showed good agreement for all stock and CCC samples. The model did not accurately predict the HCl uptake in the anodized samples compared to the other two types of aluminum. Adjusting the kinetic constants and transport properties did little to improve the prediction. X-Ray Photoelectron Spectroscopy (XPS) was used to determine that the oxide layer thickness of anodized aluminum is approximately 5,000 nm, compared to 250 nm for CCC and 50 nm for stock. XPS also revealed presence of chlorine further down in the aluminum oxide layer in anodized samples than CCC and stock samples after the samples were saturated with HCl, indicating that accounting for diffusion of HCl into the oxide layer is important for accurate prediction of HCl uptake onto anodized aluminum. Consequently, a multi-scale model was developed and tested. First, a single pore inside the anodized aluminum oxide layer was modeled and is referred to as the pore-scale model. In this model, HCl diffused through the pore and reacted with the aluminum oxide pore wall to create aluminum chloride. The sample was then saturated when the mass transfer resistance through the growing aluminum chloride layer became too large for the HCl to reach the aluminum oxide wall and continue the reaction. This pore-scale model was coupled to the reactor-scale model using a concentration-dependent diffusion coefficient, resulting in much more accurate predictions (approximately half the sum square error of the aforementioned reactor-scale model that produced good agreement for stock and CCC) for a variety of operating conditions. The amount of water vapor or relative humidity (RH) in the flow during a reactor experiment was determined to influence HCl uptake. Experiments were performed to understand the interaction of gaseous HCl with aluminum surfaces in the presence of water vapor. The results show that increasing levels of RH increased the capacity of aluminum to adsorb HCl but decreased the capacity of Nomex to uptake HCl. A series of tests were performed on individual aluminum samples after they had been saturated with a fixed concentration of HCl in dry air conditions with the goal of determining how their HCl uptake capacity changes after various treatments with water relative to the original saturation tests. HCl-saturated aluminum samples subjected to a second dry air flow at the same HCl concentration as the original test had an uptake of 23.5% of the original sample with no treatment in between. Saturated aluminum samples subjected to an air flow with a RH of 90% in between tests had an uptake of 35.6% of the original. Saturated aluminum samples submerged in distilled water for 12 hours in between tests had an uptake of 82.2% of the original sample. Previously saturated aluminum tested with HCl and a 50% RH air flow resulted in similar uptake characteristics in multiple repeated tests. The experiments show the profound effect water vapor has on HCl uptake onto aluminum surfaces. In the samples subjected to water vapor or liquid water, capillary condensation and capillary diffusion alters the transport of HCl significantly. A model was proposed that developed a relationship between RH and the coefficient of HCl diffusion in aluminum chloride. This produced an “S-shaped” curve with diffusion coefficient as a function of RH, with 45% RH represented as the point where the diffusion coefficient is halfway between no water saturation and 100% water saturation in the aluminum chloride product layer. No difference in uptake characteristics for the experiment or model were realized between 50% and 62% RH. The results from the large-scale microgravity experiment, Saffire, are discussed as they pertain to the fate of HCl throughout a spacecraft. HCl was released, both as a standalone event, and in concurrence with the burning of a structured cloth. These events only produced a small response in the far field HCl sensor, while a PMMA burn that did not produce HCl had a significantly greater response. A ground-based large-scale facility was constructed to flow acid gas at the scale and configuration realized in the Saffire experiments. A CFD model of this duct was constructed to test kinetic parameters developed in this work at a larger scale and different geometric configuration and to predict the results of the large-scale facility. The models developed in this work were used to interpret the results of the microgravity tests and lead the discussion on what further experiments and models are needed in order to predict the fate of acid gas in a spacecraft environment. To summarize, the major contributions of this work are as follows: the capacity to uptake HCl, with and without the presence of water vapor, was measured for a variety of real spacecraft surfaces. Several different models (single reactor-scale, multiscale, spacecraft-scale) were developed and with the aid of modeling, the rate of uptake for those surfaces was also predicted and validated. The kinetic parameters determined from the small-scale reactor experiments and models were used to predict large-scale and microgravity tests. Conclusions from this research will be used in the design of spacecraft vehicles and large-scale microgravity fire safety experiments. The models built by this work will aid designers in sensor placement and could be used to predict acid gas transport from fires in partial gravity, as would be seen in Lunar and Martian habitats.

fire safety↗