A SERIES OF MONTE CARLO CODES TO TRANSPORT NUCLEONS THROUGH MATTER
Monte carlo method for nucleon transport through spacecraft shielding
SEARCH · Engineering Papers
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.
Monte carlo method for nucleon transport through spacecraft shielding
Transport calculations for proton shielding
Bibliographies compiled and data evaluated on diffusivity, viscosity, transport phenomena, and turbulent flow
Quadrupole effects on gas diffusion, viscosity and thermal diffusion factors, computing collision integrals
Previous work reported from JPL has included characterization of electrode kinetics and alkali atom transport from electrodes including Mo, W, WRh(sub x), WPt(sub x)(Mn), in sodium AMTEC cells and vapor exposure cells, and Mo in potassium vapor exposure cells. These studies were generally performed in cells with small area electrodes (about 1 to 5 cm(sup 2)), and device geometry had little effect on transport. Alkali diffusion coefficients through these electrodes have been characterized, and approximate surface diffusion coefficients derived in cases of activated transport. A basic model of electrode kinetic at the alkali metal vapor/porous metal electrode/alkali beta'-alumina solid electrolyte three phase boundary has been proposed which accounts for electrochemical reaction rates with a collision frequency near the three phase boundary and tunneling from the porous electrode partially covered with adsorbed alkali metal atoms. The small electrode effect in AMTEC cells has been discussed in several papers, but quantitative investigations have described only the overall effect and the important contribution of electrolyte resistance. The quantitative characterization of transport losses in cells with large area electrodes has been limited to simulations of large area electrode effects, or characterization of transport losses from large area electrodes with significant longitudinal temperature gradients. This paper describes new investigations of electrochemical kinetics and transport, particularily with WPt(sub 3.5) electrodes, including the influence of electrode size on the mass transport loss in the AMTEC cell. These electrodes possess excellent sodium transport properties making verification of device limitations on transport much more readily attained.
The interaction energies for N2-He and N2-H2 are calculated by accurate ab initio methods. The virial coefficient and differential scattering cross section for N2-H2 are calculated; the theoretical results are compared with experimental data. The transport collision integrals for N2-H2 and N2-N2 interactions are calculated and tabulated; the results yield transport coefficients that compare well with measured data. Transport coefficients are found to be determined accurately from the interaction energies for a specific configuration of the molecule formed from the interaction partners. Comparisons with results of measurement and accurate calculations demonstrate that the transport properties of complex molecular interactions can be determined rapidly and fairly accurately from the interaction energies of simpler system using combination rules for the short-range parameters of effective interaction energies and the coefficients for the long-range forces. The coefficients for a two-parameter temperature expansion of diffusion and viscosity are tabulated for a realistic universal potential energy that is based primarily on the results of very accurate calculations of the He-He interaction energy.
Here, we present a novel spectroscopy approach to investigate impurity transport by analyzing line-radiation following high-n Rydberg transitions. While high-n Rydberg states of impurity ions are unlikely to be populated via impact excitation, they can be accessed by charge exchange (CX) reactions along the neutral beams in high-temperature plasmas. Hence, localized radiation of highly ionized impurities, free of passive contributions, can be observed at multiple wavelengths in the visible range. For the analysis and modeling of the observed Rydberg transitions, a technique for calculating effective emission coefficients is presented that can well reproduce the energy dependence seen in datasets available on the OPEN-ADAS database. By using the rate coefficients and comparing modeling results with the new high-n Rydberg CX measurements, impurity transport coefficients are determined with well-documented 2σ confidence intervals for the first time. This demonstrates that high-n Rydberg spectroscopy provides important constraints on the determination of impurity transport coefficients. By additionally considering Bolometer measurements, which provide constraints on the overall impurity emissivity and, therefore, impurity densities, error bars can be reduced even further.
A review of the important constraints on gravity wave induced diffusion of chemical tracers, heat and momentum is given. Ground-based microwave spectroscopy measurements of H2O and CO and rocket-based mass spectrometer measurements of Ar constrain the eddy diffusion coefficient for constituent transport (K sub zz) to be (1-3) x 10 to the 5th sq cm/sec in the upper mesosphere. Atomic oxygen data also limits K sub zz to a comparable value in the mesopause. From the energy balance of the upper mesosphere the eddy diffusion coefficient for heat transport (D sub H) is, at most, 6 x 10 to the 5th sq cm/sec at the mesopause and decreasing substantially with decreasing altitude. The available evidence for mean wind deceleration and the corresponding eddy diffusion coefficient for momentum stresses (D sub M) suggests that it is at least 1 x 10 to the 6th sq cm/sec in the upper mesosphere. Consequently the eddy Prandtl number for macroscopic scale lengths is greater than 3.
A review of the important constraints on gravity wave induced diffusion of chemical tracers, heat and momentum is given. Ground-based microwave spectroscopy measurements of H20 and CO and rocket-based mass spectrometer measurements of Ar constrain the eddy diffusion coefficient for constituent transport (K sub zz) to be (1-3) x 10 to the 5th sq cm/sec in the upper mesosphere. Atomic oxygen data also limits K sub zz to a comparable value in the mesopause. From the energy balance of the upper mesosphere the eddy diffusion coefficient for heat transport (D sub H) is at most, 6 x 10 to the 5th sq cm/sec at the mesopause and decreasing substantially with decreasing altitude. The available evidence for mean wind deceleration and the corresponding eddy diffusion coefficient for momentum stresses (D sub M) suggests that it is at least 1 x 10 to the 6th sq cm/sec in the upper mesosphere. Consequently the eddy Prandtl number for macroscopic scale lengths is greater than 3.
This article reviews recent molecular simulation studies of "collective" properties of aqueous electrolyte solutions, specifically free energies and activity coefficients, solubilities, nucleation rates of crystals, and transport coefficients. These are important fundamental properties for biology and geoscience, but also relevant for many technological applications. Their determination from molecular-scale calculations requires large systems and long sampling times, as well as specialized sampling algorithms. As a result, such properties have not typically been taken into account during optimization of force field parameters; thus, they provide stringent tests for the transferability and range of applicability of proposed molecular models. There has been significant progress on simulation algorithms to enable the determination of these properties with good statistical uncertainties. Comparisons of simulation results to experimental data reveal deficiencies shared by many commonly used models. Moreover, there appear to exist specific tradeoffs within existing modeling frameworks, so that good prediction of some properties is linked to poor prediction for specific other properties. For example, non-polarizable models that utilize full charges on the ions generally fail to predict accurately both activity coefficients and solubilities; the concentration dependence of viscosity and diffusivity for these models is also incorrect. Scaled-charge models improve the dynamic properties and could also perform well for solubilities, but fail in the prediction of nucleation rates. Even models that do well at room temperature for some properties generally fail to capture their experimentally observed temperature dependence. Finally, the main conclusion from the present review is that qualitatively new physics will need to be incorporated in future models of electrolyte solutions to allow description of collective properties for broad ranges of concentrations, temperatures, and solvent conditions.
In the “method of four coefficients,” electrical resistivity (ρ), Seebeck coefficient (S), Hall coefficient (RH), and Nernst coefficient (Q) of a material are measured and typically fit or modeled with theoretical expressions based on Boltzmann transport theory to glean experimental insights into features of electronic structure and/or charge carrier scattering mechanisms in materials. Although well-defined and readily available reference materials exist for validating measurements of ρ and S, none currently exists for R H or Q. We show that measurements of all four transport coefficients—ρ, S, R H , and Q—can be validated using a single reference sample, namely, the low-temperature Seebeck coefficient Standard Reference Material® (SRM) 3451 (composition Bi 2 Te 3+x ) available from the National Institute for Standards and Technology (NIST) without the need for inter-laboratory sample exchange. Here, R H and Q data for NIST SRM 3451 reported here for the temperature range 80–400 K complement the data already available for ρ and S and will therefore be of interest to researchers desiring to validate new or existing galvanomagnetic and thermomagnetic transport properties measurement systems.
Alignment charts (nomographs) designed to permit rapid calculations of the low-pressure transport properties of nonpolar gases and gas mixtures are presented. Calculations for pure gases are based on the rigorous kinetic theory of gases as applied to a realistic intermolecular force law. Mixture viscosities and conductivities are calculated from good approximations derived from rigorous theory. The properties are calculable with a precision of 2 percent or better. Accuracy depends on how well the constants characterizing the Intermolecular force law are known; if constants are derived from experimental data, results should be accurate to 5 percent or better. Force constants are tabulated for 65 gases.
The exchange current, transfer coefficient, mass-transport parameters, and electrode capacitance at the Na(g)/porous Mo/Na-Beta-double prime alumina solid electrolyte (BASE) phase boundary have been evaluated from 740 to 1220 K. The transfer coefficient exhibits a value close to 0.5 and the exchange current is dominated by collision frequency, with no significant activation energy. Since the porous Mp-electrode adopts a fairly regular microstructure on the BASE surface, the magnitude of the exchange current of mature electrodes directly depends on the actual contact zone of the porous metal film with the BASE ceramic, and decreases slightly as grain growth occurs. The exchange currents and the mass-transport parameters derived for very porous, thin Mo electrodes indicate that the charge-transfer reaction occurs at a small fraction of the interface. High-frequency limiting capacitance and resistance values due to the interface show potential dependence and a value on the order of 1 F/sq m and 0.1-1.0 Ohm-sq cm.
Recombination, coefficients vertical diffusion and transport velocities calculated for ionospheric F layer
Previous work reported from JPL has included characterization of electrode kinetics and alkali atom transport from electrodes including Mo, W, WRh x (Mn), in sodium AMTEC cells and vapor exposure cells; and Mo in potassium vapor exposure cells. These studies were generally performed in cells with small area electrodes (about 1 to 5 cm 2 ), and device geometry had little effect on transport. Alkali metal diffusion coefficients through these electrodes have been characterized, and approximate surface diffusion coefficients were derived in cases of activated transport. A basic model of electrode kinetics and transport at the alkali metal vapor/pourous metal electrode/alkali beta"-alumina solid electrolyte (BASE) three phase boundary has been proposed which accounts for electrochemical reaction rates with a collision frequency near the three phase boundary and tunneling from the porous electrode partially covered with adsorbed alkali metal atoms. The small electrode effect in AMTEC cells has been discussed in several papers, but quantitative investigations have described only the overall effect and the important contribution of electrolyte resistance. The quantitative characterization of transport losses in cells with large area electrodes has been limited to simulations of large area electrode effects, or characterization of transport losses from large area electrodes with significant longitudinal temperature gradients. This paper describes new investigations of electrochemical kinetics and transport, with four 14.4 cm 2 WPt 3.5 electrodes, including the influence of electrode size on the mass transport loss in the AMTEC cell. These electrodes exhibit very slow sintering, as well as excellent sodium transport properties, making them attractive candidates for AMTEC power conversion use. However, the facile sodium transport in WPt 3.5 electrodes makes characterization of the transport process difficult.
The methanol/stainless steel vapor chambers (flat-plate heat pipes) discussed in this paper were developed for use in spaceborne atmospheric cloud chambers. This application imposed stringent thermal and mechanical requirements on the design. Flatness, low thermal mass, vibration, and structural integrity requirements were achieved in addition to precision temperature uniformity and thermal transport. Heat transfer coefficients on the order of 0.34 to 0.40 W/sq cm -C were measured. The vapor chambers are capable of transporting 170 W-cm per cm of width in either the axial or side-to-side direction.
Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.
Two dimensional modeling has become an established technique for the simulation of the global structure of trace constituents. Such models are simpler to formulate and cheaper to operate than three dimensional general circulation models, while avoiding some of the gross simplifications of one dimensional models. Nevertheless, the parameterization of eddy fluxes required in a 2-D model is not a trivial problem. This fact has apparently led some to interpret the shortcomings of existing 2-D models as indicating that the parameterization procedure is wrong in principle. There are grounds to believe that these shortcomings result primarily from incorrect implementations of the predictions of eddy transport theory and that a properly based parameterization may provide a good basis for atmospheric modeling. The existence of these GCM-derived coefficients affords an unprecedented opportunity to test the validity of the flux-gradient parameterization. To this end, a zonally averaged (2-D) model was developed, using these coefficients in the transport parameterization. Results from this model for a number of contrived tracer experiments were compared with the parent GCM. The generally good agreement substantially validates the flus-gradient parameterization, and thus the basic principle of 2-D modeling.