Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable threshold”

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 361 records · Page 20

NASA-STD-6001B Test 1 Upward Flame Propagation; Sample Length Impact on MOC Investigation

Understanding the combustion behavior of materials in the elevated oxygen environments of habitable spacecraft is of utmost importance to crew safety and mission success. Currently, certification for unrestricted flight usage of a material with respect to flammability involves passing the Upward Flame Propagation Test of NASA-STD-6001B (Test 1). This test evaluates materials in a standardized test configuration for two failure criteria: self-extinguishment within 15 cm (6 in.) and the propensity of flame propagation by means of flaming material transfer. By the NASA standard, full-length samples are 30 cm (12 in.) in length; however, factors independent of the test method such as limited material availability or various nonstandard test configurations limit the full pretest sample lengths available for test. This paper characterizes the dependence, if any, of pretest sample length on NASA-STD-6001B Test 1 results. Testing was performed using the Maximum Oxygen Concentration (MOC) Threshold Method to obtain a data set for each sample length tested. In addition, various material types, including cloth (Nomex), foam (TA-301) and solids (Ultem), were tested to investigate potential effects of test specimen types. Though additional data needs to be generated to provide statistical confidence, preliminary findings are that use of variable sample lengths has minimal impact on NASA-STD-6001B flammability performance and MOC determination.

Harper, Susana Tapia↗

Ventilatory responses to dynamic exercise elicited by intramuscular sensors

PURPOSE: Eight subjects, aged 27.0+/-1.6 yr, performed incremental workload cycling to investigate the contribution of skeletal muscle mechano- and metaboreceptors to ventilatory control during dynamic exercise. METHODS: Each subject performed four bouts of exercise: exercise with no intervention (CON); exercise with bilateral thigh cuffs inflated to 90 mm Hg (CUFF); exercise with application of lower-body positive pressure (LBPP) to 45 torr (PP); and exercise with 90 mm Hg thigh cuff inflation and 45 torr LBPP (CUFF+PP). Ventilatory responses and pulmonary gas exchange variables were collected breath-by-breath with concomitant measurement of leg intramuscular pressure. RESULTS: Ventilation (VE) was significantly elevated from CON during PP and CUFF+PP at workloads corresponding to > or = 60% CON peak oxygen uptake (VO2peak) and during CUFF at workloads > or = 80% CON VO2peak, P < 0.05. The VO2 at which ventilatory threshold occurred was significantly reduced from CON (2.17+/-0.28 L x min(-1)) to 1.60+/-0.19 L x min(-1), 1.45+/-0.15 L x min(-1), and 1.15+/-0.11 L x min(-1) during CUFF, PP, and CUFF+PP, respectively. The slope of the linear regression describing the VE/CO2 output relationship was increased from CON by approximately 22% during CUFF, 40% during PP, and 41% during CUFF+PP. CONCLUSIONS: As intramuscular pressure was significantly elevated immediately upon application of LBPP during PP and CUFF+PP without a concomitant increase in VE, it seems unlikely that LBPP-induced increases in VE can be attributed to activation of the mechanoreflex. These findings suggest that LBPP-induced reductions in perfusion pressure and decreases in venous outflow resulting from inflation of bilateral thigh cuffs may generate a metabolite sensitive intramuscular ventilatory stimulus.

Non-NASA Center↗

Variable Shadow Screens for Imaging Optical Devices

Variable shadow screens have been proposed for reducing the apparent brightnesses of very bright light sources relative to other sources within the fields of view of diverse imaging optical devices, including video and film cameras and optical devices for imaging directly into the human eye. In other words, variable shadow screens would increase the effective dynamic ranges of such devices. Traditionally, imaging sensors are protected against excessive brightness by use of dark filters and/or reduction of iris diameters. These traditional means do not increase dynamic range; they reduce the ability to view or image dimmer features of an image because they reduce the brightness of all parts of an image by the same factor. On the other hand, a variable shadow screen would darken only the excessively bright parts of an image. For example, dim objects in a field of view that included the setting Sun or bright headlights could be seen more readily in a picture taken through a variable shadow screen than in a picture of the same scene taken through a dark filter or a narrowed iris. The figure depicts one of many potential variations of the basic concept of the variable shadow screen. The shadow screen would be a normally transparent liquid-crystal matrix placed in front of a focal-plane array of photodetectors in a charge-coupled-device video camera. The shadow screen would be placed far enough from the focal plane so as not to disrupt the focal-plane image to an unacceptable degree, yet close enough so that the out-of-focus shadows cast by the screen would still be effective in darkening the brightest parts of the image. The image detected by the photodetector array itself would be used as feedback to drive the variable shadow screen: The video output of the camera would be processed by suitable analog and/or digital electronic circuitry to generate a negative partial version of the image to be impressed on the shadow screen. The parts of the shadow screen in front of those parts of the image with brightness below a specified threshold would be left transparent; the parts of the shadow screen in front of those parts of the image where the brightness exceeded the threshold would be darkened by an amount that would increase with the excess above the threshold.

Lu, Ed↗

A Study of Precipitation Mean Particle Size with the GPM Combined Radar-Radiometer Algorithm

The Global Precipitation Measurement (GPM) mission provides quantitative estimates of the mass-weighted mean diameter, Dm, a characteristic parameter of the precipitation particle size distribution (PSD). This study investigates the three-dimensional distribution of Dm retrieved by the GPM Combined Radar-Radiometer Algorithm (CORRA) using observations from GPM’s dual-frequency precipitation radar and microwave imager. We examine Dm in individual precipitation systems in different climate regimes and investigate a six-year global climatology. The vertical structure of Dm is demonstrated with cases of deep convection, frontal rain and snow, and stratocumulus light rain. The Dm values, detectable by GPM, range from appx 0.7 mm in stratocumulus precipitation to >3.5 mm in the ice layers of intense convection. Within the constraint of the 12-dBZ detectability threshold, the smallest climatological mean Dm (~ 0.8 mm) are found in the eastern oceans, and the largest values (~ 2 mm) occur above the melting levels in convection over land in summer. The standard deviation is generally < 0.45 mm below 6 km. Climate regimes are characterized with Dm annual/seasonal variations, its convective/stratiform components, and vertical variabilities (2-10 km). The US Central Plains and Argentina are associated with the largest Dm in a deep layer. Tropical Africa has larger Dm and standard deviation than Amazon. Large convective Dm occurs at high latitudes of Eurasia and North America in summer; the southern hemisphere high latitudes have shallower systems with smaller Dm. Oceanic storm tracks in both hemispheres have relatively large Dm, particularly for convective Dm in winter.

Mei Han↗

Influence on Structural Loading of a Wave Energy Converter by Controlling Variable-Geometry Components and the Power Take-Off

Oceans are harsh environments and can impose significant loads on deployed structures. A wave energy converter (WEC) should be designed to maximize the energy absorbed while ensuring the operating wave condition does not exceed the failure limits of the device itself. Therefore, the loads endured by the support structure are a design constraint for the system. Furthermore, the WEC should be adaptable to different sea states. Herethis work uses a WEC-Sim model of a variable-geometry oscillating wave energy converter (VGOSWEC) mounted on a support structure simulated under different wave scenarios. A VGOSWEC resembles a paddle pitching about a fixed hinge perpendicular to the incoming wave fronts. The geometry of the VGOSWEC is varied by opening a series of controllable flaps on the pitching paddle when the structure experiences threshold loads. It is hypothesized that opening the flaps should result in load shedding at the base of the support structure by reducing the moments about the hinge axis. This work compares the hydrodynamic coefficients, natural periods, and response amplitude operators from completely closed to completely open configurations of the controllable flaps. This work shows that the completely open configuration can reduce the pitch and surge loads on the base of the support structure by as much as 80%. Increased loads at the structure’s natural period can be mitigated by an axial power take-off damping acting as an additional design parameter to control the loads at the WEC’s support structure.

16 TIDAL AND WAVE POWER↗

Optical Classification of Lower Amazon Waters Based on In Situ Data and Sentinel-3 Ocean and Land Color Instrument Imagery

Optical water types (OWTs) were identified from an in situ dataset of concomitant biogeochemical and optical parameters acquired in the Amazon River and its tributaries, in the Lower Amazon region, at different hydrological conditions from 2014 to 2017. A seasonal bio-optical characterization was performed. The k-means classification was applied to the in situ normalized reflectance spectra (rn(λ)), allowing the identification of four OWTs. An optical index method was also applied to the rn(λ) defining the thresholds of the OWTs. Next, level-3 Sentinel-3 Ocean and Land Color Instrument images representative of the seasonal discharge conditions were classified using the identified in situ OWTs as reference. The differences between Amazon River and clearwater tributary OWTs were dependent on the hydrological dynamics of the Amazon River, also showing a strong seasonal variability. Each OWT was associated with a specific bio-optical and biogeochemical environment assessed from the corresponding absorption coefficient values of colored dissolved organic matter (aCDOM) and particulate matter (ap), chlorophyll-a and suspended particulate matter (SPM) concentrations, and aCDOM/ap ratio. The rising water season presented a unique OWT with high SPM concentration and high relative contribution of ap to total absorption compared to the other OWTs. This bio-optical characterization of Lower Amazon River waters represents a first step for developing remote sensing inversion models adjusted to the optical complexity of this region.

54 ENVIRONMENTAL SCIENCES↗

Cloud Detection over Sunglint Regions with Observations from the Earth Polychromatic Imaging Camera (EPIC)

With the ability to observe the entire sunlit side of the Earth, EPIC data have become an important resource for studying cloud daily variability. Inaccurate cloud masking is a great source of uncertainty. One main region that is prone to error in cloud masking is the sunglint area over ocean surfaces. Cloud detection over these regions is challenging for the EPIC instrument because of its limited spectral channels. Clear sky ocean surface reflectance from visible channels over sunglint is much larger than that over the non-glint areas and can exceed reflectance from thin clouds. This paper presents an improved EPIC ocean cloud masking algorithm (Version 3). Over sunglint regions (glint angle ≤ 25°), the algorithm utilizes EPIC’s oxygen (O2) A-band ratio (764nm to 780 nm) in addition to the 780 nm reflectance observations in masking tests. Outside the sunglint regions, a dynamic reflectance threshold for the Rayleigh corrected 780 nm reflectance is applied. The thresholds are derived as a function of glint angle. When compared with co-located data from the geosynchronous Earth orbit (GEO) and the low Earth orbit (LEO) observations, the consistency of the new ocean cloud mask algorithm has increased by 4~10% and 4~6% in the glint center and granule edges respectively. The false positive rate is reduced by 10~17%. Overall global ocean cloud detection consistency increases by 2%. This algorithm, along with other improvements to the EPIC cloud masks, has been implemented in the EPIC cloud products Version 3. This algorithm will improve the cloud daily variability analysis by removing the artificial peak at local noon time in the glint center latitudes and reducing biases in the early morning and late afternoon cloud fraction over ocean surfaces.

EPIC↗

Scaler Rates from the Pierre Auger Observatory: A New Proxy of Solar Activity

The modulation of low-energy galactic cosmic rays reflects interplanetary magnetic field variations and can provide useful information on solar activity. An array of ground-surface detectors can reveal the secondary particles, which originate from the interaction of cosmic rays with the atmosphere. In this work, we present an investigation of the low-threshold rate (scaler) time series recorded in 16 yr of operation by the Pierre Auger Observatory surface detectors in Malargüe, Argentina. Through an advanced spectral analysis, we detected highly statistically significant variations in the time series with periods ranging from the decadal to the daily scale. We investigate their origin, revealing a direct connection with solar variability. Thanks to their intrinsic very low noise level, the Auger scalers allow a thorough and detailed investigation of the galactic cosmic-ray flux variations in the heliosphere at different timescales and can, therefore, be considered a new proxy of solar variability.

79 ASTRONOMY AND ASTROPHYSICS↗

Photoionization models for the wind from TW Vir

Line profiles are examined for the resonance doublets of the ions C IV, Si IV, and N V in the spectrum of the dwarf nova TW Virginis as observed by Cordova and Mason (1982). Line strengths depend on the abundances of the scattering ions and hence on the physical characteristics in the stellar wind. Results of the study show that the observed UV line ratios from TW Vir cannot be produced by photoionization by either an unattenuated single bremsstrahlung or by a single blackbody spectrum. Satisfactory values of the abundance ratios can be produced by bremsstrahlung spectra extending into the X-ray region which have low-energy cutoffs near the N IV threshold energy at 77 eV and have fluxes corresponding to an ionization parameter xi-sub-br of about 10. A finding that the flux in the soft X-ray bremsstrahlung component must exceed the EUV blackbody flux by at least an order of magnitude has implications for theories of the continuum emission from cataclysmic variables.

Kallman, T.↗

Erosive burning research

A status report is given on the results for the completed tests in a series of motor firings being carried out to measure the effects of the parameters that are considered to most strongly influence the scaling to larger rocket motor sizes of the transition to/or threshold conditions for erosive burning rate augmentation. Propellant burning rates at locations along the axis of the test motors are measured with a newly developed plasma capacitance gauge technique. The measured results are compared with erosive-burning predictions from a supporting ballistics analysis. The completed motor firings have successfully demonstrated response to the designed test variables. The trends with varying propellant burning rate, chamber pressure, and mass flow rate are consistent with existing results, but no pronounced effect of surface roughness has been observed. Rather, the influence of propellant oxidizer particle size on erosive burning is through its effect on the base, no-corssflow burning rate.

Strand, L.↗

Permeation Resistance of Chlorinated Polyethylene Against Hydrazine Fuels

The permeation resistance of chlorinated polyethylene (CPE) used in chemical protective clothing against the aerospace fuels hydrazine, monomethylhydrazine (MMH), and uns-dimethylhydrazine (UDMH) was determined by measuring breakthrough times and time-averaged vapor transmission rates using an ASTM F 739 permeation cell. Two exposure scenarios were simulated: a 2 hour (h) fuel vapor exposure, and a liquid fuel "splash" followed by a 2 h vapor exposure. To simulate internal suit pressure during operation, a positive differential pressure of 0.3 in. water (75 Pa) on the collection side of the permeation apparatus was used. Breakthrough was observed after exposure to liquid MMH, and to vapor and liquid UDMH. No breakthrough was observed after exposure to vapor and liquid hydrazine, or vapor MMH. A model was then used to calculate propellant concentrations inside a totally encapsulating chemical protective suit based on the ASTM permeation data obtained in the present study. Concentrations were calculated under conditions of fixed vapor transmission rate, variable breathing air flow rate, and variable splash exposure area. Calculations showed that the maximum allowable permeation rates of hydrazine fuels through CPE were of the order of 0.05 to 0.08 ng sq cm/min for encapsulating suits with low breathing air flow rates (of the order of 5 scfm or 140 L/min). Above these permeation rates, the 10 parts per billion (ppb) threshold limit value time - weighted average could be exceeded for chemical protective suits having a CPE torso. To evaluate suit performance at ppb level concentrations, use of a sensitive analytical method such as cation exchange high performance liquid chromatography with amperometric detection was found to be essential.

Waller, J. M.↗

Design and Calibration of Autonomous Coherent Doppler Lidar for Space Missions

Developed a new algorithm for the simulation of three dimensional homogeneous turbulent velocity fields. For typical atmospheric conditions it is impossible to produce a simulated velocity field that simultaneously satisfy a given spatial correlation and the corresponding spatial spectrum because of spectral aliasing. The new algorithms produce a turbulent velocity field which has accurate spatial correlations which is required for performance predictions from space-based systems. Developed a new algorithm for extracting the spatial statistics of the atmospheric velocity field using coherent Doppler lidar. The performance of the algorithm was compared with past methods and the new algorithm produces useful results for space-based data, which was not possible before. Developed new methods for verification of the errors in ground-based and space-based Doppler lidar wind measurements. These new methods do not require independent in situ data. This is an important issue for the verification of space-based Doppler lidar measurements of the global wind field. The performance of the new algorithm was compared with past results for both space-based and ground-based operation. The new algorithm has the best performance and is the only algorithm that performed satisfactory for spacebased operation. The performance of coherent Doppler lidar for a space missions with various scanning geometries was determined using computer simulation which contained the effects of random instrumental velocity errors, wind shear, wind variability along the range-gate and from shot-to-shot, and random variations in atmospheric aerosol backscatter over the measurement volume. The bias in the velocity estimates was small and the accuracy in the is typically less than 0.5 m/s for high signal conditions. For a large number of shot per velocity estimate, the threshold signal level for acceptable estimates is proportional to the number of shots to the minus one half power. This agrees with previous results determined for ground-based measurements. The use of multi-element optical detectors for autonomous operation of coherent Doppler lidar was shown to be a very promising technique. Optimal detector geometries were determined by computer simulation of performance: for ground-based testing with a fixed calibration target and for space-based operation using the random surface returns. The effects of refractive turbulence on ground-based calibration of coherent Doppler lidar was determined by computer simulations and compared with theoretical predictions. New techniques were required to correctly predict performance for the focused beam geometry commonly used for verification of space-based operation. An improved velocity estimator was evaluated for space-based applications were signal shot measurements are used to produce vector wind measurements. This permits more accurate measurements when the signal level is not known a priori or not available from multiple shot measurements. The average Doppler lidar signal spectrum including the effects of velocity turbulence was derived and calculated. This permits new estimation algorithms for turbulence based on spectral estimates. In situ atmospheric measurements were conducted and analyzed using an instrumented kite-platform. This work helps provide the required in situ data for verification of Doppler lidar velocity statistics.

Frehlich, Rod G.↗

Interannual Variability of Martian Global Dust Storms: Simulations with a Low-Order Model of the General Circulation

We present simulations of the interannual variability of martian global dust storms (GDSs) with a simplified low-order model (LOM) of the general circulation. The simplified model allows one to conduct computationally fast long-term simulations of the martian climate system. The LOM is constructed by Galerkin projection of a 2D (zonally averaged) general circulation model (GCM) onto a truncated set of basis functions. The resulting LOM consists of 12 coupled nonlinear ordinary differential equations describing atmospheric dynamics and dust transport within the Hadley cell. The forcing of the model is described by simplified physics based on Newtonian cooling and Rayleigh friction. The atmosphere and surface are coupled: atmospheric heating depends on the dustiness of the atmosphere, and the surface dust source depends on the strength of the atmospheric winds. Parameters of the model are tuned to fit the output of the NASA AMES GCM and the fit is generally very good. Interannual variability of GDSs is possible in the IBM, but only when stochastic forcing is added to the model. The stochastic forcing could be provided by transient weather systems or some surface process such as redistribution of the sand particles in storm generating zones on the surface. The results are sensitive to the value of the saltation threshold, which hints at a possible feedback between saltation threshold and dust storm activity. According to this hypothesis, erodable material builds up its a result of a local process, whose effect is to lower the saltation threshold until a GDS occurs. The saltation threshold adjusts its value so that dust storms are barely able to occur.

Pankine, A. A.↗

Fast and Flexible Multivariate Time Series Subsequence Search

Multivariate Time-Series (MTS) are ubiquitous, and are generated in areas as disparate as sensor recordings in aerospace systems, music and video streams, medical monitoring, and financial systems. Domain experts are often interested in searching for interesting multivariate patterns from these MTS databases which often contain several gigabytes of data. Surprisingly, research on MTS search is very limited. Most of the existing work only supports queries with the same length of data, or queries on a fixed set of variables. In this paper, we propose an efficient and flexible subsequence search framework for massive MTS databases, that, for the first time, enables querying on any subset of variables with arbitrary time delays between them. We propose two algorithms to solve this problem (1) a List Based Search (LBS) algorithm which uses sorted lists for indexing, and (2) a R*-tree Based Search (RBS) which uses Minimum Bounding Rectangles (MBR) to organize the subsequences. Both algorithms guarantee that all matching patterns within the specified thresholds will be returned (no false dismissals). The very few false alarms can be removed by a post-processing step. Since our framework is also capable of Univariate Time-Series (UTS) subsequence search, we first demonstrate the efficiency of our algorithms on several UTS datasets previously used in the literature. We follow this up with experiments using two large MTS databases from the aviation domain, each containing several millions of observations. Both these tests show that our algorithms have very high prune rates (>99%) thus needing actual disk access for only less than 1% of the observations. To the best of our knowledge, MTS subsequence search has never been attempted on datasets of the size we have used in this paper.

Bhaduri, Kanishka↗

Using an Automated Air Traffic Simulation Capability for a Parametric Study in Traffic Flow Management

Flight delays occur when demand for capacity-constrained airspace or airports exceeds predicted capacity. Demand for capacity-constrained airspace or airports can be controlled by a series of Traffic Management Initiatives (TMIs), which use departure and airborne delays, as well as pre-departure and airborne reroutes, to manage access to the constrained resources. Two systems exist in current and planned future operations to address imbalances between demand and capacity. The Collaborative Trajectory Options Program (CTOP) reduces demand to constrained resources by assigning strategic departure delay and pre-departure reroutes. Reroutes are selected from Trajectory Options Sets (TOSs) submitted by airlines. As flights approach the constrained resource, the Time-Based Flow Management System (TBFM) is used to assign tactical delay to satisfy constraints. This paper describes experiments performed to study the impact of varying levels of airline participation in CTOP via submission of TOSs on ground delay and flight time, and the impact of departure uncertainty on TBFM delays. Results suggest that as CTOP participation increases, average ground delays decrease for all airlines, but to the greatest extent for airlines participating in CTOP. A threshold in CTOP participation, which varies with the constraint capacity, is identified beyond which there is relatively little further reduction in average ground delays. Similarly, given the likely level of CTOP participation, the capacity reduction for which CTOP would be an appropriate TMI is also identified. Results also suggest that high average departure errors and high variability in departure error can make the prioritization of TBFM internal departures in TBFM metering and scheduling infeasible. Departure errors at current levels are, however, acceptable.

CTOP↗

Using an Automated Air Traffic Simulation Capability for a Parametric Study in Traffic Flow Management

Flight delays occur when demand for capacity-constrained airspace or airports exceeds predicted capacity. Demand for capacity-constrained airspace or airports can be controlled by a series of Traffic Management Initiatives (TMIs), which use departure and airborne delays, as well as pre-departure and airborne reroutes, to manage access to the constrained resources. Two systems exist in current and planned future operations to address imbalances between demand and capacity. The Collaborative Trajectory Options Program (CTOP) reduces demand to constrained resources by assigning strategic departure delay and pre-departure reroutes. Reroutes are selected from Trajectory Options Sets (TOSs) submitted by airlines. As flights approach the constrained resource, the Time-Based Flow Management System (TBFM) is used to assign tactical delay to satisfy constraints. This paper describes experiments performed to study the impact of varying levels of airline participation in CTOP via submission of TOSs on ground delay and flight time, and the impact of departure uncertainty on TBFM delays. Results suggest that as CTOP participation increases, average ground delays decrease for all airlines, but to the greatest extent for airlines participating in CTOP. A threshold in CTOP participation, which varies with the constraint capacity, is identified beyond which there is relatively little further reduction in average ground delays. Similarly, given the likely level of CTOP participation, the capacity reduction for which CTOP would be an appropriate TMI is also identified. Results also suggest that high average departure errors and high variability in departure error can make the prioritization of TBFM internal departures in TBFM metering and scheduling infeasible. Departure errors at current levels are, however, acceptable.

integrated demand management↗

The value of fusion energy to a decarbonized United States electric grid

Fusion could be a part of future decarbonized electricity systems, but it will need to compete with other technologies. In particular, pulsed tokamak plants have a unique operational mode, and evaluating which characteristics make them economically competitive can help select between design pathways. Using a capacity expansion and operations model, here we determined cost thresholds for pulsed tokamaks to reach a range of penetration levels in a future decarbonized US Eastern Interconnection. The required capital cost to reach a fusion capacity of 100 GW varied from $2,700 to $7,500 kW -1 , and the equilibrium penetration increases rapidly with decreasing cost. The value per unit power capacity depends on the variable operational cost and on the cost of its competition, particularly fission, much more than on the pulse cycle parameters. These findings can therefore provide initial cost targets for fusion more generally in the United States.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Feasibility of automated dropsize distributions from holographic data using digital image processing techniques

An automated technique for measuring particle diameters and their spatial coordinates from holographic reconstructions is being developed. Preliminary tests on actual cold-flow holograms of impinging jets indicate that a suitable discriminant algorithm consists of a Fourier-Gaussian noise filter and a contour thresholding technique. This process identifies circular as well as noncircular objects. The desired objects (in this case, circular or possibly ellipsoidal) are then selected automatically from the above set and stored with their parametric representations. From this data, dropsize distributions as a function of spatial coordinates can be generated and combustion effects due to hardware and/or physical variables studied.

Feinstein, S. P.↗