Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LAYERS”

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 433 records · Page 24

Density response of the mesospheric sodium layer to gravity wave perturbations

Lidar observations of the mesospheric sodium layer often reveal wavelike features moving through the layer. It is often assumed that these features are a layer density response to gravity waves. Chiu and Ching (1978) described the approximate form of the linear response of atmospheric layers to gravity waves. In this paper, their results are used to predict the response of the sodium layer to gravity waves. These simulations are compared with experimental observations and a good correlation is found between the two. Because of the thickness of the sodium layer and the density gradients found in it, a linear model of the layer response is not always adequate to describe gravity wave-sodium layer interactions. Inclusion of nonlinearities in the layer response is briefly discussed. Experimental data is seen to contain features consistent with the predicted nonlinearities.

Shelton, J. D.↗

Near-field turbulence properties of single- and two-stream plane mixing layers

Detailed mean-flow and turbulence measurements have been made in the near-field of two plane mixing layers in air with a maximum velocity of 21 m/s. The experimental rig enabled mixing layers of velocity ratios 0 and 0.46 to be generated simultaneously. Cases with both tripped and untripped initial boundary layers were studied. In all cases, it was found that the two-stream layer developed to the self-preserving state in a distance much shorter than the single-stream layer which followed accepted criteria for the development distance. The asymptotic levels of the turbulence quantities in the two-stream layer, and the development of the single-stream layer showed agreement with existing data. The results suggest that the two-stream mixing layer should provide a better test case for the development of turbulence models and calculation methods than the single-stream mixing layer.

Mehta, R. D.↗

An experimental study of plane mixing layer development

Mean flow and turbulence was measured in the near field of two plane mixing layers in air with a maximum velocity of 21 m/sec. The experimental rig enabled mixing layers of velocity ratios 0 and 0.46 to be generated simultaneously. Gases with both tripped and untripped initial boundary layers were studied. It was found that the two stream layer developed to the self preserving state in a distance much shorter than did the single stream layer which followed accepted criteria for the development distance. The asymptotic levels of the turbulence quantities in the two stream layer and the development of the single stream layer showed agreement with existing data. It is suggested that the two stream mixing layer should provide a better test case for the development of turbulence models and calculation methods than the single stream mixing layer.

Mehta, R. D.↗

An experimental investigation of nozzle-exit boundary layers of highly heated free jets

An experimental investigation of the effects of nozzle operating conditions on the development of nozzle-exit boundary layers of highly heated air free jets is reported in this paper. The total pressure measurements in the nozzle-exit boundary layer were obtained at a range of jet Mach numbers from 0.1 to 0.97 and jet total temperatures up to 900 K. The analysis of results shows that the nozzle-exit laminar boundary-layer development depends only on the nozzle-exit Reynolds number. For the nozzle-exit turbulent boundary layer, however, it appears that the effects of the jet total temperature on the boundary-layer integral characteristics are independent from the effects of the nozzle-exit Reynolds number. This surprizing finding has not yet been reported. Further, laminar boundary-layer profiles were compared with the Pohlhausen solution for a flat-wall converging channel and an acceptable agreement was found only for low Reynolds numbers. For turbulent boundary layers, the dependence of the shape factor on relative Mach numbers at a distance of one momentum thickness from the nozzle wall resembles Spence's prediction. Finally, the calculated total pressure loss coefficient was found to depend on the nozzle-exit Reynolds number for the laminar nozzle-exit boundary layer, while for the turbulent exit boundary layer this coefficient appears to be constant.

Lepicovsky, J.↗

The structure of a three-dimensional turbulent boundary layer

The three-dimensional turbulent boundary layer is shown to have a self-consistent two-layer asymptotic structure in the limit of large Reynolds number. In a streamline coordinate system, the streamwise velocity distribution is similar to that in two-dimensional flows, having a defect-function form in the outer layer which is adjusted to zero at the wall through an inner wall layer. An asymptotic expansion accurate to two orders is required for the cross-stream velocity which is shown to exhibit a logarithmic form in the overlap region. The inner wall-layer flow is collateral to leading order but the influence of the pressure gradient, at large but finite Reynolds numbers, is not negligible and can cause substantial skewing of the velocity profile near the wall. Conditions under which the boundary layer achieves self-similarity and the governing set of ordinary differential equations for the outer layer are derived. The calculated solution of these equations is matched asymptotically to an inner wall-layer solution and the composite profiles so formed describe the flow throughout the entire boundary layer. The effects of Reynolds number and cross-stream pressure gradient on the crossstream velocity profile are discussed and it is shown that the location of the maximum cross-stream velocity is within the overlap region.

Degani, A. T.↗

Determination of general relations for the behavior of turbulent boundary layers

An analysis has been made of a considerable amount of data for turbulent boundary layers along wings and bodies of various shapes in order to determine the fundamental variables that control the development of turbulent boundary layers. It was found that the type of velocity distribution in the boundary layer could be expressed in terms of a single parameter. This parameter was chosen as the ratio of the displacement thickness to the momentum thickness of the boundary layer. The variables that control the development of the turbulent boundary layer apparently are: (1) the ratio of the nondimensional pressure gradient, expressed in terms of the local dynamic pressure outside the boundary layer and boundary-layer thickness, to the local skin-friction coefficient and (2) the shape of the boundary layer. An empirical equation has been developed in terms of these variables that, when used with the momentum equation and the skin-friction relation, makes it possible to trace the development of the turbulent boundary layer to the separation point.

Von Doenhoff, Albert E↗

The Effect of the Inlet Mach Number and Inlet-boundary-layer Thickness on the Performance of a 23 Degree Conical-diffuser-tail-pipe Combination

An investigation was conducted to determine the effect of the inlet Mach number and entrance-boundary-layer thickness on the performance of a 23 degree 21-inch conical-diffuser - tail-pipe combination with a 2:1 area ratio. The air flows used in this investigation covered an inlet Mach number range from 0.17 to 0.89 and corresponding Reynolds numbers of 1,700,000 to 7,070,000. Results are reported for two inlet-boundary-layer thicknesses. Over the entire range of flows, the mean value of the inlet displacement thickness is about 0.034 inch for the thinner inlet boundary layer and about 0.170 inch for the case of the thicker inlet boundary layer. The performance of the diffuser - tail-pipe combination is presented together with examples of longitudinal static-pressure distribution and the results of boundary-layer pressure surveys made at six points along the diffuser wall. The results indicated a progressive diminution of the static-pressure recovery and a steady increase in the total-pressure losses as the inlet Mach number was increased for both inlet-boundary-layer thicknesses. The ratio of actual static-pressure rise to that theoretically possible was much less and the total-pressure losses were greater for the case of the thicker inlet boundary layer throughout the speed range investigated. With the thinner inlet boundary layer, flow separation occurred at the diffuser exit at all inlet Mach numbers.Unseparated flow alternating with separated flow was observed near the inlet at the higher velocities. For the case of the thicker inlet boundary layer, the origin of the separated region occurred in the vicinity of the inlet-duct-diffuser junction section at all Mach numbers.

Persh, Jerome↗

Discrete-Layer Piezoelectric Plate and Shell Models for Active Tip-Clearance Control

The objectives of this work were to develop computational tools for the analysis of active-sensory composite structures with added or embedded piezoelectric layers. The targeted application for this class of smart composite laminates and the analytical development is the accomplishment of active tip-clearance control in turbomachinery components. Two distinct theories and analytical models were developed and explored under this contract: (1) a discrete-layer plate theory and corresponding computational models, and (2) a three dimensional general discrete-layer element generated in curvilinear coordinates for modeling laminated composite piezoelectric shells. Both models were developed from the complete electromechanical constitutive relations of piezoelectric materials, and incorporate both displacements and potentials as state variables. This report describes the development and results of these models. The discrete-layer theories imply that the displacement field and electrostatic potential through-the-thickness of the laminate are described over an individual layer rather than as a smeared function over the thickness of the entire plate or shell thickness. This is especially crucial for composites with embedded piezoelectric layers, as the actuating and sensing elements within these layers are poorly represented by effective or smeared properties. Linear Lagrange interpolation polynomials were used to describe the through-thickness laminate behavior. Both analytic and finite element approximations were used in the plane or surface of the structure. In this context, theoretical developments are presented for the discrete-layer plate theory, the discrete-layer shell theory, and the formulation of an exact solution for simply-supported piezoelectric plates. Finally, evaluations and results from a number of separate examples are presented for the static and dynamic analysis of the plate geometry. Comparisons between the different approaches are provided when possible, and initial conclusions regarding the accuracy and limitations of these models are given.

Heyliger, P. R.↗

Boundary layer transition studies

A small-scale wind tunnel previously used for turbulent boundary layer experiments was modified for two sets of boundary layer transition studies. The first study concerns a laminar separation/turbulent reattachment. The pressure gradient and unit Reynolds number are the same as the fully turbulent flow of Spalart and Watmuff. Without the trip wire, a laminar layer asymptotes to a Falkner & Skan similarity solution in the FPG. Application of the APG causes the layer to separate and a highly turbulent and approximately 2D mean flow reattachment occurs downstream. In an effort to gain some physical insight into the flow processes a small impulsive disturbance was introduced at the C(sub p) minimum. The facility is totally automated and phase-averaged data are measured on a point-by-point basis using unprecedently large grids. The evolution of the disturbance has been tracked all the way into the reattachment region and beyond into the fully turbulent boundary layer. At first, the amplitude decays exponentially with streamwise distance in the APG region, where the layer remains attached, i.e. the layer is viscously stable. After separation, the rate of decay slows, and a point of minimum amplitude is reached where the contours of the wave packet exhibit dispersive characteristics. From this point, exponential growth of the amplitude of the disturbance is observed in the detached shear layer, i.e. the dominant instability mechanism is inviscid. A group of large-scale 3D vortex loops emerges in the vicinity of the reattachment. Remarkably, the second loop retains its identify far downstream in the turbulent boundary layer. The results provide a level of detail usually associated with CFD. Substantial modifications were made to the facility for the second study concerning disturbances generated by Suction Holes for laminar flow Control (LFC). The test section incorporates suction through interchangeable porous test surfaces. Detailed studies have been made using isolated holes in the impervious test plate that used to establish the Blasius base flow. The suction is perturbed harmonically and data are averaged on the basis of the phase of the disturbance, for conditions corresponding to strong suction and without suction. The technique was enhanced by using up to nine multiple probes to reduce the experimental run-time. In both cases, 3D contour surfaces in the vicinity of the hole show highly 3D TS waves which fan out in the spanwise direction forming bow-shaped waves downstream. The case without suction has proved useful for evaluating calculation methods. With suction, the perturbations on the centerline are much stronger and decay less rapidly, while the TS waves in the far field are similar to the case without suction. Downstream, the contour surfaces of the TS waves develop spanwise irregularities which eventually form into clumps. The spanwise clumping is evidence of a secondary instability that could be associated with suction vortices. Designers of porous surfaces use Goldsmith's Criterion to minimize cross-stream interactions. It is shown that partial TS wave cancellation is possible, depending on the hole spacing, disturbance frequency and free-stream velocity. New high-performance Constant Temperature Hot-Wire Anemometers were designed and built, based on a linear system theory analysis that can be extended to arbitrary order. The motivation was to achieve the highest possible frequency reponse while ensuring overall system stability. The performance is equal to or superior to commercially available instruments at about 10% of the cost. Details, such as fabrication drawings and a parts list, have been published to enable the instrument to be construced by others.

Watmuff, Jonathan H.↗

Mars Secular Obliquity Decrease and the Layered Terrain

Mars may have substantially decreased its average axial tilt over geologic time due to the waxing and waning of water ice caps through the phenomenon of climate friction (also called obliquity-oblateness feedback). Depending upon Mars' climate and internal structure, water caps of the order of 10(exp 17) - 10(exp 18) kg cycling with the obliquity oscillations could have either increased or decreased the average obliquity by possibly tens of degrees over the age of the solar system. Gravity and topography observations by the Mars Global Surveyor indicate that the south polar cap is mostly uncompensated, so that Mars may be largely rigid on the obliquity timescale. Further, Mars may be a water-rich planet, so that there is a large phase angle between insolation forcing and the size of the obliquity-driven water caps. A stiff, water-rich planet implies an obliquity decrease over the eons. Such a decrease might account for the apparent youthfulness of the polar layered terrain. The idea is that fewer volatiles were available to be cycled into and out of the terrain at high mean obliquity, because of the eveness of insolation between equator and pole, and because of small insolation variations as the obliquity oscillated; so that the movement of volatiles and dust produced thin layers or perhaps no layers at all. As the average obliquity decreased, the insolation contrast between high and low latitudes increased, plus the insolation variations over the obliquity cycle grew somewhat bigger, so that more volatiles and dust might have shuttled into and out of the polar regions, forming the observed thick layers late in Mars' history. It may also be that the average tilt has decreased to the point where the climate friction mechanism is starving itself: more and more water has gotten locked up in the polar regions, making less available for cycling with the oscillations. And the layer-forming mechanism may be starving too: not only less water, but also at low average obliquity the atmospheric pressure is lower, so that dust storms become rarer, making less dust available to be deposited in the layers. The rate of water ice sublimation and the atmospheric pressure are both expected to drop sharply for obliquities only slightly lower than the present one. Hence there are qualitative reasons for believing that (a) the mean obliquity may be about as low as it can get, and (b) the mechanism for making the layers may turn off at low average obliquity, so that the layered terrain is a transient phenomenon. Thus the formation of the layered terrain may be intimately connected to a change in the mean axial tilt.

Rubincam, D. P.↗

Turbulent Heat Transfer Characteristics in the Shear Layer of a Separated Flow

Experiments were performed to study the evolution of the heat transfer structure in a separated free shear layer region of an incompressible separated turbulent boundary layer flow behind a backward-facing step. While there is an abundance of velocity field measurements of separated flows, heat transfer measurements are rather scarce, thus limiting assessment of the heat transfer physics and its accurate modeling. The purpose of the paper is twofold: to improve an understanding of effects of flow separation on heat transfer characteristics, and to provide data for turbulence modeling and computation. The boundary layer upstream of the step was turbulent and fully developed. A constant temperature surface boundary condition was imposed upstream and downstream of the step for the heat transfer study. An internal mixing-layer like flow forms and grows from the step lip within the original boundary layer. The turbulent structure of the flow evolving downstream, however, does not switch immediately to that of a mixing layer over the entire shear layer thickness. Measurements of mean and fluctuating velocity and temperature fields indicate that the internal layer spreads gradually in the transverse direction while the outer part of the original boundary layer is effectively unperturbed. The results in this paper have not been previously reported.

Jovic, S.↗

Boundary Layer

The fundamental, practically the most important branch of the modern mechanics of a viscous fluid or a gas, is that branch which concerns itself with the study of the boundary layer. The presence of a boundary layer accounts for the origin of the resistance and lift force, the breakdown of the smooth flow about bodies, and other phenomena that are associated with the motion of a body in a real fluid. The concept of boundary layer was clearly formulated by the founder of aerodynamics, N. E. Joukowsky, in his well-known work "On the Form of Ships" published as early as 1890. In his book "Theoretical Foundations of Air Navigation," Joukowsky gave an account of the most important properties of the boundary layer and pointed out the part played by it in the production of the resistance of bodies to motion. The fundamental differential equations of the motion of a fluid in a laminar boundary layer were given by Prandtl in 1904; the first solutions of these equations date from 1907 to 1910. As regards the turbulent boundary layer, there does not exist even to this day any rigorous formulation of this problem because there is no closed system of equations for the turbulent motion of a fluid. Soviet scientists have done much toward developing a general theory of the boundary layer, and in that branch of the theory which is of greatest practical importance at the present time, namely the study of the boundary layer at large velocities of the body in a compressed gas, the efforts of the scientists of our country have borne fruit in the creation of a new theory which leaves far behind all that has been done previously in this direction. We shall herein enumerate the most important results by Soviet scientists in the development of the theory of the boundary layer.

Loitsianskii. L. G.↗

An Approximate Method for Calculation of the Laminar Boundary Layer with Suction for Bodies of Arbitrary Shape

Various ways were tried recently to decrease the friction drag of a body in a flow; they all employ influencing the boundary layer. One of them consists in keeping the boundary layer Laminar by suction; promising tests have been carried out. Since for large Reynolds numbers the friction drag of the laminar boundary layer is much lower than that of the turbulent boundary layer, a considerable saving in drag results from keeping the boundary layer laminar, even with the blower power required for suction taken into account. The boundary layer is kept laminar by suction in two ways: first, by reduction of the thickness of the boundary layer and second, by the fact that the suction changes the form of the velocity distribution so that it becomes more stable, in a manner similar to the change by a pressure drop. There by the critical Reynolds number of the boundary layer (USigma*/V) (sub crit) becomes considerably higher than for the case without suction. This latter circumstance takes full effect only if continuous suction is applied which one might visualize realized through a porous wall. Thus the suction quantities required for keeping the boundary layer laminar become so small that the suction must be regarded as a very promising auxiliary means for drag reduction.

Schlichting, H.↗

Investigation of High-Subsonic Performance Characteristics of a 12 Degree 21-Inch Conical Diffuser, Including the Effects of Change in Inlet-Boundary-Layer Thickness

Investigations were conducted of a 12 degree 21-inch conical diffuser of 2:l area ratio to determine the interrelation of boundary layer growth and performance characteristics. surveys were made of inlet and exit from, longitudinal static pressures were recorded, and velocity profiles were obtained through an inlet Reynolds number range, determined From mass flows and based on inlet diameter of 1.45 x 10(exp 6) to 7.45 x 10(exp 6) and a Mach number range of 0.11 to approximately choking. These investigations were made to two thicknesses of inlet boundary layer. The mean value, over the entire range of inlet velocities, of the displacement thickness of the thinner inlet boundary layer was approximately 0.035 inch and that of the thicker inlet boundary layer was approximately six times this value. The loss coefficient in the case of the thinner inlet boundary layer had a value between 2 to 3 percent of the inlet impact pressure over most of the air-flow range. The loss coefficient with the thicker inlet boundary layer was of the order of twice that of the thinner inlet boundary layer at low speeds and approximately three times at high speeds. In both cases the values were substantially less than those given in the literature for fully developed pipe flow. The static-pressure rise for the thinner inlet boundary layer was of the order of 95 percent of that theoretically possible over the entire speed range. For the thicker inlet boundary layer the static pressure rise, as a percentage of that theoretically possible, ranged from 82 percent at low speeds to 68 percent at high speeds.

Copp, Martin R.↗

Roles of Engineering Correlations in Hypersonic Entry Boundary Layer Transition Prediction

Efforts to design and operate hypersonic entry vehicles are constrained by many considerations that involve all aspects of an entry vehicle system. One of the more significant physical phenomenon that affect entry trajectory and thermal protection system design is the occurrence of boundary layer transition from a laminar to turbulent state. During the Space Shuttle Return To Flight activity following the loss of Columbia and her crew of seven, NASA's entry aerothermodynamics community implemented an engineering correlation based framework for the prediction of boundary layer transition on the Orbiter. The methodology for this implementation relies upon the framework of correlation techniques that have been in use for several decades. What makes the Orbiter boundary layer transition correlation implementation unique is that a statistically significant data set was acquired in multiple ground test facilities, flight data exists to assist in establishing a better correlation and the framework was founded upon state of the art chemical nonequilibrium Navier Stokes flow field simulations. The basic tenets that guided the formulation and implementation of the Orbiter Return To Flight boundary layer transition prediction capability will be reviewed as a recommended format for future empirical correlation efforts. The validity of this approach has since been demonstrated by very favorable comparison of recent entry flight testing performed with the Orbiter Discovery, which will be graphically summarized. These flight data can provide a means to validate discrete protuberance engineering correlation approaches as well as high fidelity prediction methods to higher confidence. The results of these Orbiter engineering and flight test activities only serve to reinforce the essential role that engineering correlations currently exercise in the design and operation of entry vehicles. The framework of information-related to the Orbiter empirical boundary layer transition prediction capability will be utilized to establish a fresh perspective on this role, to illustrate how quantitative statistical evaluations of empirical correlations can and should be used to assess accuracy and to discuss what the authors' perceive as a recent heightened interest in the application of high fidelity numerical modeling of boundary layer transition. Concrete results will also be developed related to empirical boundary layer transition onset correlations. This will include assessment of the discrete protuberance boundary layer transition onset data assembled for the Orbiter configuration during post-Columbia Return To Flight. Assessment of these data will conclude that momentum thickness Reynolds number based correlations have superior coefficients and uncertainty in comparison to roughness height based Reynolds numbers, aka Re(sub k) or Re(sub kk). In addition, linear regression results from roughness height Reynolds number based correlations will be evaluated, leading to a hypothesis that non-continuum effects play a role in the processes associated with incipient boundary layer transition on discrete protuberances.

Campbell, Charles H.↗

Silicon-Based Optical Modulator with Ferroelectric Layer

According to a proposal, a silicon dioxide layer in a high-speed, low-power, silicon- based electro-optical modulator would be replaced by a layer of lead zirconate titanate or other ferroelectric oxide material. The purpose of this modification is to enhance the power performance and functionality of the modulator. In its unmodified form, the particular silicon- based electro-optical modulator is of an advanced design that overcomes the speed limitation of prior silicon-based electro- optical modulators. Whereas modulation frequencies of such devices had been limited to about 20 MHz, this modulator can operate at modulation frequencies as high as 1 GHz. This modulator can be characterized as a silicon-waveguide-based metal oxide/semiconductor (MOS) capacitor phase shifter in which modulation of the index of refraction in silicon is obtained by exploiting the free-charge-carrier-plasma dispersion effect. As shown in the figure, the modulator includes an n-doped crystalline silicon slab (the silicon layer of a silicon- on-insulator wafer) and a p-doped polycrystalline silicon rib with a gate oxide layer (the aforementioned silicon dioxide layer) sandwiched between them. Under accumulation conditions, the majority charge carriers in the silicon waveguide modify the index of refraction so that a phase shift is induced in the optical mode propagating in the waveguide. The advantage of using an MOS capacitor phase shifter is that it is possible to achieve high modulation speed because there are no slow carrier-generation or -recombination processes involved in the accumulation operation. The main advantage of the proposed substitution of a ferroelectric oxide layer for the silicon dioxide layer would arise from the spontaneous polarization effect of the ferroelectric layer: This spontaneous polarization would maintain accumulation conditions in the absence of applied voltage. Consequently, once the device had been switched to a given optical state, it would remain in that state, even in the absence of applied voltage (in other words, even with power turned off). A secondary advantage is that because the ferroelectric layer would have an index of refraction larger than that of silicon dioxide, there could be some reduction of optical losses attributable to fabrication of the modulator

Sheldon, Douglas↗

Hypersonic Boundary Layer Transition Measurements Using NO2 approaches NO Photo-dissociation Tagging Velocimetry

Measurements of instantaneous and mean streamwise velocity profiles in a hypersonic laminar boundary layer as well as a boundary layer undergoing laminar-to-turbulent transition were obtained over a 10-degree half-angle wedge model. A molecular tagging velocimetry technique consisting of a NO2 approaches?NO photo-dissociation reaction and two subsequent excitations of NO was used. The measurement of the transitional boundary layer velocity profiles was made downstream of a 1-mm tall, 4-mm diameter cylindrical trip along several lines lying within a streamwise measurement plane normal to the model surface and offset 6-mm from the model centerline. For laminar and transitional boundary layer measurements, the magnitudes of streamwise velocity fluctuations are compared. In the transitional boundary layer the fluctuations were, in general, 2-4 times larger than those in the laminar boundary layer. Of particular interest were fluctuations corresponding to a height of approximately 50% of the laminar boundary layer thickness having a magnitude of nearly 30% of the mean measured velocity. For comparison, the measured fluctuations in the laminar boundary layer were approximately 5% of the mean measured velocity at the same location. For the highest 10% signal-to-noise ratio data, average single-shot uncertainties using a 1 ?Es and 50 ?Es interframe delay were ~115 m/s and ~3 m/s, respectively. By averaging single-shot measurements of the transitional boundary layer, uncertainties in mean velocity as low as 39 m/s were obtained in the wind tunnel. The wall-normal and streamwise spatial resolutions were 0.14-mm (2 pixel) and 0.82-mm (~11 pixels), respectively. These measurements were performed in the 31-inch Mach 10 Air Wind Tunnel at the NASA Langley Research Center.

Bathel, Brett F.↗

Optimization of Layer Densities for Spacecraft Multilayered Insulation Systems

Numerous tests of various multilayer insulation systems have indicated that there are optimal densities for these systems. However, the only method of calculating this optimal density was by a complex physics based algorithm developed by McIntosh. In the 1970's much data were collected on the performance of these insulation systems with many different variables analyzed. All formulas generated included number of layers and layer density as geometric variables in solving for the heat flux, none of them was in a differentiable form for a single geometric variable. It was recently discovered that by converting the equations from heat flux to thermal conductivity using Fourier's Law, the equations became functions of layer density, temperatures, and material properties only. The thickness and number of layers of the blanket were merged into a layer density. These equations were then differentiated with respect to layer density. By setting the first derivative equal to zero, and solving for the layer density, the critical layer density was determined. Taking a second derivative showed that the critical layer density is a minimum in the function and thus the optimum density for minimal heat leak, this is confirmed by plotting the original function. This method was checked and validated using test data from the Multipurpose Hydrogen Testbed which was designed using McIntosh's algorithm.

Johnson, W. L.↗