Search NASA⌕ Search

SEARCH · Search NASA

Results for “HEAT FLUX”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Shadow masks predictions in SPARC tokamak plasma-facing components using HEAT code and machine learning methods

Here, this work uses machine learning (ML) to complement HEAT (Heat flux Engineering Analysis Toolkit) by developing 3-D footprint surrogate models for fast and accurate heat load calculations in the divertor of the SPARC tokamak. The focus is on shadowed regions, or magnetic shadows, caused by the 3-D geometry of plasma-facing components (PFCs). ML classifiers are employed to create a surrogate model for HEAT generated shadow masks, predicting these shadow masks and divertor heat flux profiles based on a diverse range of equilibria and only the plasma current, safety factor(q95) at the edge, and magnetic flux angles as input parameters. The ultimate goal is to integrate the model for real-time control and future operational decisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evaluation of SPARC divertor conditions in H-mode operation using SOLPS-ITER

The predicted divertor conditions for the SPARC tokamak are calculated using SOLPS-ITER for a range of scrape-off-layer (SOL) heat flux widths λ q , input powers, and particle fueling locations. Under H-mode scenario conditions with an upstream separatrix density of 1 x 10 20 m -3 , the most conservative range of λ q extrapolations ( 0.15 mm) results in extremely high unmitigated particle and energy fluxes to the divertor, both under full field (12.2 T) and power (P SOL = 29 MW) conditions, and 2/3 field with P SOL = 10 MW. Increasing the cross-field SOL diffusivities by 2–10× reduces the magnitude of the mitigation challenge, however strategies such as impurity seeding or strike-point-sweeping will likely still be required. A combination of steady-state and time-dependent SOLPS-ITER simulations are used to map out phase space diagrams of upstream and divertor conditions. The simulations include parallel currents but neglect cross-field drifts. At low upstream density the inner and outer divertor conditions are highly asymmetric, with a large temperature difference and significant heat fluxes driven by parallel currents. The solution has sharp bifurcations with a region of hysteresis, depending on whether the initial state is at a low or high density. This behavior is observed even when the fueling location, cross-field diffusivity, and impurity level is changed, although the density window with asymmetry is reduced with increasing diffusivity. The addition of neon impurity seeding reduces the divertor heat fluxes, but also causes a drop in the upstream electron density with fixed particle throughput. This drop can be counteracted by increased main ion throughput, however too much neon results in a back transition into the asymmetric divertor regimes suggesting a need for control of both main ion and impurity seeding levels to achieve a desired divertor state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantifying heat and particle flux to primary and secondary divertors for various types of edge-localized-modes

Between 15% and 30% of edge-localized-mode (ELM) heat flux can be deposited to regions outside the main divertor including the far scrape-off-layer (SOL), private flux region (PFR), and secondary divertor inner target. Analysis shown here demonstrates that type-I and small ELM plasmas are transported to the PFR and to the secondary inner divertor, which is magnetically isolated from the outer divertor, leading to well-defined heat flux profiles and with peak values that can surpass those at the secondary outer target. Such features are consistently observed for pedestal collisionalities from $ν$$^{*}_{e}$ ~1.5 to $ν$$^{*}_{e}$ ~3.9. Heat flux profiles of the examined ELM types feature rippled structures due to bursts in the outer far-SOL region but not on the secondary inner target, causing long decay lengths in the time-averaged ELM profiles. The contribution of each ELM type to the total time-averaged power deposited to the secondary divertor has been evaluated, showing that grassy ELMs contribute ~8%, small ELMs ~67%, and type-I ELMs ~85%. These findings imply that small ELMs may yet pose a concern for future machines if some regions of the main wall are not designed to withstand significant heat and particle fluxes. Due to the low intra-ELM heat flux contribution, however, the grassy ELM regime is an attractive option for an ELMing scenario in future machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗
DOE OSTI2020

Azimuthally averaged wind and thermodynamic structures of tropical cyclones in global climate models and their sensitivity to horizontal resolution

Characteristics of tropical cyclones (TCs) in global climate models (GCMs) are known to be influenced by details of the model configurations, including horizontal resolution and parameterization schemes. Understanding model-to-model differences in TC characteristics is a prerequisite for reducing uncertainty in future TC activity projections by GCMs. This study performs a process-level examination of TC structures in eight GCM simulations that span a range of horizontal resolutions from 1° to 0.25°. A recently developed set of process-oriented diagnostics is used to examine the azimuthally averaged wind and thermodynamic structures of the GCM-simulated TCs. Results indicate that the inner-core wind structures of simulated TCs are more strongly constrained by the horizontal resolutions of the models than are the thermodynamic structures of those TCs. As expected, the structures of TC circulations become more realistic with smaller horizontal grid spacing, such that the radii of maximum wind (RMW) become smaller, and the maximum vertical velocities occur off the center. However, the RMWs are still too large, especially at higher intensities, and there are rising motions occurring at the storm centers, inconsistently with observations. The distributions of precipitation, moisture, radiative and surface turbulent heat fluxes around TCs are diverse, even across models with similar horizontal resolutions. At the same horizontal resolution, models that produce greater rainfall in the inner-core regions tend to simulate stronger TCs. Further, when TCs are weak, the radial gradient of net column radiative flux convergence is comparable to that of surface turbulent heat fluxes, emphasizing the importance of cloud-radiative feedbacks during the early developmental phases of TCs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Hypervelocity impact in stellar media: heat shielding, shock fronts and ablation clouds

The Frontier Plasma Science experiments proposed by our team will model the processes occurring during spacecraft atmospheric entries and meteorite planetary collisions by inserting material targets into the Scrape off Layer (SOL) and edge plasma of a tokamak. This project will investigate important plasma-material processes during high-enthalpy atmospheric entries, including heat flux generation and material ablation. Of particular interest is to explore heat and particle flux detachment in front of the sample using concepts and techniques developed in the latest studies of divertor detachment. Exploration missions to the Solar System’s gaseous giants and hyperbolic re-entries into the Earth’s atmosphere require spacecraft that can withstand high velocity (>10 km/s), high heat fluxes (>10 MW/m2), and corresponding enthalpies. Ablative materials have been used as thermal shields to protect the spacecraft from severe heating during entry. However, developing high enthalpy ablating materials is challenging due to the lack of adequate ground testing facilities. Our team stated that the performance of candidate shielding materials could be assessed in a laboratory environment with multiple in-situ diagnostics, which was enabled by the progress in the tokamak research. Modern machines feature relatively long discharges (~10 s) with well-controlled stable plasma conditions at the edge where the heat flux and the flow speed are similar to those experienced during atmospheric entries. The experiments were carried out at the DIII-D National Fusion Facility in San Diego, CA, operated by General Atomics (GA) for the DoE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Medium-Scale Methanol Pool Fire Model Validation

In this work, medium scale (30 cm diameter) methanol pool fires were simulated using the latest fire modeling suite implemented in Sierra/Fuego, a low Mach number multiphysics reacting flow code. The sensitivity of model outputs to various model parameters was studied with the objective of providing model validation. This work also assesses model performance relative to other recently published large eddy simulations (LES) of the same validation case. Two pool surface boundary conditions were simulated. The first was a prescribed fuel mass flux and the second used an algorithm to predict mass flux based on a mass and energy balance at the fuel surface. Gray gas radiation model parameters (absorption coefficients and gas radiation sources) were varied to assess radiant heat losses to the surroundings and pool surface. The radiation model was calibrated by comparing the simulated radiant fraction of the plume to experimental data. The effects of mesh resolution were also quantified starting with a grid resolution representative of engineering type fire calculations and then uniformly refining that mesh in the plume region. Simulation data were compared to experimental data collected at the University of Waterloo and the National Institute of Standards and Technology (NIST). Validation data included plume temperature, radial and axial velocities, velocity temperature turbulent correlations, velocity velocity turbulent correlations, radiant and convective heat fluxes to the pool surface, and plume radiant fraction. Additional analyses were performed in the pool boundary layer to assess simulated flame anchoring and the effect on convective heat fluxes. This work assesses the capability of the latest Fuego physics and chemistry model suite and provides additional insight into pool fire modeling for nonluminous, non-sooting flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

DIII-D High Field Side Lower Hybrid Current Drive Experiment Overview

In preparation for high field side lower hybrid current drive (HFS LHCD) experiments in DIII-D, the HFS LHCD was to be commissioned and physics experiments commence once the system operated up to 300 kW for 0.5 s. The initial physics experiments sought to characterize coupling, wave propagation, and driven current measurements. In HFS LHCD first campaign, the maximum power was limited to <200 kW due to waveguide pressure leaks limiting the number of available modules and power per module. Here we summarize commissioning progress and initial physics observations. The HFS LHCD coupler was optimized for high qmin, DIII-D discharges where efficient off-axis current at r/a~0.6-0.8 is desired. The coupler n|| spectrum is peaked at 2.7 and is predicted to generate ~0.14 MA/MW coupled for 1.6 MW injected. In preparation, the HFS scrape-off layer density profile was characterized and found to have steeper profiles and lower fluctuation levels than the low field side. Furthermore, the HFS SOL density profile can be accurately predicted using global plasma quantities using machine learning. Thus far, one module has injected ~100 kW for 0.5 s with <5% reflected power. Nonthermal electrons have been observed on lower frequency channels of the electron cyclotron emission radiometer correlated with the LH power indicating core wave absorption. To avoid 30R neutral beam heat flux, a split launcher is proposed to avoid the high heat flux region while maintaining power spectrum and directivity.

Wukitch, Stephen J. [Massachusetts Inst. of Techno↗
DOE OSTI2025

CFD Simulations of Molten Salts Within Plasma-Facing Components with Twisted Tape Inserts Under One-Sided–Heated Conditions

For liquid immersion blankets (LIBs), the heat removal within different plasma-facing components (PFCs) requires substantial heat removal capabilities that could benefit from the usage of heat transfer enhancements. Of consideration are twisted tape insert heat transfer enhancements that are ribbon-like metal pieces inserted within cooling channels that increase fluid mixing, resulting in higher heat transfer and pressure drop. Unlike traditional applications, PFCs such as divertor coolant tubes experience a significantly higher heat flux on one side versus the other or entirely one-sided heating. These conditions are not well studied for molten salts where the higher Prandtl numbers characteristic of the fluids are significantly higher than water, gases, and liquid metals. As an effort to study this knowledge gap, large eddy simulations were done with Prandtl numbers representative of FLiBe at proposed operating conditions. In this study, the effects of applying a surface heat flux to one side of a tube versus the entire tube with a twisted tape insert were investigated for the heat transfer coefficient and effects on the temperature profile. The findings of the study suggest that for Prandtl number fluids above unity, the Nusselt number was largely unaffected by the change in boundary condition. Both boundary conditions showed good agreement with the Manglik and Bergles Nusselt number correlation for tight-fitting twisted tapes. Finally, the outcomes of this study will be useful for future design studies of LIBs to provide fundamental guidance implementing twisted tape inserts and their performance in one-sided–heated conditions.

molten salts↗
DOE OSTI2021

Strong influence of convective heat transfer efficiency on the cooling benefits of green roof irrigation

Abstract By enhancing evapotranspiration (ET), green roofs provide cooling benefits for the urban environment and are recognized as a promising heat mitigation strategy. The evaporative cooling effects of green roofs strongly depend on the soil moisture conditions and thus irrigation may be needed to sustain the cooling benefits. It has been shown that the magnitude of cooling benefits offered by green roof irrigation varies spatially, but its controlling factors remain elusive. In this study, we combine a surface energy balance (SEB) model with global simulations generated by an improved Earth System Model to illustrate the key factors controlling the cooling benefits of green roof irrigation. We employ a simple irrigation scheme, which is only active when there is no ice in soil layers and when the soil moisture is below field capacity. As a result, most of the irrigation water leaves the green roof system via ET. We find that the magnitude and also the spatial variability of the cooling benefits of green roof irrigation are controlled by the irrigation amount, and a surface energy redistribution factor that encodes the efficiencies of different SEB components in transferring heat. Further analysis indicates that the enhancement of latent heat flux due to irrigation is largely balanced by the reduction of sensible heat flux on green roofs. Therefore, the amount of irrigation needed per unit decrease of green roof surface temperature is mainly controlled by the convective heat transfer efficiency. A lower convective heat transfer efficiency (e.g. under a lower wind speed) helps reduce the amount of irrigation needed per unit decrease of green roof surface temperature. This study highlights the importance of SEB in constraining the cooling benefits of green roof irrigation and provides valuable guidance for urban planners and policy makers in terms of heat mitigation and sustainable water management.

Wang, Linying (ORCID:0000000308769220)↗
DOE OSTI2021

A Case for Electron-Astrophysics

Abstract The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the research theme of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. In addition, plasma electrons often play an important role for the spatial transfer of thermal energy due to the high heat flux associated with their velocity distribution. The regulation of this electron heat flux is likewise not understood. By focussing on these and other fundamental electron processes, the research theme of electron-astrophysics links outstanding science questions of great importance to the fields of space physics, astrophysics, and laboratory plasma physics. In this White Paper, submitted to ESA in response to the Voyage 2050 call, we review a selection of these outstanding questions, discuss their importance, and present a roadmap for answering them through novel space-mission concepts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Development of experimental and computational frameworks to predict subcooled flow boiling in the LANL Isotope Production Facility

Cooling is crucial to maintain the integrity of target systems in isotope production facilities. At Los Alamos National Laboratory (LANL)’s Isotope Production Facility (IPF), multiple encapsulated targets are stacked and irradiated in tandem with a 100 MeV, ~250μA proton beam. To facilitate effective heat removal, these stacked targets are separated and cooled via a series of water channels. At these beam currents, this high-energy proton beam heats the target system, likely initiating subcooled flow boiling in the cooling channels. However, in-beam monitoring of the IPF target system is not possible due to the extreme radiation environment, and the necessarily significant shielding. To better understand high-power target performance, we developed ex-situ experimental and computational frameworks to predict the behavior of subcooled flow boiling at IPF. Subcooled flow boiling experiments on Inconel 625 samples under IPF conditions (2 bar pressure, 10 GPM flow rate (i.e., 2249 kg/m 2 /s), 85 K subcooling) revealed that IPF's average operating power is at the early stage of boiling with a heat transfer coefficient of 48,000 W/m 2 /s. The proposed modeling framework enables us to predict a complete boiling curve, i.e., single-phase heat transfer, onset of nucleate boiling, two-phase heat transfer, and critical heat flux (CHF), with specification of input boiling parameters up to intermediate heat flux levels. The estimated CHF under IPF conditions is 5.2 MW/m 2 . Experimental data under reduced conditions (2 bar pressure, 1.5 GPM flow rate (i.e., 337 kg/m 2 /s), 45 K subcooling) served as validation cases for the computational modeling. This computational model can be further extended to more complicated systems replicating the real IPF configuration, for instance, to study void distribution as a function of the incident proton beam profile and coolant velocity profile of multiple cooling channels. Finally, the proposed experimental and computational frameworks provide a means to better understand cooling systems in the isotope production facilities at different accelerators, where in-beam monitoring of the cooling process is not available.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

A Machine Learning-Assisted Framework to Control Thermally Anisotropic Building Envelopes in Residential Buildings

To curb the energy consumption of buildings and their related CO2 emissions, Oak Ridge National Laboratory (ORNL) has developed the thermally anisotropic building envelope (TABE) —a multi-layer design comprising insulation materials and metal foils connected to thermal loops. In this study, a machine learning-assisted framework was developed to control the TABE in residential buildings to reduce the computation load for future optimal rule-based control and application. First, a 2D finite element model was established in COMSOL to calculate the hourly heat flux through exterior walls installed with the TABE. Then, TABE wall heat fluxes were simulated for various given indoor and outdoor boundary conditions, as well as thermal loops fluid temperatures and flow rates. Since the finite element simulations are computationally expensive, an artificial neural network (ANN) was then trained to use as a proxy of the finite element (COMSOL) modeling. Finally, the trained ANN model was coupled with the EnergyPlus model to predict the energy consumption of a US Department of Energy prototype single-family house installed with the TABE. An optimal simple rule-based control was determined from predefined rules for a case study. The results demonstrate that the developed machine learning–assisted framework can reduce 99.9% of the computation time while efficiently managing residential building energy for installed TABE walls.

Shen, Zhenglai↗
DOE OSTI2026

Thermodynamics-guided machine learning model for predicting convective boundary layer height and its multi-site applicability

Accurate estimation of convective boundary layer height (CBLH) is vital for weather, climate, and air quality modeling. Machine learning (ML) shows promise in CBLH prediction, but input parameter selection often lacks physical grounding, limiting generalizability. This study introduces a novel ML framework for CBLH prediction, integrating thermodynamic constraints and the diurnal CBLH cycle as an implicit physical guide. Boundary layer growth is modeled as driven by surface heat fluxes and atmospheric heat absorption represented with the low tropospheric stability, using the diurnal cycle as input and output. TPOT and AutoKeras are employed to select optimal models, validated against Doppler lidar-derived CBLH data, achieving an R 2 of 0.84 across untrained years. Comparisons of eddy covariance (ECOR) and energy balance Bowen ratio (EBBR) flux measurements show the same prediction capability. Models trained on the ARM SGP C1 site with ECOR data and tested at E37 and E39 yield R 2 values of 0.79 and 0.81, respectively, demonstrating their adaptability. The ML model trained with all sites' data slightly enhances the performance compared with ML models trained over single-site data. The interquartile range for predicted CBLH is consistently narrower than that for DL-derived CBLH, reflecting lower variability in predicted CBLH compared to DL-derived CBLH, which is influenced by additional factors, which are not well represented with the model inputs. The model's generalizability across multiple sites at the ARM SGP site demonstrates its potential for transfer to greater distances, offering a scalable approach for enhancing boundary layer parameterization in atmospheric models.

Chu, Yufei [Stony Brook Univ., NY (United States)]↗
DOE OSTI2021

Observation of fully detached divertor integrated with improved core confinement for tokamak fusion plasmas

Integration of divertor detachment with a high-performance (β N ~ 3, β p > 2, H 98 ~ 1.5) core plasma has been demonstrated in DIII-D high-β p (poloidal beta) plasmas associated with a sustained core internal transport barrier (ITB) and an H-mode edge transport barrier (ETB). Such good core-edge integration has been achieved for both neon and nitrogen seeding, for both favorable and unfavorable B-field directions, independently from the impurity puffing locations, though these variations play important roles on divertor characteristics. Compared to the standard H-mode plasmas, the high-β p plasma exhibits a much wider window of detachment compatible with high confinement core. Fully detached divertor plasmas with low plasma temperature (Te < 5 eV), low particle flux, and low heat flux across the entire divertor target plate were obtained by using nitrogen seeding. This detached high-β p plasma is compatible with a newly developed detachment control system which can help optimize the nitrogen gas flow rate. Several features, i.e., the high edge safety factor in the high-β p scenario, impurity injection, closed divertor and reduced heating power requirement due to the high confinement, facilitate the achievement of full divertor detachment at lower density. Instead of degrading global performance, the divertor detachment facilitates the access to an even stronger ITB at large radius with a relatively weak ETB through self-organized synergy between ITB and ETB, leading to sustained high confinement. The strengthening of the large-radius ITB compensates for the ETB degradation associated with divertor detachment. In addition, a weak ETB naturally has smaller edge localized modes (ELMs). In particular, with neon injection, a long-period no-ELM H-mode phase has been achieved simultaneously with high-performance core and partially detached divertor plasmas. Furthermore, these results demonstrate the possibility of integrating excellent core plasma performance with an effective divertor solution, an essential step toward steady-state operation of reactor-grade plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Real-time gas cooling of flowing liquid lithium limiter for the EAST

A novel continuously flowing liquid lithium limiter (FLiLi) which employs an in-vessel electro-magnetic pump to drive liquid Li flowing on the surface of limiter has been successfully designed and tested in the experimental advanced superconducting tokamak (EAST) device in 2014. In order to better control the surface temperature, an upgraded design and the real-time gas cooling were performed in 2016. Two kinds of cooling gas, helium (He) and argon (Ar) were tested on the FLiLi system prior to the experiment. It was found that the cooling rates of He and Ar at 2.5 MPa are 34.2 °C/min, 10.2 °C/min or 31.7 °C/kL and 28.1 °C/kL, respectively. The cooling performance of He is more effective than Ar. Experimental results show that the real-time He cooling could effectively contain the FLiLi surface temperature increase and prevent the strong passive evaporation of lithium during plasma discharges. Limiter heat flux up to 0.2 MW/m2 was removed during ohmic discharge with 2.5 MPa He gas. Finally, Because of the upgraded design and real-time gas cooling, the surface of FLiLi was not damaged by heat flux after entire plasma discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Thermal and magnetic evolution of a crystallizing basal magma ocean in Earth's mantle

We present the thermochemical evolution of a downward crystallizing basal magma ocean (BMO) overlying the liquid outer core and probe its capability to dissipate enough power to generate and sustain an early dynamo. Furthermore, a total of 61 out of 112 scenarios for a BMO with imposed, present-day heat flux, Q BMO , values of 15, 18, and 21 TW and radiogenic heat, Q r , values of 4, 8, and 12 TW fully crystallized during the age of the Earth. Most of these models are energetically capable of inducing magnetic activity for the first 1.5 Gyrs, at least, with durations extending to 2.5 Gyrs; with final core-mantle boundary (CMB) temperatures of 4400 ±500 K-well within current best estimates for inferred temperatures. None of the models with Q BMO = 12 TW achieved a fully crystallized state, which may reflect a lower bound on the present-day heat flux across the CMB. BMO-powered dynamos exhibit strong dependence on the partition coefficient of iron into the liquid layer and its associated melting-point depression for a lower mantle composition at near-CMB conditions-parameters which are poorly constrained to date. Nonetheless, we show that a crystallizing BMO is a plausible mechanism to sustain an early magnetic field.

58 GEOSCIENCES↗