Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A flux-limited treatment for the conductive evaporation of spherical interstellar gas clouds

In this work, we present and analyze a new analytic solution for the saturated (flux-limited) thermal evaporation of a spherical cloud. This work is distinguished from earlier analytic studies by allowing the thermal conductivity to change continuously from a diffusive to a saturated form, in a manner usually employed only in numerical calculations. This closed form solution will be of interest as a computational benchmark. Using our calculated temperature profiles and mass-loss rates, we model the thermal evaporation of such a cloud under typical interstellar medium (ISM) conditions, with some restrictions. We examine the ionization structure of the cloud-ISM interface and evaluate column densities of carbon, nitrogen, oxygen, neon, and silicon ions toward the cloud. In accord with other investigations, we find that ionization equilibrium is far from satisfied under the assumed conditions. Since the inclusion of saturation effects in the heat flux narrows the thermal interface relative to its classical structure, we also find that saturation effects tend to lower predicted column densities.

Dalton, William W.↗
DOE OSTI2023

Analysis of ducted fuel injection at high-pressure transcritical conditions using large-eddy simulations

Ducted fuel injection (DFI) is a proposed fuel injection concept for achieving substantial reductions in emissions. In this concept, the fuel is injected through a coannular duct, resulting in increased fuel-air mixing and minimized formation of soot and other unwanted combustion products. Apart from comprehensive experimental investigations on DFI, so far computational studies have been limited to single-point Reynolds-averaged Navier Stokes simulations. Therefore, the objective of this work is to complement these studies by performing large-eddy simulations using a diffuse-interface method to examine the physical mechanisms and combustion processes of DFI, specifically focusing on the mixing process and the effect of fuel-ducting on combustion and pollutant emissions. To this end, finite-rate chemistry simulations are performed of the DFI configuration corresponding to the Engine Combustion Network Spray A injector at transcritical conditions (n-dodecane fuel, 60 bar pressure and 1000 K temperature chamber conditions). A two-equation soot model is employed for the qualitative analysis of soot emissions. Direct comparisons of averaged and instantaneous flow field results with the Spray A configuration are performed to assess the effect of DFI on the first- and second-stage ignition and soot formation. Compared to the free-spray condition, the results show that the DFI case exhibits a combination of (i) increased mass flow rate and entrained air, (ii) larger pressure drop magnitude and flow velocity, and (iii) a closer-to-stoichiometric mixture composition (both globally and locally), each of which is conjectured to contribute toward reduced soot production.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1998

Experiment 1: Orbital Processing of High-Quality Zn-Alloyed CdTe Compound Semiconductors

Zinc-alloyed cadmium telluride (Zn:CdTe) crystals were grown in unit gravity and in microgravity for comparative analysis. Two crystals were grown on USML-1/STS-50 in 1992, and two additional crystals were grown on USML-2/STS-73 in 1995. The Crystal Growth Furnace (CGF) in the seeded Bridgman-Stockbarger crystal growth geometry was utilized on both missions. Crystals grown on USML-1/STS-50 were found to have solidified with partial wall contact due to the near-absence of the hydrostatic pressure in microgravity, a residual g-vector that was not axial, and the nonwetting sample/ampoule wetting conditions. Crystals grown on USML-2/STS-73 included: a sample/ampoule identical to the USML-1/STS-73 sample/ampoule, with the addition of a restraining spring to simulate hydrostatic pressure internally, and a tapered ampoule which accomplished 2.2 cm of crystal growth without wall contact. Infrared transmission of all ground and flight samples was found to be 63-66 percent, very close to the theoretical 66 percent, suggesting good stoichiometric control. Infrared microscopy confirmed that the primary precipitates were Te and their size (1-10 microns) and density suggested that the flight- and ground-base samples experienced similar thermal histories. Longitudinal macrosegregation, calculated using scaling analysis, was predicted to be low. Nearly diffusion controlled growth was achieved even in unit gravity and macrosegregation data could be fit with a diffusion controlled model. Radial segregation was monitored and was found to vary with fraction solidified, particularly through the shoulder region, where the sample cross section was varying significantly. It was also disturbed in the flight samples in regions where asymmetric wall contact was noted. In regions where a steady state was established, the radial segregation was invariant within our experimental measurement error. Flight samples were found to be much higher in structural perfection than samples processed in unit gravity under identical growth conditions. In regions where solidification had occurred without wall contact, the free surfaces evidenced virtually no twinning, although twins appeared in the flight samples in regions of wall contact and were pervasive in the ground samples. These results were confirmed using optical microscopy and synchrotron x-ray white beam topography. Full-width half-maximum rocking curve widths (FWHM), recorded in arc-seconds, were significantly reduced from 20 a-s (1 g) to 9 a-s (micro-g) for the best regions of the crystals. The 9 a-s (FWHM) rocking curve value in the unconfined flight samples equals the best value reported terrestrially for this material. The ground samples exhibited a fully developed (111)[110] dislocation mosaic structure, whereas dislocations within the flight samples were discrete and no mosaic structure was evident. The defect density was quantitatively reduced from 75,000 (1 g) to 800 (micro-g) +/-50 percent. Dislocation etch pit density results were confirmed using transmission synchrotron white beam and monochromated beam topography. The low defect density is thought to have resulted from the near absence of hydrostatic pressure, which allowed the molten boule to solidify with little or no wall contact. This minimized the transfer of hoop stresses during solidification and post-solidification processing.

Larson, D. J., Jr.↗
DOE OSTI2026

Multi-physics Preconditioning for Thermally Activated Batteries

Thermal batteries, also known as molten-salt batteries, are single-use reserve power systems activated by pyrotechnic heat generation, which transitions the solid electrolyte into a molten state. The simulation of these batteries relies on multiphysics modeling to evaluate performance and behavior under various conditions. This paper presents advancements in scalable preconditioning strategies for the Thermally Activated Battery Simulator (TABS) tool, enabling efficient solutions to the coupled electrochemical systems that dominate computational costs in thermal battery simulations. We propose a hierarchical block Gauss-Seidel preconditioner implemented through the Teko package in Trilinos, which effectively addresses the challenges posed by tightly coupled physics, including charge transport, porous flow, and species diffusion. The preconditioner leverages scalable subblock solvers, including smoothed aggregation algebraic multigrid (SA-AMG) methods and domain-decomposition techniques, to achieve robust convergence and parallel scalability. Strong and weak scaling studies demonstrate the solver’s ability to handle problem sizes up to 51.3 million degrees of freedom on 2048 processors, achieving near sub-second setup and solve times for the end-to-end electrochemical solve. These advancements significantly improve the computational efficiency and turnaround time of thermal battery simulations, paving the way for higher-resolution models and enabling the transition from 2D axisymmetric to full 3D simulations.

25 ENERGY STORAGE↗
DOE OSTI2026

Dataset_for_Conserved_macromolecular_architecture_of_Poplar_secondary_cell_walls_revealed_by_ssNMR_and_atomistic_modeling

This dataset contains solid-state 13C NMR data and atomistic molecular dynamics simulation files supporting the study of nanoscale secondary cell wall architecture across 13 genetically diverse Populus trichocarpa genotypes grown under uniform greenhouse conditions in 13C-enriched CO2 atmospheres (~89% 13C enrichment).The dataset contains two collections of solid-state 13C NMR data. (1) 200 MHz data (Bruker Avance III HD, 4 mm HX probe, 10 kHz MAS): raw Bruker TopSpin experiment folders and DMFIT-exported ascii spectra for selective and non-selective 1D 13C-13C spin diffusion experiments (3000 ms mixing) used to quantify inter-polymer spatial proximities, and short-mixing (1 ms) reference spectra used for polymeric abundance quantification by spectral deconvolution. (2) 600 MHz data (Bruker Avance III, 1.6 mm PhoenixNMR HXY probe, 30 kHz MAS): raw Bruker TopSpin experiment folders containing 2D CORD, 2D CP-INADEQUATE, and 13C/1H relaxation (T1, T1rho) experiments for all 13 genotypes, with processed Excel workbooks per experiment type. Molecular dynamics simulation code, coordinate files, and analysis scripts (NAMD/CHARMM/Python) for six atomistic cell wall models are included. Summarized ssNMR data are compiled into a single excel file and subjected to statistical analysis. Multivariate analysis code (PCA, Pearson correlation) and summary data are provided as excel worksheets and Jupyter notebooks (Python 3).

09 BIOMASS FUELS↗
DOE OSTI2019

Highly resolved Eulerian simulations of fuel spray transients in single and multi-hole injectors: Nozzle flow and near-exit dynamics

In high pressure fuel injectors, needle opening and closing transients cause complex off-design fluid dynamics behaviors that profoundly impact the spray and mixture formation processes. These dynamics are completely different from what is known to occur in steady state conditions. In this study, diesel spray transients have been investigated in single-hole and 3-hole nozzles, encompassing internal and external nozzle flow and including needle motion, performing highly resolved (2.5 μm) computational fluid dynamics (CFD) simulations. We focused on end-of-injection (EOI) and start-of-injection (SOI) processes, in order to provide insights in to the physics. The liquid fuel, vapor and gas species are modeled with a single-fluid multiphase mixture approach, with diffuse interface, and with large eddy simulations (LES) of the turbulence. Occurrence of phase change due to cavitation is accounted for, and the spray dispersion is described with a turbulent dispersion model. Detailed needle motion data and orifice internal surface are available from xray synchrotron source measurements carried out at Argonne National Laboratory, and shared through the Engine Combustion Network (ECN) community. Simulations are compared against x-ray phase contrast imaging and radiography of the internal and near-exit flow, in addition to optical microscopy data of the near-exit sprays. Simulation results are found to agree well with available experimental data, and are able to realistically capture local and global features. The simulations allow to gain insight into the physics of gas ingestion and dribbles at EOI, for different hole diameters, operating conditions and number of holes. At SOI, timing of liquid appearance out of the injector and spray tip penetration are adequately predicted, by using the EOI flow field as in-nozzle initialization, and by prescribing the measured tip needle displacement with an informed effective valve opening point inferred from the x-ray observations. Lastly, the variation of spreading angle over time is also discussed in detail for the multi-hole case, including hole-to-hole variations. Due to real geometry features and asymmetric needle motion with eccentric components, it is found that the three holes exhibit swirling flows of increasing intensity as the lift decreases, causing the near cone angle to open and spread, in a quasi-hollow cone structure. These features are not observed in axial single-hole injectors because of their relative simplicity and intrinsic symmetry.

42 ENGINEERING↗
DOE OSTI2021

A NIR fluorescent smart probe for imaging tumor hypoxia

Abstract Background Tumor hypoxia is a characteristic of paramount importance due to low oxygenation levels in tissue negatively correlating with resistance to traditional therapies. The ability to noninvasively identify such could provide for personalized treatment(s) and enhance survival rates. Accordingly, we recently developed an NIR fluorescent hypoxia‐sensitive smart probe ( NO 2 ‐Rosol ) for identifying hypoxia via selectively imaging nitroreductase (NTR) activity, which could correlate to oxygen deprivation levels in cells, thereby serving as a proxy. We demonstrated proof of concept by subjecting a glioblastoma (GBM) cell line to extreme stress by evaluating such under radiobiological hypoxic ( p O 2 ≤ ~0.5%) conditions, which is a far cry from representative levels for hypoxia for brain glioma ( p O 2 = ~1.7%) which fluctuate little from physiological hypoxic ( p O 2 = 1.0‐3.0%) conditions. Aim We aimed to evaluate the robustness, suitability, and feasibility of NO 2 ‐Rosol for imaging hypoxia in vitro and in vivo via assessing NTR activity in diverse GBM models under relevant oxygenation levels ( p O 2 = 2.0%) within physiological hypoxic conditions that mimic oxygenation levels in GBM tumor tissue in the brain. Methods We evaluated multiple GBM cell lines to determine their relative sensitivity to oxygenation levels via measuring carbonic anhydrase IX (CAIX) levels, which is a surrogate marker for indirectly identifying hypoxia by reporting on oxygen deprivation levels and upregulated NTR activity. We evaluated for hypoxia via measuring NTR activity when employing NO 2 ‐Rosol in in vitro and tumor hypoxia imaging studies in vivo. Results The GBM39 cell line demonstrated the highest CAIX expression under hypoxic conditions representing that of GBM in the brain. NO 2 ‐Rosol displayed an 8‐fold fluorescence enhancement when evaluated in GBM39 cells ( p O 2 = 2.0%), thereby establishing its robustness and suitability for imaging hypoxia under relevant physiological conditions. We demonstrated the feasibility of NO 2 ‐Rosol to afford tumor hypoxia imaging in vivo via it demonstrating a tumor‐to‐background of 5 upon (i) diffusion throughout, (ii) bioreductive activation by NTR activity in, and (iii) retention within, GBM39 tumor tissue. Conclusion We established the robustness, suitability, and feasibility of NO 2 ‐Rosol for imaging hypoxia under relevant oxygenation levels in vitro and in vivo via assessing NTR activity in GBM39 models.

Hettie, Kenneth S.↗
DOE OSTI2013

An Efficient Algorithm for Mapping Imaging Data to 3D Unstructured Grids in Computational Biomechanics

Geometries for organ scale and multiscale simulations of organ function are now routinely derived from imaging data. However, medical images may also contain spatially heterogeneous information other than geometry that are relevant to such simulations either as initial conditions or in the form of model parameters. In this manuscript, we present an algorithm for the efficient and robust mapping of such data to imaging based unstructured polyhedral grids in parallel. We then illustrate the application of our mapping algorithm to three different mapping problems: 1) the mapping of MRI diffusion tensor data to an unstuctured ventricular grid; 2) the mapping of serial cyro-section histology data to an unstructured mouse brain grid; and 3) the mapping of CT-derived volumetric strain data to an unstructured multiscale lung grid. Execution times and parallel performance are reported for each case.

Einstein, Daniel R.↗
DOE OSTI2024

3D printed carbon fiber reinforced carbon as an energy efficient alternative to graphite for EFAS tooling

As Electric Field Assisted Sintering (EFAS) gains more industrial acceptance and use, it becomes more important to develop more efficient means to implement this technology. To this aim, 3D printed continuous carbon fiber reinforced carbon (CCC) was manufactured and fabricated into tooling for EFAS systems as an alternative to traditional graphite tooling. The impact of fiber orientation on the thermal and electrical properties of the CCC was characterized. Sample material was sintered in Tokai G535 graphite tooling, under common processing conditions and compared with CCC tooling. There was nearly 50 % energy savings compared to graphite while maintaining equivalent sample density and microstructure plus keeping ram temperatures 39 % cooler. This is due to spatial control of generated heat and thermal diffusivity within the molds, by means of fiber orientation anisotropy. Finite element modeling of the tooling design supported the experimental results as well as displays the effect of optimization of this 3D printed CCC material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-order partitioned spectral deferred correction solvers for multiphysics problems

We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1], [2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems with moderate coupling strength. They include a stiff linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1], [2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Ab Initio Thermodynamic Model for Magnesium Carbonates and Hydrates

An ab initio thermodynamic framework for predicting properties of hydrated magnesium carbonate minerals has been developed using density-functional theory linked to macroscopic thermodynamics through the experimental chemical potentials for MgO, water, and CO2. Including semiempirical dispersion via the Grimme method and small corrections to the generalized gradient approximation of Perdew, Burke, and Ernzerhof for the heat of formation yields a model with quantitative agreement for the benchmark minerals brucite, magnesite, nesquehonite, and hydromagnesite. The model shows how small differences in experimental conditions determine whether nesquehonite, hydromagnesite, or magnesite is the result of laboratory synthesis from carbonation of brucite, and what transformations are expected to occur on geological time scales. Because of the reliance on parameter-free first principles methods, the model is reliably extensible to experimental conditions not readily accessible to experiment and to any mineral composition for which the structure is known or can be hypothesized, including structures containing defects, substitutions, or transitional structures during solid state transformations induced by temperature changes or processes such as water, CO2, or O2 diffusion. Demonstrated applications of the ab initio thermodynamic framework include an independent means to evaluate differences in thermodynamic data for lansfordite, predicting the properties of Mg analogs of Ca-based hydrated carbonates monohydrocalcite and ikaite which have not been observed in nature, and an estimation of the thermodynamics of barringtonite from the stoichiometry and a single experimental observation.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2026

Towards a theoretical understanding of excitonic properties of phthalocyanine thin films. I. Low-temperature exciton absorption spectra

Phthalocyanine (Pc)-based molecular thin films have emerged in recent years as a promising class of organic semiconductor materials for optoelectronic applications owing to their long exciton coherence length and fast exciton diffusion. However, the dependence of their exciton properties on the dimensionality and thermodynamic conditions, presence of metal ions, and effects of chemical modifications to the PC systems is not yet fully understood. As a first step towards a more comprehensive theoretical understanding of the excitonic properties of Pc thin films, we model their low-temperature exciton absorption spectra by employing the Frenkel Hamiltonian. The latter is derived from quantum-chemical estimates of site energies and exciton–exciton couplings. The predicted exciton absorption spectra of octabutoxy phthalocyanine (H 2 OBPc) is found to be strongly dependent on the dimensionality of the model as well as the distance cutoff for including the monomer–monomer exciton coupling. We also caution that the widely used dipole–dipole approximation could substantially overestimate the excitonic coupling between different monomers compared to a more accurate evaluation using the respective transition densities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling Isotope Separation in Electrochemical Lithium Deposition

Naturally occurring Li consists of two stable isotopes, 6 Li with an abundance of about 7.5%, and 7 Li making up the remainder with 92.5%. The development of a 6 Li enrichment technique, in terms of technical reliability and environmental safety to reach 6 Li future requirements, represents a key step in the roadmap for nuclear fusion energy supply worldwide. This paper uses finite element analysis-based models to simulate electrochemical Li isotope separation, which is an attractive method in terms of simplicity, safety, and scalability. In the model, we quantitatively analyze how different electrochemical factors including thermodynamics, charge-transfer kinetics, and diffusivities affect the separation process (separation factor), together with cell parameters, such as cell length and current density. The maximum separation factor of 1.128 could be obtained with the cell under the optimal thermodynamic, kinetic, and diffusive conditions, which is among the highest separation factors ever reported. Furthermore, these results will assist in designing the actual isotope separation setup with large separation factor and appropriate timing for sample collection.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2017

High-Throughput Study of Diffusion and Phase Transformation Kinetics of Magnesium-Based Systems for Automotive Cast Magnesium Alloys

The objective of the proposed study is to establish a scientific foundation on kinetic modeling of diffusion, phase precipitation, and casting/solidification, in order to accelerate the design and optimization of cast magnesium (Mg) alloys for weight reduction of U.S. automotive fleet. The team has performed the following tasks: 1) study diffusion kinetics of various Mg-containing binary systems using high-throughput diffusion multiples to establish reliable diffusivity and mobility databases for the Mg-aluminum (Al)-zinc (Zn)-tin (Sn)-calcium (Ca)-strontium (Sr)-manganese (Mn) systems; 2) study the precipitation kinetics (nucleation, growth and coarsening) using both innovative dual-anneal diffusion multiples and cast model alloys to provide large amounts of kinetic data (including interfacial energy) and microstructure atlases to enable implementation of the Kampmann-Wagner numerical model to simulate phase transformation kinetics of non-spherical/non-cuboidal precipitates in Mg alloys; 3) implement a micromodel to take into account back diffusion in the solid phase in order to predict microstructure and microsegregation in multicomponent Mg alloys during dendritic solidification especially under high pressure die-casting (HPDC) conditions; and, 4) widely disseminate the data, knowledge and information using the Materials Genome Initiative infrastructure (http://www.mgidata.org) as well as publications and digital data sharing to enable researchers to identify new pathways/routes to better cast Mg alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2016

HST emission line galaxies at z ∼ 2: Comparing physical properties of Lyman alpha and optical emission line selected galaxies

We compare the physical and morphological properties of z ∼ 2 Lyα emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ∼ 2 optical emission line selected galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5<logM/M{sub ⊙}<10.5), size (0.5 < R < 3.0 kpc), and star formation rate (∼1<SFR<100 M{sub ⊙} yr{sup −1}). Remarkably, a comparison of the most commonly used physical and morphological parameters—stellar mass, half-light radius, UV slope, SFR, ellipticity, nearest neighbor distance, star formation surface density, specific SFR, [O iii] luminosity, and [O iii] equivalent width—reveals no statistically significant differences between the populations. This suggests that the processes and conditions which regulate the escape of Lyα from a z ∼ 2 star-forming galaxy do not depend on these quantities. In particular, the lack of dependence on the UV slope suggests that Lyα emission is not being significantly modulated by diffuse dust in the interstellar medium. We develop a simple model of Lyα emission that connects LAEs to all high-redshift star-forming galaxies where the escape of Lyα depends on the sightline through the galaxy. Using this model, we find that mean solid angle for Lyα escape is Ω{sub Lyα}=2.4±0.8 steradians; this value is consistent with those calculated from other studies.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1979

Energetic solar particle behaviour in the magnetosphere

The behavior of energetic solar flare particles in the magnetosphere is discussed. In the absence of magnetospheric motion, the problem of particle transport can be treated as simple propagation of charged particles in a stationary magnetic field configuration using, for instance, trajectory calculations in model fields. This single particle approach is the basis for the determination of intensity and anisotropy structures over the polar caps and in the geomagnetic tail from different interplanetary conditions. Particle transport on closed field lines is in addition strongly affected by resonant interaction processes as pitch angle scattering and radial diffusion.

Scholer, M.↗
NASA NTRS1980

Baffle aperture design study of hollow cathode equipped ion thrusters

A simple theoretical model which can be used as an aid in the design of the baffle aperture region of a hollow cathode equipped ion thruster was developed. An analysis of the ion and electron currents in both the main and cathode discharge chambers is presented. From this analysis a model of current flow through the aperture, which is required as an input to the design model, was developed. This model was verified experimentally. The dominant force driving electrons through the aperture was the force due to the electrical potential gradient. The diffusion process was modeled according to the Bolm diffusion theory. A number of simplifications were made to limit the amount of detailed plasma information required as input to the model to facilitate the use of the model in thruster design. This simplified model gave remarkably consistant results with experimental results obtained with a given thruster geometry over substantial changes in operating conditions. The model was uncertain to about a factor of two for different thruster cathode region geometries. The design usefulness was limited by this factor of two uncertainty and by the accuracy to which the plasma parameters required as inputs to the model were specified.

Brophy, J. R., Jr.↗