Open Access. Powered by Scholars. Published by Universities.®

Mechanical Engineering Commons

Open Access. Powered by Scholars. Published by Universities.®

Mechanical & Aerospace Engineering Theses & Dissertations

Discipline
Keyword
Publication Year

Articles 391 - 420 of 551

Full-Text Articles in Mechanical Engineering

Navier-Stokes Simulation Of Quasi-Axisymmetric And Three-Dimensional Supersonic Vortex Breakdown, Hamdy A. Kandil Apr 1993

Navier-Stokes Simulation Of Quasi-Axisymmetric And Three-Dimensional Supersonic Vortex Breakdown, Hamdy A. Kandil

Mechanical & Aerospace Engineering Theses & Dissertations

Computational simulation of supersonic vortex breakdown is considered for internal and external flow applications. The interaction of a supersonic swirling flow with a shock wave in bounded and unbounded domains is studied. The problem is formulated using the unsteady, compressible, full Navier-Stokes equations which are solved using an implicit, flux-difference splitting, finite-volume scheme. Solutions are obtained for quasi-axisymmetric and three-dimensional flows. The quasi-axisymmetric solutions are obtained by forcing the components of the flowfield vector to be equal on two axial planes, which are in close proximity to each other. For the flow in a bounded domain, a supersonic swirling flow …


Effects Of Radiative Interactions In Molecular Gases Under Local And Nonlocal Thermodynamic Equilibrium Conditions, Manoj Kumar Jha Apr 1993

Effects Of Radiative Interactions In Molecular Gases Under Local And Nonlocal Thermodynamic Equilibrium Conditions, Manoj Kumar Jha

Mechanical & Aerospace Engineering Theses & Dissertations

Basic formulations, analyses, and numerical procedures are presented to investigate radiative heat interactions in diatomic and polyatomic gases under local and nonlocal thermodynamic equilibrium conditions. Essential governing equations are presented for both gray and nongray gases. Information is provided on absorption models, relaxation times, and transfer equations. Radiative flux equations are developed which are applicable under local and nonlocal thermodynamic equilibrium conditions.

The problem is solved for fully developed laminar incompressible flows between two parallel plates under the boundary condition of a uniform surface heat flux. For specific applications, three diatomic and three polyatomic gases are considered. The results are …


An Analysis Of The Response Of A Buckled Beam With Fixed Ends To An Harmonic Base Excitation Using The Finite Element Method, Peter John O'Donoghue Apr 1993

An Analysis Of The Response Of A Buckled Beam With Fixed Ends To An Harmonic Base Excitation Using The Finite Element Method, Peter John O'Donoghue

Mechanical & Aerospace Engineering Theses & Dissertations

A finite element formulation is developed to analyze the large deflection vibration of a buckled isotropic beam. The beam is buckled mechanically by an applied end displacement, and the dynamic loading is derived from the forced periodic base motion. The principle of virtual work is used to derive the system equations of motion. A Newton-Raphson iteration routine modified for the case of an applied end displacement is used with the incremental equations of motion to determine static deflection. The 0l time numerical integration algorithm is applied to the system of equations of motion in physical coordinates to determine dynamic response …


Theoretical Results Supporting The Use Of Passive Damping As Augmentation To The Active Control Of Flexible Structures, Joseph Vincent Harrell Jan 1993

Theoretical Results Supporting The Use Of Passive Damping As Augmentation To The Active Control Of Flexible Structures, Joseph Vincent Harrell

Mechanical & Aerospace Engineering Theses & Dissertations

One challenge of modern control technology is how to control a flexible structure with accuracy, speed, and economy of effort. Controlling a structure with many degrees of freedom by purely active means implies the implementation of inordinate sensors and actuators and creates the need for numerous calculations that must be done instantly. Experiments have shown that practical structures under active control alone can suffer instabilities due to modal vibrations beyond the bandwidth of the active controller. Furthermore, if there is a high degree of model uncertainty, instabilities can be produced by inputs of modal vibrations not occurring in the system …


A Fracture Mechanics Analysis Of A 38 Inch Diameter Hot Pressed Beryllium Telescope Mirror, Phillip L. Brown Oct 1992

A Fracture Mechanics Analysis Of A 38 Inch Diameter Hot Pressed Beryllium Telescope Mirror, Phillip L. Brown

Mechanical & Aerospace Engineering Theses & Dissertations

A 38 inch diameter hot pressed beryllium telescope mirror was analyzed for structural integrity using finite-element modeling and a fracture mechanics analysis program, ZIP3D, developed for three-dimensional (3-0) solid and cracked bodies. The calculation of stress and mixed-mode fracture mechanics parameters for the mirror critical stress area was the objective of the mirror analysis.

Finite-element modeling for ZIP3D solution was performed in two stages. The first stage modeling modeled the mirror configuration on a global basis for calculation of accurate nodal displacements. The second stage modeling reduced the modeled mirror configuration to a set of localized models for calculation of …


Analysis Of Heat Pump Cycles Using Thermodynamic Equations Of State, Stephen C. Ericson Oct 1992

Analysis Of Heat Pump Cycles Using Thermodynamic Equations Of State, Stephen C. Ericson

Mechanical & Aerospace Engineering Theses & Dissertations

In this study, thermal performance of heat pump cycles is analyzed by employing several equations of state. The analysis presented here is particularly useful for applications where refrigerant properties need to be calculated from a thermodynamic equation of state. Such a situation, for example, would arise during a long-term numerical simulation of a solar-assisted heat pump system. In these cases, the direct use of thermodynamic property tables of refrigerants would not be possible in a computational procedure. The approach of using equations of state to predict thermodynamic properties is also useful for applications where new refrigerants or blends of refrigerants …


Transient Simulation Of A Trombe Wall Passive Heating System Coupled To A Thermally Stratified Room, Ali Saleh Kheireddine Oct 1992

Transient Simulation Of A Trombe Wall Passive Heating System Coupled To A Thermally Stratified Room, Ali Saleh Kheireddine

Mechanical & Aerospace Engineering Theses & Dissertations

Thermal performance of a Trombe-wall space heating system has been studied numerically for one-day and multiple-day periods. This study analyzes the system performance by simultaneously considering the coupled nature of the channel convection, wall heat conduction and transient room energy balance. A computer program has been written and implemented to assess the effects of system parameters, weather conditions and peak temperature control strategies on the system performance.

The system performance is influenced by such system parameters as the storage wall thickness, channel width and the temperature stratification in the room. Increasing the wall thickness decreases the temperature swings considerably, but …


Laminar And Turbulent Natural Convection Heat Transfer In Trombe Wall Channels, Tony D. T. Chen Jul 1992

Laminar And Turbulent Natural Convection Heat Transfer In Trombe Wall Channels, Tony D. T. Chen

Mechanical & Aerospace Engineering Theses & Dissertations

The natural convective heat transfer and air movement in a Trombe wall solar passive system has been studied analytically and numerically. Three Trombe wall channel geometries including the parallel channel with axial inlet and exit, parallel channel with side vents and Trombe wall channel coupled to the room have been considered. Several models representing these Trombe wall geometries have been formulated. For the parallel channel with axial inlet and exit geometry, a momentum-integral method has been used to solve parabolic governing equations for two-dimensional laminar flow. This formulation leads to a second order ordinary differential equation for pressure defect in …


Runge-Kutta Upwind Multigrid Multi-Block Three-Dimensional Thin Layer Navier-Stokes Solver, Frank E. Cannizzaro Jul 1992

Runge-Kutta Upwind Multigrid Multi-Block Three-Dimensional Thin Layer Navier-Stokes Solver, Frank E. Cannizzaro

Mechanical & Aerospace Engineering Theses & Dissertations

A state-of-the-art computer code has been developed that incorporates a modified Runge-Kutta time integration scheme, Upwind numerical techniques, Multigrid acceleration, and Multi-block capabilities (RUMM). A three-dimensional thin-layer formulation of the Navier-Stokes equations is employed. For turbulent flow cases, the Baldwin-Lomax algebraic turbulence model is used. Two different upwind techniques are available, van Leer's flux-vector splitting and Roe's flux-difference splitting. Full approximation multigrid plus implicit residual and corrector smoothing were implemented to enhance the rate of convergence. Multi-block capabilities were developed to provide geometric flexibility. This feature allows the developed computer code to accommodate any grid topology or grid configuration with …


Experimental Investigation Of Crossflow Instability Waves In Rotating Disk Flow, Alan E. Blanchard Jul 1992

Experimental Investigation Of Crossflow Instability Waves In Rotating Disk Flow, Alan E. Blanchard

Mechanical & Aerospace Engineering Theses & Dissertations

For the past 40 years, the crossflow (or inflectional) instability associated with rotating disk flow has been generally assumed to be a pure vortex that is fixed with respect to the disk; however, the terms "wave" and "vortex" have also been used interchangeably throughout the computational and experimental literature to describe this instability. Through the use of both hot-wire surveys and smoke-wire flow visualizations in the present investigation, the existence of crossflow instability waves prior to vortex formation has been shown. The results of the hot-wire surveys and smoke-wire flow visualizations for both a "clean" disk case and a disk …


Dynamics And Control Of A Five Degree-Of-Freedom Magnetic Suspension System, Anwar Mohammed Haj Jul 1992

Dynamics And Control Of A Five Degree-Of-Freedom Magnetic Suspension System, Anwar Mohammed Haj

Mechanical & Aerospace Engineering Theses & Dissertations

A large-gap magnetic suspension system with five degrees-of-freedom is presented. The system is multi-input multi-output with coupling between degrees-of-freedom. Simulation was performed on this multi degree-of-freedom system in order to control each degree-of-freedom separately. Two types of controllers are considered by adding white noise to a single degree-of-freedom system in order to test their behavior and determine which is the best choice for the system. The responses of the system are produced in continuous and discrete time where a sample interval and delay time was introduced. Using these responses, a comparison between each degree-of-freedom was made and the maximum value …


Shape Sensitivity Analysis And Optimization Of Skeletal Structures And Geometrically Nonlinear Solids, Ching-Hung Chuang Apr 1992

Shape Sensitivity Analysis And Optimization Of Skeletal Structures And Geometrically Nonlinear Solids, Ching-Hung Chuang

Mechanical & Aerospace Engineering Theses & Dissertations

Formulations and computational schemes for shape design sensitivity analysis and optimization have been developed for both skeletal structures and geometrically nonlinear elastic solids. The continuum approach, which is based on the weak variational form of the governing differential equation and the concept of the material derivative, plays a central role in such a development.

In the first part of this work, the eigenvalue and eigenvector sensitivity equations for skeletal structures are derived with respect to configuration variables of joint and support locations. This derivation is done by the domain method as well as the boundary method. The discrete approach for …


Unsteady Euler And Navier-Stokes Computations Around Oscillating Delta Wing Including Dynamics, Ahmed Abd-El-Bar Ahmed Salman Apr 1992

Unsteady Euler And Navier-Stokes Computations Around Oscillating Delta Wing Including Dynamics, Ahmed Abd-El-Bar Ahmed Salman

Mechanical & Aerospace Engineering Theses & Dissertations

Unsteady flows around rigid or flexible delta wings with and without oscillating leading-edge flaps are considered. These unsteady flow problems are categorized under two classes of problems. In the first class, the wing motion is prescribed a priori and in the second class, the wing motion is obtained as a part of the solution. The formulation of the first class includes either the unsteady Euler or unsteady Navier-Stokes equations for the fluid dynamics and the unsteady linearized Navier-displacement (ND) equations for the grid deformation.

The problem of unsteady transonic flow past a bicircular-arc airfoil undergoing prescribed thickening-thinning oscillation is studied …


Control Of Low-Speed Turbulent Separated Flow Over A Backward-Facing Ramp, John C. Lin Apr 1992

Control Of Low-Speed Turbulent Separated Flow Over A Backward-Facing Ramp, John C. Lin

Mechanical & Aerospace Engineering Theses & Dissertations

The relative performance and flow phenomena associated with several devices for controlling turbulent separated flow were investigated at low speeds. Relative performance of the devices was examined for flow over a curved, backward-facing ramp in a wind tunnel, and the flow phenomena were examined in a water tunnel using dye-flow visualization. Surface static pressure measurements and oil-flow visualization results from the wind tunnel tests indicated that transverse grooves, longitudinal grooves, submerged vortex generators, vortex generator jets (VGJ’s), Viets’ fluidic flappers, elongated arches +a + a (positive angle of attack), and large-eddy breakup devices (LEBU’s) +a + a placed near the …


Characterization Of Trailing Line Vortices, Karla M. Mossi Apr 1992

Characterization Of Trailing Line Vortices, Karla M. Mossi

Mechanical & Aerospace Engineering Theses & Dissertations

Mean velocity surveys were taken for axial or trailing line vortices produced behind a bi-wing vortex generator in a low-speed wind tunnel. Data were taken primarily with a seven-hole probe at a speed of 30m/s; but vortex steadiness was studied using smoke flow visualization over a speed range between 10 and 40 m/s. Various generator center body configurations were studied, with particular interest directed toward how a centerline jet and an asymmetric, unsteady circumferential flow influenced the mean velocity profiles. This study found that a centerline jet did not produce a well-behaved two-cell vortex over the parametric range available in …


Development Of A Model Updating Technique For Controls-Structures Integrated Design Optimization, Angela Georges Apr 1992

Development Of A Model Updating Technique For Controls-Structures Integrated Design Optimization, Angela Georges

Mechanical & Aerospace Engineering Theses & Dissertations

A model updating technique is developed with the intent of providing an efficient means of improving the accuracy of the modelling of the vibrational qualities of a large, flexible structure as it undergoes an integrated controls-structures design optimization. The proposed technique utilizes eigenvalue/eigenvector sensitivity information, two different model reduction methods which assign weights of importance to the vibrational modes of the structure, and modal switching based upon the modal weights determined by the model reduction methods. Issues associated with the numerical implementation of the proposed technique, such as discontinuities induced by modal switching, are also addressed.

The proposed model updating …


Eigenvalue Reanalysis Techniques Using Msc/Nastran, Christopher B. Barton Apr 1992

Eigenvalue Reanalysis Techniques Using Msc/Nastran, Christopher B. Barton

Mechanical & Aerospace Engineering Theses & Dissertations

Reanalysis is a numerical analysis technique that uses results from the total analysis of a design to compute the new response at a reduced computational cost. While there are many techniques available for reanalysis, the first order perturbation technique is investigated in this study. Because of its wide use in industry, NASTRAN finite element analysis software is used for total analysis, and as a basis for error. Since error seems rather unpredictable for this technique, it is sought, in this thesis, to find a reliable index for error so that a designer can know when the technique can be trusted. …


Preliminary Structural Design Of A Lunar Transfer Vehicle Aerobrake, Lance B. Bush Apr 1992

Preliminary Structural Design Of A Lunar Transfer Vehicle Aerobrake, Lance B. Bush

Mechanical & Aerospace Engineering Theses & Dissertations

An aerobrake concept for a Lunar transfer vehicle was weight optimized through the use of the Taguchi design method, structural finite element analyses and structural sizing routines. Six design parameters were chosen to represent the aerobrake structural configuration. The design parameters included honeycomb core thickness, diameter-to-depth ratio, shape, material, number of concentric ring frames, and number of radial frames. Each parameter was assigned three levels. The minimum weight aerobrake configuration resulting from the study was approximately half the weight of the average of all twenty seven experimental configurations. The parameters having the most significant impact on the aerobrake structural weight …


Modal Analysis Of Gyroscopically Coupled Sound-Structure Interaction Problems, Vijay B. Bokil Apr 1992

Modal Analysis Of Gyroscopically Coupled Sound-Structure Interaction Problems, Vijay B. Bokil

Mechanical & Aerospace Engineering Theses & Dissertations

In many practical situations we encounter sound-structure interaction between the fluid in a cavity and the flexible wall on the boundary. The aircraft cabin interacting with its fuselage and double-leaf partition walls interacting with the enclosed cavity are some of the typical examples. The vibratory behavior of the structure and the acoustical pressure inside the cavity are affected by this interaction.

In most of the cases, the governing differential equations for the coupled system can be shown to involve gyroscopic coupling. A special method has been applied to decouple the equations in terms of pairs of generalized coordinates and the …


A Computational Aerodynamic Design Optimization Method Using Sensitivity Analysis, Mohamed E. Eleshaky Apr 1992

A Computational Aerodynamic Design Optimization Method Using Sensitivity Analysis, Mohamed E. Eleshaky

Mechanical & Aerospace Engineering Theses & Dissertations

A new and efficient procedure for aerodynamic shape optimization is presented. The salient lineaments of this procedure are: (1) using a discrete sensitivity analysis approach to determine analytically the aerodynamic sensitivity coefficients; (2) obtaining the flowfield solution either by a computational fluid dynamics (CFD) analysis or, alternatively, by a flowfield extrapolation method which is based on a truncated Taylor's series; (3) defining the aerodynamic shape in such a way that it is not restricted to any class of surfaces and the optimizer automatically shapes the aerodynamic configuration to any arbitrary geometry; and (4) requiring no expertise other than that needed …


Explicit Multistage Schemes For The Solution Of The Three-Dimensional Compressible Euler/Navier-Stokes Equations, Alaa Ahmed Elmiligui Jan 1992

Explicit Multistage Schemes For The Solution Of The Three-Dimensional Compressible Euler/Navier-Stokes Equations, Alaa Ahmed Elmiligui

Mechanical & Aerospace Engineering Theses & Dissertations

The objective of this study was to develop a high-resolution-explicit-multi-block numerical algorithm, suitable for efficient computation of the three-dimensional, time-dependent Euler and Navier-Stokes equations. The resulting algorithm has employed a finite volume approach, using MUSCL-type differencing to obtain state variables at cell interface. Variable interpolations were written in the $\kappa$-scheme formulation. Inviscid fluxes were calculated via Roe's flux-difference splitting, and van Leer's flux-vector splitting techniques, which are considered state of the art. The viscous terms were discretized using a second-order, central-difference operator.

Two classes of explicit time integration has been investigated for solving the compressible inviscid/viscous flow problems--two-stage predictor-corrector schemes, …


Approaches To Control Of The Large Angle Magnetic Suspension Test Fixture, Mehran Ghofrani Jan 1992

Approaches To Control Of The Large Angle Magnetic Suspension Test Fixture, Mehran Ghofrani

Mechanical & Aerospace Engineering Theses & Dissertations

The Large Angle Magnetic Suspension Test Fixture is a five degree-of-freedom system, developed and built at NASA Langley Research Center. It is intended for study of control techniques in magnetic suspension systems with large angular capabilities. In this study, steps have been taken to prove the system in practice, using the existing hardware. A classical control approach, using dual phase advance compensators, is applied in simulation and hardware. A single decoupled degree-of-freedom of the system is stabilized and controlled in simulation. The procedure is then employed for all five degrees-of-freedom. The design and implementation of an analog and a digital …


Experimental Study Of Transition In A Supersonic Free Shear Layer, Rudolph A. King Oct 1991

Experimental Study Of Transition In A Supersonic Free Shear Layer, Rudolph A. King

Mechanical & Aerospace Engineering Theses & Dissertations

The transition behavior of a free shear layer above a cavity with high and low levels of freestream acoustic disturbances has been investigated at Mach 3.5 in the Supersonic Low-Disturbance Pilot Tunnel located at NASA Langley Research Center. Schlieren photography, mean pitot pressure data, and hot-wire surveys were used to detect transition locations. The transition Reynolds numbers obtained were between ReTr = 3.6 x 105 and 5.3 x 105 in good agreement with conventional ("noisy") wind tunnel results. The study indicated that the effects of a decrease in freestream noise on ReTr were minimal but, as …


Finite Element Frequency Domain Solution Of Nonlinear Panel Flutter With Temperature Effects And Fatigue Life Analysis, David Yongxiang Xue Oct 1991

Finite Element Frequency Domain Solution Of Nonlinear Panel Flutter With Temperature Effects And Fatigue Life Analysis, David Yongxiang Xue

Mechanical & Aerospace Engineering Theses & Dissertations

A frequency domain solution method for nonlinear panel flutter with thermal effects using a consistent finite element formulation has been developed. The von Karman nonlinear strain-displacement relation is used to account for large deflections, the quasi-steady first-order piston theory is employed for aerodynamic loading and the quasi-steady thermal stress theory is applied for the thermal stresses with a given change of the temperature distribution, ΔΤ (x, y, z). The equation of motion under a combined thermal-aerodynamic loading can be mathematically separated into two equations and then solved in sequence: (1) thermal-aerodynamic postbuckling and (2) limit-cycle oscillation. The Newton-Raphson iteration technique …


Conceptual Design Of Dual Cam Mechanisms Using Graph Theory, Raymond Lee Mcclendon Jul 1991

Conceptual Design Of Dual Cam Mechanisms Using Graph Theory, Raymond Lee Mcclendon

Mechanical & Aerospace Engineering Theses & Dissertations

This paper presents the theory of graphs and how they can be applied to the conceptual design of dual-cam mechanisms. The graph theory was used to validate the results of previous authors (Hain and Pryor). Several applications of dual cam mechanisms are reviewed.

Algorithms using a process of bit mapping were used to developed all possible non-isomorphic graphs for various classes of mechanisms. A classification method was developed using common occurring substructures found in the graphs. This classification scheme was used for conversion of graphs into schematics. Once the results of Hain were verified, the investigation turned to mechanisms with …


Dynamic Test/Analysis Correlation Using Reduced Analytical Models, Paul E. Mcgowan Jul 1991

Dynamic Test/Analysis Correlation Using Reduced Analytical Models, Paul E. Mcgowan

Mechanical & Aerospace Engineering Theses & Dissertations

Test/analysis correlation is an important aspect of the verification of analysis models which are used to predict on-orbit response characteristics of large space structures. This thesis summarizes the study of using reduced analysis models for performing dynamic test/analysis correlation. The reduced test-analysis model (TAM) has the same number of degrees-of-freedom (DOF) as the number of test measurements. This allows the direct comparison of modal parameters between ground vibration test results and finite element model (FEM) analysis results. The TAM can also be used for pre-test analysis and to aid in the post-test mathematical model updating process.

Herein, four methods for …


Investigation Of Radiative Interactions In Supersonic Chemically Reacting Channel Flows, Almuttil M. Thomas Jul 1991

Investigation Of Radiative Interactions In Supersonic Chemically Reacting Channel Flows, Almuttil M. Thomas

Mechanical & Aerospace Engineering Theses & Dissertations

Analyses and numerical procedures are presented to investigate the radiative interactions of absorbing emitting species in chemically reacting supersonic flow in various ducts. The two-dimensional time dependent Navier-Stokes equations in conjunction with radiative flux equation are used to investigate supersonic flows undergoing finite rate chemical reaction in a hydrogen air system. The specific problem considered is the flow of premixed radiating gas between parallel plates. Specific attention is directed in investigating the radiative contribution of H2O, OH and NO under realistic physical and flow conditions. Results are presented for the radiative flux obtained for different gases and for various combination …


Integrated System Identification And Adaptive State Estimation For Control Of Flexible Space Structures, Chung-Wen Chen Jul 1991

Integrated System Identification And Adaptive State Estimation For Control Of Flexible Space Structures, Chung-Wen Chen

Mechanical & Aerospace Engineering Theses & Dissertations

Accurate state information is crucial for control of flexible space structures in which the state feedback strategy is used. The performance of a state estimator relies on accurate knowledge about both the system and its disturbances, which are represented by system model and noise covariances respectively. For flexible space structures, due to their great flexibility, obtaining good models from ground testing is not possible. In addition, the characteristics of the systems in operation may vary due to temperature gradient, reorientation, and deterioration of material, etc. Moreover, the disturbances during operation are usually not known. Therefore, adaptive methods for system identification …


An Experimental And Analytical Investigation Of Wing Flutter On A Trial Rotor V/Stol Aircraft, David L. Soistmann Jul 1991

An Experimental And Analytical Investigation Of Wing Flutter On A Trial Rotor V/Stol Aircraft, David L. Soistmann

Mechanical & Aerospace Engineering Theses & Dissertations

A wind tunnel model based on a trail rotor vehicle concept was tested for wing flutter in the subsonic flow regime. The wind tunnel model was flutter tested in the forward flight mode only. From this test it was determined that this conceptual vehicle may be susceptible to wing flutter. However, the flutter encountered was not catastrophic but was of the limited amplitude type. The trailing rotor blades were found to have a positive aerodynamic damping effect on flutter. The flutter mechanism for each flutter point was due to a coupling of the first wing bending and first wing torsion …


Simulation And Measurement Of The Radiant Field Produced By Quartz Heating Systems, Travis Lee Turner Jul 1991

Simulation And Measurement Of The Radiant Field Produced By Quartz Heating Systems, Travis Lee Turner

Mechanical & Aerospace Engineering Theses & Dissertations

A method is developed for predicting the radiant heat flux distribution produced by tungsten-filament tubular quartz heating systems with reflectors. The method is an application of Monte Carlo simulation, which takes the form of a ray tracing or random walk scheme. The method is applied to four systems of increasing complexity including; a single unreflected lamp, a single lamp with a flat reflector, a single lamp with a parabolic reflector, and up to six lamps in a six-lamp contoured reflector heating unit. The general methodology of the Monte Carlo method is discussed and applied to the simulation of thermal radiation. …