Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional thermal state”

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 685 records · Page 38

Elemental Fractionation During Rapid Accretion of the Moon Triggered by a Giant Impact

Recently, Ida et al. made an N-body simulation of lunar accretion from a protolunar disk formed by a giant impact. One of their important conclusions is that the accretion time of the Moon is as short as one month. Such rapid accretion is a necessary consequence of the high surface density of a lunar mass disk accreting just beyond the Roche limit (about 3Re); the Safronov accretion time (a few days) is even shorter. The energy of accretion always exceeds the gravitational binding energy of newly arriving matter. Hence, without an energy sink, the accreting body is thermally unstable. For the Earth and other planets, radiation acts as the sink. However, in such a short accretion time, the Moon cannot radiate the accretional energy. Even radiating at a silicate cloudtop temperature of roughly 2000 K, it would take more than 100 yr to radiatively cool the Moon. The plausible alternative heat sinks are heat capacity, latent heat of vaporization, and thermal escape of the gas to space (i.e., hydrodynamic blowoff). The latter becomes plausible for the Moon because the scale height at 2000 K (about 300 km) is a significant fraction of the lunar radius. The early stages of lunar (or "lunatesimal") growth release relatively little energy and can occur simply by heating the material, especially if the accreting material is originally cold. However, the material is unlikely to be cold, because the disk itself is hot and cooling time is long, while the lunar accretion time iss very short. Therefore, the moon is likely to accrete condensed material just after it condenses. Accordingly, the newly accreted material will be on the verge of vaporization and will have very little heat capacity to spare. The immediate heat sink is the latent heat of vaporization. Most of the vapor will escape from the moon, because the thermal energy in the gas can be used to drive escape. However, vaporization is generally incomplete. the latent heat of vaporization exceeds the energy of accretion. Viewed globally, the accretional energy is about half the energy required to vaporize the entire Moon. Thus to first approximation, half of the Moon-forming material can be vaporized and lost during accretion. During this process, we would expect preferential loss of relatively volatile elements. Escape will retard the rate of accretion. To test these ideas, we computed detailed models of the thermal state of the Moon during accretion. We pay special attention to the structure of the silicate atmosphere and its loss rate by calculating the chemical species at equilibrium. We used the PHEQ program which includes 12 elements (H,O,C,Mg,Si,Fe,Ca, Al, Na,Ti, and N.) and 272 compounds (including ionic compounds). Because of the large heats of vaporization and ionization, the adiabatic atmosphere is nearly isothermal and massive escape is expected. The pressure of the atmosphere is determined by the balance between vaporization of a accreting material and escape. If the accretion time is one month, a 0.3 bar atmosphere is expected. Elemental fractionation depends strongly on the temperature of the accreting material. The initial temperature of the material can be estimated from the condition of gravitational instability in the protolunar disk. As shown by Ida et al, accretion starts when gravitational instability occurs when more than 99% of the material condenses. At this point, all of Ca, Al, Si, Mg, and Fe, and 95% of Na (probably K also), are in condensed phases. If the moon is formed from the accretion of such material, volatile elements such as Na, and K are retained by the moon only early in accretion. At later times, K and Na are lost and a fraction of the MG, Si and Fe is lost. However, refractory elements such as Ca and Al are retained and so achieve a mild degree (factor 2) of superabundance.

Abe, Y.↗

Numerical Simulations of Crystal Growth of an Alloy Under Microgravity Conditions

The directional solidification of a dilute binary alloy (Bi-1.0 at.%Sn)is investigated. Results are obtained at a gravity level of I pg. Computations are performed in two dimensions with a fixed, non-uniform grid. The simulation involves a solution of the species concentration equation (modified to account for solute rejection at the interface) and energy equation (modified to account for phase-change) for both the solid and liquid phases, in addition to the constitutive equations for describing convective flow in the melt. The effects of conductive heat transfer in the ampoule and in a capillary tube in the sample are included. To gauge the effects of including this growth capillary tube in the apparatus, simulations both with and without the capillary tube are presented and compared. Fully transient simulations have been performed; no simplifying steady-state approximations are used, however, the influence of solute on the melting temperature at the interface is not included. Both thermal and solutal convective cells are seen to form. Convective velocities are significantly damped inside the capillary, causing less segregation due to convection. As solidification proceeds beyond the capillary tube, longitudinal segregation arises as a result of the change in cross-sectional area of solidifying material. The magnitudes of the velocities in this cell increase significantly once the solid/liquid front passes beyond the end of the capillary tube; this causes a corresponding increase in the level of radial solute segregation in the solidified material.

Simpson, James E.↗

Space processing applications payload equipment study. Volume 2E: Commercial equipment utility

Examination of commercial equipment technologies revealed that the functional performance requirements of space processing equipment could generally be met by state-of-the-art design practices. Thus, an apparatus could be evolved from a standard item or derived by custom design using present technologies. About 15 percent of the equipment needed has no analogous commercial base of derivation and requires special development. This equipment is involved primarily with contactless heating and position control. The derivation of payloads using commercial equipment sources provides a broad and potentially cost-effective base upon which to draw. The derivation of payload equipment from commercial technologies poses other issues beyond that of the identifiable functional performance, but preliminary results on testing of selected equipment testing appear quite favorable. During this phase of the SPA study, several aspects of commercial equipment utility were assessed and considered. These included safety, packaging and structural, power conditioning (electrical/electronic), thermal and materials of construction.

Smith, A. G.↗

Advanced Modular Power System Electronics Enclosure Thermal Testing

An analysis was set up to model the temperature of the advanced modular power system (AMPS) power distribution cards when installed within the electronics enclosure case. The analysis was used to determine the steady-state temperature distribution of the cards within the case. To verify the analysis, an experiment was set up and conducted to simulate the operation of the cards within the enclosure. Four tests were conducted. The tests varied the position of the cold plate and evaluated the use of a thermal compound to reduce the contact resistance between the joints within the thermal path between the cards and the cold plate. Three of the four cases examined showed very good agreement between the analysis and the experiment with a less than 1-percent variation in the predicated temperatures determined through the analysis and the experimentally derived temperatures. In the remaining case, the difference between the analysis and experiment was approximately 12 percent. Both the experiment and analysis showed that the modular power conditioning cards can be maintained within their desired maximum operating temperature range of 40 to 45 °C through thermal conduction to a cold plate when operating with their estimated maximum heat output of 16 W per card.

Colozza, Anthony J.↗

Nonlinear Evolution of Instabilities in a Laminar Separation Bubble at a Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

Hypersonic flow↗

Nonlinear Evolution of Instability in a Laminar Separation Bubble at Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

High speed flow↗

Laser Forming of Sheet Metal for In-Space Manufacturing Applications

As the in-space economy matures, the migration of typically Earth-bound manufacturing methods, such as metal forming and joining, up to space will be required to meet the increasing demand for refueling and supply depots, habitats, laboratories, lunar surface structures, and sites for deep-space exploration shipbuilding. Only in this way can the magnitude of individual structures begin to exceed what can be launched, deployed, docked, or inflated and allow true economies of scale to be unlocked in space. However, typical terrestrial methods of forming and shaping sheet metal require massive machinery with the ability to deflect reaction forces to a large mass in a fixed reference frame (Earth). To both limit up-mass required for manufacture and to mitigate the issues of reaction forces in the relative reference frames of space, a non-contact, laser-based method of forming metal sheet is proposed. NASA Marshall Space Flight Center (MSFC), in cooperation with DARPA and the University of Florida (UF), have begun investigation of non-contact methods of laser forming sheet metal in thermal vacuum (TVAC) conditions to great success. To date, several grades of aluminum, stainless steel, and titanium coupons have been formed via laser-induced internal stresses in these metals at both ambient atmosphere and vacuum conditions. Temperatures in these tests ranged from -120 to 60 C, and though the total energy required to induce strain does vary with workpiece temperature and environmental pressure, the method remains predictable and controllable for the materials investigated. It is believed that this method will revolutionize manufacturing in-space as it not only continues to build upon the capabilities of the laser (sensing, marking, cutting, drilling, joining, and now forming) as an all-in-one tool, but also reduces the need to counteract forming forces as would be required in mechanical bending. In theory, all required stress to induce bending is generated through thermal gradient and strain within the workpiece, meaning that the need for reaction control and fuel expense on-orbit will be limited. Archimedes once stated, “Give me a lever long enough and a fulcrum on which to place it, and I shall move the world”. With a laser as our lever and our workpiece its own fulcrum, we can move the world’s manufacturing into space.

Laser↗

Bond Line Thickness Estimation in Composite Structures Using Multiple Inspection Techniques

Imaging and other nondestructive evaluation techniques are commonly used for material characterization and defect recognition in safety critical aerospace applications, with data fusion providing the framework for uncertainty quantification in these contexts. Most commonly, forward physics-based modeling predicts the response conditioned on material properties and defect assumptions, and probabilistic methods are used to infer the hidden state of subject of the inspection from a combination of prior information, likelihoods, and inspection data. In this paper Bayesian methods are used to estimate bond thickness in lap joints comprised of aluminum adherends using a combination of infrared thermography, digital radiography, and ultrasound. The accuracy of the fused estimates are validated against data generated from synthetic specimens, and by comparison against high resolution X-ray computed tomography inspections of built specimens.

thermal nondestructive evaluation↗

Results from the vehicle charging and potential experiment on STS-3

The Vehicle Charging and Potential Experiment formed part of the scientific experiments comprising the NASA Office of Space Sciences-1 pallet flown at an altitude of 250 km and an orbital inclination of 40 deg on the third flight of the Space Shuttle Orbiter in March 1982. Details of the objectives, the instrumentation of the experiment, and some preliminary results from the individual instruments are presented. The experiment studied passive vehicle charging using thermal plasma probes, and it is concluded that apart from moderate v x B (L)-induced electrical potential resulting from the large size and restricted plasma contact area, the Space Shuttle Orbiter behaved similarly to smaller spacecraft orbiting in the same ionospheric region. However, an unexpected passive orbital result was the detection of enhanced plasma density in the vicinity of the Orbiter under restricted orbital conditions, apparently composed of contaminant molecular species and with an electron component having a temperature considerably higher than the ambient ionosphere electrons. The active vehicle charging experiments were performed by using a 100 mA/1 keV electron beam emitted out of the payload bay. Charging measurements made during short electron beam pulse emissions provide an upper limit of about 1 mF for the Orbiter capacitance. Steady-state charging measurements showed that for the vehicle attitudes covered by the mission, the induced electrical potential was normally only a few volts, except during some nighttime conditions when it rose above 40 V.

Banks, P. M.↗

Fatigue life prediction modeling for turbine hot section materials

A major objective of the fatigue and fracture efforts under the Hot Section Technology (HOST) program was to significantly improve the analytic life prediction tools used by the aeronautical gas turbine engine industry. This was achieved in the areas of high-temperature thermal and mechanical fatigue of bare and coated high-temperature superalloys. The cyclic crack initiation and propagation resistance of nominally isotropic polycrystalline and highly anisotropic single crystal alloys were addressed. Life prediction modeling efforts were devoted to creep-fatigue interaction, oxidation, coatings interactions, multiaxiality of stress-strain states, mean stress effects, cumulative damage, and thermomechanical fatigue. The fatigue crack initiation life models developed to date include the Cyclic Damage Accumulation (CDA) and the Total Strain Version of Strainrange Partitioning (TS-SRP) for nominally isotropic materials, and the Tensile Hysteretic Energy Model for anisotropic superalloys. A fatigue model is being developed based upon the concepts of Path-Independent Integrals (PII) for describing cyclic crack growth under complex nonlinear response at the crack tip due to thermomechanical loading conditions. A micromechanistic oxidation crack extension model was derived. The models are described and discussed.

Halford, G. R.↗

Fatigue life prediction modeling for turbine hot section materials

A major objective of the fatigue and fracture efforts under the NASA Hot Section Technology (HOST) program was to significantly improve the analytic life prediction tools used by the aeronautical gas turbine engine industry. This was achieved in the areas of high-temperature thermal and mechanical fatigue of bare and coated high-temperature superalloys. The cyclic crack initiation and propagation resistance of nominally isotropic polycrystalline and highly anisotropic single crystal alloys were addressed. Life prediction modeling efforts were devoted to creep-fatigue interaction, oxidation, coatings interactions, multiaxiality of stress-strain states, mean stress effects, cumulative damage, and thermomechanical fatigue. The fatigue crack initiation life models developed to date include the Cyclic Damage Accumulation (CDA) and the Total Strain Version of Strainrange Partitioning (TS-SRP) for nominally isotropic materials, and the Tensile Hysteretic Energy Model for anisotropic superalloys. A fatigue model is being developed based upon the concepts of Path-Independent Integrals (PII) for describing cyclic crack growth under complex nonlinear response at the crack tip due to thermomechanical loading conditions. A micromechanistic oxidation crack extension model was derived. The models are described and discussed.

G R Halford↗

Numerical simulation of heat and mass transport during space crystal growth with MEPHISTO

The MEPHISTO space experiments are collaborative United States and French investigations aimed at understanding the fundamentals of crystal growth. Microgravity experiments were conducted aboard the USMP-1 and -2 missions on STS-52 and 62 in October 1992 and March 1994 respectively. MEPHISTO is a French designed and built Bridgman type furnace which uses the Seebeck technique to monitor the solid/liquid interface temperature and Peltier pulsing to mark the location and shape of the solid/liquid interface. In this paper the Bridgman growth of Sn-Bi and Bi-Sn under terrestrial and microgravity conditions is modeled using the finite element code, FIDAP*. The numerical model considers fully coupled heat and mass transport, fluid motion and solid/liquid phase changes in the crystal growth process. The primary goals of this work are: to provide a quantitative study of the thermal buoyancy-induced convection in the melt for the two flight experiments; to compare the vertical and horizontal growth configurations and systematically evaluate the effects of various gravity levels on the solute segregation. Numerical results of the vertical and horizontal Bridgman growth configurations are presented.

Yao, Minwu↗

Durability, Performance, and Emission of Diesel Engines Using Carbon Fiber Piston and Liner

This report summarizes the research conducted by NC State University in investigating the durability, performance and emission of a carbon fiber piston and liner in our single cylinder research Diesel engine. Both the piston and liner were supplied to NC State University by NASA LaRC and manufactured by C-CAT under a separate contract to NASA LaRC. The carbon-carbon material used to manufacture the piston and liner has significantly lower thermal conductivity, coefficient of thermal expansion, and superior strength characteristics at elevated temperatures when compared to conventional piston materials such as aluminum. The results of the carbon-carbon fiber piston testing were compared to a baseline configuration, which used a conventional aluminum piston in a steel liner. The parameters measured were the brake specific fuel consumption, ignition delay, frictional horsepower, volumetric efficiency, and durability characteristics of the two pistons. Testing was performed using a naturally aspirated Labeco Direct Injection single cylinder diesel engine. Two test cases were performed over a range of loads and speeds. The fixed test condition between the aluminum and carbon-carbon piston configurations was the brake mean effective pressure. The measured data was the fuel consumption rate, volumetric efficiency, load, speed, cylinder pressure, needle lift, and exhaust gas temperature. The cylinder pressure, and fuel consumption, exhaust gas temperature, and needle lift were recorded using a National Instruments DAQ board and a PC. All test cases used Diesel no. 2 for fuel.

Afify, E. M.↗

Spatial and Temporal Varying Thresholds for Cloud Detection in Satellite Imagery

A new cloud detection technique has been developed and applied to both geostationary and polar orbiting satellite imagery having channels in the thermal infrared and short wave infrared spectral regions. The bispectral composite threshold (BCT) technique uses only the 11 micron and 3.9 micron channels, and composite imagery generated from these channels, in a four-step cloud detection procedure to produce a binary cloud mask at single pixel resolution. A unique aspect of this algorithm is the use of 20-day composites of the 11 micron and the 11 - 3.9 micron channel difference imagery to represent spatially and temporally varying clear-sky thresholds for the bispectral cloud tests. The BCT cloud detection algorithm has been applied to GOES and MODIS data over the continental United States over the last three years with good success. The resulting products have been validated against "truth" datasets (generated by the manual determination of the sky conditions from available satellite imagery) for various seasons from the 2003-2005 periods. The day and night algorithm has been shown to determine the correct sky conditions 80-90% of the time (on average) over land and ocean areas. Only a small variation in algorithm performance occurs between day-night, land-ocean, and between seasons. The algorithm performs least well. during he winter season with only 80% of the sky conditions determined correctly. The algorithm was found to under-determine clouds at night and during times of low sun angle (in geostationary satellite data) and tends to over-determine the presence of clouds during the day, particularly in the summertime. Since the spectral tests use only the short- and long-wave channels common to most multispectral scanners; the application of the BCT technique to a variety of satellite sensors including SEVERI should be straightforward and produce similar performance results.

Jedlovec, Gary↗

A Perspective on Multi-Source Energy Harvesting Technology Using Piezoelectric and Thermoelectric Materials

Energy harvesting using various locally available energy sources such as vibration energy, heat, sound, or magnetic field have become attractive topics for supplying power to modular electronic devices making them run independently in extreme environments. In this paper, we will be discussing the perspectives on thermoelectric (TE) and piezoelectric materials and devices, and then the concept of multi-source energy harvester using piezoelectric and thermoelectric devices and integration of them into a reliable and independent power source. TE materials having low thermal conductivity and high figure-of-merit ( zT ) are developed to convert even a small temperature gradient efficiently into electrical energy with the state-of-the-art conversion efficiency of ~15% and output power of ~56 W from single device. The piezoelectric device architecture is configured using high performance piezoelectric ceramics (Cu-Mn-PIN-PMN-PT). These ceramics exhibit high piezoelectric coefficient with high mechanical quality factor and low dielectric loss factor. Using these piezoelectric materials, power density as high a 2 mW/cm 2 is demonstrated in 1-1.5 g vibration environments. The piezoelectric device is attached on the surface of TE module to capture both the vibration and thermal energy sources to realize dual mode energy harvester. The multi-energy transfer strategy opens opportunities for a future generation of wireless and modular electronic devices. These devices would be useful in powering wearable electronic devices, micro sensor chargers, etc. in extreme environmental conditions using body heat/thermal sources and induced motion/vibrations.

Piezoelectric↗

On the internal structure of the major satellites of the outer planets

The structures and thermal evolutions of the large icy-satellites of the outer solar system are considered. It is shown (for bodies comparable in size and mass to the Galilean satellites, having sizeable mass fraction of H2O, and with meteoritic abundances of radioactive materials contained within their silicate fractions) that the crust of solid ice over a liquid mantle predicted by conductive heat-transfer calculations is unstable to large-scale solid-state convection. For appropriate material parameters, convective heat-transfer rates are sufficient to freeze a large liquid mantle on a time scale that is short compared to the lifetime of the body. It is also concluded that the ice layer is convecting at the present time. A reevaluation of previous work, using improved values for material parameters and boundary conditions, reverses earlier conclusions and implies a rigid outer crust with resulting long-term stability of surface features to creep deformation. The combination of a rigid crust with active internal convection presents the additional possibility of surface features that are produced and maintained by dynamic internal processes.

Reynolds, R. T.↗

Thermal and electric properties of Nd(1.85)Ce(0.15)CuO(4-y) and Pr(1.85)Ce(0.15)CuO(4-y)

Electric resistivity, magnetic susceptibility, thermoelectric power, and Hall coefficient of Nd(1.85)Ce(0.15)CuO(4-y) and Pr(1.85)Ce(0.15)CuO(4-y) whose onset temperature of the superconductivity are 24 and 23 K were measured. Experimental results show many interesting features. In particular, the Hall coefficients are negative and relatively flat as a function of temperature. However, the temperature dependence of the thermoelectric power (TEP) for these two samples shows the positive sign for both samples in contrast to the previous results. Moreover, TEP for both samples remains flat in the normal state below 250 K, but decreases rapidly above 250 K. TEP of only Pr(1.85)Ce(0.15)CuO(4-y) shows a peak near 50 K. Finally, onset temperatures of sudden drop of TEP are higher than those of resistance drop. The physical properties of these samples produced at different conditions such as different heat treatment temperatures, atmospheres were also measured. TEP and resistance measurement show that oxygen deficiency is essential to produce better superconducting samples. Correlation between TEP and superconductivity for these different samples are discussed.

Lim, Z. S.↗

Bridgman crystal growth

The objective of this theoretical research effort was to improve the understanding of the growth of Pb(x)Sn(1-x)Te and especially how crystal quality could be improved utilizing the microgravity environment of space. All theoretical growths are done using the vertical Bridgman method. It is believed that improved single crystal yields can be achieved by systematically identifying and studying system parameters both theoretically and experimentally. A computational model was developed to study and eventually optimize the growth process. The model is primarily concerned with the prediction of the thermal field, although mass transfer in the melt and the state of stress in the crystal were of considerable interest. The evolution is presented of the computer simulation and some of the important results obtained. Diffusion controlled growth was first studied since it represented a relatively simple, but nontheless realistic situation. In fact, results from this analysis prompted a study of the triple junction region where the melt, crystal, and ampoule wall meet. Since microgravity applications were sought because of the low level of fluid movement, the effect of gravitational field strength on the thermal and concentration field was also of interest. A study of the strength of coriolis acceleration on the growth process during space flight was deemed necessary since it would surely produce asymmetries in the flow field if strong enough. Finally, thermosolutal convection in a steady microgravity field for thermally stable conditions and both stable and unstable solutal conditions was simulated.

Carlson, Frederick↗