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.
263-page report addresses what advanced electrical and electronic spacecraft technologies can be applied in commercial airplanes. Report discusses a study in which technologies used in the Space Shuttle were evaluated for their potential in commercial transports. Describes new technologies, airplanes, tradeoffs and methods of evaluation. Concludes that major beneficiary would be an advanced wide-body transport (500 passenger).
A new method for in-situ Sb doping of CdTe that uses a modified vapor transport deposition system is described. This modification enables control of the Sb concentration with a pyrolysis stage to enhance the doping efficiency. CdTe:Sb films under different deposition conditions are characterized by SEM, XRD, and CV measurements for determining morphology, crystal structure, and hole concentration. Variations of the Sb dopant heater and pyrolyzer temperatures do not affect the CdTe morphology and crystal structure. However, CV measurements show that a higher dopant heater or pyrolyzer temperature leads to higher hole concentration. In this study, CdTe: Sb films achieve a hole concentration of 1016 cm-3 and 10% doping efficiency when the dopant heater is 600C and the pyrolyzer temperature is 1100C. This demonstrates a path to produce high hole concentration polycrystalline CdTe film with a low concentration of dopant-induced defect.
Griffin, a MOOSE-based reactor multiphysics code jointly developed by Idaho National Laboratory and Argonne National Laboratory under the DOE Office of Nuclear Energy’s NEAMS program, has pursued the development of an online multigroup cross section generation capability for a few years to enable high-fidelity, problem-dependent neutronics analyses of advanced thermal reactors. Recent advancements in Griffin’s online multigroup cross section generation capability have significantly improved the accuracy, robustness, and efficiency of self-shielding calculations for both prismatic and pebble-bed TRISO-fueled reactor applications. Key developments include a unified fuel self-shielding method applicable to both TRISO and annular compact/spherical shell fuel zone geometries; an advanced Dancoff Category-based Equivalence Theory using a bell function for non-fuel resonance treatment, achieving more than an order-of-magnitude speedup compared to the Tone method; an on-the-fly multigroup equivalence approach to mitigate group condensation errors; and a streaming correction method for pebble-bed homogenization. A proof-of-concept demonstration of on-the-fly group condensation with consistent P0 transport correction was also achieved. The method reproduced direct fine-group solutions with excellent accuracy (eigenvalue errors within 10 pcm and pin-power differences within 0.5%), but due to performance limitations of the current fixed-source solver, improvements to solver efficiency will be addressed in future work. Verification tests were performed on graphite-moderated TRISO-fueled two-dimensional core benchmark problems representing gas-cooled microreactors, heat pipe-cooled microreactors, gas-cooled pebble-bed reactors, and fluoride salt-cooled high-temperature reactors. Across all cases, Griffin showed excellent agreement with Serpent2 continuous energy Monte Carlo solutions: eigenvalue errors within 200 pcm, pin-power root-mean-square errors within 2%, and control rod and drum worth errors less than 2%. It should be noted that, for the benchmark problem, cross section generation contributed less than 3% of the total simulation times. These results demonstrate that Griffin’s online cross section generation capability delivers accurate and efficient reactor physics solutions across a wide spectrum of TRISO-fueled advanced reactor designs. With further improvements to the fine-group fixed-source solver and planned extensions to depletion, transients, and coupled neutron–gamma transport, Griffin will be well-positioned to become a powerful and comprehensive tool for advanced reactor analysis.
A description is given of the computational codes that are to be used in studies of the heating and aerodynamic performance of entry probes for the scientific exploration of the planets Jupiter, Saturn, and Uranus. These codes include a time-asymptotic numerical method which takes into account radiative transport and equilibrium chemistry for the direct solution of the radiative, inviscid flow over an axisymmetric or two-dimensional body. The method considers both supersonic and subsonic flow in the downstream region. Another code represents a laminar, viscous-shock-layer solution.
Neutrino quantum kinetics is a rapidly evolving field in computational astrophysics, with a primary focus on collective neutrino oscillations in core-collapse supernovae and postmerger phases of binary neutron star mergers. In recent years, there has been considerable debate concerning resolution dependence in numerical simulations. In this paper, we conduct a comprehensive resolution study in both angular- and spatial directions by using two independent schemes of quantum kinetic neutrino transport: finite volume and pseudospectral methods. We complement our discussion by linear stability analysis including inhomogeneous modes. Our result suggests that decreasing spatial resolutions underestimates the growth of flavor instability, and then leads to wrong asymptotic states of flavor conversions, which potentially has a critical impact on astrophysical consequences. We further delve into numerical results of low resolution simulations, that reveals the underlying mechanism responsible for numerical artifacts caused by insufficient resolutions. Finally, this study settles the debate on requirements of resolutions and serves as a guideline for numerical modeling of quantum kinetic neutrino transport.
The United States Department of Energy (DOE) Office of Nuclear Energy is applying knowledge and understanding in the areas of systems level engineering, analysis, and decision making to better inform the waste management pathways for U.S. spent nuclear fuel (SNF). Currently there are 93 operating and 23 shutdown commercial nuclear reactors in the United States. SNF at most of these locations is being stored in spent fuel pools, dry storage, or both [1, 2]. This paper initially goes over the basics of the agent-based simulation tool known as the Next Generation System Analysis Model (NGSAM) that has the capability to model interaction and movements of individual components or groups like casks, railcars etc. The next section covers some of the methods implemented in the Java Transportation Operations Model (JTOM) to facilitate the movement of assemblies, casks, railcars etc. and the various intermodal transfer options. This section also covers the current cask loading and intermodal transfer times implemented in NGSAM. The last section goes over the proposed values for intermodal transportation, and cask transfer/loading operations.
In this work, we present a method for accelerating discrete ordinates radiative transfer calculations for radiative transfer. Our method works with nonlinear positivity fixes, in contrast to most acceleration schemes. The method is based on the dynamic mode decomposition (DMD) and using a sequence of rank-one updates to compute the singular value decomposition needed for DMD. Using a sequential method allows us to automatically determine the number of solution vectors to include in the DMD acceleration. We present results for slab geometry discrete ordinates calculations with the standard temperature linearization. Compared with positive source iteration, our results demonstrate that our acceleration method reduces the number of transport sweeps required to solve the problem by a factor of about 3 on a standard diffusive Marshak wave problem, a factor of several thousand on a cooling problem where the effective scattering ratio approaches unity, and a factor of 20 improvement in a realistic, multimaterial radiating shock problem.
Presented herein is a method of including density fluctuations in the equations of turbulent transport. Results of a numerical analysis indicate that the method may be used to predict heat transfer for the case of near-critical para-hydrogen in turbulent upflow inside vertical tubes. Wall temperatures, heat transfer coefficients, and velocities obtained by coupling the equations of turbulent momentum and heat transfer with a perturbed equation of state show good agreement with experiment for inlet reduced pressures of 1.28-5.83.
An efficient method to account for the chemically frozen thermodynamic and transport properties of air in three-dimensional Navier-Stokes calculations has been demonstrated. This approach uses an explicitly specified equation of state (EOS) so that the fluid pressure, temperature, and transport properties are directly related to the flow variables. The method is efficient since no subiterations are required to deduce the pressure and temperature from the flux variables and is modular by allowing different equations of state to be easily supplied to the code. The flexibility of the EOS approach is shown by its implementation into a high-order total variation diminishing upwinding scheme as well as a standard central-differencing scheme. The EOS approach is then demonstrated by computing the hypersonic flow through the corner region of two mutually perpendicular flat plates using both the upwind and central-differencing schemes.
This presentation discusses recent work being conducted by the National Aeronautics and Space Administration (NASA) at Marshall Space Flight Center (MSFC) to evaluate various methods that could be employed to provide for planetary protection of those solar system bodies that are candidates for extraterrestrial life, thus preventing contamination of such bodies. MSFC is presently involved in the development phase of the Europa Lander De-Orbital Stage (DOS) braking motor. In order to prevent bio-contamination of this Jovian satellite, three paths are currently being considered. The first is (1) Bio-Reduction of those microscopic organisms in or on the vehicle (in this case a solid rocket motor (SRM)) that might otherwise be transported during the mission. Possible methods being investigated include heat sterilization, application or incorporation of biocide materials, and irradiation. While each method can be made to work, effects on the SRM's components (propellant, liner, insulation, etc.) could well prove deleterious. A second path would be use of (2) Bio-Barrier material(s). So long as such barrier(s) can maintain their integrity, planetary protection should be afforded. Under the harsh conditions encountered during extended spaceflight (vacuum, temperature extremes, radiation), however, such barrier(s) could well experience a breach. Finally, a third path would be to perform (3) Pyrotechnic Sterilization of the SRM during its end-of-mission phase. Multiple pyrotechnic units would be triggered to ensure activation of such an event and provide for a final sterilization before vehicle impact. In light of Europa's stringent bio-reduction targets, the final and best choice to minimize risk will probably be some combination of the above.
The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)
A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.
In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.
Because of its simplicity and the ability to produce a stable, slow-propagating crack, the Double-Torsion (DT) method has been used widely for investigating the critical and subcritical propagation of a slow-propagating tensile (mode-I) crack. However, to determine the complex relationship between the crack velocity $\mathcal{v_c}$ vs. the strain energy release rate $\mathscr{G}$ (or the stress intensity factor K) from laboratory measurements, several corrections must be made to account for the impact of sample and crack geometry. Particularly, DT test typically produces a crack with a curved edge profile instead of a straight line, causing the local $\mathcal{v_c}$ and $\mathscr{G}$ vary along the crack front. The experimentally measured $\mathcal{v_c}$ and $\mathscr{G}$ data merely reflect collective, averaged behavior of the crack. This makes inversion for the intrinsic, “true” crack growth kinetics necessary, based upon the knowledge of the crack geometry. Simple and effective correction methods have been proposed and validated for the slow, chemical-reaction-controlled part (Region I) of the $\mathcal{v_c}$-$\mathscr{G}$ curve. However, reliable methods for the highly nonlinear, transport-dominated part (Region II) and its sudden transition to the dynamic propagation part (Region III) are still lacking. Here we propose a method for determining the intrinsic $\mathcal{v_c}$-$\mathscr{G}$ relationship across all three Regions based upon DT test data, using a simple model function and its numerical inversion. The performance of this approach is examined and demonstrated using both synthetic and laboratory data for subcritical crack growth in soda lime glass.
Volume I of this report presents a new method for synthesizing hybrid systems directly from design requirements, and applies the method to design of a hybrid system for longitudinal control of transport aircraft. The resulting system satisfies general requirement for safety and effectiveness specified a priori, enabling formal validation to be achieved. Volume II contains seven appendices intended to make the report accessible to readers with backgrounds in human factors, fli ght dynamics and control. and formal logic. Major design goals are (1) system desi g n integrity based on proof of correctness at the design level, (2), significant simplification and cost reduction in system development and certification, and (3) improved operational efficiency, with significant alleviation of human-factors problems encountered by pilots in current transport aircraft. This report provides for the first time a firm technical basis for criteria governing design and certification of avionic systems for transport aircraft. It should be of primary interest to designers of next-generation avionic systems.
Volume I of this report presents a new method for synthesizing hybrid systems directly from desi gn requirements, and applies the method to design of a hybrid system for longitudinal control of transport aircraft. The resulting system satisfies general requirement for safety and effectiveness specified a priori, enabling formal validation to be achieved. Volume II contains seven appendices intended to make the report accessible to readers with backgrounds in human factors, flight dynamics and control, and formal logic. Major design goals are (1) system design integrity based on proof of correctness at the design level, (2) significant simplification and cost reduction in system development and certification, and (3) improved operational efficiency, with significant alleviation of human-factors problems encountered by pilots in current transport aircraft. This report provides for the first time a firm technical basis for criteria governing design and certification of avionic systems for transport aircraft. It should be of primary interest to designers of next-generation avionic systems.
The Parker Solar Probe and the incoming Solar Orbiter mission will provide measurements of solar energetic particle (SEP) events at close heliocentric distances from the Sun. Up to present, the largest data set of SEP events in the inner heliosphere are the observations by the two Helios spacecraft. Aims.We re-visit a sample of 15 solar relativistic electron events measured by the Helios mission with the goal of better characterizing the injection histories of solar energetic particles and their interplanetary transport conditions at heliocentric distances <1AU. Methods. The measurements provided by the E6 instrument on board Helios provide us with the electron directional distributionsin eight di erent sectors that we use to infer the detailed evolution of the electron pitch-angle distributions. The results of a Monte Carlo interplanetary transport model, combined with a full inversion procedure, were used to fit the observed directional intensities in the 300–800 keV nominal energy channel. Unlike previous studies, we have considered both the energy and angular responses of the detector. This method allowed us to infer the electron release time profile at the source and determine the electron interplanetary transport conditions. Results. We discuss the duration of the release time profiles and the values of the radial mean free path, and compare them with thevalues reported previously in the literature using earlier approaches. Five of the events show short injection histories (<30 min) at the Sun and ten events show long-lasting (>30 min) injections. The values of mean free path range from 0.02AU to 0.27AU.Conclusions. The inferred injection histories match with the radio and soft X-ray emissions found in literature.We find no dependence of the radial mean free path on the radial distance. In addition, we find no apparent relation between the strength of interplanetary scattering and the size of the solar particle release.
Beyond its fascinating chemistry as the first element in the Periodic Table, hydrogen is of high societal importance in energy technologies and of growing importance in energy-efficient computing. In energy, hydrogen has reemerged as a potential solution to long-term energy storage and as a carbon-free input for materials manufacturing. Its utilization and production rely on the availability of proton-conducting electrolytes and mixed proton–electron conductors for the components in fuel cells and electrolyzers. In computing, proton mediation of electronic properties has garnered attention for electrochemically controlled energy-efficient neuromorphic computing. Incorporation of substitutional and interstitial hydride ions in oxides, though only recently established, enables tuning of electronic and magnetic properties, inviting a range of possible exotic applications. This article addresses common themes in the fundamental science of hydrogen incorporation and transport in oxides as relevant to pressing technological needs. The content covers (1) lattice (or bulk) mechanisms of hydrogen transport, primarily addressing proton transport, but also touching on hydride ion transport; (2) interfacial transport; (3) exploitation of extreme external drivers to achieve unusual response; and (4) advances in methods to probe the hydrogen environment and transport pathway. The snapshot of research activities in the field of hydrogen-laden materials described here underscores exciting recent breakthroughs, remaining open questions, and breathtaking experimental tools now available for unveiling the nature of hydrogen in solid-state matter.