Search NASA⌕ Search

SEARCH · Search NASA

Results for “Effective diffusion coefficient”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Analytical prediction of moisture absorption/desorption in resin matrix composites exposed to aircraft environments

The moisture absorption/desorption behavior of resin matrix composites was mathematically modeled by classical diffusion theory using an effective diffusion coefficient. Good agreement was found between calculated moisture content and published data for T300/5208 graphite fiber reinforced epoxy matrix composite. Weather Bureau data for Langley Air Force Base and Norfolk, Va., were used to calculate the amount of moisture a T300/5208 composite panel would contain if exposed outdoors. Results obtained by using average monthly weather data for several high aircraft usage locations around the world suggest that, except for desert areas, geographical locations should have only minimal effect on the moisture absorption level reached in composites. Solar radiation data together with cloud and wind information were included in the analysis to estimate an effective temperature of the composite panel during ground exposure.

Unnam, J.↗
DOE OSTI2023

Hydrogen Diffusion Coefficient Measures on Thin Film Uranium Oxide

Thin films of uranium oxides, putatively UO 2 and U 3 O 8 , were deposited on palladium/silver (75/25) foil discs (10 μm thick, 10 mm diameter). Methodology for sealing these foils, applying a hydrogen pressure (~ 1 atm) to the oxide side, and measuring the pressure on the permeate side (opposite side of foil) is reported. The experimental apparatus is held at 100 °C to speed diffusion, aided by the pressure differential across the foil (~1 atm (750 torr) on oxide side, initially low vacuum on permeate side, ~1x 10 -4 torr). Early results indicate an effective diffusion coefficient of roughly 7.66 x 10 -17 +/- 2.22 x 10 -17 cm 2 /sec for UO 2 . These values are in line with expectations and prior measures relative to large lag times for system baseline (substrate only). Evaluation of the technique for U 3 O 8 thin films suggests further development will be needed to extend the technique to more oxidized films (U 3 O 8 , UO 3 ). Characterization of the foils (thickness and speciation) post diffusion experiments will be carried out by SIGMA (Eric Tegtmeier and Andy Richards) in FY24 which will sharpen the uncertainty quantification for UO 2 diffusion coefficients, aiding the nucleation model for DRACO.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Effects of thermal-solutal convection on temperature and solutal fields under various gravitational orientations

Semiconductor crystals such as Hg(1-x)CD(x)Te grown by unidirectional solidification Bridgmann method have shown compositional segregations in both the axial and radial directions. Due to the wide separation between the liquidus and the solidus of its pseudobinary phase diagram, there is a diffusion layer of higher HgTe content built up in the melt near the melt-solid interface which gives a solute concentration gradient in the axial direction. The value of effective diffusion coefficient calculated from fitting of the data to 1D model varies with Hg(1-x)Cd(x)Te growth conditions. This indicates that the growth condition of the Hg(1-x)Cd(x)Te is not purely diffusion controlled. Because of the higher thermal conductivity in the melt than that in the crystal in the growth system, there is a thermal leakage through the fused silica crucible wall near the melt-solid interface. This gives a thermal gradient in the radial direction. Hart, and Thorpe, Hutt and Soulsby have shown that under such conditon a fluid will become convectively unstable as a result of different diffusitivities of temperature and solute. It is quite important to understand the effects of this thermosolute convection on the compositonal segregation in both axial and radial directions in the unidirectionally solidified crystals under various gravitational directions. To reach this goal, we start with a simplified problem to study the effects of thermal-solutal convection on the temperature and solutal fields under various gravitional orientations. We begin by reviewing model governing equations.

Wang, Jai-Ching↗
DOE OSTI2025

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗
NASA NTRS1976

Modeling of three-dimensional mixing and reacting ducted flows

A computer code, based upon a finite element solution algorithm, was developed to solve the governing equations for three-dimensional, reacting boundary region, and constant area ducted flow fields. Effective diffusion coefficients are employed to allow analyses of turbulent, transitional or laminar flows. The code was used to investigate mixing and reacting hydrogen jets injected from multiple orifices, transverse and parallel to a supersonic air stream. Computational results provide a three-dimensional description of velocity, temperature, and species-concentration fields downstream of injection. Experimental data for eight cases covering different injection conditions and geometries were modeled using mixing length theory (MLT). These results were used as a baseline for examining the relative merits of other mixing models. Calculations were made using a two-equation turbulence model (k+d) and comparisons were made between experiment and mixing length theory predictions. The k+d model shows only a slight improvement in predictive capability over MLT. Results of an examination of the effect of tensorial transport coefficients on mass and momentum field distribution are also presented. Solutions demonstrating the ability of the code to model ducted flows and parallel strut injection are presented and discussed.

Zelazny, S. W.↗
DOE OSTI2023

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

HTO and selenate diffusion through compacted Na-, Na–Ca-, and Ca-montmorillonite

Radionuclide transport in smectite clay barrier systems used for nuclear waste disposal is controlled by diffusion, with adsorption significantly retarding transport rates. While a relatively minor component of spent nuclear fuel, 79 Se is a major driver of the safety case for spent fuel disposal due to its long half-life (3.3×10 5 yr) and its low adsorption to clay (K D < 10 L/kg), thus a thorough understanding of Se diffusion through clay is critical for understanding the long-term safety of spent fuel disposal systems. Through-diffusion experiments with tritiated water (HTO, conservative tracer) and Se(VI) were conducted with a well-characterized, purified montmorillonite source clay (SWy-2) under a constant ionic strength (0.1 M) and three different electrolyte compositions: Na + , Ca 2+ , and a Na + -Ca 2+ mixture at pH 6.5 in order to probe the effects of electrolyte composition and interlayer cation composition on clay microstructure, Se(VI) aqueous speciation, and ultimately diffusion. Further, the results were modeled using a reactive transport modeling approach to determine values of porosity (ε), D e (effective diffusion coefficient), and K D (distribution coefficient for adsorption). HTO diffusive flux was higher in Ca-montmorillonite (D e =1.68×10 -10 m 2 s -1 ) compared to Na-montmorillonite (De=7.83×10 -11 m 2 s -1 ). This increase in flux is likely due to a greater degree of clay layer stacking in the presence of Ca 2+ compared to Na + , which leads to larger inter-particle pores. Overall, the Se(VI) flux was much lower than the HTO flux due to anion exclusion, with Se(VI) flux following the order Ca (D e = 1.03×10 -11 m 2 s -1 ) > Na–Ca (D e = 2.12×10 -12 m 2 s -1 ) > Na (D e = 1.28×10 -12 m 2 s -1 ). These differences in Se(VI) flux are due to a combination of factors, including (1) larger accessible porosity in Ca-montmorillonite due to clay layer stacking and smaller electrostatic effects compared to Na-montmorillonite, (2) larger accessible porosity for neutral-charge CaSeO4 species which makes up 32% of aqueous Se(VI) in the pure Ca system, and (3) possibly higher Se(VI) adsorption for Ca-montmorillonite. Through a combination of experimental and modeling work, this study highlights the compounding effects that electrolyte and counterion compositions can have on radionuclide transport through clay. Diffusion models that neglect these effects are not transferable from laboratory experimental conditions to in situ repository conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A time-dependent chloride diffusion model for predicting initial corrosion time of reinforced concrete with slag addition

The effect of granulated blast furnace slag (GBFS) addition on the threshold chloride concentration (TCC) for rebar corrosion was investigated. A modified diffusion model, coupled with a time-dependent effective diffusion coefficient and surface chloride concentration, was proposed to predict the chloride profile. The corrosion initiation time was estimated based on the model predictions and the measured TCC. The results indicate that adding GBFS decreases the TCC by lowering the pH value of the pore solution. The evolution of corrosion potential and current density is found to obey a 3-parameter Weibull distribution. MIP tests show that adding GBFS contributes to refinement of pore structure by decreasing the fraction of large capillary pores. The time-dependent model exhibits good predictive strength and helps understand how GBFS addition delays the corrosion initiation by retarding the chloride diffusion, though a lower TCC is obtained.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Hydrogen and deuterium permeation in Hastelloy N

Hastelloy N was chosen as the fluoride salt-contacting structural material for the Molten Salt Reactor Experiment due to its excellent compatibility with the fuel salt FLiBe. FLiBe is currently investigated for several advanced fusion and fission reactor concepts where tritium generation in the FLiBe is anticipated. Knowledge of hydrogen transport properties through Hastelloy N is important to understand how tritium would permeate through this material and result in an unintentional release. In this study, the hydrogen and deuterium permeability, diffusivity, and solubility were measured from 500 to 700 °C at a primary-side pressure of 10 kPa in a well-characterized sample of Hastelloy N. The prepared polycrystalline Hastelloy N had C and O impurities present on the surface. These impurities were investigated using Auger Emission Spectroscopy and Ar depth profiling. Here, the adventitious C was removed upon the first Ar sputter cycle whereas O persisted deeper into the sample. For permeation experiments, applied deuterium pressures ranged from 13 Pa to 100 kPa and deuterium transport remained in the diffusion-limited regime (J ∝ P 0.5 ) throughout the pressure range examined. Two methods are employed to measure the effective hydrogen and deuterium diffusivity: rise and decline. The decline method produced improved statistical model fits for calculating the effective diffusion coefficient compared to the rise method. The resultant transport properties compared well to published values for other nickel alloys.

08 HYDROGEN↗
DOE OSTI2025

Experimental Characterization of Hydrogen Diffusion in Shale Rocks for Geologic Storage Applications

As global energy systems undergo a transition to cleaner alternatives, geologic hydrogen storage has emerged as a promising solution for large-scale energy storage. A critical factor in determining the feasibility of this approach is the effectiveness of caprock formations, such as shale, in preventing hydrogen migration. This study investigates the diffusion behavior of hydrogen through shale to assess its suitability as a caprock for geologic hydrogen storage. Using a novel double-seal core holder design and a through-diffusion apparatus, hydrogen diffusion was measured through shale rock from the Eagle Ford and Wolfcamp Formations under dry conditions. These measurements were complemented by microstructural and mineralogical analyses using low-pressure nitrogen adsorption and X-ray diffraction. The effective diffusion coefficient of hydrogen in these shale caprocks ranged from 2.51 × 10 –8 to 9.85 × 10 –8 m 2 /s. Notably, we observed that the diffusion behavior was more related to the pore network structure and could not be attributed to differences in the total pore volume between shale types alone. Here, to further understand the role of pore network complexity, a fractal pore model was developed to correlate tortuosity with the fractal dimension of the pore structure (a measure of pore network complexity). The proposed model closely matched tortuosity values obtained from diffusion experiments, outperforming existing theoretical tortuosity–porosity correlations. These findings provide key quantitative parameters needed to assess the feasibility of geologic hydrogen storage as well as insights that can be applied to hydrogen storage in a range of geologic formations.

08 HYDROGEN↗
DOE OSTI2021

Impact of sodium vanadium oxide (NaV 3 O 8 , NVO) material synthesis conditions on charge storage mechanism in Zn-ion aqueous batteries

The electrochemical charge storage of sodium vanadate (NaV 3 O 8 or NVO) cathodes in aqueous Zn-ion batteries has been hypothesized to be influenced by the inclusion of structural water for facilitating ion transfer in the material. Materials properties considered important (morphology, crystallite and particle size, surface area) are systematically studied herein through investigation of two NVO materials, NaV 3 O 8 ·0.34H 2 O [NVO(300)] and NaV 3 O 8 ·0.05H 2 O [NVO(500)], with different water content, acicular morphologies with different size and surface area achieved via post-synthesis heat treatment. The electrochemistry of the two materials was evaluated in aqueous Zn-ion cells with 2 M ZnSO 4 electrolyte using cyclic voltammetry, galvanostatic cycling, and rate capability testing. Furthermore, the thinner NVO(300) nanobelts (0.13 μm) demonstrate greater specific capacities and higher effective diffusion coefficients relative to the thicker NVO(500) nanorods. Notably however, while cells containing NVO(500) deliver lower specific capacity, they demonstrate enhanced capacity retention with cycling. The structural changes accompanying oxidation and reduction are elucidated via ex situ X-ray diffraction, transmission electron microscopy, and operando V K-edge X-ray absorption spectroscopy (XAS), where NVO material properties are shown to influence the ion insertion. Operando XAS verified that electron transfer corresponds directly to change in vanadium oxidation state, affirming vanadium redox as the governing electrochemical process.

25 ENERGY STORAGE↗
DOE OSTI2024

Resistive drift wave turbulence and anomalous transport of multi-species plasma

Anomalous transport of multi-species plasma is considered with the generalized Hasegawa–Wakatani model. It is shown that the transport of all plasma species is described by fractional diffusion equations with the same effective diffusion coefficient. Strongly enhanced perturbations of heavy impurity density are found in long-living plasma flow vortices.

Physics↗
DOE OSTI2021

Atomic dynamics of metallic glass melts La 50 Ni 15 Al 35 and Ce 70 Cu 19 Al 11 studied by quasielastic neutron scattering

By employing quasielastic neutron scattering, we studied the atomic-scale relaxation dynamics and transport mechanism of La 50 Ni 15 Al 35 and Ce 70 Cu 19 Al 11 metallic glass melts in the temperature range of > 200K above their liquidus temperatures. The results show that both liquids exhibit stretched exponential relaxation and Arrhenius-type temperature dependence of the effective diffusion coefficient. The La 50 Ni 15 Al 35 melt exhibits an activation energy of 0.545 ± 0.008 eV and a stretching exponent ~0.77 to 0.86 in the studied temperature range; no change of activation energy, as suggested in previous reports, associated with liquid-liquid phase transition was observed. In contrast, the Ce 70 Cu 19 Al 11 melt exhibits larger diffusivity with a much smaller activation energy of 0.201±0.003eV and a smaller stretching exponent ~0.51 to 0.60, suggestive of more heterogeneous dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Electrochemistry of the NaI-AlBr 3 Molten Salt System: A Redox-Active, Low-Temperature Molten Salt Electrolyte

NaI-AlBr 3 is a very appealing low melting temperature (<100 °C), salt system for use as an electrochemically-active electrolyte. This system was investigated for its electrochemical and physical properties with focus to energy storage considerations. A simple phase diagram was generated; at >100 °C, lower NaI concentrations had two partially miscible liquid phases, while higher NaI concentrations had solid particles. Considering the fully molten regime, electrical conductivities were evaluated over 5–25 mol% NaI and 110 °C–140 °C. Conductivities of 6.8–38.9 mS cm −1 were observed, increasing with temperature and NaI concentration. Effective diffusion coefficients of the I − /I 3 − redox species were found to decrease with both increasing NaI concentration and increasing applied potential. Regardless, oxidation current density at 3.6 V vs Na/Na + was observed to increase with increasing NaI concentration over 5–25 mol%. Finally, the critical interface between the molten salt electrolyte and electrode materials was found to significantly affect reaction kinetics. When carbon was used instead of tungsten, an adsorbed species, most likely I 2 , blocked surface sites and significantly decreased current densities at high potentials. This study shows the NaI-AlBr 3 system offers an attractive, low-temperature molten salt electrolyte that could be useful to many applied systems, though composition and electrode material must be considered.

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

Hydraulic Properties Data Package for the E-Area Soils, Cementitious Materials, and Waste Zones–Update

Hydraulic property estimates for the soils, cementitious materials, and waste zones associated with the E-Area low-level radioactive waste disposal units have been provided to support the Performance Assessments (PA) for the E-Area Low-Level Waste Facility (LLWF). Nominal or “best estimate” hydraulic property values for use in the deterministic modeling are provided along with representations of the hydraulic property value uncertainty for use in sensitivity and uncertainty modeling. The hydraulic properties provided for each of the E-Area materials include: porosity ( η ), dry bulk density ( ρ b ), particle density (ρ p ), saturated hydraulic conductivity ( K sat ), characteristic curves (suction head, saturation, and relative permeability), and effective diffusion coefficient ( D e ). A representation of the uncertainty associated with each property, except for the characteristic curves, is provided for each material, except for the E-Area waste zones. These nominal parameter values for each of the E-Area soils, cementitious materials, and waste zones, and, where indicated, parameter uncertainty representations, are based upon the following in order of priority: Site-specific field data, Site-specific laboratory data, Similarity to material with site-specific field or laboratory data, and Literature data. Additionally, a methodology to represent long-term concrete degradation is provided. Finally, because much of the nominal hydraulic property values and uncertainty representations for the E-Area soils, cementitious materials, and waste zones are based on similarity to other materials or literature data, a methodology to prioritize additional work to better define these values and representations is outlined. Prioritization should be based on the importance of the material and/or property to the results of deterministic, sensitivity, and uncertainty modeling. This prioritization should be established through a process of sensitivity modeling. This report supersedes the previous hydraulic properties data package. Revision 1 contains clarifications, updates and corrections to text, figures and tables in a number of places in the report along with some minor editing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1976

The radial variation of solar flare proton anisotropies observed in deep space on Pioneers 10 and 11

The anisotropy of solar-flare protons (greater than 3 MeV) was evaluated at the time of maximum for each of 24 events measured by Pioneers 10 and 11 and was observed to decrease gradually between 1 and 4.7 AU. Since diffusive streaming is related directly to the effective diffusion coefficient parallel to the interplanetary spiral magnetic field, these observations provide a measure of the radial dependence of interplanetary diffusion parameters. Analysis of these results suggests that the spiral field increases gradually with heliocentric distance. An expression is given for the radial dependence of that field.

Mccarthy, J.↗