SEARCH · Search NASA

Results for “concentration dependent properties”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Influence of aluminum source and Ni/Al ratio in a batch stirred tank reactor on the structure, morphology, and electrochemical performance of Ni-rich NMA cathodes

Here, the structural, morphological, and electrochemical performance of Ni-rich LiNi 0.9 Mn 0.05 Al 0.05 O 2 (955NMA) and LiNi 0.85 Mn 0.05 Al 0.1 O 2 (85,510) cathodes strongly depends on the properties of their hydroxide precursors. Ni-Mn-Al hydroxide precursors were synthesized through controlled co-precipitation in a batch stirred tank reactor, where pH, reaction time, metal-ion feed rate, aluminum source, and aluminum concentration were systematically varied to tailor particle morphology, phase composition, and dopant distribution. Two aluminum sources, aluminum nitrate and sodium aluminate produced two distinct hydroxide precursors NMA(OH) 2 -1 and NMA(OH) 2 -2, which were lithiated to form LiNMA1 (Li 0.992 [Ni 0.905 Mn 0.049 Al 0.046 ]O 2 ) and LiNMA2 (Li 0.990 [Ni 0.850 Mn 0.047 Al 0.103 ]O 2 ). Structural and compositional analyses revealed that aluminum incorporation and phase formation in Ni–Mn–Al hydroxides are governed by local supersaturation and interfacial growth kinetics. Rapid dilute aluminum addition produced aluminum-free β-phase hydroxides, intermediate conditions generated mixed α/β phases, whereas slow concentrated dosing enabled uniform aluminum incorporation and stabilization of the β-phase structure. LiNMA1 delivers a high initial discharge capacity of 223 mAhg −1 but significant capacity fading with 67% retention after 100 cycles, associated with structural instability. In contrast, LiNMA2 delivers a lower initial capacity 172 mAhg −1 yet excellent cycling stability 91% retention, attributed to improved TM–O framework stability and reduced cation disorder.

Capacity↗

DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗

NASA NTRS1994

Stratospheric aerosol acidity, density, and refractive index deduced from SAGE 2 and NMC temperature data

Water vapor concentrations obtained by the Stratospheric Aerosol and Gas Experiment 2 (SAGE 2) and collocated temperatures provided by the National Meteorological Center (NMC) from 1986 to 1990 are used to deduce seasonally and zonally averaged acidity, density, and refractive index of stratospheric aerosols. It is found that the weight percentage of sulfuric acid in the aerosols increases from about 60 just above the tropopause to about 86 at 35 km. The density increases from about 1.55 to 1.85 g/cu cm between the same altitude limits. Some seasonal variations of composition and density are evident at high latitudes. The refractive indices at 1.02, 0.694, and 0.532 micrometers increase, respectively, from about 1.425, 1.430, and 1.435 just above the tropopause to about 1.445, 1.455, and 1.458 at altitudes above 27 km, depending on the season and latitude. The aerosol properties presented can be used in models to study the effectiveness of heterogeneous chemistry, the mass loading of stratospheric aerosols, and the extinction and backscatter of aerosols at different wavelengths. Computed aerosol surface areas, rate coefficients for the heterogeneous reaction ClONO2 + H2O yields HOCl + HNO3 and aerosol mass concentrations before and after the Pinatubo eruption in June 1991 are shown as sample applications.

Yue, G. K.↗

DOE OSTI2022

High Performance Computing to Quantify the Evolution of Microscopic Concentration Gradients During Flash Processing

During the Flash process, the cross section of a plain-carbon or a low-alloy steel is austenitized through rapid heating and transformed on rapid cooling to a predominantly martensite + bainite structure with small amounts of retained austenite. Unlike conventional heat treating, homogeneity is intentionally avoided during Flash processing of steels. The Flash process assembly consists of a pair of rolls that transfer the steel sheets through the heating and cooling stage of the thermal cycle. The initial microstructure of the steel consists of ferrite (body-centered cubic iron) + carbide ((Fe,X)mCn) mixture. The heating rate through the peak temperature is a function of temperature and reaches a peak of about 300-400°C/s and the cooling rate has a maximum value of 3,000-4,000°C/s. The on-heating phase transformations include carbide dissolution, austenite (face-centered cubic iron) nucleation and growth, and diffusion of carbon and other substitutional elements in the steel. The on-cooling phase transformations include formation of martensite (body-centered tetragonal phase containing supersaturated solute) and bainite (ferrite plates with or without fine carbides). In this project, the focus is on Fe-C-Cr steels that are currently Flash processed for armor applications. The modeling effort proposed here will help optimize the Flash thermal cycle for these low alloy steels to achieve the target performance, which is an ongoing effort at SFP Works. A significant feature of Flash processed Fe-C-Cr steels is the presence of scatter in the through-thickness in the sheet. The variability in hardness results from a variability in the bainite + martensite microstructure that is sensitive to the local chemical concentration of C and Cr. Such a chemical inhomogeneity is intentionally obtained in the Flash process. Although such a microstructural gradient is presumably responsible for the exceptional properties of the Flash processed steel, it is very important to quantify the gradients as a function of Flash variabilities in processing parameters and the input microstructure. Understanding the mechanistic pathway that leads to microstructural gradients could be ground-breaking and instrumental for achieving better process control and optimized microstructural state to meet application-specific strength-ductility requirements. Since the final microstructure depends on setting up precise solute concentration gradients through a rapid heating process, and transforming these regions into various phases, it is important to understand how small changes in steel chemistry, input microstructure (carbide size and distribution), and process variables (Flash thermal cycle) will impact the solute concentration gradients.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Ab-initio informed cluster dynamics simulation of self- and Xe diffusivity in uranium mononitride under irradiation

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuels considered for light water reactor and advanced reactor designs, as it presents significant advantages such as high uranium density (better economics) and high thermal conductivity and melting point (increased safety). Self- and fission gas diffusivities need to be better understood, given that they influence key fuel performance phenomena like swelling and fission gas release. Recently, radiation enhanced diffusivity was investigated in UN by means of cluster dynamics simulations relying on empirical potential-based parameterizations, the reliability of which highly depends on the interatomic potential accuracy. Here, in this work, we refine this approach by determining, using ab-initio calculations, the properties of defect clusters containing vacancies, self-interstitials and Xe impurities. We also consider larger clusters than previous studies. The obtained dataset (formation enthalpies, entropies, and migration barriers) is used to parameterize a cluster dynamics model of mobile clusters, and to calculate the defect cluster concentrations under irradiation. This gives us access to the radiation enhanced self- and fission gas diffusivities. Although the resulting diffusivities are close to the values reported in the literature, we find important qualitative differences in the diffusion mechanisms. Capturing the correct mechanisms is crucial to properly describe the chemistry and fission rate dependence of the model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Quasilinear dispersion in electronic band structure and high Seebeck coefficient in CuFeS 2 -based thermoelectric materials

The earth-abundant natural mineral chalcopyrite CuFeS 2 is a potential n-type thermoelectric material because of its large Seebeck coefficient at high carrier concentrations. For a long time, the large Seebeck coefficient of CuFeS 2 has been attributed to a large electron effective mass, but the reasons for this and the unusual carrier concentration dependent behavior have rarely been discussed. Here, we systematically investigated the special transport behavior of CuFeS 2 and found the classical parabolic band model to be inadequate in explaining it. Our experimental and theoretical studies indicate that there are two flat electronic pockets at the and Z points of the Brillouin zone near the conduction band edge of CuFeS 2 that dominate the charge transport. These electronic pockets result from nonparabolic quasilinearly dispersing bands that give rise to a linear wave vector dependent energy (E ~ k) and a carrier density dependent effective mass (m* ~ m 0 n 1/3 ). Such a strong carrier concentration dependent carrier effective mass results in the high Seebeck coefficient of CuFeS 2 compound under a large carrier density. Here, the work demonstrates that quasilinearly dispersing bands can give strongly enhanced Seebeck coefficient, and could be useful in optimizing the properties of thermoelectric materials.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Evaluating Transfer and Pumping of Slurries from Pulsed Jet Mixed Vessels

Due to gravity, solids in slurries will settle if density differences between the solids and liquid are positive (i.e., particle has a negative buoyant force) unless rheological properties and flow conditions are adequate to overcome the gravitational effects. The rate of settling depends on the force balance of the particle, which includes the surface forces associated with fluid rheology. Given the same fluid and solid properties, particles that are larger and denser tend to settle faster. When pumping slurry into a vessel at concentrations precluding hindered settling with insufficient mixing, particle and density distributions can result in preferential settling, creating stratification in the solids concentration within the vessel. For vessels with transfer line inlets located in the lower portion of the tank, the stratified solids concentration may be detrimental to the transfer system performance. Elevated concentrations of solids in the slurry entrained at the inlet to the transfer line can result in the effective viscosity or slurry bulk density exceeding the design limits of the pump. These conditions could result in plugging of the transfer line or onset of cavitation of the pumps because of excessive pressure drop. These conditions can be exacerbated with periodic inlet conditions existing at the transfer line inlet. Periodic conditions can result when vessel mixing is intermittent such as with pulsed jet mixers (PJM). The transfer line inlet conditions are impacted by the periodic nature of the PJM operations with respect to suspension of solids and their transport to the inlet of the transfer line. A scaling approach is presented, and corresponding test requirements are developed for assessing the prevention of plugging the pipeline. Line plugging mechanisms are addressed that do not include plugging due to steady-state high-density slurry entering the transfer line and reducing the net positive suction head available (NPSHA) at the pump inlet to below that required for pump operation. Items considered include the transition to reduced relative flow velocities, such that the critical pipe velocity for solids deposition, Ucd, is not maintained, and segregation of heavy solids during the transport. The recommended requirements to prevent plugging include: • Limits for viscosity and density for entrained slurry to prevent the pressure drop in the pipeline from exceeding pump capacity. • Limits for viscosity and density for entrained slurry to prevent the net positive suction head available (NPSHA) from falling below the net positive suction head required (NPSHR) for operating the pump. • Transfer line velocity and flow rate requirements to maintain solids in suspension, while avoiding line plugging that results from deposition of solids within the transfer line. This paper describes the development of the scaling and testing requirements to verify that proposed approaches for transfer and pump out are appropriately developed

periodic flow, mixing, pipeline plugging, scaling,↗

NASA NTRS1972

A molecular theory of elastomer deformation and rupture.

The mechanical properties of elastomers, including rupture and its time dependence, can be semiquantitatively predicted from nine molecular parameters which are characteristic for a given species, from the initial molecular weight of a sample before cross-linking and from the effective chain concentration (which must still be determined for each vulcanizate). Only two empirical quantities are involved - the Plazek retardation function for entanglement slippage and a related constant which serves to locate this function on the time scale. The theory leads to a new method of estimating fatigue lifetimes from short-time data.

Landel, R. F.↗

NASA NTRS1973

Supersonic turbulent boundary-layer flows with tangential slot injection.

An implicit finite-difference method has been developed for the solution of the compressible boundary-layer equations. This method is applied to tangential slot injection into supersonic turbulent boundary-layer flows. In addition, the effects induced by the interaction between the boundary-layer displacement thickness and the external pressure field are considered. Three different eddy viscosity models have been used to specify the turbulent momentum exchange. One model depends on the species concentration profile, and the species conservation equation has been included in the system of governing partial differential equations. For air injected into air, the freestream and injected gases are treated as separate species which have common fluid properties. Calculations were made and results were compared with experimental data at freestream Mach numbers of 2.4 and 6.0 and with results of another finite-difference method. Good agreement was obtained for the reduction of wall skin friction with slot injection.

Miner, E. W.↗

NASA NTRS1977

The moon's permanent magnetic field - A cratered-shell model

The possibility is considered that the properties of the larger-scale lunar magnetic field may depend in a rather straightforward manner on the geometrical properties of the distribution of larger craters. The relevant effects of cratering are analyzed in the framework of a model of a spherically symmetric shell magnetized to some specific intensity by a concentric dipolar magnetic field. The remanent dipole moment of the model is determined by calculating the dipole moment of the material removed to form craters and them computing the remanent dipole moment for three mutually orthogonal orientations of the magnetizing dipole. The resulting absolute values for the remanent dipole moment are found to be comparable to the upper limit of 10 to the -19th power G-cu cm estimated from Apollo subsatellite measurements. It is concluded that for a crust with an empirically reasonable mean magnetization, the craters can account for a dipole moment as great as that indicated by the estimated upper limit.

Weiss, H.↗

NASA NTRS1995

Single-drop reactive extraction/extractive reaction with forced convective diffusion and interphase mass transfer

An algorithm has been developed for the forced convective diffusion-reaction problem for convection inside and outside a droplet by a recirculating flow field hydrodynamically coupled at the droplet interface with an external flow field that at infinity becomes a uniform streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet or reactions can take place in both phases. The algorithm has been implemented and results are shown here for the case of no reaction and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Kleinman, Leonid S.↗

NASA NTRS1999

Exploiting the Temperature/Concentration Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

The objective of this new research project is to demonstrate by experiment, supplemented by mathematical modeling and physical property measurement, that the effects of buoyancy driven convection can be largely eliminated in ground-based experiments, and further reduced in flight, by applying a new technique. That technique exploits the dependence of magnetic susceptibility on composition or temperature. It is emphasized at the outset that the phenomenon to be exploited is fundamentally and practically different from the magnetic damping of convection in conducting liquids that has been the subject of much prior research. The concept suggesting this research is that all materials, even non-conductors, when placed in a magnetic field gradient, experience a force. Of particular interest here are paramagnetic and diamagnetic materials, classes which embrace the "model alloys", such as succinonitrile-acetone, that have been used by others investigating the fundamentals of solidification. Such alloys will exhibit a dependence of susceptibility on composition. The consequence is that, with a properly oriented field (gradient) a force will arise that can be made to be equal to, but opposite, the buoyancy force arising from concentration (or temperature) gradients. In this way convection can be stilled. The role of convection in determining the microstructure, and thereby properties, of materials is well known. Elimination of that convection has both scientific and technological consequences. Our knowledge of diffusive phenomena in solidification, phenomena normally hidden by the dominance of convection, is enhanced if we can study solidification of quiescent liquids. Furthermore, the microstructure, microchemistry and properties of materials (thereby practical value) are affected by the convection occurring during their solidification. Hitherto the method of choice for elimination of convection has been experimentation in microgravity. However, even in low Earth orbit, residual convection has effects. That residual convection arises from acceleration (drag on the spacecraft), displacement from the center of mass or transients in the gravitational field (g-jitter). There is therefore a need for both further reducing buoyancy driven flow in flight and allowing the simulation of microgravity during ground based experiments. Previous investigations, the research project description, theory behind the study and experimental methods as well as plots of magnetic fields and forces are presented.

Evans, J. W.↗

DOE OSTI2024

Soft matter approach for creating novel protein hydrogels using fractal whey protein assemblies as building blocks

In this study, mesoscopic sized fractal assembly (FA) particles were prepared using whey protein isolate; the coldset gelation properties of FA particles were investigated in-depth. Two types of FA particles (FA -62 and FA -90) with different mean sizes were synthesized through controlled thermal treatment of whey protein solutions at two concentrations (62 g/L and 90 g/L). Particle characteristics e.g. hydrodynamic radius, ζ-potential and surface hydrophobicity were dependent on pH and structure of FA. Transmission electron microscopy (TEM) observation confirmed the fractal morphology and small-angle X-ray scattering (SAXS) analysis suggested an internal fractal structure for the obtained FA particles. Eight cold-set FA protein gels (2 % w/v) were manufactured by controlling two gelling factors at two levels: pH (5.8 and 7.0) and Ca 2+ content (5 mM and 10 mM). Rheological characteristics in the large amplitude oscillatory shear regime revealed that pH 7.0 gels were softer and elastic while pH 5.8 gels were harder and brittle. Rheology Pipkin diagrams demonstrated that the strain softening/stiffening and the shear-thinning/thickening behaviors may be fine-tuned by manipulating the key gelation factors: e.g. FA structure, pH, and Ca 2+ . The entrainment speed-dependent friction coefficient curves showed that at an intermediate velocity regime (6-250 mm s -1 ), FA -90 particles induced hydrogels had superior lubrication effect compared to FA -62 gels. Further, this work demonstrated a food structure design approach regarding tuning texture and lubrication properties of protein gels without changing protein content and protein composition. The optimized protein hydrogels may be used as texturizer for "cleaner label" food formulas and/or as delivery system for carrying micronutrients.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Arctic Cloud‐Base Ice Precipitation Properties Retrieved Using Bayesian Inference

Abstract Cloud‐climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long‐lived, and nearly continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This Arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at Arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud‐base properties a more appealing target for inference and evaluation of model simulations. Here I describe an inverse model for the estimation of cloud base ice precipitation properties over Utqiagvik, North Slope of Alaska, using the synthesis of ground‐based radar and lidar measurements. By leveraging a Markov Chain Monte Carlo algorithm as the core of the inverse model, a wide range of particle size distributions are sampled, and different combinations of ice habit models are examined, both of which are typically fixed in other retrieval methods. Results show intriguing links between different cloud base thermodynamic and ice precipitation properties. Apparent ice number concentration enhancements at temperatures of −5 and −15°C suggest possible secondary ice production (SIP). The analysis alludes to an overestimation of SIP occurrence and intensity, especially in studies relying only on radar or lidar measurements. Finally, reflectivity‐dependent ice precipitation rate and ice water content parameterizations are presented.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Characterization of Arsenic and Selenium in Coal Fly Ash to Improve Evaluations for Disposal and Reuse Potential (Final Technical Report)

Coal fly ash is a high volume waste material that is discarded in landfills and surface water impoundments across the U.S. and is also widely recycled for a variety of applications. The leaching of potential of contaminants of concern, such as arsenic (As) and selenium (Se), is often the driver of risk assessments for coal ash disposal and reuse. The extent of leachable As and Se depends on several factors related to environmental conditions and fly ash characteristics. Previous studies employed various methods to delineate the concentration, chemical form, and distribution of As and Se in fly ash materials. However, few studies have attempted to directly correlate these properties to mobilization parameters relevant to disposal and reuse. Instead, the coal residuals industries often rely upon standardized leaching protocols that can be laborious or involve hazardous chemicals. The goals of the project were to: 1) Develop and evaluate a characterization protocol that can be used to screen fly ash samples for leachability of As and Se; 2) Characterize As, Se, and associated constituents of fly ash particles at multiple length scales (nanometer to micrometer) to determine if elemental associations differ as a function of the resolution of characterization; and 3) Establish a predictive model for the chemical composition of coal ash produced annually at major U.S. coal fired power facilities on 50-year national coal supply records. For the first objective, we performed leaching experiments with 52 fly ash samples collected from 15 different U.S. power plants and representing coal feedstocks from the three major domestic coal regions. For this work, we assessed the mobilization potential of As and Se in fly ash based on standardized leaching protocols and performed multivariate and lasso regression analyses to explore correlations of leachable As and Se contents with characteristics such as major element contents, loss on ignition (LOI) and pH. The results of regression models indicated that major elements (Fe, Ca, Al) for a wide range of fly ashes can serve as predictor variables for the leaching potential of As, but not for Se. LOI and pH were not important predictive variables in the models. Both regression approaches resulted in relatively strong fits for leachable As (correlation coefficient R 2 = 0.78 for both models) compared to models for leachable Se (R 2 = 0.49). Overall, these results suggest that correlation models combined with on-site elemental analysis with portable analyzers may enable a screening method for leachable As in coal ash. For the second objective, we utilized nanoscale 2-D imaging (30-50 nm spot size) with the Hard X-ray Nanoprobe (HXN) in combination with microprobe X-ray capabilities (~5 µm resolution) to determine As and Se elemental associations in fly ash particles. Speciation of As and Se was also measured at the nano- to microscale with X-ray absorption spectroscopy. The enhanced resolution of HXN showed As and Se that were diffusely located around or comingled with Ca- and Fe-rich particles. The results also showed nanoparticles of Se attached to the surface of fly ash grains. Overall, a comparison of As and Se species across scales highlights the heterogeneity and complexity of chemical associations for these trace elements of concern in coal fly ash. For the final objective, we developed a predictive model for major element composition of coal ash in reserve at disposal sites of major U.S. coal fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973-2022 and was trained on coal ash composition data showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe 2 O 3 ), as coal-fired power stations changed their source of coal over this time frame. Our approach enables an estimation of coal ash chemical composition that is stored in waste impoundments at individual power stations. Such information can help delineate the regional market potential for material applications that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗

NASA NTRS1983

The 3D inelastic analysis methods for hot section components

The objective of this research is to develop an analytical tool capable of economically evaluating the cyclic time dependent plasticity which occurs in hot section engine components in areas of strain concentration resulting from the combination of both mechanical and thermal stresses. The techniques developed must be capable of accommodating large excursions in temperatures with the associated variations in material properties including plasticity and creep. The overall objective of this proposed program is to develop advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate and yet more cost effective analysis of combustors, turbine blades, and vanes. The approach will be to develop four different theories, one linear and three higher order with increasing complexities including embedded singularities.

Dame, L. T.↗

NASA NTRS2004

Chemical Reaction and Flow Modeling in Fullerene and Nanotube Production

The development of processes to produce fullerenes and carbon nanotubes has largely been empirical. Fullerenes were first discovered in the soot produced by laser ablation of graphite [1]and then in the soot of electric arc evaporated carbon. Techniques and conditions for producing larger and larger quantities of fullerenes depended mainly on trial and error empirical variations of these processes, with attempts to scale them up by using larger electrodes and targets and higher power. Various concepts of how fullerenes and carbon nanotubes were formed were put forth, but very little was done based on chemical kinetics of the reactions. This was mainly due to the complex mixture of species and complex nature of conditions in the reactors. Temperatures in the reactors varied from several thousand degrees Kelvin down to near room temperature. There are hundreds of species possible, ranging from atomic carbon to large clusters of carbonaceous soot, and metallic catalyst atoms to metal clusters, to complexes of metals and carbon. Most of the chemical kinetics of the reactions and the thermodynamic properties of clusters and complexes have only been approximated. In addition, flow conditions in the reactors are transient or unsteady, and three dimensional, with steep spatial gradients of temperature and species concentrations. All these factors make computational simulations of reactors very complex and challenging. This article addresses the development of the chemical reaction involved in fullerene production and extends this to production of carbon nanotubes by the laser ablation/oven process and by the electric arc evaporation process. In addition, the high-pressure carbon monoxide (HiPco) process is discussed. The article is in several parts. The first one addresses the thermochemical aspects of modeling; and considers the development of chemical rate equations, estimates of reaction rates, and thermodynamic properties where they are available. The second part addresses modeling of the arc process for fullerene and carbon nanotube production using O-D, 1-D and 2-D fluid flow models. The third part addresses simulations of the pulsed laser ablation process using time-dependent techniques in 2-D, and a steady state 2-D simulation of a continuous laser ablation process. The fourth part addresses steady state modeling in O-D and 2-D of the HiPco process. In each of the simulations, there is a variety of simplifications that are made that enable one to concentrate on one aspect or another of the process. There are simplifications that can be made to the chemical reaction models , e.g. reduction in number of species by lumping some of them together in a representative species. Other simulations are carried out by eliminating the chemistry altogether in order to concentrate on the fluid dynamics. When solving problems with a large number of species in more than one spatial dimension, it is almost imperative that the problem be decoupled by solving for the fluid dynamics to find the fluid motion and temperature history of "particles" of fluid moving through a reactor. Then one can solve the chemical rate equations with complex chemistry following the temperature and pressure history. One difficulty is that often mixing with an ambient gas is involved. Therefore, one needs to take dilution and mixing into account. This changes the ratio of carbon species to background gas. Commercially available codes may have no provision for including dilution as part of the input. One must the write special solvers for including dilution in decoupled problems. The article addresses both ful1erene production and single-walled carbon nanotube (SWNT) production. There are at least two schemes or concepts of SWNT growth. This article will only address growth in the gas phase by carbon and catalyst cluster growth and SW T formation by the addition of carbon. There are other models that conceive of SWNT growth as a phase separation process from clusters me up carbon and metal catalyst, with the carbon precipitating from the cluster as it cools. We will not deal with that concept in this article. Further research is needed to determine the rates at which these composite clusters form, evaporate, and segregate.

Scott, Carl D.↗

DOE OSTI2021

First-principles calculation of lattice distortions in four single phase high entropy alloys with experimental validation

Exceptional properties of high entropy alloys (HEAs) are attributed to the disordered random solid solution of multiple alloying elements. Despite numerous studies, the fundamental understanding at the electronic and atomic levels is still missing. We report a comparative study of four fcc HEAs of NiFeCoCr and NiFeCoCrX (X = Mn, Cu, or Pd) based on ab initio calculations using large supercells with 500 atoms in equal composition. After fully optimizing their structures using the VASP package, their electronic structure, interatomic bonding, partial charge distribution, and mechanical properties are calculated and compared, revealing the intricate interdependence among them. A novel parameter based on the quantum mechanical metric for internal cohesion, the total bond order density (TBOD), is used to interpret the calculated properties. The highest TBOD is found in Cantor alloy NiFeCoCrMn but lower in NiFeCoCrPd. The atomic radii vary depending on their local chemical environment. The resulting lattice distortions is validated experimentally in NiFeCoCrMn and NiFeCoCrPd. Moreover, modeling of Cu and Pd clustering in the supercell shows they have lower total energy in agreement with the observation of Cu-enhanced nano-participates in NiFeCoCrCu and Pd-induced concentration wave in NiFeCoCrPd HEAs.

36 MATERIALS SCIENCE↗