Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW THEORY”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

High-Rate, Selective Electrosynthesis of Cyclohexanone Oxime via In Situ Generation and Release of Hydroxylamine on Bismuth

Oxime compounds are key industrial intermediates for nylon precursors and commodity chemicals. However, conventional routes rely on multistep reactions and hydroxylamine (NH 2 OH) salts, raising significant safety and sustainability concerns. Although electrosynthesis offers an alternative, oxime formation on d-block transition metals suffers from poor selectivity, as nitrogen oxyanion intermediates bind strongly to the surface and are readily over-reduced to ammonia. Here, we report morphology-controlled p-block bismuth rhombic dodecahedra (Bi RDs) that promote in situ NH 2 OH generation and its desorption into the electrolyte, enabling an electrochemical-chemical decoupled route for cyclohexanone oxime (CHO) synthesis. Bi RDs deliver nearly 100% Faradaic efficiency (FE) at −0.5 V vs. RHE and a yield of 1.4 mmol h –1 cm –2 at −0.9 V vs. RHE in an H-cell, while maintaining a CHO selectivity of nearly 100% at 100 mA cm –2 in a flow cell. Under identical conditions, d-block electrodes (Cu, Pd, Ag) show FE below 30%. Density functional theory calculations reveal that Bi 6p orbital-derived surface states weaken intermediate binding and facilitate NH 2 OH desorption, suppressing over-reduction. Kinetic analysis, post-addition trapping experiments, and in situ ATR-FTIR and Raman spectroscopy suggest the following reaction mechanism: NH 2 OH is selectively generated at the electrode surface, released as a freely diffusing intermediate, and undergoes homogeneous condensation with cyclohexanone in the bulk electrolyte, bypassing the surface-confined Langmuir–Hinshelwood pathway. These findings demonstrate that regulating intermediate desorption through p-block orbital chemistry provides a general strategy for achieving high selectivity in electro-organic nitrogen synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts

Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Pressure–strain interaction. I. On compression, deformation, and implications for Pi-D

The pressure–strain interaction describes the rate per unit volume that energy is converted between bulk flow and thermal energy in neutral fluids or plasmas. The term has been written as a sum of the pressure dilatation and the collisionless analog of viscous heating referred to as Pi–D , which isolates the power density due to compressible and incompressible effects, respectively. It has been shown that Pi–D can be negative, which makes its identification as collisionless viscous heating troubling. We argue that an alternate decomposition of pressure–strain interaction can be useful for interpreting the underlying physics. Since Pi–D contains both normal deformation and shear deformation, we propose grouping the normal deformation with the pressure dilatation to describe the power density due to converging/diverging flows, with the balance describing the power density purely due to shear deformation. We then develop a kinetic theory interpretation of compression, normal deformation, and shear deformation. We use the results to determine the physical mechanisms that can make Pi–D negative. We argue that both decompositions can be useful for the study of energy conversion in weakly collisional or collisionless fluids and plasmas, and implications are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dynamics of active liquid interfaces

Controlling interfaces of phase-separating fluid mixtures is key to the creation of diverse functional soft materials. Traditionally, this is accomplished with surface-modifying chemical agents. Using experiment and theory, we studied how mechanical activity shapes soft interfaces that separate an active and a passive fluid. Chaotic flows in the active fluid give rise to giant interfacial fluctuations and noninertial propagating active waves. At high activities, stresses disrupt interface continuity and drive droplet generation, producing an emulsion-like active state composed of finite-sized droplets. When in contact with a solid boundary, active interfaces exhibit nonequilibrium wetting transitions, in which the fluid climbs the wall against gravity. These results demonstrate the promise of mechanically driven interfaces for creating a new class of soft active matter.

Science & Technology - Other Topics↗
DOE OSTI2025

Construction of 3D MHD pressure drop correlation and flow characterization in the contraction region of a fusion blanket manifold

Inlet and outlet manifolds are typical components of liquid metal (LM) blanket designs of a fusion power reactor to be used to distribute the LM flow into breeding channels and collect it at the exit of the blanket. High pressure loss in the magnetohydrodynamic (MHD) flows featuring abrupt geometrical changes is one of the main feasibility issues of such designs. Recently, optimization studies were conducted to construct 3D MHD pressure drop correlations for a LM flow in an electrically insulating manifold with gradual expansion. Here, the 3D computational approach developed in that study is applied to the outlet manifold featuring gradual contraction. A systematic analysis was performed with a total number of 135 flow cases computed with COMSOL Multiphysics for Hartmann numbers 1000 < Ha < 10,000, Reynolds numbers 100 < Re < 12,000, and contraction angles 45° < θ < 75° for a fixed contraction ratio of 4. The effects of Ha, Re and θ on the flow recirculation, development length and the total pressure drop were carefully examined. A linear regression analysis was used to determine the power rule of pressure drop coefficient k related to Ha and Re, demonstrating a good match with the Ludford layer theory. Eventually, a correlation for the 3D MHD pressure drop coefficient was constructed as a function of Ha, Re and θ. Further, the results were compared against the inlet manifold. It was found that the flow in the inlet manifold exhibits larger recirculation zones. In the investigated range of Ha, Re and θ, the pressure drop coefficient k of the LM MHD flow in the gradual contraction is only slightly lower (< 8 %) than that in the gradual expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Invariant Forms of Dissolution Fingers

Dissolution of fractured and porous media introduces a positive feedback between fluid flow and reactant transport, leading to the emergence of pronounced, fingerlike channels. We investigate the formation of these structures using a microfluidic Hele-Shaw cell with a soluble bottom. Our experiments show that the shape of dissolution fingers is invariant and reveals itself over time as the fingers extend into the system. Here, by combining reactive-transport theory and conformal mapping techniques we derive these invariant forms. We relate these results to natural dissolution fingers in karst landscapes, and illustrate how to determine the groundwater flow rate responsible for their formation based on the finger shape.

Flows in porous media↗
DOE OSTI2024

Dark matter raining on DUNE and other large volume detectors

Direct detection is a powerful means of searching for particle physics evidence of dark matter (DM) heavier than about a GeV with o(kiloton) volume, low-threshold detectors. In many scenarios, some fraction of the DM may be boosted to large velocities enhancing and generally modifying possible detection signatures. We investigate the scenario where 100% of the DM is boosted at the Earth due to new attractive long-range forces. This leads to two main improvements in detection capabilities: (1) the large boost allows for detectable signatures of DM well below a GeV at large-volume neutrino detectors, such as DUNE, Super-K, Hyper-K, and JUNO, as possible DM detectors, and (2) the flux at the Earth’s surface is enhanced by a focusing effect. In addition, the model leads to a significant anisotropy in the signal with the DM flowing dominantly vertically at the Earth’s surface instead of the typical approximately isotropic DM signal. We develop the theory behind this model and also calculate realistic constraints using a detailed GENIE simulation of the signal inside detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Ligand-coordination effects on the selective hydrogenation of acetylene in single-site Pd-ligand supported catalysts

The selective hydrogenation of acetylene to ethylene is a critical step in the synthesis of polyethylenes. Achieving high conversion to ethylene without over-hydrogenation to ethane is a challenge that requires control of the transition metal site, which we achieve through a ligand-coordinated supported catalyst (LCSC) strategy. Using Pd catalysts coordinated to 1,10-phenanthroline-5,6-dione (PDO) ligands on CeO2 supports, we have discovered that the reaction selectivity depends strongly on the ligand:metal ratio with higher selectivity when more ligand is present in the catalyst. Catalyst structure was examined by extended X-ray absorption fine structure spectroscopy, transmission electron microscopy, and CO adsorption, which indicate single-atom character of the Pd. The ligand:metal ratio is determined by Xray photoelectron spectroscopy measurements and correlated with hydrogenation reactions under steady-state flow conditions to examine trends in hydrogenation activity and selectivity. Those trends can be better understood by density functional theory calculations that indicate hydrogen binding on the ligand to guide reaction selectivity toward the desired ethylene hydrogenation product. Here these results demonstrate the importance of considering the dynamic character of LCSCs and inform the design of future single-site heterogeneous catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Adsorption of Ar into zeolite Al-MFI (NH 4 )

The crystal structure of anhydrous Al-MFI (NH 4 ) containing adsorbed Ar has been determined and refined using synchrotron X-ray powder diffraction data taken at 90 K, and optimized using density functional theory techniques. Six highly occupied Ar sites almost completely fill the pore volume of the zeolite. Changing the gas flow from Ar to He at 90 K decreases the Ar occupancies of all six sites, but two decrease more than the others. Warming the sample from 90 to 295 K in Ar flow results in further decreases in site occupancies, but five of the original six sites persist.

Materials Science↗
DOE OSTI2023

Revealing topology in metals using experimental protocols inspired by K-theory

Topological metals are conducting materials with gapless band structures and nontrivial edge-localized resonances. Their discovery has proven elusive because traditional topological classification methods require band gaps to define topological robustness. Inspired by recent theoretical developments that leverage techniques from the field of C*-algebras to identify topological metals, here, we directly observe topological phenomena in gapless acoustic crystals and realize a general experimental technique to demonstrate their topology. Specifically, we not only observe robust boundary-localized states in a topological acoustic metal, but also re-interpret a composite operator—mathematically derived from the K-theory of the problem—as a new Hamiltonian whose physical implementation allows us to directly observe a topological spectral flow and measure the topological invariants. Our observations and experimental protocols may offer insights for discovering topological behaviour across a wide array of artificial and natural materials that lack bulk band gaps.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Toward a least-effort principle for evaluating prices of elements as indicators of sustainability

In this article, we explore simplifying mechanisms in materials economics for conceptionally organizing the flows of chemical elements into society. We use the concept of rank, which plays a crucial role in information theory where power laws are involved, and is a useful measure for seeking extremal principles. Borrowing from an analogy to Zipf’s law leads to a power-law relationship between abundance and rank that—pressing the Zipf’s law analogy further—may suggest that the price paid for elemental materials optimizes efficient utilization. Factors such as weight, volume, and moles in the earth's crust or in global mining production of elements are used in our analysis. For this research, 42 elements with price and production data for the year 2010 (a relatively typical year) were considered. We observe that the role of the volume from global mining production dominates in materials economics and indicate that elements such as Te and He should be sold at prices at leas ~ 100–1000 times more than their present price for efficient utilization.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2026

Process-oriented evaluation of quasi-stationary Rossby waves and their impact on surface air temperature extremes in dynamical downscaling over North America

Quasi-stationary Rossby waves are a crucial component of the general circulation and play a significant role in regional water and energy cycles, as well as in extreme events. However, process-oriented evaluation for Rossby waves is rarely performed for dynamical downscaling simulations. To close this gap, we evaluate three classes of dynamical downscaling approaches, with a focus on quasi-stationary Rossby waves and their impact on surface air temperature over North America during Northern Hemisphere summer. The three classes of models differ in the way large-scale forcing is provided: a limited-area model (LAM) constrained only by lateral boundary conditions, represented by RegCM4 from the North American branch of the Coordinated Regional Downscaling Experiment (NA-CORDEX), a LAM with spectral nudging to maintain consistency in large-scale dynamics with the forcing data, represented by the Weather Research and Forecasting (WRF) model simulation in NA-CORDEX, and a global variable-resolution model with smoothly varying grid spacings, represented by the Community Atmosphere Model version 5.4, with the Model for Prediction Across Scales (MPAS) as its dynamical core (CAM-MPAS). With no constraints on the atmospheric dynamics, CAM-MPAS exhibits several mean biases in the upper-level circulations over the Pacific Coast region: a weaker subtropical jet, a northward-shifted mid-latitude jet, and an overestimated southerly flow. With the lateral boundary constraint alone, RegCM4 also exhibits weaker jets and overestimated southerly winds off the West Coast. Rossby ray theory reveals that those wind biases direct incoming Rossby waves northward. The erroneously routed Rossby waves distort the relationship between the accumulation of wave activity over the US West Coast and surface temperature anomalies over the Southern Great Plains, which emerges approximately 4 d after the convergence of wave-activity flux in the ERA-Interim reanalysis. Furthermore, the response of heatwaves to the extreme wave activity flux is not reproduced by the two models, a serious drawback as a dynamical downscaling framework is expected to connect large-scale forcing to local-scale phenomena. The WRF model employing spectral nudging is largely free from the aforementioned problems. A pair of sensitivity simulations suggests that spectral nudging is the key to improving the dynamics of quasi-stationary Rossby waves and their impact on surface air temperature. Our results also demonstrate the effectiveness of Rossby wave diagnostics that allow for realistic background flows for assessing the credibility of dynamical downscaling over North America, where incoming Rossby waves propagate through complex circulation patterns before traveling across the continent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Reexamination of hydrodynamic phonon transport in thin graphite

The recent hydrodynamic phonon transport theory for graphitic materials has been supported by the measurements of the second sound at temperatures up to about 100 K. When boundary scattering becomes comparable to momentum-conserving normal phonon scattering processes that are responsible for phonon hydrodynamics, Poiseuille phonon flow phenomena can emerge to give rise to unique size-dependent thermal conductivity in thin graphite. Here, we examine the thickness range for the Poiseuille phonon flow to become observable in thin graphite with the use of both deviational Monte Carlo simulation of the Peierls-Boltzmann transport equation and four-probe thermal transport measurements. As the basal-plane thermal conductivity calculated by prior first-principles theories saturates to the graphite value when the thickness is increased to five graphene layers, the phonon dispersion of graphite is used in the current calculations of thin graphite of micrometer thickness and a 23-layer thick ultrathin graphite (UTG) sample. The calculations show that diffuse surface scattering by surface defects can lead to Poiseuille phonon flow at 50 K in thin graphite with the thickness close to several micrometers but not in the 65 μm thin graphite and 23-layer UTG, where phonon scattering by the top and bottom surfaces become, respectively, much less and more frequent than the normal processes. In addition, the calculation results with the bulk graphite dispersion and diffuse surface scattering show decreased basal-plane thermal conductivity with decreasing thickness, opposite to recent thermocouple measurements of thin graphite samples. In comparison, the calculation results reveal that partially diffuse surface defect scattering can yield the four-probe measurement results of UTG samples, which are prepared here with an improved process to minimize surface contamination.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structure Formation Through Magnetohydrodynamical Instabilities in Primordial Disks

The shear flow instabilities under the presence of magnetic fields in the primordial disk can greatly facilitate the formation of density structures that serve as seeds prior to the onset of the gravitational Jeans instability. We evaluate the effects of the Parker, magnetorotational and kinematic dynamo instabilities by comparing the properties of these instabilities. We calculate the mass spectra of coagulated density structures by the above mechanism in the radial direction for an axisymmetric magnetohydrodynamic (MHD) torus equilibrium and power density profile models. Our local three-dimensional MHD simulation indicates that the coupling of the Parker and magnetorotational instability creates spiral arms and gas blobs in an accretion disk, reinforcing the theory and model. Such a mechanism for the early structure formation may be tested in a laboratory. The recent progress in experiments involving shear flows in rotating tokamak, field reversed configuration (FRC) and laser plasmas may become a key element to advance in nonlinear studies.

Noguchi, Koichi↗