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.
Oral presentation on NETL research delivered at the 2022 Resource Sustainability Meeting. This effort concerns advancing a concept for early kick detection through the development of synthetic (simulation and experiment) multiphase acoustic data corresponding to well characterized kick events. The simulation activities are the focus of the talk.
Disclosed herein are multiphase metal anodes useful in non-aqueous batteries. The anodes include at least one active metal and at least one conductive metal.
A solid multiphase polymer blend material comprising: (i) a polyphenolic substance having a molecular weight of at least 500 g/mol; and (ii) a polyester having a molecular weight of at least 500 g/mol; wherein at least a portion of the polyphenolic substance is covalently bonded directly or through a linking moiety to the polyester. Methods for producing the blend material are also described, e.g., homogeneously melt blending a mixture comprising components (i) and (ii) under conditions resulting in covalent attachment of at least a portion of the polyphenolic substance directly or through a linking moiety to the polyester. Methods for producing objects made of the blend material by melt extrusion are also described.
There is a significant growth in the field of decommissioning of nuclear facilities due the fact that for a lot of them their end of life is nearing. In Taiwan, the government plans to steadily reduce the dependence on nuclear power and moves gradually towards a nuclear-free homeland. Therefore, the development of decommissioning technology becomes the significant issue for the next coming years. The decommissioning of a nuclear power plant is a long process and practical work including radiological surveys, decontamination, dismantling, and the subsequent waste management. Chemical decontamination is an important task in the decommissioning of nuclear power plants and is usually performed by circulating the chemical solution in selected systems or full systems in order to reduce dose levels. However, the decontamination of segmented metal parts is also considered to reduce the contamination to such levels that they may be disposed of at a lower classification or be unconditionally released, and it helps to achieve the goal of waste minimization and provide a more economical cost for disposal management. Chemical decontamination technology is used to remove the contaminated oxides deposited on the surface of metals by dissolving them with chemical mechanisms, such as chelation, acidic or alkaline dissolution, redox reaction and so on. Carbon steel is the major material used in systems and pipelines in the nuclear power plant, thus it is also commonly selected as the target to be decontaminated in the decommissioning process. In this study, the effectiveness of chemical decontamination is evaluated by using multiphase treatment process including the reduction and the oxidation steps. In the step of redox reaction, the oxidation state of the substance is very important. Especially metals of some oxidation states show solubility in water, which has a decisive influence on the chemical decontamination. The experiment of decontamination of metal scraps is carried out by the immersion method within a chemical solution. The relative operation parameters such as temperature, pH, and concentration of the chemical solution are recorded to demonstrate the effect of the decontamination process. (authors)
This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two-phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection-dominated to diffusion-dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.
The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.
In this work we present a comprehensive investigation of state-of-the-art closure models employed to represent interphase momentum, thermal, and work exchange between the gas and particulate phases for Euler–Lagrange (EL) simulations in shock-driven flows. A complete list of closures for the force, torque, heat transfer, and work exchange models is provided. In particular, the present work includes a stochastic closure for the particle-to-particle variation in the quasi-steady force and a deterministic closure for particle-to-particle variation in the added mass force in an EL framework. These variations arise due to the presence of neighboring particles and particle–particle interactions. To investigate the importance of each closure term, we carry out fully three-dimensional simulations for a planar shock propagating over a random bed of inert particles. The primary goal is to evaluate the role of each closure term on the gas dynamic features (such as transmitted and reflected shock locations) and particle curtain features (such as upstream and downstream curtain locations). To this end, thirteen cases are considered, with each case progressively including a closure model with the goal to identify and quantify its contribution to the simulated dynamics. We show that the volume fraction dependence of the mean force models plays an important role in generating wave-like instabilities that lead to concentration bands. In addition, fluctuations in quasi-steady and added mass forces primarily decrease the internal instabilities that tend to enhance local volume fraction variations. Particle rotation is primarily due to inter-particle collisions, is generally weak, and does not play an important role in the translational dynamics for the present configuration. Inter-phase heat transfer has a strong effect on gas phase temperature, slows down the transmitted and reflected shocks, and decreases the width of the curtain. Furthermore, the absence of a work-coupling model fails to conserve the total energy, greatly under-predicts the gas temperature which in turn affects the particle dynamics.
In the present work, a general formulation is proposed to implement the contact angle boundary conditions for the second-order Phase-Field models, which is applicable to N-phase (N ≥ 2) moving contact line problems. To remedy the issue of mass change due to the contact angle boundary condition, a source term or Lagrange multiplier is added to the original second-order Phase-Field models, which is determined by the consistent and conservative volume distribution algorithm so that the summation of the order parameters and the consistency of reduction are not influenced. To physically couple the proposed formulation to the hydrodynamics, especially for large-density-ratio problems, the consistent formulation is employed. The reduction-consistent conservative Allen-Cahn models are chosen as examples to illustrate the application of the proposed formulation. The numerical scheme that preserves the consistency and conservation of the proposed formulation is employed to demonstrate its effectiveness. Results produced by the proposed formulation are in good agreement with the exact and/or asymptotic solutions. The proposed method captures complex dynamics of moving contact line problems having large density ratios.
MPLBM-UT is a specialized lattice-Boltzmann library that makes running single- and two-phase flow simulations in porous media accessible to everyone. We provide a suite of tools to pre-process computational domains for simulation, to set up custom boundary conditions, to run simulations, to post-process simulation outputs, and to visualize simulation results and data. All of these tools are easily accessible to users through the mplbm_utils Python package included in and automatically installed with MPLBM-UT. The high-performance, highly parallel library Palabos is used as the solver backend. MPLBM-UT is easily deployed in a variety of systems, from laptops to supercomputer clusters. MPLBM-UT also features multiple examples and benchmark templates that allow for fast prototyping of different porous media problems. We also provide an interface for reading in different file types and downloading domains from the Digital Rocks Portal to perform simulations.
The age of the nearest particle pair is introduced as the difference between the current time and the most recent time when the nearest particle pair was formed. The evolution equation for the age-included nearest pair distribution function is derived. With the assumption of random destruction of the nearest particle pairs, the evolution equation predicts the exponential probability distribution of the ages of the nearest particle pairs. Particle-resolved numerical simulations with moving particles are performed to verify this prediction. The equation is then used to derive the evolution equation for the particle–fluid–particle (PFP) stress, which is known to be related to hyperbolicity of the two-fluid equations. It is found that the relaxation time of the age probability distribution is also the relaxation time for the PFP stress. Guided by the closure terms in the PFP stress evolution equation, we study kinematics of the nearest particle pairs in the particle-resolved simulations for flows caused by sedimentation of the particles with initially isotropic and homogeneous particle distributions. At the steady states, the particle Reynolds numbers are around 20. Anisotropy and inhomogeneity of particle distributions are seen to develop in these flows. The mean distances to the nearest particles and evolution of the distribution of the Voronoi cell volumes are studied. We also found the PFP stress is closely related to the changes in these inter-particle scale quantities.
The Rothman-Keller color-gradient (CG) lattice Boltzmann method is a popular method to simulate two-phase flow because of its ability to deal with fluids with large viscosity contrasts and a wide range of interfacial tensions. Here, two fluids are labeled red and blue, and the gradient in the color difference is used to compute the effect of interfacial tension. It is well known that finite-difference errors in the color-gradient calculation lead to anisotropy of interfacial tension and errors such as spurious currents. Here, we investigate the accuracy of the CG calculation for interfaces between fluids with several radii of curvature and find that the standard CG calculations lead to significant inaccuracy. Specifically, we observe significant anisotropy of the color gradient of order 7% for high curvature of an interface such as when a pinchout occurs. We derive a second order accurate color gradient and find that the diagonal nearest neighbors can be weighted differently than in the usual color-gradient calculation such that anisotropy is minimized to a fraction of a percent. The optimal weights that minimize anisotropy for the smallest radius of curvature interface are found to be w = (0.298, 0.284, 0.275) for diagonal nearest neighbors for the cases of the interface smoothing parameter β = (0.5, 0.7, 0.99), somewhat higher than the w = 0.25 value derived by Leclaire et al. [Leclaire, Reggio, and Trepanier, Computers and Fluids 48, 98 (2011)] based on obtaining isotropic errors to second order. We find that use of these optimal w values yields over a factor of 10 decrease in anisotropy and over a factor of 30 decrease in mean anisotropy relative to using the standard w = 1 value. And we find a factor of about 2 decrease in the anisotropic error and up to factor 15 decrease in mean anisotropic error relative to the choice of w = 0.25 for small radius of curvature interfaces. The improved CG calculations will allow the method to be more reliably applied to studies of phenomenology and pore scale processes such as viscous and capillary fingering, and droplet formation where surface-tension isotropy of narrow fingers and small droplets plays a crucial role in correctly capturing phenomenology. We present an example illustrating how different phenomena can be captured using the improved color-gradient method. Namely, we present simulations of a wetting fluid invading a fluid filled pipe where the viscosity ratio of fluids is unity in which droplets form at the transition to fingering using the improved CG calculations that are not captured using the standard CG calculations. We present an explanation of why this is so which relates to anisotropy of the surface tension, which inhibits the pinchouts needed to form droplets.
A general model is presented for mathematically describing the physics and chemistry of biomass pyrolysis in a fluidized bed reactor.
This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.
Abstract Mechanical properties of traditional engineering materials are typically coupled to each other, presenting a challenge to practitioners with multi‐dimensional material property requirements. In this work, continuous, independent control over multiple mechanical properties is demonstrated in composite materials realized using additive manufacturing. For the first time, composites additively manufactured from rigid plastic, soft elastomer, and liquid constituents are experimentally characterized, demonstrating materials which span four orders of magnitude in modulus and two orders of magnitude in toughness. By forming analytical mappings between relative concentrations of constituents at the microscale and resulting macroscale material properties, inverse material design is enabled; the method is showcased by printing artifacts with prescribed toughness and elasticity distributions. The properties of these composites are placed in the context of biological tissues, showing they have promise as mechanically plausible tissue mimics.
Not provided.