Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffuse interface methods”

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.

319 records · Page 18

Immiscible phase incorporation during directional solidification of hypermonotectics

Solidification processes in immiscible samples were investigated by directly observing the events taking place at the solid-liquid interface during directional solidification. Visualization of these events was made possible through the use of a transparent metal analog system and a temperature gradient stage assembly fitted to an optical microscope. The immiscible transparent analog system utilized was the succinonitrile-glycerol system. This system has been shown to exhibit the same morphological transitions as observed in metallic alloys of monotectic composition. Both monotectic and hypermonotectic composition samples were directionally solidified in order to gain an improved understanding of the manner in which the excess hypermonotectic liquid is incorporated into the solidifying structure. The processing conditions utilized prevented sedimentation of the excess hypermonotectic liquid by directionally solidifying the samples in very thin (13 microns), horizontally oriented cells. High thermal gradient to growth rate ratios (G/R) were used in an effort to prevent constitutional supercooling and the subsequent formation of L(sub 2) droplets in advance of the solidification front during the growth of fibrous composite structures. Results demonstrated that hypermonotectic composites could be produced in samples up to two weight percent off of the monotectic composition by using a G/R ratio greater than or equal to 4.6 x 10(exp 4) C(s)/mm(sup 2) to avoid constitutional supercooling. For hypermonotectic samples processed with G/R ratios below 4.6 x 10(exp 4) C(s)/mm(sup 2), constitutional supercooling occurred and resulted in slight interfacial instability. For these samples, two methods of incorporation of the hypermonotectic liquid were observed and are reported. The correlation between the phase spacing, lambda, and the growth rate, R, was examined and was found to obey a relationship generally associated with a diffusion controlled coupled growth process. For samples with compositions ranging from the monotectic composition up to 2 percent off of the monotectic composition, data indicated that the square of the phase spacing (lambda) varied linearly with the inverse of the growth rate (R).

Andrews, J. Barry↗

Cooling rate and dendrite spacing control in direct metal deposition printed Cu-Fe alloys

Over the last several decades, additive manufacturing (AM) has been primarily used for rapid prototyping or to create novel geometries that would be difficult or impossible to create by normal manufacturing methods. More recently, research has been focused on expanding the list of materials that can be made through additive manufacturing, opening a greater range of material properties for this manufacturing method. Due to the unusual conditions during AM, including the high cooling rates and voxel by voxel method of production, AM parts often have anisotropic material microstructures and properties. In this investigation, the laser power, composition, and nozzle head speed during direct metal deposition of copper-iron alloys was varied to understand how the grain structure within the printed parts could be changed and controlled. The resulting dendrite spacing was measured and compared to calculated cooling rates from Gaussian beams on flat plates under similar material and laser properties, which resulted in a cooling rate to dendrite spacing relationship following an inverse square root, as is found in other dendritic systems [Young and Kerkwood, Metall. Trans. A 6, 197–205 (1975)]. Thus, it is demonstrated that in the Cu-Fe system, dendrite spacing can be controlled through manipulation of printing parameters.

36 MATERIALS SCIENCE↗

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗

Rapid Simulation of Unconventional Reservoirs by Multidomain Multiresolution Modeling Based on the Diffusive Time of Flight

Summary The fast marching method (FMM)-based rapid flow simulation has been shown to accelerate simulation efficiency by orders of magnitude by transforming 3D simulation to equivalent 1D simulation using the concept of the “diffusive time of flight” (DTOF). However, the 1D transformation does not directly apply to multiwell problems. In this paper, we propose a novel DTOF-based multidomain multiresolution discretization scheme to accelerate multiwell simulation of unconventional reservoirs. Our method formulates multiwell simulation problems based on the DTOF which displays the pressure front propagation in unconventional reservoirs. The DTOF contours are used to partition the reservoir into local and shared domains. A local domain is where the flow is dominated by a single well, and the shared domain is where the fluid flow is influenced by multiple wells. The DTOF contours expand independently in local domains and interfere in the shared domain. After the partitioning, each domain is discretized using a multiresolution scheme whereby the original 3D fine mesh is preserved near the wells to account for detailed physics including gravity, and the rest of the domain is discretized into 1D mesh based on the DTOF contours to alleviate the simulation workload. The power and efficacy of our approach are demonstrated using synthetic and field-scale simulation models with different degrees of geologic and well-completion complexity. The simulation results, number of active cells, and computation time for the proposed discretization scheme are compared with the original high-fidelity 3D model for each case. The results show that the proposed method is suitable for multiwell simulation problems in unconventional reservoirs and can accelerate flow simulations by orders of magnitude with minimal loss of accuracy. The novelty of this work is the creation of DTOF-derived multiresolution discretization with local and shared domains to simplify and accelerate the calculation of subsurface flow problems, especially in unconventional reservoirs. Our workflow can be easily interfaced with commercial simulators, making it suitable for large-scale field applications.

Engineering↗

Phase-field modeling of diffusion bonding in 316H stainless steel: Impact of processing conditions on grain morphology and bonding quality

A novel multi-phase, multi-component phase‐field model is presented to study the diffusion bonding of 316H stainless steel. Combined with targeted experimental investigations, this model simulates the bond-growth process and predicts the bonding quality. Unlike previous models, our approach captures the simultaneous evolution of voids and grain structures, while quantifying bonding quality using defined bonding ratio. A comprehensive analysis of bond process control is performed by changing temperature, pressure and surface roughness observing the resulting bond structure, which is consistent with experimental observations and analytical predictions. Temperature is determined to be the dominant factor, with the transition from a flat to a robust bond occurring between 1000 °C and 1050 °C. At the ideal bonding temperature of 1050 °C, a surface roughness exceeding 0.6 μm or an applied stress below 4 MPa results in poor bonding quality. Beyond this, higher pressures and smoother surfaces reduce void size, accelerate void shrinkage, and lead to improved bond integrity. This diffuse-interface model can be extended to other material systems if supplied with appropriate thermodynamic and kinetic data. In conclusion, this makes it an effective modeling platform for optimizing high-temperature diffusion bonding and developing reliable bonded components such as compact heat exchangers.

Diffusion bonding↗

Profile Evolution Simulation in Etching Systems Using Level Set Methods

Semiconductor device profiles are determined by the characteristics of both etching and deposition processes. In particular, a highly anisotropic etch is required to achieve vertical sidewalls. However, etching is comprised of both anisotropic and isotropic components, due to ion and neutral fluxes, respectively. In Ar/Cl2 plasmas, for example, neutral chlorine reacts with the Si surfaces to form silicon chlorides. These compounds are then removed by the impinging ion fluxes. Hence the directionality of the ions (and thus the ion angular distribution function, or IAD), as well as the relative fluxes of neutrals and ions determines the amount of undercutting. One method of modeling device profile evolution is to simulate the moving solid-gas interface between the semiconductor and the plasma as a string of nodes. The velocity of each node is calculated and then the nodes are advanced accordingly. Although this technique appears to be relatively straightforward, extensive looping schemes are required at the profile corners. An alternate method is to use level set theory, which involves embedding the location of the interface in a field variable. The normal speed is calculated at each mesh point, and the field variable is updated. The profile comers are more accurately modeled as the need for looping algorithms is eliminated. The model we have developed is a 2-D Level Set Profile Evolution Simulation (LSPES). The LSPES calculates etch rates of a substrate in low pressure plasmas due to the incident ion and neutral fluxes. For a Si substrate in an Ar/C12 gas mixture, for example, the predictions of the LSPES are identical to those from a string evolution model for high neutral fluxes and two different ion angular distributions.(2) In the figure shown, the relative neutral to ion flux in the bulk plasma is 100 to 1. For a moderately isotropic ion angular distribution function as shown in the cases in the left hand column, both the LSPES (top row) and rude's string method (bottom row) predict tapered profiles. The LSPES uses an AD with a FWHM = 13.5 degrees, and rude's model uses a ratio of sheath voltage to ion temperature of 50. The more anisotropic IADs produce profiles with more vertical sidewalls and a wider bottom surface, as shown in the right hand column. Here, the LSPES has an AD with a FWHM = 2 degrees, and rude's model uses a ratio of 500 of sheath voltage to ion temperature. The agreement between the LSPES and rude's model is excellent in both cases. We will present etching profiles generated by the LSPES, including calculations of the re-emitted fluxes of both neutrals and ions off of the profile walls. In addition, we will show the effect of geometric structures (overhangs, etc.) on the etching profiles. Other physical aspects, such as surface diffusion, will also be included in the model.

Hwang, Helen H.↗

Symmetry in Seeded Metal Nanocrystal Growth

Symmetry underpins the structure and function of the world around us and is also captured in modern nanomaterials, where changing the symmetry of a nanocrystal or the interparticle spacing and orientation of nanocrystal building blocks in a superlattice can give new function. However, the synthesis and assembly of nanocrystals have been limited largely to simple compositions and structures. It remains a grand challenge to achieve nanocrystals with compositional and structural complexity while maintaining the monodispersity required for their use. This Account will illustrate through recent examples that seeded methods enable the synthesis of compositionally and structurally complex multimetallic crystals with defined and predictable symmetries for applications in plasmonics and catalysis. This outcome arises because the barrier for heterogeneous nucleation (i.e., seeded) is lower than that of homogeneous nucleation, where seeds can serve as preferential sites for the growth of complex structures and crystal phases. Our analysis begins by considering metal overgrowth from single-crystalline seeds of different shapes and symmetries, where the kinetics of adatom addition to seeds relative to their diffusion across seeds accounts for the expressed nanocrystal shapes. These results are then compared to overgrowth from seeds with different internal structures (i.e., planar defects), where the relationships between nanocrystal size and volumetric strain energy and surface energy are discussed. A major finding from this analysis is that often the underlying symmetry of seeds can be predictably transferred to the final crystals during overgrowth processes. Consequences of this finding are the predictable syntheses of crystals with different hierarchies akin to snowcrystals as well as nanocrystals with complex compositions (e.g., quaternary nanoparticles). Yet, there are subtle aspects to seeded growth that pave a path toward examples where nanocrystal symmetry has been reduced compared to the original seeds in a controlled manner. As we found, both the concentrations of metal precursors and capping agents can impact whether symmetry is transferred or reduced during overgrowth. Examples from our laboratory will be placed in context to other reported strategies for symmetry breaking. As will be argued, understanding what conditions favor symmetry preservation versus symmetry reduction during seeded crystal growth is central to accessing next-generation crystal forms. The Account concludes by outlining synthetic challenges associated with forming nanoscale heterostructures with precise 3-D placement of different materials within a given nanocrystal as well as facet control within different material domain and interface engineering. Furthermore, we envision meeting these challenges through regioselective and chemoselective seeded syntheses for which a foundation is outlined herein.

36 MATERIALS SCIENCE↗

Improving stability of ALD ZrN thin film coatings over U-Mo dispersion fuel

Atomic layer deposition (ALD) of ZrN is a candidate technology for coating U-Mo dispersion fuel as a diffusion barrier coating. During the early development stages of the coating the ZrN deposition showed to be mechanically unstable and ultimately resulted in spalling. Based on experimental outcomes it was found that instabillity can be eliminated through the introduction of a thin amorphous (a) Al 2 O 3 interlayer coating deposited in between the ~1-μm thick ZrN and the U-Mo. substrate. To elucidate these findings simulations were performed with density functional theory (DFT) to measure work of adhesions at different interfaces while finite element modeling (FEM) was performed to measure the residual stress distribution. DFT indicated that the ZrN coating deposited over (a)-Al 2 O 3 is ~2.5 times stronger when compared to direct deposition over U or UO 2 substrate. While calculations from FEM recognized; (1) large stress concentrations can originate from the irregularly distributed native surface oxides (UO 2 ) over U-Mo and (2) stress concentrations can be reduced if those surface UO 2 can be modified to a uniform layer (achieved after application of a 8 nm (a)-Al 2 O 3 interlayer). These results provided explanations and confirmed the role played by interlayer in eliminating the original mechanical instability of the ZrN.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Over-Expression, Purification and Crystallization of Human Dihydrolipoamide Dehydrogenase

Dehydrolipoamide dehydrogenase (E3; dihydrolipoan-tide:NAD+ oxidoreductase, EC 1.8.1.4) is a common catalytic component found in pyruvate dehydrogenase complex, alpha-ketoglutarate dehydrogenase complex, and branched-chain cc-keto acid dehydrogenase complex. E3 is also a component (referred to as L protein) of the glycine cleavage system in bacterial metabolism (2). Active E3 forms a homodimer with four distinctive subdomain structures (FAD binding, NAD+ binding, central and interface domains) with non-covalently but tightly bound FAD in the holoenzyme. Deduced amino acids from cloned full-length human E3 gene showed a total of 509 amino acids with a leader sequence (N-terminal 35 amino acids) that is excised (mature form) during transportation of expressed E3 into mitochondria membrane. So far, three-dimensional structure of human E3 has not been reported. Our effort to achieve the elucidation of the X-ray crystal structure of human E3 will be presented. Recombinant pPROEX-1 expression vector (from GIBCO BRL Life Technologies) having the human E3 gene without leader sequence was constructed by Polymerase Chain Reaction (PCR) and subsequent ligation, and cloned in E.coli XL1-Blue by transformation. Since pPROEX-1 vector has an internal His-tag (six histidine peptide) located at the upstream region of a multicloning site, one-step affinity purification of E3 using nickelnitriloacetic acid (Ni-NTA) agarose resin, which has a strong affinity to His-tag, was feasible. Also a seven-amino-acid spacer peptide and a recombinant tobacco etch virus protease recognition site (seven amino acids peptide) found between His-tag and first amino acid of expressed E3 facilitated the cleavage of His-tag from E3 after the affinity purification. By IPTG induction, ca. 15 mg of human E3 (mature form) was obtained from 1L LB culture with overnight incubation at 25C. Over 98% of purity of E3 from one-step Ni-NTA agarose affinity purification was confirmed by SDS-PAGE analysis. For crystallization, E3 samples were prepared with and without His-tag. To minimize the aggregation of E3, apo- and holo- forms of E3s were tested, as well as a mutated E3. Dynamic light scattering measurements revealed that the E3 preparations without His-tag and substrate are highly monodispersive with regard to homodimers. Consequent crystallization trials of this E3 preparation led to single crystals of E3 grown by the vapor diffusion method. Crystals were obtained within a few days from solution containing poly (ethylene glycol) monomethyl ether 5000 as a precipitant. Autoindexing and integration of the X-ray diffraction data showed that E3 crystals belong to an orthorhombic system with unit cell parameters a-- 123. 1, b= 165.3 and c=214.3A. Further optimization of protein preparation and crystallization experiments for the structural determination will be discussed.

Hong, Y. S.↗

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗

Single-Molecule Detection of the Encounter and Productive Electron Transfer Complexes of a Photosynthetic Reaction Center

Small, diffusible redox proteins play an essential role in electron transfer (ET) in respiration and photosynthesis, sustaining life on Earth by shuttling electrons between membrane-bound complexes via finely tuned and reversible interactions. Ensemble kinetic studies show transient ET complexes form in two distinct stages: an “encounter” complex largely mediated by electrostatic interactions, which subsequently, through subtle reorganization of the binding interface, forms a “productive” ET complex stabilized by additional hydrophobic interactions around the redox-active cofactors. Here, using single-molecule force spectroscopy (SMFS) we dissected the transient ET complexes formed between the photosynthetic reaction center-light harvesting complex 1 (RC-LH1) of Rhodobacter sphaeroides and its native electron donor cytochrome c 2 (cyt c 2 ). Importantly, SMFS resolves the distribution of interaction forces into low (~150 pN) and high (~330 pN) components, with the former more susceptible to salt concentration and to alteration of key charged residues on the RC. Thus, the low force component is suggested to reflect the contribution of electrostatic interactions in forming the initial encounter complex, whereas the high force component reflects the additional stabilization provided by hydrophobic interactions to the productive ET complex. Employing molecular dynamics simulations, we resolve five intermediate states that comprise the encounter, productive ET and leaving complexes, predicting a weak interaction between cyt c 2 and the LH1 ring near the RC-L subunit that could lie along the exit path for oxidized cyt c 2 . The multimodal nature of the interactions of ET complexes captured here may have wider implications for ET in all domains of life.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chain Dynamics in Magnetorheological Suspensions

Magnetorheological (MR) suspensions are composed of colloidal particles which acquire dipole moments when subjected to an external magnetic field. At sufficient field strengths and concentrations, the dipolar particles rapidly aggregate to form long chains. Subsequent lateral cross-linking of the dipolar chains is responsible for a rapid liquid-to-solid-like rheological transition. The unique, magnetically-activated rheological properties of MR suspensions make them ideal for interfacing mechanical systems to electronic controls. Additionally, the ability to experimentally probe colloidal suspensions interacting through tunable anisotropic potentials is of fundamental interest. Our current experimental work has focused on understanding the fluctuations of dipolar chains. It has been proposed by Halsey and Toor (HT) that the strong Landau-Peierls thermal fluctuations of dipolar chains could be responsible for long-range attractions between chains. Such interactions will govern the long-time relaxation of MR suspensions. We have synthesized monodisperse neutrally buoyant MR suspensions by density matching stabilized ferrofluid emulsion droplets with D2O. This allows us to probe the dynamics of the dipolar chains using light scattering without gravitational, interfacial, and polydispersity effects to resolve the short-wavelength dynamics of the dipolar chains. We used diffusing wave spectroscopy to measure these dynamics. The particle displacements at short times that show an independence to the field strength, but at long times exhibit a constrained, sub-diffusive motion that slows as the dipole strength is increased. The experiments are in good qualitative agreement with Brownian dynamics simulations of dipolar chains. Although there have been several important and detailed studies of the structure and interactions in MR suspensions, there has not been conclusive evidence that supports or contradicts the HT model prediction that long-range interactions exist between fluctuating chains of dipolar particles. Resolving this issue would contribute greatly to the understanding of these interesting and important materials. We have begun to test the predictions of the HT model by both examining the dynamics of individual chains and by measuring the forces between dipolar chains directly to accurately and quantitatively assess the interactions that they experience. To do so, we employ optical trapping techniques and video-microscopy to manipulate and observe our samples on the microscopic level. With these techniques, it is possible to observe chains that are fluctuating freely in three-dimensions, independent of interfacial effects. More importantly, we are able to controllably observe the interactions of two chains at various separations to measure the force-distance profile. The techniques also allow us to study the mechanical properties of individual chains and chain clusters. Our work to this point has focused on reversibly-formed dipolar chains due to field induced dipoles where the combination of this chaining, the dipolar forces, and the hydrodynamic interactions that dictate the rheology of the suspensions. One can envision, however, many situations where optical, electronic, or rheological behavior may be optimized with magneto-responsive anisotropic particles. Chains of polarizable particles may have the best properties as they can coil and flex in the absence of a field and stiffen and orient when a field is applied. We have recently demonstrated a synthesis of stable, permanent paramagnetic chains by both covalently and physically linking paramagnetic colloidal particles. The method employed allows us to create monodisperse chains of controlled length. We observed the stability, field-alignment, and rigidity of this new class of materials. The chains may exhibit unique rheological properties in an applied magnetic field over isotropic suspensions of paramagnetic particles. They are also useful rheological models as bead-spring systems. These chains form the basis for our current experiments with optical traps.

Gast, A. P.↗

MAGNETICALLY SUSPENDED CANNED ROTOR PUMPS FOR THE INTEGRAL MOLTEN SALT REACTOR

This project developed a molten salt pump design for small modular reactors, concentrated solar, and Gen IV nuclear reactors with a magnetically levitated rotor that can operate at temperatures of up to 700 °C. This eliminates the need for rotating seals and roller element bearings which require maintenance and are prone to failure. The pump design is also more compact than existing molten salt pump designs which require long shafts to thermally isolate the electric motor. The aims of the project are to: 1) design and test the molten salt pump, 2) develop and test a high temperature position sensor needed for magnetic bearing control, and 3) develop and test high temperature coil fabrication methods for electric motor and magnetic bearing fabrication. To develop the high temperature coils, conductor materials were analyzed for oxidation, diffusion, cladding, and resistivity. Of the candidate conductor materials that met the requirements, silver was chosen. The long-term impact of radiation on silver resistivity was analyzed by developing a model of a generic molten salt fast reactor to model the conductor irradiation in worst case scenarios. The analysis showed that increased resistivity due to transmutation of silver decreases exponentially with the reactor reflector thickness. Different electrical insulation materials were tested and the final high temperature coils developed were able to withstand conductor to ground voltages up to 1800 V before insulation breakdown. Magnetic bearings require high-precision high-speed measurements of the rotor position to operate. Current commercial position sensors can only operate up to 550 °C so a novel position sensor was developed that can operate up to 800 °C. These sensors were also designed to measure the position of the metallic rotor through a thin metallic corrosion barrier that protects that coils, insulation, and magnetic cores so that the salt containment barrier does not have any penetrations. The high-temperature sensor was tested and had comparable performance to commercially available inductive position sensors. The sensor design and signal processing has led to several invention disclosures. The magnetically levitated molten salt pump design was completed and analyzed. The chloride salt test loop at Oak Ridge National Laboratory was chosen to test the pump and the pump hydraulics and power were designed for the test loop flow and pressure requirements. Classes of motor designs that do not utilize permanent magnetic were considered for the high-temperature motor. After extensive analysis of their performance, a synchronous reluctance motor was chosen for the pump. The design was refined and analyzed to optimize performance. The magnetic bearings were designed to meet the force and frequency requirements necessary to levitate the rotor and reject disturbances from rotor imbalance, motor forces, and hydraulic forces. Compliant mounts for the position sensors were designed that would not damage the ceramic sensor components and maintain sensor alignment over a wide temperature range. A thermal hydraulic analysis of the molten salt coolant flow used to cool the rotor and stator was performed to ensure that the maximum temperature in the rotor and stator will not exceed 725 °C. The corrosion barrier and casing were designed to maintain the salt fluid boundary by using compliant features to absorb large expansion mismatches due to different coefficients of thermal expansion (CTE). Finally, the salt drain and fill systems were designed along with the pump support structure, integration with the test loop, and thermal insulation and heating. Corrosion barrier fabrication experiments were performed to assess using laser cladding to apply a sub-millimeter corrosion barrier to the surface of a Fe-Co-V magnetic core material. The cladding method showed minimal mixing between materials at their interface and withstood temperature cycling without damage to the cladding. The technology commercialization is focused on the high-temperature position sensors and there has been significant interest from the nuclear industry in utilizing these sensors along with some interest from the aerospace industry. The project also led to a research collaboration with a fusion energy industrial partner to study barriers to scaling the pump design from the kW size to the MW size. This research collaboration will also study methods to create high-temperature coils that have conductor to ground breakdown voltages above 5 kV. The successful development of high-temperature electromagnetic coils and high-temperature position sensors greatly increases the temperature limits for electromagnetic devices that don’t rely on permanent magnets. This includes electric motors, magnetic bearings, linear actuators, rotary position and velocity sensing, actuated valves, and generators to name a few. This will fill a current need in molten salt reactors, concentrated solar, and fusion energy for these critical peripheral devices needed for practical reactor designs that do not currently have commercial solutions.

Hines, J Wesley↗