Vapor chambers with hydrophobic and biphilic evaporators in moderate to high heat flux applications
Not Available
SEARCH · Engineering Papers
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.
Not Available
Not Available
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The Hanford Waste Treatment and Immobilization Plant (WTP) Low Activity Waste (LAW) vitrification facility will generate an aqueous condensate recycle stream, LAW Off-Gas Condensate, from the off-gas system. The baseline plan for disposition of this stream is to send it to the WTP Pretreatment Facility, where it will be blended with LAW, concentrated by evaporation and recycled to the LAW vitrification facility again. Alternate disposition of this stream would eliminate recycling of problematic components, and would enable de-coupled operation of the LAW melter and the Pretreatment Facilities. Eliminating this stream from recycling within WTP would also decrease the LAW vitrification mission duration and quantity of canistered glass waste forms. This LAW Off-Gas Condensate stream contains components that are volatile at melter temperatures and are problematic for the glass waste form. Because this stream recycles within WTP, these components accumulate in the Condensate stream, exacerbating their impact on the number of LAW glass containers that must be produced. Approximately 32% of the sodium in Supplemental LAW comes from glass formers used to make the extra glass to dilute the halides to be within acceptable concentration ranges in the LAW glass. Diverting the stream reduces the halides in the recycled Condensate and is a key outcome of this work. Additionally, under possible scenarios where the LAW vitrification facility commences operation prior to the WTP Pretreatment facility, identifying a disposition path becomes vitally important. This task examines the impact of potential future disposition of this stream in the Hanford tank farms, and investigates auxiliary evaporation to enable another disposition path. Unless an auxiliary evaporator is used, returning the stream to the tank farms would require evaporation in the 242-A evaporator. This stream is expected to be unusual because it will be very high in corrosive species that are volatile in the melter (chloride, fluoride, sulfur), will have high ammonia, and will contain carryover particulates of glass-former chemicals. These species have potential to cause corrosion of tanks and equipment, precipitation of solids, release of ammonia gas vapors, and scale in the tank farm evaporator. Routing this stream to the tank farms does not permanently divert it from recycling into the WTP, only temporarily stores it prior to reprocessing. Testing is normally performed to demonstrate acceptable conditions and limits for these compounds in wastes sent to the tank farms. The primary parameter of this phase of the test program was measuring the formation of solids during evaporation in order to assess the compatibility of the stream with the evaporator and transfer and storage equipment. The origin of this LAW Off-Gas Condensate stream will be the liquids from the Submerged Bed Scrubber (SBS) and the Wet Electrostatic Precipitator (WESP) from the LAW facility melter offgas system. The stream is expected to be a dilute salt solution with near neutral pH, and will likely contain some insoluble solids from melter carryover. The soluble components are expected to be mostly sodium and ammonium salts of nitrate, chloride, and fluoride. This stream has not been generated yet, and, thus, the composition will not be available until the WTP begins operation, but a simulant has been produced based on models, calculations, and comparison with pilot-scale tests. This report discusses results of evaporation testing of the simulant. Two conditions were tested, one with the simulant at near neutral pH, and a second at alkaline pH. The neutral pH test is comparable to the conditions in the Hanford Effluent Treatment Facility (ETF) evaporator, although that evaporator operates at near atmospheric pressure and tests were done under vacuum. For the alkaline test, the target pH was based on the tank farm corrosion control program requirements, and the test protocol and equipment was comparable to that used for routine evaluation of feed compatibility studies for the 242-A evaporator. One of the radionuclides that is volatile in the melter and expected to be in high concentration in this LAW Off-Gas Condensate stream is Technetium-99 (99Tc). Technetium will not be removed from the aqueous waste in the Hanford WTP, and will primarily end up immobilized in the LAW glass by repeated recycle of the off-gas condensate into the LAW melter. Other radionuclides that are also expected to be in appreciable concentrations in the LAW Off-Gas Condensate are 129I, 90Sr, 137Cs, and 241Am. The concentrations of these radionuclides in this stream will be much lower than in the LAW, but they will still be higher than limits for some of the other disposition pathways currently available. At this time, these scoping tests did not evaluate the partitioning of the radionuclides to the evaporator condensate, since ample data are available separately from other experience in the DOE complex. Results from the evaporation testing show that the neutral SBS simulant first forms turbidity at ~7.5X concentration, while the alkaline-adjusted simulant became turbid at ~3X concentration. The major solid in both cases was Kogarkoite, Na3FSO4. Sodium and lithium fluorides were also detected. Minimal solids were formed in the evaporator bottoms until a substantial fraction of liquid was removed, indicating that evaporation could minimize storage volume issues. Achievable concentration factors without significant insoluble solids were 17X at alkaline pH, and 23X at neutral pH. In both runs, significant ammonia carried over and was captured in the condenser with the water condensate. Results also indicate that with low insoluble solids formation in the initial testing at neutral pH, the use of Reverse Osmosis is a potential alternate method for concentrating the solution, although an evaluation is needed to identify equipment that can tolerate insoluble solids. Most of the ammonia remains in the evaporator bottoms during the neutral pH evaporation, but partitions to the condensate during alkaline evaporation. Disposition of both streams needs to consider the management of ammonia vapor and its release. Since this is an initial phase of testing, additional tasks related to evaporation methods are expected to be identified for development. These tasks likely include evaluation and testing of composition variability testing and evaluations, corrosion and erosion testing, slurry storage and immobilization investigations, and evaporator condensate disposition.
Evaporation of Cu or Sn from liquid iron alloys containing C and S was experimentally investigated. The initial C concentration, [pct C] 0 , in the liquid alloy was varied from zero to C saturation, and the evaporation temperature was varied from 1513 K to 1773 K (1240 °C to 1500 °C). Along with the report by one of the present authors, the evaporation mechanism of Cu and Sn from liquid Fe-C-S alloy is proposed, after a modification from the previous mechanism. It was proposed that Cu and Sn evaporate as Cu(g) and Sn(g) and also evaporate as CuS(g) and SnS(g), which are more volatile species. Therefore, availability of S in the alloy affects the overall evaporation rate of Cu and Sn. At the same time, C in the alloy also forms volatile carbosulfides CS(g) and CS 2 (g), thereby competing with Cu and Sn. Moreover, C increases the activity coefficients of Cu, Sn, and S. This increases the thermodynamic driving force for the formation of CuS(g) and SnS(g). Therefore, increasing [pct C] partly accelerates the evaporation rate of Cu and Sn by increasing the activity coefficient but partly decelerates the evaporation rate by lowering the available S content. S partly accelerates the evaporation rate by increasing the available S for the sulfide gas species but partly decelerates the evaporation rate due to the surface poisoning effect. Increasing the reaction temperature increases the overall evaporation rate. All these facts were taken into account in order to develop an evaporation rate model. This model was extended from the present authors’ previous one by taking into account (1) CS(g), S(g), and CS{ 2 (g) (therefore, the following species were considered as dominant evaporating species: Cu(g), CuS(g), Sn(g), SnS(g), S(g), CS(g), and CS 2 (g)); (2) the effect of C and temperature on the activity coefficients of Cu, Sn, and S; (3) the effect of C and temperature on the density of the liquid alloy; and (4) the effect of temperature on the S adsorption coefficient. This revised evaporation model was used in order to explain the experimental data, and it showed good agreement. In particular, it was found that the temperature showed a significant effect on the evaporation rate, and the effect of temperature and C content on the activity coefficients of Cu, Sn, and S also significantly affected the evaporation rate. The chemical reaction rate constant of the individual evaporation reaction (k i } R ) and residual rate constant (k i r ) could be obtained as a function of temperature. The activation energy of each evaporation reaction was derived and discussed. The evaporation rate model can be applied in order to predict the content of Cu and Sn remaining in liquid iron under various conditions of temperature and [pct C].
Interfacial solar evaporation using three-dimensional evaporator materials has shown promise to achieve high evaporation rates exceeding the photothermal limit for water treatment and resource recovery processes. However, challenges remain in optimizing the material geometry to balance the vertical capillary uptake of water and the rate of evaporation to achieve both stable and high evaporation rates. Delignified wood can serve as a highly porous and hydrophilic substrate material to achieve high evaporation rates. In this work, we designed a suspended biochar-coated delignified wood rod evaporator and investigated the influence of its geometric parameters on its evaporation performance. The height of the evaporator exposed to air for evaporation was observed to be an important parameter in governing the evaporation rate, as shorter evaporators limited the available surface area for evaporation, but taller evaporators could not achieve sufficient rates of capillary uptake to maintain stable evaporation rates over 24 h under 1-sun illumination (1 kW m –2 ). In conclusion, the exposed height of the evaporator was optimized to achieve an average hourly evaporation rate of 3.05 ± 0.23 kg m –2 h –1 , which could be maintained in saline solutions containing up to 5 wt % NaCl.
During the last few years lot of research has been conducted to investigate and predict the Crud deposition in PWRs. One of the common findings of various experimental and computational studies is that thickness or amount of crud deposits is directly related to evaporative flux maps. Evaporation driven deposition profile on the surfaces takes different forms, i.e. from ring to disk shapes during initial stages and in the form of porous to impervious continuous layers covering the hot fuel surface at the later stages. The amount and thickness of these deposits although are directly related to evaporative flux, the relationship is much complex and can be dependent upon several other factors. A popular problem relating evaporative flux to particulate deposition has been well studied in the form of drying of sessile colloidal drops. It was proposed that, this phenomenon is due to a geometrical constraint: the free surface, constrained by a pinned contact line squeezes the fluid outwards to compensate for evaporative losses and thereby transporting particles to the contact line. Later on a quantitative relationship between evaporative flux and deposition rate was proposed and validated by the same group, which showed power-law relationship. Deposition studies conducted with controlled concentration of particles with subcooled nucleate boiling also showed a power-law relationship rather than a simpler linear relationship whereas physics of evaporation would tend to concentrate particles at linear rate. In a typical subcooled nucleate boiling system, bubble growth and departure time is of the order 1-2 ms, measured experimentally. Diffusion time constant for dilute spherical metallic colloids such as crud particulates, generally present as particulate in process piping, range from 0.0025-2.5 sec in water for a size range 0.1-1.0 μm. Colloidal behavioral characteristics during the growth of bubbles can be computed with escape time constant at collision, which is dependent on effective Hamaker constant. For minimum effective Hamaker constant of Hematite is -2.5 kT and diffusion time constant τ = 0.0025 sec, escape time constant is computed as Kramer rule, τexp[A/kT] magnitude is 0.03 sec or 30 ms. This clearly indicates that evaporation will cause particles to come together and they cannot escape at a fast enough rate which would lead them to form agglomerates. Therefore both effects of evaporative flux in conjunction with the agglomeration kinetics must be considered for modeling the crud deposition if the evaporative flux is known. In this paper, we propose a phenomenological model of colloidal transport and aggregation with evaporative flux. The dynamic of particles interactions while drying of sessile colloidal drops is developed by coupling general transport equation with Smoluchowshki Coagulation equation. Crud deposition on the nuclear fuel has been a serious problem with severe consequences in nuclear power plants. In several studies it has been consistently reported that crud deposition is directly related to evaporative flux on the nuclear fuel surfaces. However, it has not been established how to quantify the deposition rates of crud type particulates if the evaporative flux is known. This complex multiphysics problem is similar to coffee-ring effect, where the governing physics is understood more clearly. In this work, a simplified geometry is used to model the effects of evaporative flux acting at a boundary on the colloidal concentration. The effect of this condition on change in concentration is shown numerically. Further, agglomeration is suspected to play a significant role in high evaporative flux such as bubble nucleation. The coagulation kinetics is modeled using uniform particle distribution and the effective total concentration of colloids including original particles and agglomerates is studied. Once this model is validated with a simple separate effect experiments, it can be integrated to large scale PWR models of evaluating evaporative flux and chemistry models developed by Electrical Power Research Institute (EPRI). (authors)
The rate of 2H evaporator scale solids dissolution in dilute nitric acid has been experimentally evaluated under laboratory conditions in the SRNL shielded cells. The 2H scale sample used for the dissolution study came from the bottom of the evaporator cone section and the wall section of the evaporator cone. The accumulation rate of aluminum and silicon, assumed to be the two principal elemental constituents of the 2H evaporator scale aluminosilicate mineral, were monitored in solution. Aluminum and silicon concentration changes, with heating time at a constant oven temperature of 90 deg C, were used to ascertain the extent of dissolution of the 2H evaporator scale mineral. The 2H evaporator scale solids, assumed to be composed of mostly aluminosilicate mineral, readily dissolves in 1.5 and 1.25 M dilute nitric acid solutions yielding principal elemental components of aluminum and silicon in solution. The 2H scale dissolution rate constant, based on aluminum accumulation in 1.5 and 1.25 M dilute nitric acid solution are, respectively, 9.21E-04 ± 6.39E-04 min{sup -1} and 1.07E-03 ± 7.51E-05 min{sup -1}. Silicon accumulation rate in solution does track the aluminum accumulation profile during the first few minutes of scale dissolution. It however diverges towards the end of the scale dissolution. This divergence therefore means the aluminum-to-silicon ratio in the first phase of the scale dissolution (non-steady state conditions) is different from the ratio towards the end of the scale dissolution. Possible causes of this change in silicon accumulation in solution as the scale dissolution progresses may include silicon precipitation from solution or the 2H evaporator scale is a heterogeneous mixture of aluminosilicate minerals with several impurities. The average half-life for the decomposition of the 2H evaporator scale mineral in 1.5 M nitric acid is 12.5 hours, while the half-life for the decomposition of the 2H evaporator scale in 1.25 M nitric acid is 10.8 hours. Based on averaging the two half-lives from the 2H scale acid dissolution in 1.25 and 1.5 M nitric acid solutions, a reasonable half-live for the dissolution of 2H scales in dilute nitric acid is 11.7 ± 1.3 hours. The plant operational time for chemically cleaning (soaking) the 2H evaporator with dilute nitric acid is 32 hours. It therefore may require about 3 half-lives or less to completely dissolve most of the scales in the Evaporator pot which come into contact with the dilute nitric acid solution. On a mass basis, the Al-to-Si ratio for the scale dissolution in 1.5 M nitric acid averaged 1.30 ± 0.20 and averaged 1.18 ± 0.10 for the 2H scale dissolution in 1.25 M nitric acid. These aluminum-to-silicon ratios are in fairly good agreement with ratios from previous studies. Therefore, there is still more aluminum in the 2H evaporator scales than silicon which implies that there are no significant changes in scale properties which will exclude nitric acid as a viable protic solvent for aluminosilicate scale buildup dissolution from the 2H evaporator. Overall, the monitoring of the scale decomposition reaction in 1.25 and 1.5 M nitric acid may be better ascertained through the determination of aluminum concentration in solution than monitoring silicon in solution. Silicon solution chemistry may lead to partial precipitating of silicon with time as the scale and acid solution is heated.
In this study, we demonstrate a 3-dimensional graphene oxide (3D GO) stalk that operates near the capillary wicking limit to achieve an evaporation flux of 34.7 kg m -2 h -1 under 1 sun conditions (1 kW/m 2 ). This flux represents nearly a 100 times enhancement over a conventional solar evaporation pond. Interfacial solar evaporation traditionally uses 2D evaporators to vaporize water using sunlight, but their low evaporative water flux limits their practical applicability for desalination. Some recent studies using 3D evaporators demonstrate potential for more efficient water transfer, but the flux improvement has been marginal because of a low evaporation area index (EAI), which is defined as the ratio of the total evaporative surface area to the projected ground area. By using a 3D GO stalk with an ultrahigh EAI of 70, we achieved nearly a 20-fold enhancement over a 2D GO evaporator. The 3D GO stalk also exhibited additional advantages including omnidirectional sunlight utilization, a high evaporation flux under dark conditions from more efficient utilization of ambient heating, a dramatic increase of the evaporation rate by introducing wind, and scaling resistance in evaporating brines with a salt content of up to 17.5 wt %. This performance makes the 3D GO stalk well suited for the development of a low-cost, reduced footprint technology for zero liquid discharge in brine management applications.
Knowing how actual evaporation is down-regulated from potential evaporation during periods with soil moisture deficits is one of the greatest challenges towards computing evaporation everywhere on a regular basis. We propose the hypothesis that vegetated landscapes transmit information on soil moisture deficits through its down-regulation of evaporation, which in turn affects the humidification and growth of the planetary boundary layer. To test this hypothesis, here, we examined how the evaporative fraction, defined as the ratio between actual and potential evaporation, corresponds with the soil moisture stress index, defined as the relative humidity with a power law dependence on vapor pressure deficit times a parameter, β. We tested the parameterized soil moisture stress index with direct eddy covariance measurements of actual evaporation and computations of potential evaporation based on meteorological conditions, averaged on monthly time steps. The analysis was conducted on a dataset obtained from 144 FLUXNET sites. These sites spanned much of the world's climates and biomes and contained over 5900 months of observations. Observations of the evaporative fraction fit the model of the soil moisture stress index best for semi-arid and arid ecosystems, and least for the humid and wet ecosystems. Consequently, a significant relationship between measurements of the evaporative fraction and soil moisture stress index held for about one-half of the population of sites. Under this condition, the median value of the parameter, β, was 1.41. To investigate the mechanism of this empirical soil moisture stress index, we diagnosed it with a coupled evaporation-planetary boundary layer model. The soil moisture stress index is strongly related to surface resistance, as defined by inverting Penman-Monteith equation. Consequently, this index provides an independent estimate of surface resistance based on easy to measure mean monthly weather conditions like relative humidity and temperature. Thereby, this soil moisture stress index has potential to be applied to weather, climate and biogeochemical models and the interpretation of satellite derived evaporation products, like the one provided by the ECOSTRESS mission.
Extreme daily evaporation from reservoir surfaces can lead to significant short-term water losses, affecting water quality, water supply, and reservoir operation strategies. Historical trends in daily reservoir evaporation events have eluded the scientific and operational communities, largely due to a lack of long-term, consistent data record. This study quantifies trends in extreme daily reservoir evaporation events at 165 major reservoirs located in the western U.S. Here, we use the place-based energy balance and aerodynamic Daily Lake Evaporation Model (DLEM) driven by multiple meteorological data products (RTMA, gridMET, Daymet) to estimate daily evaporation rates at these reservoirs from 1981 to 2022. The results—while are based on different meteorological forcing datasets—consistently indicate that the California, Lower Colorado, and Rio Grande hydrologic regions are more prone to higher daily evaporation extremes. Compared to the rest of western U.S, these three regions also experience a more pronounced increasing trend in the annual maximum daily evaporation rate, at about 0.3 mm day -1 decade -1 during 1981-2022. The results show that heat and dryness are the main drivers to the increasing trend of extreme evaporation, while extreme wind speed is the primary contributor to exceptionally high daily evaporation events across all regions. This phenomenon is particularly prominent in the arid Lower Colorado region, but less significant in the cold and humid Pacific Northwest region. We also find that the correlation between extreme wind speed and extreme evaporation degrades as the time scale increases from daily, to monthly and seasonal. Our findings have strong implications for the pattern and distribution of extreme evaporation events at the western U.S. reservoirs, and illustrate how various drivers influence extreme evaporation across different time scales.
Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.
This paper describes thermal performance of a loop heat pipe (LHP) with two evaporators and two condensers in ambient testing. Each evaporator has an outer diameter of 15mm and a length of 76mm, and has an integral compensation chamber (CC). An aluminum mass of 500 grams is attached to each evaporator to simulate the instrument mass. A thermal electric cooler (TEC) is installed on each CC to provide heating as well as cooling for CC temperature control. A flow regulator is installed in the condenser section to prevent vapor from going back to the evaporators in the event that one of condenser is fully utilized. Ammonia was used ad the working fluid. Tests conducted included start-up, power cycle, heat load sharing, sink temperature cycle, operating temperature control with TECs, and capillary limit tests. Experimental data showed that the loop could start with a heat load of less than 1OW even with added thermal masses. The loop operated stably with even and uneven evaporator heat loads, and even and uneven condenser sink temperatures. The operating temperature could be controlled within +/-0.5K of the set point temperature using either or both TECs, and the required TEC control heater power was less than 2W under most test conditions. Heat load sharing between the two evaporators was also successfully demonstrated. The loop had a heat transport capability of 120W to 140W, and could recover from a dry-out when the heat load was reduced. The 500-gram aluminum mass on each evaporator had a negligible effect on the loop operation. Existing LHPs servicing the orbiting spacecraft have a single evaporator with an outer diameter of about 25mm. Important performance characteristics demonstrated by this LHP included: 1) Operation of an LHP with 15mm diameter evaporators; 2) Robustness and reliability of an LHP with multiple evaporators and multiple condensers under various test conditions; 3) Heat load sharing among LHP evaporators; 4) Effectiveness of TECs in controlling the LHP operating temperature; and 5) Effectiveness of the flow regulator in preventing vapor from going back the evaporators.
Secondary organic aerosols (SOA) dominate atmospheric organic aerosols that affect climate, air quality, and health. Recent studies indicate that, contrary to previously held assumptions, at low relative humidity (RH) these particles are semi-solid and evaporate orders of magnitude slower than expected. Elevated relative humidity has the potential to affect significantly formation, properties, and atmospheric evolution of SOA particles. Here we present a study of the effect of RH on the room-temperature evaporation kinetics of SOA particles formed by ozonolysis of α-pinene and limonene. Experiments were carried out on SOA particles generated, evaporated, and aged at 0%, 50% and 90% RH. We find that in all cases evaporation begins with a relatively fast phase, during which 30% to 70% of the particle mass evaporates in 2 hours, followed by a much slower evaporation rate. Evaporation kinetics at 0% and 50% RH are nearly the same, while at 90% RH a slightly larger fraction evaporates. In all cases, aging the particles prior to inducing evaporation reduces the evaporative losses, with aging at elevated RH leading to more significant effect. In all cases, SOA evaporation is nearly size-independent, providing direct evidence that oligomers play a crucial role in determining the evaporation kinetics.
Evaporation experiments were conducted for SiO2 in three different conditions: in equilibrium, in vacuum, and in hydrogen gas. Evaporation rate in vacuum is about two orders of magnitude smaller than that in equilibrium, which is consistent with previous works. The rate in hydrogen gas changes depending on hydrogen pressure. The rate at 10 exp -7 bar of hydrogen pressure is as small as that of free evaporation, but at 10 exp -5 bar of hydrogen pressure it is larger than that in equilibrium. In equilibrium and in vacuum, the evaporation rate is limited by decomposition of SiO2 on the crystal surface, but it is limited by a diffusion process for evaporation in hydrogen gas. Therefore, evaporation rate of minerals in the solar nebula can be shown neither by that in equilibrium nor by that in vacuum. The maximum temperature of the solar nebula at the midplane at 2-3 AU where chondrites are believed to have originated is calculated to be as low as 150 K, 1500 K, or in between them. The temperature is, in any case, not high enough for total evaporation of the interstellar materials. Therefore, evaporation of interstellar materials is one of the most important processes for the origin and fractionation of solid materials. The fundamental process of evaporation of minerals has been intensively studied for these several years. Those experiments were carried out either in equilibrium or in vacuum; however, evaporation in the solar nebula is in hydrogen (and much smaller amount of helium) gas. In order to investigate evaporation rate and compositional (including isotopic) fractionation during evaporation, vaporization experiments for various minerals in various conditions are conducted. At first, SiO2 was adopted for a starting material, because thermochemical data and its nature of congruent vaporization are well known. Experiments were carried out in a vacuum furnace system.