Search NASA⌕ Search

SEARCH · Search NASA

Results for “Direct simulation Monte Carlo (DSMC)”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Estimating Drag and Heating Coefficients for Hollow Reentry Objects in Transitional Flow Using DSMC

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These methods are typically limited to convex, blunt objects (such as spheres) and are applied to other objects such as boxes and cylinders using multiplicative “shape factors” to account for the different behavior. Previous literature has analyzed the aerodynamic and aerothermodynamic properties of flow around sharp-edged objects like boxes and cylinders in transitional flow, though only those objects with solid external boundaries. However, many reentry objects we have encountered in real spacecraft have been hollow (i.e., with the potential to allow flow through them). We present here preliminary results from analyses performed using the NASA Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) on hollow cylinders and boxes (with varying wall thickness-diameter ratio).

Marichalar, Jeremiah J.↗
NASA NTRS2004

Remote Sensing of Icy Galilean Moon Surface and Atmospheric Composition Using Low Energy (1 eV-4 keV) Neutral Atom Imaging

We describe a low energy neutral atom imager suitable for composition measurements Europa and other icy Galilean moons in the Jovian magnetosphere. This instrument employs conversion surface technology and is sensitive to either neutrals converted to negative ions, neutrals converted to positive ions and the positive ions themselves depending on the power supply. On a mission such as the Jupiter Icy Moons Orbiter (JIMO), two back-to-back sensors would be flown with separate power supplies fitted to the neutral atom and iodneutral atom sides. This will allow both remote imaging of 1 eV < E < 4 keV neutrals from icy moon surfaces and atmospheres, and in situ measurements of ions at similar energies in the moon ionospheres and Jovian magnetospheric plasma. The instrument provides composition measurements of the neutrals and ions that enter the spectrometer with a mass resolution dependent on the time-of-flight subsystem and capable of resolving molecules. The lower energy neutrals, up to tens of eV, arise from atoms and molecules sputtered off the moon surfaces and out of the moon atmospheres by impacts of more energetic (keV to MeV) ions from the magnetosphere. Direct Simulation Monte Carlo (DSMC) models are used to convert measured neutral abundances to compositional distributions of primary and trace species in the sputtered surfaces and atmospheres. The escaping neutrals can also be detected as ions after photo- or plasma-ionization and pickup. Higher energy, keV neutrals come from charge exchange of magnetospheric ions in the moon atmospheres and provide information on atmospheric structure. At the jovicentric orbits of the icy moons the presence of toroidal gas clouds, as detected at Europa's orbit, provide M e r opportunities to analyze both the composition of neutrals and ions originating from the moon surfaces, and the characteristics of magnetospheric ions interacting with neutral cloud material. Charge exchange of low energy ions near the moons, and directional distributions of the resultant neutrals, allow indirect global mapping of magnetic field structures around the moons. Temporal variation of the magnetic structures can be linked to induced magnetic fields associated with subsurface oceans.

Collier, M. R.↗
NASA NTRS2001

Chemical and Flowfield Modeling for Enhanced Analysis of Contamination Experiments

This paper describes the application of a new Direct Simulation Monte Carlo (DSMC) code, the Molecular Beam Simulator (MBS), which is designed to analyze laboratory scale molecular beam-surface (and crossed-beam) experiments. The MBS is primarily intended to model experiments associated with spacecraft contamination effects, but it can also be used to simulate a variety of surface chemistry and reactive flow measurements. The MBS code is fully three-dimensional, includes a wide-range of chemical processes, and can model one or multiple pulsed (non-steady) sources. As an example application of the MBS code, a fast, pulsed, oxygen atom-surface experiment which examines the chemistry behind erosion of graphite by oxygen atoms is analyzed. Unsteady DSMC simulations show that experimental observations of excited molecular states after the pulse has hit the surface are consistent with two distinct chemical mechanisms: a direct one where the excited molecules are formed on the surface, and a two-step mechanism where ground state molecules formed on the surface are collisionally excited after they leave the surface by trailing oxygen atoms in the pulse. Further DSMC calculations suggest experiments which can distinguish between these mechanisms.

Braunstein, Matthew↗
NASA NTRS1993

Simulated rarefied entry of the Galileo probe into the atmosphere of Jupiter

Flow properties and aerodynamics are computed with a direct simulation Monte Carlo (DSMC) method for rarefied entry of the Galileo Probe into the atmosphere of Jupiter. Accurate predictions of vehicle drag coefficients are needed in order to assess atmospheric properties from the onboard Atmosphere Structure Experiment where highly-sensitive accelerometers will measure the drag force to within 10-6 barr during the initial entry phase at high altitudes. The corresponding flow rarefraction extends from the free molecule limit to the near continuum transition regime (Re less than 1000). Simulation results indicate that C(sub D) varies from 2.1 at the free molecule limit down to 1.6 at Re(infinity) = 1,000. Temperatures, densities, and internal energies throughout the flow field were also computed at each altitude ranging from 735 km to 353 km above the 1 barr level in the Jovian atmosphere. Surface heating and temperatures of the probe were computed directly in the DSMC code by assuming radiative equilibrium. Material response was re-asssessed accurately during entry by accounting for conductivity, heat capacity, and pyrolysis which led to surface material mass efflux several times that of the freestream mass influx. The simulation also accounted for the quantum nature of the rotational energy mode of the dominant atmospheric species H2 through partial internal excitation in the freestream gas.

Hass, Brian L.↗
NASA NTRS2002

Simulations of Ground and Space-Based Oxygen Atom Experiments

Fast, pulsed atomic oxygen sources are a key tool in ground-based investigations of spacecraft contamination and surface erosion effects. These technically challenging ground-based studies provide a before and after picture of materials under low-earth-orbit (LEO) conditions. It would be of great interest to track in real time the pulsed flux from the source to the surface sample target and beyond in order to characterize the population of atoms and molecules that actually impact the surface and those that make it downstream to any coincident detectors. We have performed simulations in order to provide such detailed descriptions of these ground-based measurements and to provide an assessment of their correspondence to the actual LEO environment. Where possible we also make comparisons to measured fluxes and erosion yields. To perform the calculations we use a detailed description of a measurement beam and surface geometry based on the W, pulsed apparatus at Montana State University. In this system, a short pulse (on the order of 10 microseconds) of an O/O2 beam impacts a flat sample about 40 cm downstream and slightly displaced &om the beam s central axis. Past this target, at the end of the beam axis is a quadrupole mass spectrometer that measures the relative in situ flux of 0102 to give an overall normalized erosion yield. In our simulations we use the Direct Simulation Monte Carlo (DSMC) method, and track individual atoms within the atomic oxygen pulse. DSMC techniques are typically used to model rarefied (few collision) gas-flows which occur at altitudes above approximately 110 kilometers. These techniques are well suited for the conditions here, and multi-collision effects that can only be treated by this or a similar technique are included. This simulation includes collisions with the surface and among gas atoms that have scattered from the surface. The simulation also includes descriptions of the velocity spread and spatial profiles of the O/O2 beam obtained from separate measurements. These computations use basic engineering models for the gas-gas and gas-surface scattering and focus on the influence of multi-collision effects. These simulations characterize many important quantities of interest including the actual flux of atoms that reach the surface, the energy distribution of this flux, as well as the direction of the velocity of the flux that strikes the surface. These quantities are important in characterizing the conditions which give rise to measured surface erosion. The calculations also yield time- snapshots of the pulse as it impacts and flows around the surface. These snapshots reveal the local environment of gas near the surface for the duration of the pulse. We are also able to compute the flux of molecules that travel downstream and reach the spectrometer, and we characterize their velocity distribution. The number of atoms that reach the spectrometer can in fact be influenced by the presence of the surface due to gas-gas collisions from atoms scattered h m the surface, and it will generally be less than that with the surface absent. This amounts to an overall normalization factor in computing erosion yields. We discuss these quantities and their relationship to the gas-surf$ce interaction parameters. We have also performed similar calculations corresponding to conditions (number densities, temperatures, and velocities) of low-earth orbit. The steady-state nature and lower overall flux of the actual space environment give rise to differences in the nature of the gas-impacts on the surface from those of the ground-based measurements using a pulsed source.

Minton, T. K.↗
NASA NTRS1993

Particle kinetic simulation of high altitude hypervelocity flight

In this grant period, the focus has been on enhancement and application of the direct simulation Monte Carlo (DSMC) particle method for computing hypersonic flows of re-entry vehicles. Enhancement efforts dealt with modeling gas-gas interactions for thermal non-equilibrium relaxation processes and gas-surface interactions for prediction of vehicle surface temperatures. Both are important for application to problems of engineering interest. The code was employed in a parametric study to improve future applications, and in simulations of aeropass maneuvers in support of the Magellan mission. Detailed comparisons between continuum models for internal energy relaxation and DSMC models reveals that several discrepancies exist. These include definitions of relaxation parameters and the methodologies for implementing them in DSMC codes. These issues were clarified and all differences were rectified in a paper (Appendix A) submitted to Physics of Fluids A, featuring several key figures in the DSMC community as co-authors and B. Haas as first author. This material will be presented at the Fluid Dynamics meeting of the American Physical Society on November 21, 1993. The aerodynamics of space vehicles in highly rarefied flows are very sensitive to the vehicle surface temperatures. Rather than require prescribed temperature estimates for spacecraft as is typically done in DSMC methods, a new technique was developed which couples the dynamic surface heat transfer characteristics into the DSMC flow simulation code to compute surface temperatures directly. This model, when applied to thin planar bodies such as solar panels, was described in AIAA Paper No. 93-2765 (Appendix B) and was presented at the Thermophysics Conference in July 1993. The paper has been submitted to the Journal of Thermophysics and Heat Transfer. Application of the DSMC method to problems of practical interest requires a trade off between solution accuracy and computational expense and limitations. A parametric study was performed and reported in AIAA Paper No. 93-2806 (Appendix C) which assessed the accuracy penalties associated with simulations of varying grid resolution and flow domain size. The paper was also presented at the Thermophysics Conference and will be submitted to the journal shortly. Finally, the DSMC code was employed to assess the pitch, yaw, and roll aerodynamics of the Magellan spacecraft during entry into the Venus atmosphere at off-design attitudes. This work was in support of the Magellan aerobraking maneuver of May 25-Aug. 3, 1993. Furthermore, analysis of the roll characteristics of the configuration with canted solar panels was performed in support of the proposed 'Windmill' experiment. Results were reported in AIAA Paper No. 93-3676 (Appendix D) presented at the Atmospheric Flight Mechanics Conference in August 1993, and were submitted to Journal of Spacecraft and Rockets.

Haas, Brian L.↗
NASA NTRSDate not supplied

Simulation of etch pit formation through active sites in carbon fiber micro-structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation as a result of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gasses such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗
NASA NTRS1997

Predictive Modeling in Plasma Reactor and Process Design

Research continues toward the improvement and increased understanding of high-density plasma tools. Such reactor systems are lauded for their independent control of ion flux and energy enabling high etch rates with low ion damage and for their improved ion velocity anisotropy resulting from thin collisionless sheaths and low neutral pressures. Still, with the transition to 300 mm processing, achieving etch uniformity and high etch rates concurrently may be a formidable task for such large diameter wafers for which computational modeling can play an important role in successful reactor and process design. The inductively coupled plasma (ICP) reactor is the focus of the present investigation. The present work attempts to understand the fundamental physical phenomena of such systems through computational modeling. Simulations will be presented using both computational fluid dynamics (CFD) techniques and the direct simulation Monte Carlo (DSMC) method for argon and chlorine discharges. ICP reactors generally operate at pressures on the order of 1 to 10 mTorr. At such low pressures, rarefaction can be significant to the degree that the constitutive relations used in typical CFD techniques become invalid and a particle simulation must be employed. This work will assess the extent to which CFD can be applied and evaluate the degree to which accuracy is lost in prediction of the phenomenon of interest; i.e., etch rate. If the CFD approach is found reasonably accurate and bench-marked with DSMC and experimental results, it has the potential to serve as a design tool due to the rapid time relative to DSMC. The continuum CFD simulation solves the governing equations for plasma flow using a finite difference technique with an implicit Gauss-Seidel Line Relaxation method for time marching toward a converged solution. The equation set consists of mass conservation for each species, separate energy equations for the electrons and heavy species, and momentum equations for the gas. The sheath is modeled by imposing the Bohm velocity to the ions near the walls. The DSMC method simulates each constituent of the gas as a separate species which would be analogous in CFD to employing separate species mass, momentum, and energy equations. All particles including electrons are moved and allowed to collide with one another with the stipulation that the electrons remain tied to the ions consistent with the concept of ambipolar diffusion. The velocities of the electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. These benefits come at a price in terms of computational time and memory. The DSMC and CFD are made as consistent as possible by using similar chemistry and power deposition models. Although the comparison of CFD and DSMC is interesting, the main goal of this work is the increased understanding of high-density plasma flowfields that can then direct improvements in both techniques. This work is unique in the level of the physical models employed in both the DSMC and CFD for high-density plasma reactor applications. For example, the electrons are simulated in the present DSMC work which has not been done before for low temperature plasma processing problems. In the CFD approach, for the first time, the charged particle transport (discharge physics) has been self-consistently coupled to the gas flow and heat transfer.

Hash, D. B.↗
NASA NTRS1984

Direct simulation of transitional flow for hypersonic reentry conditions

This paper presents results of flowfield calculations for typical hypersonic reentry conditions encountered by the nose region of the Space Shuttle Orbiter. Most of the transitional flow regime is covered by the altitude range of 150 to 92 km. Calculations were made with the Direct Simulation Monte Carlo (DSMC) method that accounts for translational, rotational, vibrational, and chemical nonequilibrium effects. Comparison of the DSMC heating results with both Shuttle flight data and continuum predictions showed good agreement at the lowest altitude considered. However, as the altitude increased, the continuum predictions, which did not include slip effects, departed rapidly from the DSMC results by overpredicting both heating and drag. The results demonstrate the effects of rarefaction on the shock and the shock layer, along with the extent of the slip and temperature jump at the surface. Also, the sensitivity of the flow structure to the gas-surface interaction model, thermal accommodation, and surface catalysis are studied.

Moss, J. N.↗
NASA NTRS1985

Monte Carlo simulations in support of the Shuttle upper atmospheric mass spectrometer experiment

This paper presents the results of calculations obtained with a Direct Simulation Monte Carlo (DSMC) method that describes both the external flow about the nose region of the Shuttle Orbiter and the internal flow in an inlet tube that provides the initial path between the shock-processed gases and a mass spectrometer mounted inside the Orbiter. A dedicated, three-dimensional version of the DSMC was developed during this study for the internal flow simulation along with the procedure for interfacing the external and internal flows. The calculations span the 140 to 95 km altitude range, which includes most of the hypersonic transitional flow encountered during reentry. All calculations were for a multicomponent gas mixture consisting of five chemical species while simulating the effects of transitional, rotational, vibrational, and chemical nonequilibrium. The results of the calculations show that within the entry region of the inlet tube where the gas has equilibrated with the sidewall tube temperature, the pressure is substantially less than the pressure at the external surface. This pressure correction for the entry region is significant for all conditions investigated and increases with altitude. The results highlight the structure of both the external and internal flow. Results of parametric studies show the sensitivity of the equilibrated inlet sidewall pressure to mass flow rates, gas-surface reflection model, tube side-wall temperature variations, and surface recombinations.

Moss, J. N.↗
NASA NTRS1986

Direct simulation of hypersonic flows over blunt slender bodies

Results of a numerical study of low-density hypersonic flow about cylindrically blunted wedges and spherically blunted cones with body half angles of 0, 5, and 10 deg are presented. Most of the transitional flow regime encountered during entry between the free molecule and continuum regimes is simulated for a reentry velocity of 7.5 km/s by including freestream conditions of 70 to 100 km. The bodies are at zero angle of incidence and have diffuse and finite catalytic surfaces. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method. The numerical simulations show that noncontinuum effects such as surface temperature jump, and velocity slip are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of 96 km for the two-dimensional configurations. Comparisons between the DSMC and continuum viscous shock-layer calculations highlight the significant difference in flowfield structure predicted by the two methods.

Moss, J. N.↗
NASA NTRS1986

Hypersonic low-density solutions of the Navier-Stokes equations with chemical nonequilibrium and multicomponent surface slip

Solutions of the Navier-Stokes equations with chemical nonequilibrium and multicomponent surface slip are presented along the stagnation streamline under low-density hypersonic flight conditions. The conditions analyzed are those encountered by the nose region of the Space Shuttle Orbiter during reentry. A detailed comparison of the Navier-Stokes (NS) results is made with the viscous shock-layer (VSL) and Direct Simulation Monte Carlo (DSMC) predictions. With the inclusion of surface-slip boundary conditions in NS calculations, the surface heat transfer and other flow field quantities adjacent to the surface are predicted favorably with the DSMC calculations from 75 km to 115 km in altitude. Therefore, the practical range for the applicability of Navier-Stokes solutions is much wider than previously thought. This is appealing because the continuum (NS and VSL) methods are commonly used to solve the fluid flow problems and are less demanding in terms of computer resource requirements than the noncontinuum (DSMC) methods. The NS solutions agree well with the VSL results for altitudes less than 92 km. An assessment is made of the frozen flow approximation employed in the VSL calculations.

Gupta, R. N.↗
NASA NTRS1986

Direct simulation of hypersonic transitional flow

A review of recent calculations obtained with Bird's direct simulation Monte Carlo (DSMC) method are presented for the transitional flowfield encountered at reentry conditions. Consequently, the emphasis is to simulate the real-gas effects resulting from a highly energetic nonequilibrium flow. The DSMC calculations for both wide-angle and slender bodies are compared with experimental and continuum calculations. The wide angle body calculations simulate the windward nose of the Shuttle Orbiter and projected aeroassisted orbital transfer vehicles. For the slender, blunted configurations, both two-dimensional and axisymmetric results are presented. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant.

Moss, J. N.↗
NASA NTRS1987

Nonequilibrium effects for hypersonic transitional flows

Presented are the results of numerical simulations of hypersonic flow about blunt cones and hemispherical nose configurations for reentry velocities of 7.5 and 10 km/s. Cone half angles 0, 5, and 10 deg are considered at zero angle of incidence; however, the focus is for the 5 deg cone. The body size and altitude ranges considered (70 to 110 km) are such that the flow is in the transitional regime. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method of Bird. The DSMC results are compared with those obtained with viscous shock-layer and Navier-Stokes methods. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant. The current calculations show that the binary scaling similitude provides a means of correlating the blunt body surface quantities in the hypersonic, transitional regime. Furthermore, for the higher velocity entry conditions, the results highlight some of the concerns in the application of multitemperature continuum formulations, particularly the use of some proposed functional relations for the chemical rate constants under thermodynamic nonequilibrium conditions.

Moss, James N.↗
NASA NTRS1987

Direct simulation of stagnation streamline flow for hypersonic reentry

The direct simulation Monte Carlo (DSMC) method is applied to simulate one-dimensional flow along the stagnation streamline. The freestream conditions considered are those encountered by the nose region of the Space Shuttle Orbiter during the hypersonic reentry. The range of altitudes (75 to 92 km) considered in the present calculations covers continuum to the less rarefied portion of the transition flow regime. The calculations account for thermal as well as the chemical nonequilibrium effects. The attention is focused on the flow structure along the stagnation streamline for different specified shock locations at a given altitude. The effect of shock location on the stagnation-point gas composition and heat transfer is analyzed for specific altitude; consequently, these results are appropriate for hemispherical bodies with different nose radii. Finally, comparison is made between the present DSMC results and the viscous shock-layer calculations, which show good agreement at the lowest altitude.

Dogra, Virendra K.↗
NASA NTRS1986

Numerical simulations of rarefied reentry flows

A review of recent calculations obtained with Bird's direct simulation Monte Carlo (DSMC) method is presented for the transitional flowfield encountered at reentry conditions. Consequently, the emphasis is the real-gas effects resulting from a highly energetic nonequilibrium flow. The DSMC calculations for both wide-angle and slender bodies are compared with continuum calculations. The wide-angle body calculations simulate the windward nose of the Shuttle Orbiter and projected aeroassisted orbital transfer vehicles. For the slender, blunted configurations, both two-dimensional and axisymmetric results are presented. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant.

Moss, James N.↗
NASA NTRS1987

Direct simulation of high-vorticity gas flows

The computational limitations associated with the molecular dynamics (MD) method and the direct simulation Monte Carlo (DSMC) method are reviewed in the context of the computation of dilute gas flows with high vorticity. It is concluded that the MD method is generally limited to the dense gas case in which the molecular diameter is one-tenth or more of the mean free path. It is shown that the cell size in DSMC calculations should be small in comparison with the mean free path, and that this may be facilitated by a new subcell procedure for the selection of collision partners.

Bird, G. A.↗
NASA NTRS1987

Direct simulation of aerothermal loads for an aeroassist flight experiment vehicle

Results of a numerical study using the direct simulation Monte Carlo (DSMC) method are presented for the hypersonic flow about an elliptically blunted cone. The flow conditions are those for a proposed Aeroassist Flight Experiment (AFE) vehicle. The altitude range considered is that from 130 to 90 km which encompasses most of the transitional flow regime for the AFE vehicle, that is, the region bounded by free molecular and continuum flow. Freestream velocities of 9.9 to 7.5 km/sec are considered. The numerical simulations show that noncontinuum effects are evident for all cases considered. The onset of chemical dissociation occurs at a simulated altitude of about 130 km. Results presented highlight the thermal and chemical nonequilibrium nature of the flowfield and the impact of these effects on the surface heating and body drag. A calculation which included the additional effects of ionization and thermal radiation demonstrates that the inclusion of such efects would not significantly alter the surface quantities calculated in the present study. The radiative heating is negligible when compared with the convective heating, and the same would be true for the other conditions considered.

Dogra, Virendra K.↗