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.
Detailed descriptions are presented of the input data for each of the MSFC computer codes applied to the analysis of a realistic nuclear propelled vehicle. The analytical techniques employed include cross section data, preparation, one and two dimensional discrete ordinates transport, point kernel, and single scatter methods.
Structural shape optimization couples with other discipline optimization in the design of complex engineering systems. For instance, the wing structural weight and elastic deformations couple to aerodynamic loads and aircraft performance through drag. This coupling makes structural shape optimization a subtask in the overall vehicle synthesis. Decomposition methods for optimization and sensitivity analysis allow the specialized disciplinary methods to be used while the disciplines are temporarily decoupled, after which the interdisciplinary couplings are restored at the system level. Application of decomposition methods to structures-aerodynamics coupling in aircraft is outlined and illustrated with a numerical example of a transport aircraft. It is concluded that these methods may integrate structural and aerodynamic shape optimizations with the unified objective of the maximum aircraft performance.
A new solution method has recently been developed to analyze the thermal and structural response of polymeric composites during chemical decomposition (Sullivan and Salamon, 1992a). This method couples the equations of energy, mass transport and momentum and solves these equations simultaneously at each time step. In this paper, the couples solution method is applied to analyze the response of a glass phenolic composite during its chemical decomposition. The numerical solution to this particular problem was originally performed by Henderson, et al. (1985) and the actual response was measured by Henderson and Hagen (1985) and Ramamurthy (1988). In addition to further verification of the coupled solution method, the present solution offers additional insight into certain physical events which occur during chemical decomposition.
The use of electrothermal vaporization (ETV) has not been readily applied to nuclear safeguards or pre inspection check (PIC) sample analysis, even though the technique is ideally suited for low level (pg) isotopic quantification due to its high analyte transport efficiency. Here demonstrate the ETV method for PIC safeguard analysis. The ETV system is an automated, commercial of the shelf (COTS) technology, with the capability to introduce solid samples to a mass spectrometer platform. The system has been shown to simultaneously determine uranium assay, and the major isotopic ratios for uranium and plutonium. The ETV unit has demonstrated analysis up to 50 samples a day autonomously. The utilization of this ETV methodology provide a secondary option to the current PIC analysis and offers high throughput, enhanced sensitivity, and broader isotopic detection platform.
The primary objective of this effort was to develop a method for suppressing convective flows during the growth of mercurous chloride crystals by vapor transport in closed tubes to levels approaching those obtained in the microgravity environment. Mercurous chloride was chosen because it is a technologically interesting acoustical optical material whose optical properties are believed to be affected by convective flows. Since the Grashof number scales as the cube of the smallest dimension in the flow system, reduction of the size scale can be extremely effective in reducing unwanted convective flows. However, since materials of practical interest must be grown at least on the cm scale, reduction of the overall growth system is not feasible. But if the region just above the growing crystal could be restricted to a few mm, considerable reduction in flow velocity would result. By suspending an effusive barrier in the growth ampoule just above the growth interface, it should be possible to reduce the convective velocity in this vicinity to levels approaching flows in microgravity. If successful, this growth technique will offer a screening test for proposed space experiments that involve vapor transport to see if reduction of convection will result in improved material and will set a new standard against which the improvements obtained in microgravity may be judged. In addition, it may provide an improved method for preparing materials on Earth whose growth is affected adversely by convection. If the properties of this material can be improved there is a potential commercial interest from Brimrose Inc., who has agreed to fabricate and test devices from the crystals we have grown. This report describes the development of the growth facility, the purification processes developed for preparing the starting material, and the results from growth experiments with and without the effusive baffle. Mercurous chloride turned out to be a more difficult material to deal with than originally anticipated. At growth temperatures, it is extremely sensitive to practically any impurity which causes it to form oxychlorides and/or to decompose into elemental mercury and bichloride of mercury. We were unable to find a suitable method for protecting the magnetic material used to suspend the effusion barrier from the attack of mercurous chloride vapor. Although we were successful in growing single crystals of mercurous chloride without the effusion baffle, they exhibited severe microcracking which we attribute to wall-induced thermal stresses. This leads us to believe that uncontrolled convection may not be the most important problem in the development of this material and a new growth process was attempted that eliminates the wall-induced stress. Unfortunately, the grant ran out before this new method could be adequately tested.
The Asian Marginal Seas are interconnected by a number of narrow straits, such as the Makassar Strait connecting the Pacific Ocean with the Indian Ocean, the Luzon Strait connecting the South China Sea with the Pacific Ocean, and the Korea/Tsushima Strait connecting the East China Sea with the Japan/East Sea. Here we propose a method, the combination of the "geostrophic control" formula of Garrett and Toulany (1982) and the "hydraulic control" theory of Whitehead et al. (1974), allowing the use of satellite-observed sea-surface-height (SSH) and ocean-bottom-pressure (OBP) data for estimating interbasin transport. The new method also allows separating the interbasin transport into surface and bottom fluxes that play an important role in maintaining the mass balance of the regional oceans. Comparison with model results demonstrates that the combined method can estimate the seasonal variability of the strait transports and is significantly better than the method of using SSH or OBP alone.
MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.
Over the solar-activity cycle, there are extended periods where significant discrepancies occur between the spacecraft-observed total (unsigned) open magnetic flux and that determined from coronal models. In this article, the total open heliospheric magnetic flux is computed using two different methods and then compared with results obtained from insitu interplanetary magnetic-field observations. The first method uses two different types of photospheric magnetic-field maps as input to the Wang–Sheeley–Arge (WSA) model: i) traditional Carrington or diachronic maps, and ii) Air Force Data Assimilative Photospheric Flux Transport model synchronic maps. The second method uses observationally derived helium and extreme-ultraviolet coronal-hole maps overlaid on the same magnetic-field maps in order to compute total open magnetic flux. The diachronic and synchronic maps are both constructed using magnetograms from the same source, namely the National Solar Observatory Kitt Peak Vacuum Telescope and Vector Spectromagnetograph. The results of this work show that the total open flux obtained from observationally derived coronal holes agrees remarkably well with that derived from WSA, especially near solar minimum. This suggests that, on average, coronal models capture well the observed large-scale coronal-hole structure over most of the solar cycle. Both methods show considerable deviations from total open flux deduced from spacecraft data, especially near solar maximum, pointing to something other than poorly determined coronal-hole area specification as the source of these discrepancies.
The back end of the fuel cycle focuses on the interim storage, transportation, and final disposition of the spent nuclear fuel from nuclear reactors. Commercial light water nuclear power stations across the United States operate with fuel assemblies being irradiated for up to 6 years (planned) in the reactor pressure vessel. After their planned irradiation, the fuel assemblies are moved to a spent fuel pool within the facility complex. After the spent nuclear fuel is removed from the fuel pool, it is typically inserted into a welded metal canister that can be transferred between overpacks for storage, transportation and possibly disposal. Most dry storage systems being used by industry today use dual purpose canisters (DPCs), designed for use in storage and transportation overpacks, but not specifically designed for disposal. Triple purpose-canisters, designed for disposal in addition to storage and transportation, have also been researched. Traditional manufacturing methods for spent fuel canisters involve fusion welding along the length or circumference of the canister which results in high tensile residual stresses in the joint weld zone (WZ) and heat affected zone (HAZ). In this paper, spent fuel canister designs were printed by wire arc additive manufacturing (AM) using the 316L SS welding wire to demonstrate: 1. Feasibility of spent nuclear fuel canister fabrication using this advanced manufacturing method, 2. Dynamic response of the additive manufacture canister design when subjected to the federally mandated Normal Conditions of Transport (NCT) and Hypothetical Accident Conditions (HAC) physical tests for Type B packages. This paper will focus on the canister printing design and structural tests. The AM 3D printed design, regulatory testing, and post-test evaluation of the canister tested to the 10 CFR 71.71 and 10 CFR 71.73 requirements will be presented. One AM canister design was subjected to the penetration, free drop, and puncture test. Before and after the dynamic structural tests, the AM canister design was scanned with a handheld scanner to capture a 3D CAD geometry to compare to the 3D printed canister design in the deformed shape. The scanned geometry was sectioned in areas with deformation and the cross-section profile was measured to determine accurate and repeated results of the deformed shape of the AM canister design.
The overarching goal of this project, which unfolded over a period of seventeen years, including two no-cost extensions in the last two years, was to extend and deepen our ability to use time-dependent density functional theory and related theoretical methods for the description of the transport of charge, spin, or energy in quantum electronic systems at the nanoscale. The output of our efforts is contained in 52 papers, listed below, which we have published since 2006. There have also been many presentations at international conferences, workshops, colloquia and seminars. This body of work has been carried out with the essential input of 8 postdoctoral students (who were supported by the grant) and several external senior collaborators and their postdocs.
In virtue of adding one elliptic equation that must be implicitly solved, the 'barely implicit correction' method presented removes the stringent sidestep limit imposed by sound speed in the explicit methods. The barely implicit method is presently combined with a flux-corrected transport algorithm, which has been rendered fourth-order by the appropriate subtraction of corrected fluxes, in order to accurately represent the sharp gradients in subsonic flows. The multigrid MGRID technique (DeVore, 1984) is used to efficiently solve the elliptic pressure equation.
A numerical solution of the turbulent mass transport equation utilizing the concept of eddy diffusivity is presented as an efficient method of investigating turbulent mass transport in boundary layer type flows. A FORTRAN computer program is used to study the two-dimensional diffusion of ammonia, from a line source on the surface, into a turbulent boundary layer over a flat plate. The results of the numerical solution are compared with experimental data to verify the results of the solution. Several other solutions to diffusion problems are presented to illustrate the versatility of the computer program and to provide some insight into the problem of mass diffusion as a whole.
In this study, we developed a new method for in situ Sb doping of CdTe thin films combining vapor transport deposition with a Group V pyrolyzer to address Sb doping concentration and doping efficiency. The Sb doped CdSeTe (CdSeTe:Sb) films were deposited in solar cell structures under variations of Sb dopant source heater, vapor pyrolyzer temperature, and Cd vapor excess. Results indicate that although these parameters do not affect the CdTe morphology or crystal structure, they critically influence doping efficiency and trap concentration. Capacitance-voltage measurements show that a higher dopant heater (TD) or pyrolyzer (TP) temperature leads to higher net carrier concentration, achieving a net carrier concentration of 1016 cm-3 and 20% doping efficiency with a TD/TP combination of 600 degrees C/1100 degrees C. By tuning the Cd/Sb flux ratio during CdSeTe:Sb deposition, the lowest defect concentration is achieved at Cd/Sb of 1.4:1, which produced the best VOC CdSeTe:Sb cell. This demonstrates a path to produce high net carrier concentration polycrystalline CdTe thin film with a low concentration of dopant-induced defects.
A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.
There is an opportunity to save energy and reduce operational expenses when choosing a suitable desalination method aided by computational modeling. Existing models are not conducive to generalized comparisons between different desalination methods. Therefore, the work in this study developed metamodels for six desalination methods, grouped them into thermal and molecular transport families, and validated their predictive performance within 9% difference from published data. This validated framework allowed comparisons of desalination methods at their prescribed ranges of operational conditions that they were designed for. These conditions specify feed salinity ranges of 1.6 to 2.4 g/kg for Capacitive Deionization and Reverse Osmosis (RO), 2.8 to 4.2 g/kg for Electrodialysis, 28 to 42 g/kg for Thermovapor Compression and Humidification-Dehumidification, and 37 to 55 g/kg for Multi-Effect Distillation (MED). Despite different operational conditions, all models exhibit non-linear, positive correlation between energy consumption and system size in response to feed salinity and production rate. The framework is also employed in a cross-comparative analysis between MED and RO whose results suggest that energy intensity for MED is an order of magnitude greater than RO for the same operational conditions, but actual operational costs are comparable. Overall, the framework is ready for deployment in case studies of actual desalination plants.
A modified Bhatnagar-Gross-Krook (BGK) collision method is used for the derivation of reactive collision terms to allow for complete higher-moment approximations in fluid transport theory. The general reactive collision terms for the 8-, 10- and 13-moment approximation for binary reactions are derived. Then, reactive collision terms for specific chemistry are derived with the assumptions: (1) all reactants and products have 13-moment distribution functions; (2) simple, classical assumptions provide sufficient models for reaction dynamics; and (3) the reaction rate is energy-independent. The derived specific formulas reflect two extremes in reactive collision dynamics: direct and indirect reaction mechanisms. Finally, considerations for correct use of the reactive collision terms are discussed in the context of space plasma environments.
The primary mode of spent nuclear fuel transportation within the United States will be by railcar. One such system is the Atlas railcar, which is designed to transport 17 different spent nuclear fuel cask systems, including bare fuel systems and canister fuel systems. In the latter configuration, multipurpose canisters containing spent nuclear fuel may be placed within an overpack for storage, or within a cask for transportation. Compared to bare fuel systems, canister fuel systems have additional degrees of freedom for motion during transportation, because clearance between the cask and canister allows for some motion of the canister to occur relative to the cask. This work investigates the effect of canister motion on the shock and vibration imparted to the spent nuclear fuel within. Structural dynamic analyses have been conducted to identify the effects of canister to cask clearance, presence and type of dunnage, and loading direction and frequency. This modeling study calculates anticipated cask motion, canister motion, and spent nuclear fuel structural dynamic response to normal conditions of transportation railcar motion using finite element analysis methods that were developed to model the rail segment of the ENSA/DOE (Equipos Nucleares S.A., U.S. Department of Energy) multimodal transportation test of 2017.
A variable-complexity strategy of combining simple and detailed analysis methods is presented for the design optimization of a high-speed civil transport (HSCT) wing. Two sets of results are shown: the aerodynamic design of the wing using algebraic weight equations for structural considerations, and optimization results of the internal wing structure for a fixed wing configuration. We show example results indicating that using simple analysis methods alone for the calculation of a critical constraint can allow an optimizer to exploit weaknesses in the analysis. The structural optimization results provide a valuable check for the weight equations used in the aerodynamic design. In addition, these results confirm the need for using simple, algebraic models in conjunction with more detailed analysis methods. A strategy of interlaced aerodynanic-structural design is proposed.