Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Methods in transport”

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.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods

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.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

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.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Inversion-based correction of Double-Torsion (DT) subcritical crack growth tests for crack profile geometry

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.

42 ENGINEERING↗
DOE OSTI2024

Hydrogen in energy and information sciences

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.

08 HYDROGEN↗
DOE OSTI2024

A mineral precipitation model based on the volume of fluid method

A novel volume of fluid method is presented for mineral precipitation coupled with fluid flow and reactive transport. The approach describes the fluid-solid interface as a smooth transitional region, which is designed to provide the same precipitation rate and viscous drag force as a sharp interface. Specifically, the governing equation of mineral precipitation is discretized by an upwind scheme, and a rigorous effective viscosity model is derived around the interface. The model is validated against analytical solutions for mineral precipitation in channel and ring-shaped structures. It also compares well with interface tracking simulations of advection-diffusion-reaction problems. Here, the methodology is finally employed to model mineral precipitation in fracture networks, which is challenging due to the low porosity and complex geometry. Compared to other approaches, the proposed model has a concise algorithm and contains no free parameters. In the modeling, only the pore space requires meshing, which improves the computational effciency especially for low-porosity media.

58 GEOSCIENCES↗
DOE OSTI2020

Accelerating the solution of the S N equations with highly anisotropic scattering using the Fokker-Planck approximation

The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accelerating such problems because these methods do not correct higher order Legendre-moments of angular flux. In this paper, we explore the idea of using the Fokker-Planck approximation as a preconditioner to accelerate the convergence of solution of forward-peaked transport problems. Finally, we observe that such preconditioning can successfully reduce the iteration count by up to four orders of magnitude, and solver runtime by up to three orders of magnitude when compared to diffusion synthetic acceleration.

42 ENGINEERING↗
DOE OSTI2020

The origin, transport, and evolution of ash in engine particulate filters

Engine particulate filters have been widely applied across the world to control engine exhaust particulate matter (or particulate number) emissions. With increasing vehicle mileage, ash accumulation deteriorates vehicle fuel economy and complicates on-board control. Extending filter service life with ash loading has significant economic and environmental impacts. Many studies have been conducted in characterizing ash accumulation and evaluating its impacts on filter performance. However, comprehensive reviews covering all the key issues in the field are rather rare. This paper reviews the extensive prior research on filter ash, and not only summarizes the experimental observations but also elucidates the fundamental mechanisms. The review covers the areas of ash origin, accumulation, transport, evolution, and artificial acceleration methods. The previously reported data of ash properties is compiled and analyzed. Furthermore, the advantages and disadvantages of ash acceleration approaches are also discussed in detail. Based on the cumulative understanding, a few potential ways to improve ash management are discussed in this paper. In summary, the present work systematically reviews the previous observations and understanding of ash aging in particulate filters and identifies areas that need further research, which can be useful guidance for future studies.

42 ENGINEERING↗
DOE OSTI2023

Using micro-XRF to characterize chloride ingress through cold joints in 3D printed concrete

Digital fabrication methods with concrete have been rapidly developing, with many problems related to component production and material control being solved in recent years. These processes produce inherently layered cementitious components that are anisotropic, and in many cases, produces a weak interface between layers, which are generally referred to as cold joints. While material strength at these interfaces has been well studied in recent years, durability has received less attention, even though cold joints can function as channels for aggressive agents, such as chlorides. This work presents a method using micro-X-ray fluorescence (μXRF) to image chloride ingress into layer interfaces of 3D printed fine-grained concrete specimens produced with varying layer deposition time intervals, and also compares it to neutron imaging of moisture uptake. The results show that cold joints formed after a 1 day time interval are highly susceptible to chloride ingress, and that curing conditions play a major role in how quickly interfacial transport can take place. The μXRF method is also shown to be useful for study of transport of chlorides in cold joints, due to its spatial resolution and direct analysis of an aggressive species of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Probing and Understanding the Spatial and Energy Distributions of Plasmonic Hot Carriers via Single-Molecule Quantum Transport

The generation of hot carriers in plasmonic nanostructures offers new transformative opportunities for a range of applications, including solar energy harvesting, photochemistry, and photodetection. To better understand such phenomena, this project aims to develop an experimental approach to study and understand hot carrier distributions in photoexcited plasmonic nanostructures. The developed experimental approach, combining scanning probe-based single-molecule quantum transport measurements with plasmonic excitation methods, leverages nanoscopic single molecules as quantum energy filters to selectively access hot carriers in a variety of plasmonic nanostructures. This report summarizes the key outcomes, impacts, and lessons learned in the first two years of the project. Over the past two years, important progress was made in achieving the project's goals, including the development of the scanning probe-based experimental platform, fabrication of plasmonic nanostructures, and design and synthesis of the desired molecular candidates. Insights from this research will pave the path toward the development of technologies that will improve the way we harvest solar energy and drive chemical reactions and ultimately help transition the U.S. to a clean and sustainable society.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Evaluation of Attila and MCNP computational methods for dose and exposure estimation

Radiation transport calculations are often used to estimate dose or exposure to components and personnel surrounding a radiation source. The sources for these calculations are decaying radionuclides within various nuclear materials. Historically, dose calculations use MCNP (Monte Carlo N-Particle) transport code as the primary particle transport tool without a secondary computational tool to validate the results from the MCNP simulations [1]. The goal of this study is to make an independent check of the Monte Carlo solution from MCNP6 Version 6.2.1 with the discrete ordinates solution from Attila 10.2.0 Beta 3. As an example problem for this study, water-filled, stainless-steel vessels, modeled with an unstructured mesh (UM) with both MCNP and Attila [2], are exposed to 252Cf and 60Co point sources. This report also includes a discussion of the limitations of unstructured mesh in a MCNP calculation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Depletion capabilities in the OpenMC Monte Carlo particle transport code

A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of k eff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate that it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗
DOE OSTI2025

Software For Advanced Large-scale Analysis Of Magnetic Confinement For Numerical Design, Engineering & Research (salamander)

As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems

Simon, Pierre-Clement [Idaho National Laboratory (↗
DOE OSTI2023

Methane and methane-carbon dioxide activated synergystic biomass gasification for hydrogen rich syngas production

In one aspect, the disclosure relates to methods for biomass gasification to produce sustainable and renewable alternatives to fossil fuels including, but not limited to syngas having a high H2 content. The method can produce an H2/CO ratio close to 2:1, which is desirable for further chemical or transportation fuel synthesis. In another aspect, the methods disclosed herein have high yields and make use of agricultural and industrial waste (e.g., hardwood pellets and grain stovers) as starting materials. In a further aspect, the methods disclosed herein can produce useful byproducts including, but not limited to, carbon nanofibers (CNF). This abstract is intended as a scanning tool for purposes of searching in the particular art and is not intended to be limiting of the present disclosure.

Hu, Jianli↗
DOE OSTI2022

Real-Time Highly Resolved Spatial-Temporal Vehicle Energy Consumption Estimation Using Machine Learning and Probe Data

Real-time highly resolved spatial-temporal vehicle energy consumption is a key missing dimension in transportation data. Most roadway link-level vehicle energy consumption data are estimated using average annual daily traffic measures derived from the Highway Performance Monitoring System; however, this method does not reflect day-to-day energy consumption fluctuations. As transportation planners and operators are becoming more environmentally attentive, they need accurate real-time link-level vehicle energy consumption data to assess energy and emissions; to incentivize energy-efficient routing; and to estimate energy impact caused by congestion, major events, and severe weather. This paper presents a computational workflow to automate the estimation of time-resolved vehicle energy consumption for each link in a road network of interest using vehicle probe speed and count data in conjunction with machine learning methods in real time. The real-time pipeline can deliver energy estimates within a couple seconds on query to its interface. The proposed method was evaluated on the transportation network of the metropolitan area of Chattanooga, Tennessee. The volume estimation results were validated with ground truth traffic volume data collected in the field. To demonstrate the effectiveness of the proposed method, the energy consumption pipeline was applied to real-world data to quantify road transportation-related energy reduction because of mitigation policies to slow the spread of COVID-19 and to measure energy loss resulting from congestion.

Severino, Joseph↗
DOE OSTI2023

What is the long-term demand for liquid hydrocarbon fuels and feedstocks?

Liquid hydrocarbons made from crude oil serve many functions: (1) a dense, easy-to-store, easy-to-transport energy source, (2) a method for daily-to-seasonal energy storage, (3) a chemical feedstock, (4) a chemical reducing agent and (5) a method to enhance high-temperature heat transfer in many furnaces and industrial processes. Liquid hydrocarbons can be produced and used without increasing atmospheric carbon dioxide levels if made from non-fossil feedstocks such as carbon dioxide or biomass. Understanding future liquid hydrocarbon demand is the starting point in assessing the viability of such options. Our assessment is that U.S. demand for liquid hydrocarbons is unlikely to go below the equivalent of 10 million barrels per day of crude oil but could be as high as 20 million barrels per day. The costs to replace liquid hydrocarbons increases rapidly at lower liquid hydrocarbon consumption rates. Hydrocarbon biofuels from cellulosic feedstocks can meet such demands if produced with large external inputs of heat and hydrogen. Options based on more limited feedstocks (starches, plant oils, sugars, etc.) can’t meet such demands.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

The collisional particle-in-cell method for the Vlasov–Maxwell–Landau equations

We introduce an extension of the particle-in-cell method that captures the Landau collisional effects in the Vlasov–Maxwell–Landau equations. The method arises from a regularisation of the variational formulation of the Landau equation, leading to a discretisation of the collision operator that conserves mass, charge, momentum and energy, while increasing the (regularised) entropy. The collisional effects appear as a fully deterministic effective force, thus the method does not require any transport–collision splitting. The scheme can be used in arbitrary dimension, and for a general interaction, including the Coulomb case. We validate the scheme on scenarios such as the Landau damping, the two-stream instability and the Weibel instability, demonstrating its effectiveness in the numerical simulation of plasma.

Bailo, Rafael (ORCID:0000000180183799)↗
DOE OSTI2021

Microscale Mechanical‐Chemical Modeling of Granular Salt: Insights for Creep

Abstract Microscale numerical modeling is potentially an effective approach for understanding salt creep, but it is quite challenging because natural salt grain boundaries can have complex and evolving geometry, and because contacts between these deformable salt mineral grains are dynamic as a result of compaction, chemical reaction, fluid flow, and heat transfer. In this study, we have overcome these challenges and developed a new microscale mechanical‐chemical (MC) model to analyze creep of salt at the microscale, accounting for coupled deformation, dynamic contacts, and chemical reaction in granular systems with realistic geometric representations. The MC model was realized by linking a new microscale mechanical code based on the numerical manifold method (NMM) to a reactive transport code named Crunch. We simulated the processes of reorganization of the salt grains, microfracturing, and pressure solution that contribute to creep at larger scales. Based on this first quantitative microscale model, we found that sharp corners of mineral grains can dominate the contact dynamics, microfracturing, and pressure solution when salt is compacted, thus governing the structural changes and porosity loss of the system. We found that pressure solution, which preferentially dissolves sharp corners and edges, can lead to relatively high porosity loss in the system, thus playing an important role in the creep of salt. Our analysis shows that the dynamic changes of salt granular systems involving grain relocation and pressure solution can occur repeatedly and continuously, thus contributing to longer‐term creep of salt at larger scales.

58 GEOSCIENCES↗