Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Extrapolation”

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 1,189 records · Page 66

Assessment of Antarctic Ice-Sheet Mass Balance Estimates: 1992 - 2009

Published mass balance estimates for the Antarctic Ice Sheet (AIS) lie between approximately +50 to -250 Gt/year for 1992 to 2009, which span a range equivalent to 15% of the annual mass input and 0.8 mm/year Sea Level Equivalent (SLE). Two estimates from radar-altimeter measurements of elevation change by European Remote-sensing Satellites (ERS) (+28 and -31 Gt/year) lie in the upper part, whereas estimates from the Input-minus-Output Method (IOM) and the Gravity Recovery and Climate Experiment (GRACE) lie in the lower part (-40 to -246 Gt/year). We compare the various estimates, discuss the methodology used, and critically assess the results. Although recent reports of large and accelerating rates of mass loss from GRACE=based studies cite agreement with IOM results, our evaluation does not support that conclusion. We find that the extrapolation used in the published IOM estimates for the 15 % of the periphery for which discharge velocities are not observed gives twice the rate of discharge per unit of associated ice-sheet area than the 85% faster-moving parts. Our calculations show that the published extrapolation overestimates the ice discharge by 282 Gt/yr compared to our assumption that the slower moving areas have 70% as much discharge per area as the faster moving parts. Also, published data on the time-series of discharge velocities and accumulation/precipitation do not support mass output increases or input decreases with time, respectively. Our modified IOM estimate, using the 70% discharge assumption and substituting input from a field-data compilation for input from an atmospheric model over 6% of area, gives a loss of only 13 Gt/year (versus 136 Gt/year) for the period around 2000. Two ERS-based estimates, our modified IOM, and a GRACE-based estimate for observations within 1992 to 2005 lie in a narrowed range of +27 to - 40 Gt/year, which is about 3% of the annual mass input and only 0.2 mm/year SLE. Our preferred estimate for 1992-2001 is - 47 Gt/year for West Antarctica, + 16 Gt/year for East Antarctica, and -31 Gt/year overall (+0.1 mm/year SLE), not including part of the Antarctic Peninsula (1.07 % of the AIS area)

Zwally, H. Jay↗

Semiempirical Theories of the Affinities of Negative Atomic Ions

The determination of the electron affinities of negative atomic ions by means of direct experimental investigation is limited. To supplement the meager experimental results, several semiempirical theories have been advanced. One commonly used technique involves extrapolating the electron affinities along the isoelectronic sequences, The most recent of these extrapolations Is studied by extending the method to Include one more member of the isoelectronic sequence, When the results show that this extension does not increase the accuracy of the calculations, several possible explanations for this situation are explored. A different approach to the problem is suggested by the regularities appearing in the electron affinities. Noting that the regular linear pattern that exists for the ionization potentials of the p electrons as a function of Z, repeats itself for different degrees of ionization q, the slopes and intercepts of these curves are extrapolated to the case of the negative Ion. The method is placed on a theoretical basis by calculating the Slater parameters as functions of q and n, the number of equivalent p-electrons. These functions are no more than quadratic in q and n. The electron affinities are calculated by extending the linear relations that exist for the neutral atoms and positive ions to the negative ions. The extrapolated. slopes are apparently correct, but the intercepts must be slightly altered to agree with experiment. For this purpose one or two experimental affinities (depending on the extrapolation method) are used in each of the two short periods. The two extrapolation methods used are: (A) an isoelectronic sequence extrapolation of the linear pattern as such; (B) the same extrapolation of a linearization of this pattern (configuration centers) combined with an extrapolation of the other terms of the ground configurations. The latter method Is preferable, since it requires only experimental point for each period. The results agree within experimental error with all data, except with the most recent value of C, which lies 10% lower.

Edie, John W.↗

Long-Period Tidal Variations in the Length of Day

A new model of long-period tidal variations in length of day is developed. The model comprises 80 spectral lines with periods between 18.6 years and 4.7 days, and it consistently includes effects of mantle anelasticity and dynamic ocean tides for all lines. The anelastic properties followWahr and Bergen; experimental confirmation for their results now exists at the fortnightly period, but there remains uncertainty when extrapolating to the longest periods. The ocean modeling builds on recent work with the fortnightly constituent, which suggests that oceanic tidal angular momentum can be reliably predicted at these periods without data assimilation. This is a critical property when modeling most long-period tides, for which little observational data exist. Dynamic ocean effects are quite pronounced at shortest periods as out-of-phase rotation components become nearly as large as in-phase components. The model is tested against a 20 year time series of space geodetic measurements of length of day. The current international standard model is shown to leave significant residual tidal energy, and the new model is found to mostly eliminate that energy, with especially large variance reduction for constituents Sa, Ssa, Mf, and Mt.

LOD tidal oscillations↗

Experimental Aging and Lifetime Prediction in Grid Applications for Large-Format Commercial Li-Ion Batteries

Due to the growth of electric vehicle and stationary energy storage markets, the production and use of lithium-ion batteries has grown exponentially in recent years. For many of these applications, large-format lithium-ion batteries are being utilized, as large cells have less inactive material relative to their energy capacity and require fewer electrical connections to assemble into packs. And especially for stationary energy storage systems, where energy delivered is the only revenue source, the economics of these battery systems is highly dependent on cell lifetime. However, testing of large-format lithium-ion batteries is time consuming and requires high current channels and large testing chambers, making information on the performance of commercial, large-format lithium-ion batteries hard to come by. Here, accelerated aging test data from four commercial large-format lithium-ion batteries is reported. These batteries span both NMC-Gr and LFP-Gr cell chemistries, pouch and prismatic formats, and a range of cell designs with varying power capabilities. Accelerated aging test results are analyzed to examine both cell performance, in terms of efficiency and thermal response under load, as well as cell lifetime. Cell thermal response is characterized by measuring temperature during cycle aging, which is used to calculated a normalized thermal resistance value that may help estimate both cell cooling needs or to help extrapolate aging test results to different thermal environments. Cell lifetime is evaluated qualitatively, considering simply the average calendar and cycle life across a range of conditions, as well as quantitatively, using statistical modeling and machine-learning methods to identify predictive aging models from the accelerated aging data. These predictive aging models are then used to investigate cell sensitivities to stressors, such as cycling temperature, voltage window, and C-rate, as well as to predict cell lifetime in various stationary storage applications. Results from this work show that cell lifetime and sensitivity to aging conditions varies substantially across commercial cells, necessitating testing for specific cell formats to make quantitative lifetime predictions. That being said, all commercial cells tested here are predicted to reach at least 10-year lifetimes for stationary storage applications. Based on the aging test results and modeling, some cells are expected to be relatively insensitive to temperature and use-case, making them suited for simple use cases with little or no thermal management and simple controls, while the lifetime of other cells could be extended to 20+ years if operated with thermal management and degradation-aware controls.

battery↗

Is Radioactive Tank Waste a Bingham Fluid?

Slurries and sludges across the United States Department of Energy (DOE) complex rank among the most rheologically interesting. Their composition is heterogeneous, spanning a very broad range of particle sizes, densities, and interparticle forces. All exhibit shear thinning, some have yield stresses, and many are thixotropic. Despite the variety, these complex fluids are often represented using the historic Bingham fluid model, which fits higher shear rate data to a simple straight line. The intercept provides a yield stress, which has been a key design parameter in construction of large-scale waste processing facilities. However, many radioactive wastes are simply not Bingham fluids, which extrapolate but poorly across low to intermediate shear rates that are characteristic of typical processing conditions. Indeed, processing shear rates as high as 200 1/s, which has been a typical minimum shear rate used in fitting, are seldom encountered in nuclear waste processing. Therefore, more realistic rheological models are necessary to accurately predict waste processing performance. PNNL recently re-evaluated the rheology of reconstituted Hanford REDOX sludge waste against a wide variety of rheological models including the Bingham, Cross, Cross with yield stress, Carreau, biviscous, Herschel-Bulkley (which includes a power law dependence), Casson, and Gay models. They found that all of the models provided a closer fit than the Bingham model and that the biviscous model and Cross with yield stress model were convincing. However, reconstituted Hanford REDOX sludge waste is but one type of DOE waste and a direct contrast and comparison of these three models against undiluted, unmixed tank waste (actual not simulant) has not been performed previously. Therefore, the purpose of this paper is to evaluate the rheology of actual tank waste with these more accurate rheological models. In this paper, we evaluate select rheological data for slurry samples from Hanford’s AZ-101, AZ-102, and SY-101 waste tanks. In each of these cases, we find that Cross’ model with yield stress and the biviscous model significantly outperform the Bingham fluid model. Furthermore, the AZ-101 data also shows that the shear stress peak that startup significantly exceeds the Bingham yield stress, which is commonly observed in the incipient moments of rheological measurements on simulants. In contrast, Cross’ model may empirically accommodate an initial spike in shear stress at modest shear rates. These are important observations because computational and analytical fluid dynamics simulations rely on rheological constitutive models for accurately and conservatively predict waste processing performance. These findings suggest the need for better rheological modeling of and validation against radioactive waste.

Rheology, Radioactive Waste↗

Combining data from lidar and in situ instruments to characterize the vertical structure of aerosol optical properties

Over the last decade, the quantification of tropospheric aerosol abundance, composition and radiative impacts has become an important research endeavor. For the most part, the interest in tropospheric aerosols is derived from questions related to the global and local (instantaneous) radiative forcing of climate due to these aerosols. One approach is to study local forcing under well-defined conditions, and to extrapolate such results to global scales. To estimate local aerosol forcing, appropriate radiative transfer models can be employed (e.g., the Fu-Liou radiative transfer code, [Fu and Liou, 1993]). In general, such models require information on derived aerosol properties [Toon, 1994]; namely the aerosol optical depth, single-scattering albedo, and asymmetry factor (phase function), all of which appear in the equations of radiative transfer. In this paper, we report on a method that utilizes lidar data and in situ aerosol size distribution measurements to deduce the vertical structure of the aerosol complex index of refraction in the near IR, thus identifying the aerosol type. Together with aerosol size distributions obtained in situ, the aerosol refractive index can be used to calculate the necessary derived aerosol properties. The data analyzed here were collected during NASA's PEM West-B (Pacific Exploratory Mission) experiment, which took place in February/March 1994. The platform for the measurements was the NASA DC-8 aircraft. The primary goal of the PEM West missions [Browell et al., 1996] was the assessment of potential anthropogenic perturbations of the chemistry in the Pacific Basin troposphere. For this purpose the timing of PEM West-B corresponded to the seasonal peak in transport from the Asian continent into the Pacific basin [Merrill et al., in press]. This period normally occurs during Northern Hemisphere spring, when the Japan jet is well developed.

Redemann, J.↗

Empirical estimation of densities in NaCl-KCl-UCl 3 and NaCl-KCl-YCl 3 molten salts using Redlich-Kister expansion

Densities of molten KCl-NaCl-UCl 3 and KCl-NaCl-YCl 3 ternary systems have been estimated using a multidimensional Redlich-Kister model. Temperature and composition dependent Redlich-Kister functions have been used to generate binary interaction parameters in the outlined ternary salt systems. These binary interactions have been used in the extrapolation to ternary system densities. The results of the density extrapolations by Muggianu interpolation scheme provide agreement within 2–3% for the NaCl-KCl-YCl 3 liquids and 11% in NaCl-KCl-UCl 3 liquids compared to the available experimental data. Modeling NaCl-KCl-UCl 3 molten phase density with a ternary interaction parameter improved the agreement within 4%. Thermophysical modeling used in this study has shown promising results for use in other material properties, such as viscosity, thermal conductivity, and heat capacity of the molten salts. Lastly, the outlined modeling method applied in these specific molten salt ternaries can be used for quaternary or higher multicomponent molten salt systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A critical review of the history of fabricating monolithic U-Mo fuel plates

The fabrication of monolithic uranium-molybdenum alloy fuels, specifically those developed for high-performance research and test reactors, began in the early 2000s. The primary fuel form consists of uranium alloyed with 10-wt% molybdenum in a thin foil coated with zirconium and encapsulated in aluminum-6061 cladding. Over the years, the process has evolved with different types of casting, heat treatments, rolling schedules, and cladding applications. This review examines the history of these fabrication processes and their impact on microstructure and fuel-swelling performance. Even though various fabrication methods were used, we found little correlation between fabrication variation and fuel swelling. This insensitive relationship between fabrication variation and fuel swelling is primarily due to inhomogeneous microstructures that formed during casting and grain refinement that occurred during rolling. We conclude that the fabrication processes we examined produced similar microstructures, indicating that the fuel microstructure is somewhat insensitive to the fabrication parameters evaluated. However, the relatively small amount of historic data, such as those for grain sizes, limited this analysis. More recently fabricated materials, such as those from ongoing irradiation experiment, Mini-Plate-1 and Mini-Plate-2, were also excluded from this analysis and are intended to be reviewed separately. The findings, that fuel microstructure is somewhat insensitive to the fabrication parameters, do not imply that any fabrication method is acceptable, given the uncertainties in data and fuel-swelling observations. For example, only arc melting and vacuum induction melting casting processes were previously explored in the historic fabrication efforts. In conclusion, the findings should not be extrapolated to other casting processes.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Association of Cl with C 2 H 2 by unified variable-reaction-coordinate and reaction-path variational transition-state theory

Significance A key issue in chemical kinetics is advancing the theoretical framework to handle reactions beyond the domain of textbook transition-state theory by including––for example––anharmonicity, barrierless transition states, transition states in series, and the effect of conformational flexibility on equilibrium constants. Other key issues are validating affordable electronic structure methods for direct dynamics and the direct calculation of high-pressure limiting rate constants, which are often obtainable experimentally only by extrapolation. The present article addresses all these issues for a prototype radical–molecule association reaction and thereby demonstrates how to combine improved theoretical methods to provide rate constants in cases where experimental data are uncertain or missing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparison of Spectroscopic Techniques for Determining the Peroxide Value of 19 Classes of Naturally Aged, Plant-Based Edible Oils

The peroxide value of edible oils is a measure of the degree of oxidation, which directly relates to the freshness of the oil sample. Several studies previously reported in the literature have paired various spectroscopic techniques with multivariate analyses to rapidly determine peroxide values using field portable and process instrumentation; those efforts presented “best-case scenarios” with oils from narrowly defined training and test sets. The purpose of this paper is to evaluate the use of near- and mid-infrared absorption and Raman scattering spectroscopies on oil samples from different oil classes, including seasonal and vendor variations, to determine which measurement technique or combination thereof is best for predicting peroxide values. Following peroxide value assays of each oil class using an established titration-based method, global and global-subset calibration models were constructed from spectroscopic data collected on the 19 oil classes used in this study. Spectra from each optical technique were used to create partial least squares regression calibration models to predict the peroxide value of unknown oil samples. A global peroxide value model based on near-infrared (8 mm optical path length) oil measurements produced the lowest RMSEP (4.9), followed by 24 mm optical path length near infrared (5.1), Raman (6.9) and 50 µm optical path length mid-infrared (7.3). However, it was determined that the Raman RMSEP resulted from chance correlations. Global peroxide value models based on low-level fusion of the NIR (8 and 24 mm optical path length) data and all infrared data produced the same RMSEP of 5.1. Global subset models, based on any of the spectroscopies and olive oil training sets from any class (pure, extra light, extra virgin), all failed to extrapolate to the non-olive oils. However, the near-infrared global subset model built on extra virgin olive oil could extrapolate to test samples from other olive oil classes. This work demonstrates the difficulty of developing a truly global method for determining peroxide value of oils.

36 MATERIALS SCIENCE↗

Theory of interstellar medium diagnostics

The theoretical interpretation of observed interplanetary resonance luminescence patterns is used as one of the must promising methods to determine the state of the local interstellar medium (LISM). However, these methods lead to discrepant results that would be hard to understand in the framework of any physical LISM scenario. Assuming that the observational data are reliable, two possibilities which could help to resolve these discrepancies are discussed: (1) the current modeling of resonance luminescence patterns is unsatisfactory and has to be improved, and (2) the extrapolated interstellar parameters are not indicative of the unperturbed LISM state, but rather designate an intermediate state attained in the outer regions of the solar system. It is shown that a quantitative treatment of the neutral gas-plasma interaction effects in the interface between the heliospheric and the interstellar plasmas is of major importance for the correct understanding of the whole complex.

Fahr, H. J.↗

Measuring Neutrino Oscillations in the NOvA and CHIPS Detectors.

NOvA is a long baseline neutrino experiment which uses muon neutrinos produced by the NuMI beam at the Fermi National Accelerator Laboratory to measure muon neutrino disappearance and electron neutrino appearance, among other analyses. The Far/Near Extrapolation analysis framework (FNEX) was used to measure muon neutrino disappearance using the events observed by NOvA while the NuMI beam was in Forward Horn Current configuration (i.e. this analysis uses neutrino data, not antineutrino data). The values of the parameters measured in this analysis are sin 2 θ 23 - 0.499$^{+0.098}_{-0.069}$ and Δm$^2_{32}$ - 2.454$^{+0.108}_{0.119}$ × 10 -3 eV, which is consistent with the results of NOvA’s main muon neutrino disappearance analysis using the experiment’s usual framework, known as CAF. The CHIPS detector is a water Cherenkov detector currently in construction. It will also use the NuMI beam to perform muon neutrino disappearance and electron neutrino appearance analyses. The necessary water clarity for a detector with a diameter of 25 m and a height of 6 m is discussed, as well as the results of testing water filtration systems. The filtration that will be used for the detector is a 10 µm carbon block filter and a 0.5 µm filter in series, with no reverse osmosis or ultraviolet sterilisation being used. This provides an attenuation length of 16.63 ± 0.25 m for light of 532 nm wavelength (green), so light with a wavelength of 405 nm (violet) would have an attenuation length of 133±2 m. The greatest possible distance between a PMT and any point in the detector’s volume is approximately 26 m, so 133 m is sufficient for the detector to measure Cherenkov light produced anywhere in the detector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region

The project Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Rocky Mountain Carbon Capture and Storage (RMCCS) project investigated multiple geologic formations and characterized a local site on the Colorado Plateau for future CCS opportunities. The RMCCS project focused on the Cretaceous Dakota, Jurassic Entrada, and Pennsylvanian Weber Sandstones, the three largest regional formations. All formations in this project are potential CO2 storage resources for future power plants, natural gas processing plants, cement plants, and oil shale development projects. The area adjacent to Craig, Colorado, (Sand Wash Basin) was the area selected for detailed geologic characterization on the RMCCS project. The basin was selected in part because the geology can be extrapolated to other sites on the Colorado Plateau. Field mapping and seismic surveys were conducted to identify and evaluate the basin's structural configuration. A 9,745-foot deep characterization well was drilled to collect 131 feet of core and a suite of geophysical well log data. Petrophysical tests on samples of core were used to calibrate geophysical log data, which can be used to obtain storage resource estimates and evaluate associated uncertainty as well as simulate the hydrologic behavior of injected CO2. A detailed analysis of the primary formations (Dakota, Entrada and Weber sandstones) yielded a more accurate CO2 storage resource assessment for these formations within the Colorado Plateau; RMCCS estimates indicate a total CO2 storage resource of more than 38,000 million metric tons. The characterization of the Sand Wash Basin (2-D seismic surveys, multiple well logs and lithological, petrophysical and geochemical analyses) allowed for a detailed 3-D model to be constructed. The model served as the framework for analyses ranging from CO2 storage resource, injectivity, and subsurface flow to uncertainty estimates to evaluation of risk.

2-D seismic↗

A method of limit point calculation in finite element structural analysis

An approach is presented for the calculation of limit points for structures described by discrete coordinates, and whose governing equations derive from finite element concepts. The nonlinear load-displacement path of the imperfect structure is first traced by use of a direct iteration scheme and the determinant of the governing algebraic equations is calculated at each solution point. The limit point is then established by extrapolation and imposition of the condition of zero slope of the plot of load vs. determinant. Three problems are solved in illustration of the approach and in comparison with alternative procedures and test data.

Gallagher, R. H.↗

Stellar occultation measurements of atmospheric ozone and chlorine from OAO3

The Princeton University telescope on the Orbiting Astronomical Observatory 3 was used for occultation measurements of upper-atmospheric species. The target star, Beta Cen, and observing date were selected so that the line of sight grazed earth's limb at local midnight. Observations at 2580, 2825, 2997, 3100, and 3428 A produced ozone absorption data of high statistical accuracy between 48-km and 114-km tangent altitude. Near-exponential density profiles are obtained up to 85 km. Near 100 km, a significant density excess above the extrapolated low-altitude profile is observed. Near-equatorial observations at altitudes up to approximately 75 km agree well with theoretical calculations. Ozone densities measured at 6.5 to 12.5 deg S latitude above approximately 75 km exceed the theoretical calculations and are not understood at the present time. A search for the 1188-A chlorine line in the 1187-A O2 absorption feature was carried out at tangent altitudes of 96-116 km. An upper limit of 3 ppb to the atomic chlorine mixing ratio at 106 km is found. This upper limit exceeds the expected value for atomic chlorine by a factor of approximately 2.5.

Riegler, G. R.↗

The 10 January 1986 disconnection event in Comet Halley

The disconnection event (DE) in the plasma tail of Comet Halley on January 9-12, 1986 is examined. The distances between the comet head and the disconnected tail are measured for a series of images for that time period and then extrapolated to the nucleus to determine the disconnection time: January 9.60 +/- 0.2 days. The approximate solar-wind conditions at the time of the DE were obtained by corotation of IMP-8 satellite data in earth orbit to Comet Halley. At the time of the DE, Comet Halley is inferred to have been close to a magnetic-sector boundary and a high-speed stream-compression region. However, a heliographic latitude separation of 22 degrees (between the comet and IMP-8) and gaps in the IMP-8 data render a more definitive statement about the linkage of the DE to external conditions quite difficult. It is not possible to resolve the effects of magnetic changes associated with the sector boundary and plasma pressure in the compression region.

Niedner, Malcolm B., Jr.↗

X-38 Experimental Aerothermodynamics

The X-38 program seeks to demonstrate an autonomously returned orbital test flight vehicle to support the development of an operational Crew Return Vehicle for the International Space Station. The test flight, anticipated in 2002, is intended to demonstrate the entire mission profile of returning Space Station crew members safely back to earth in the event of medical or mechanical emergency. Integral to the formulation of the X-38 flight data book and the design of the thermal protection system, the aerothermodynamic environment is being defined through a synergistic combination of ground based testing and computational fluid dynamics. This report provides an overview of the hypersonic aerothermodynamic wind tunnel program conducted at the NASA Langley Research Center in support of the X-38 development. Global and discrete surface heat transfer force and moment, surface streamline patterns, and shock shapes were measured on scaled models of the proposed X-38 configuration in different test gases at Mach 6, 10 and 20. The test parametrics include angle of attack from 0 to 50 degs, unit Reynolds numbers from 0.3 x 10 (exp 6) to 16 x 10 (exp 6)/ ft, rudder deflections of 0, 2, and 5 deg. and body flap deflections from 0 to 30 deg. Results from hypersonic aerodynamic screening studies that were conducted as the configuration evolved to the present shape at, presented. Heavy gas simulation tests have indicated that the primary real gas effects on X-38 aerodynamics at trim conditions are expected to favorably influence flap effectiveness. Comparisons of the experimental heating and force and moment data to prediction and the current aerodynamic data book are highlighted. The effects of discrete roughness elements on boundary layer transition were investigated at Mach 6 and the development of a transition correlation for the X-38 vehicle is described. Extrapolation of ground based heating measurements to flight radiation equilibrium wall temperatures at Mach 6 and 10 were made and generally compared to within 50 deg F of flight prediction.

Horvath, Thomas J.↗

Evaporation Loss of Light Elements as a Function of Cooling Rate: Logarithmic Law

Knowledge about the evaporation loss of light elements is important to our understanding of chondrule formation processes. The evaporative loss of light elements (such as B and Li) as a function of cooling rate is of special interest because recent investigations of the distribution of Li, Be and B in meteoritic chondrules have revealed that Li varies by 25 times, and B and Be varies by about 10 times. Therefore, if we can extrapolate and interpolate with confidence the evaporation loss of B and Li (and other light elements such as K, Na) at a wide range of cooling rates of interest based upon limited experimental data, we would be able to assess the full range of scenarios relating to chondrule formation processes. Here, we propose that evaporation loss of light elements as a function of cooling rate should obey the logarithmic law.

Xiong, Yong-Liang↗