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.
The Hilbert-Smith conjecture claims that if a compact group G acts freely on a manifold, then it is a Lie group. For a finite-dimensional orbit space a reduction of the Hilbert-Smith conjecture to certain other problems in geometric topology is presented; in these the key problem is the existence of an essential sequence of lens spaces of increasing dimension. Bibliography: 52 titles. (paper)
This study investigates the extrusion processes of deuterium and protium using ANSYS-Polyflow. The geometries and computational fluid dynamics (CFD) settings closely replicate the experimental setups and data acquired from the extruder experiments at Oak Ridge National Laboratory (ORNL) for validation purposes. We explore the impacts of (1) slip versus non-slip boundary conditions and (2) the use of constant, temperature-, and shear rate–dependent viscosities, concluding that the implementation of non-slip wall boundary conditions combined with shear rate–dependent viscosity produced more accurate predictions. The simulations achieved excellent agreement with the experimental data, with relative differences of only 5% for deuterium, and 3% to 6% for protium. This is the first time that experimental extrusion data at ORNL have been accurately predicted through high-fidelity CFD modeling. In conclusion, the advancements offer valuable insights and a foundational modeling tool for optimizing pellet injectors for ITER and other future reactor-scale devices.
Not Available
Morse-Smale (MS) complexes have been gaining popularity as a tool for feature-driven data analysis and visualization. However, the quality of their geometric embedding and the sole dependence on the input scalar field data can limit their applicability when expressing application-dependent features. In this paper we introduce a new combinatorial technique to compute an MS complex that conforms to both an input scalar field and an additional, prior segmentation of the domain. The segmentation constrains the MS complex computation guaranteeing that boundaries in the segmentation are captured as separatrices of the MS complex. We demonstrate the utility and versatility of our approach with two applications. First, we use streamline integration to determine numerically computed basins/mountains and use the resulting segmentation as an input to our algorithm. This strategy enables the incorporation of prior flow path knowledge, effectively resulting in an MS complex that is as geometrically accurate as the employed numerical integration. Our second use case is motivated by the observation that often the data itself does not explicitly contain features known to be present by a domain expert. We introduce edit operations for MS complexes so that a user can directly modify their features while maintaining all the advantages of a robust topology-based representation.
We present an algorithm to compute an approximation of the generalized Voronoi diagram (GVD) on arbitrary collections of 2D or 3D geometric objects. In particular, we focus on datasets with closely spaced objects; GVD approximation is expensive and sometimes intractable on these datasets using previous algorithms. With our approach, the GVD can be computed using commodity hardware even on datasets with many, extremely tightly packed objects. Our approach is to subdivide the space with an octree that is represented with an adjacency structure. We then use a novel adaptive distance transform to compute the distance function on octree vertices. The computed distance field is sampled more densely in areas of close object spacing, enabling robust and parallelizable GVD surface generation. We demonstrate our method on a variety of data and show example applications of the GVD in 2D and 3D.
Boiler tubes in steam power plants experience tube blockages due to exfoliation of oxide grown on the inner side of the tubes. In extreme cases, significant tube blockages can lead to forced power plant outages. It is thus desired to predict through modeling the amount of tube blockage in order to inform power plant operators of possible forced outages. SpalLoop solves for the stress-strain equations in an axisymmetric geometry, tracking the stress/strain evolution during boiler operation including outages for the entire boiler tube length. At each operational outage, i.e., temperature excursions down to room temperature, the amount of exfoliated area for the entire tube loop is estimated the amount of tube blockage is predicted based assumed blockage geometry and site. The SpaLLoop code contains modules developed for oxide growth, stress analysis, tube loop geometry, blockage area by taking into account the following phenomena and features, (a) Plant operation schedule with periodic alternate full-load and partial-load regimes and shut-downs, i.e., temperature excursions from high-load to room temperature, (b) axisymmetric formulation for cylindrical tubes, (c) oxide growth in a temperature gradient with multiple oxide layers, (d) geometry of a boiler tube with a single tube loop or two tube loops, (e) temperature variation along the tube length based on hot gas temperature distribution outside the tube and inlet steam temperature, (f) non-uniform oxide growth along the tube length according to the local steam tube temperature, (g) exfoliated area module: at each operational outage considered, the amount of exfoliated area and exfoliated volume along the tube is estimated, (h) blockage module: at each operational outage considered, the exfoliated volume/mass for each tube loop is estimated from which the amount of tube blockage is predicted based on given blockage geometry (length, location, and geometry). The computer program is written in FORTRAN90. Its modular structure was sought for allowing the best flexibility in updating the program by implementing new constitutive equations due to availability of new material property data and/or new physical phenomena.
Adsorbents developed for the recovery of uranium from seawater display poor selectivity over other transition metals present in the ocean, with vanadium particularly problematic. To improve selectivity, an indispensable step is the positive identification of metal binding environments following actual seawater deployment. In this work we apply x-ray absorption fine structure (XAFS) spectroscopy to directly investigate the vanadium binding environment on seawater-deployed polyamidoxime adsorbents. Comparison of the x-ray absorption near edge spectra (XANES) reveal marked similarities to recently a reported non-oxido vanadium (V) structure formed upon binding with cyclic imidedioxime, a byproduct of generating amidoxime functionalities. Density functional theory (DFT) calculations provided a series of putative vanadium binding environments for both vanadium (IV) and vanadium (V) oxidation states, and with both amidoxime and cyclic imidedioxime. Fits of the extended XAFS (EXAFS) data confirmed vanadium (V) is bound exclusively by the cyclic imidedioxime moiety in a 1:2 metal:ligand fashion, though a modest structural distortion is also observed compared to crystal structure data and computationally optimized geometries which is attributed to morphology effects from the polymer graft chain and the absence of crystal packing interactions. These results demonstrate that improved selectivity for uranium over vanadium can be achieved by suppressing the formation of cyclic imidedioxime during preparation of polyamidoxime adsorbents for seawater uranium recovery.
A geometric cell of the Monte Carlo N-Particle (MCNP)1 transport code is traditionally created by using Boolean operators on defined surfaces. This constructive solid geometry (CSG) capability has been available in the MCNP code since its beginning. However, a CSG model approach is limited when it comes to constructing a representative geometry for a complex model in its ability to capture a correct model representation. Starting with the version 6.0, the MCNP code has the ability of embedding an unstructured mesh (UM) model into a CSG cell to create a hybrid geometry [1]. The MCNP UM feature provides the flexibility of defining very complex geometries because computer aided design (CAD) and mesh generation software packages can be utilized to construct UM models for MCNP simulations.
The small length scales of the dissipative processes of physical viscosity and heat conduction are typically not resolved in the numerical simulation of high Reynolds number flows in the discrete geometry of computational grids. Historically, the simulations of flows with shocks and/or turbulence have relied on solving the Euler equations with dissipative regularization. In this paper, we begin by reviewing the regularization strategies used in shock wave calculations in both a Lagrangian and an Eulerian framework. We exhibit the essential similarities with Large Eddy Simulation models of turbulence, namely that almost all of these depend on the square of the size of the computational cell. In our principal result, we justify that dependence by deriving the evolution equations for a finite-sized volume of fluid. Those evolution equations, termed finite scale Navier-Stokes (FSNS), contain dissipative terms similar to the artificial viscosity first proposed by von Neumann and Richtmyer. Here, we describe the properties of FSNS, provide a physical interpretation of the dissipative terms and show the connection to recent concepts in fluid dynamics, including inviscid dissipation and bi-velocity hydrodynamics.
Here, we present a comparison of the numerics and simulation results for two "full" Stokes ice sheet models, FELIX-S (Leng et al. 2012) and Elmer/Ice. The models are applied to the Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d). For the diagnostic experiment (P75D) the two models give similar results (< 2 % difference with respect to along-flow velocities) when using identical geometries and computational meshes, which we interpret as an indication of inherent consistencies and similarities between the two models. For the standard (Stnd), P75S, and P75R prognostic experiments, we find that FELIX-S (Elmer/Ice) grounding lines are relatively more retreated (advanced), results that are consistent with minor differences observed in the diagnostic experiment results and that we show to be due to different choices in the implementation of basal boundary conditions in the two models. While we are not able to argue for the relative favorability of either implementation, we do show that these differences decrease with increasing horizontal (i.e., both along- and across-flow) grid resolution and that grounding-line positions for FELIX-S and Elmer/Ice converge to within the estimated truncation error for Elmer/Ice. Stokes model solutions are often treated as an accuracy metric in model intercomparison experiments, but computational cost may not always allow for the use of model resolution within the regime of asymptotic convergence. In this case, we propose that an alternative estimate for the uncertainty in the grounding-line position is the span of grounding-line positions predicted by multiple Stokes models.
FASTER computer program modified for tracking neutrons through reactor slabs penetrated by helical voided ducts
Recent trends in turbine aerodynamics are summarized. Areas discussed include cooled turbine aerodynamics, high work factor turbines, pneumatic variable geometry, and computer analyses.
Brief description of some of the recent trends in turbine aerodynamics research. Areas covered include cooled turbine aerodynamics, high work-factor turbines, pneumatic variable geometry, and computer analysis. Most of this research is shown to have been directed primarily toward allowing the turbine to accept efficiently higher inlet temperatures and blade loading.
Studies to evaluate advanced technology relative to chemical upper stages and orbit-to-orbit stages are reported. The work described includes: development of LH2/LOX stage data, development of data to indicate stage sensitivity to engine tolerance, modified thermal routines to accommodate storable propellants, added stage geometries to computer program for monopropellant configurations, determination of the relative gain obtainable through improvement of stage mass fraction, future propulsion concepts, effect of ultrahigh chamber-pressure increases, and relative gains obtainable through improved mass fraction.
A digital computer program for evaluating the weight and costs of advanced transport designs was developed. The resultant program, intended for use at the preliminary design level, incorporates both batch mode and interactive graphics run capability. The basis of the weight and cost estimation method developed is a unique way of predicting the physical design of each detail part of a vehicle structure at a time when only configuration concept drawings are available. In addition, the technique relies on methods to predict the precise manufacturing processes and the associated material required to produce each detail part. Weight data are generated in four areas of the program. Overall vehicle system weights are derived on a statistical basis as part of the vehicle sizing process. Theoretical weights, actual weights, and the weight of the raw material to be purchased are derived as part of the structural synthesis and part definition processes based on the computed part geometry.
Transonic flow fields about wing geometries are computed using an Euler/Navier-Stokes approach in which the flow field is divided into several zones. The grid zones immediately adjacent to the wing surface are suitably clustered and solved with the Navier-Stokes equations. Grid zones removed from the wing are less finely clustered and are solved with the Euler equations. Wind tunnel wall effects are easily and accurately modeled with the new grid-zoning algorithm because the wind tunnel grid is constructed as an exact subset of the corresponding free-air grid. Solutions are obtained that are in good agreement with experiment, including cases with significant wind tunnel wall effects and shock-induced separation on the upper wing surface.
The geometry and computer programming for mathematical models of a one-dimensional conventional photoconverter, a one-dimensional vertical junction photoconverter, a three-dimensional conventinal photoconverter, and a three-dimensional vertical junction solar cell are discussed.