Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transitional flow”

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 253 records · Page 14

Evidence of short chains in liquid sulfur

High energy x-ray pair distribution function measurements show the average coordination number of the first shell in liquid sulfur is 1.86 ± 0.04 across the λ-transition, not precisely 2.0 as widely accepted. This indicates that upon melting, liquid sulfur does not comprise solely of S 8 rings but also possesses a significant number of short chains. Intensities of the pre-peak and first diffraction peak of the x-ray structure factor and third peak height of the pair distribution function all show deviations at the λ-transition temperature T λ , associated with the break-up of S 8 rings and the start of oligomer polymerization. A significant number of non-bonded or loosely bonded “interstitial atoms,” with an average coordination number of 0.20 ± 0.005, are also observed in the so-called “forbidden zone” between the first and second shells upon melting. The number of interstitial atoms is found to decrease to a minimum at the λ-transition, but the majority persist into the high temperature polymerized liquid. Furthermore, the existence of short chains and nearby interstitial atoms represent the two main factors required to initiate the S 8 -ring to chain transition, as proposed by recent molecular dynamics simulations.

Chemical bonding↗

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optical transient grating pumped X-ray diffraction microscopy for studying mesoscale structural dynamics

Abstract A fundamental understanding of materials’ structural dynamics, with fine spatial and temporal control, underpins future developments in electronic and quantum materials. Here, we introduce an optical transient grating pump and focused X-ray diffraction probe technique (TGXD) to examine the structural evolution of materials excited by modulated light with a precisely controlled spatial profile. This method adds spatial resolution and direct structural sensitivity to the established utility of a sinusoidal transient-grating excitation. We demonstrate TGXD using two thin-film samples: epitaxial BiFeO 3 , which exhibits a photoinduced strain (structural grating) with an amplitude proportional to the optical fluence, and FeRh, which undergoes a magnetostructural phase transformation. In BiFeO 3 , structural relaxation is location independent, and the strain persists on the order of microseconds, consistent with the optical excitation of long-lived charge carriers. The strain profile of the structural grating in FeRh, in comparison, deviates from the sinusoidal excitation and exhibits both higher-order spatial frequencies and a location-dependent relaxation. The focused X-ray probe provides spatial resolution within the engineered optical excitation profile, resolving the spatiotemporal flow of heat through FeRh locally heated above the phase transition temperature. TGXD successfully characterizes mesoscopic energy transport in functional materials without relying on a specific transport model.

36 MATERIALS SCIENCE↗

Impact of plasma density/collisionality on divertor heat flux width

Both ASDEX-Upgrade (AUG) data and the generalized HD (GHD) model showed that the scrape-off width broadens as the density/collisionality increases [1, 2]. A series of BOUT++ transport simulations are performed to study the physics of the scaling characteristics of the divertor heat flux width vs density/collisionality via a plasma density scan with either fixed pressure profile or fixed temperature profile inside separatrix. Additionally, the simulations show that even in the drift dominated regime, the divertor heat flux width can be broadened due to the transition of the SOL residence time from the parallel particle flow time to the enhanced parallel conduction time as the collisionality/density increases as posited in the GHD model. In addition, the heat flux width is found to be proportional to the square root of ion mass for low collisionality while it has a weakly dependence on ion mass for high collisionality. Furthermore, our simulations show that as the density increases, the radial electric field (E r ) well shallows, which potentially weakens E r × B flow shear stabilization of turbulence at high density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large-eddy simulations of turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines rotating at low tip speed ratios

Turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines (VAWTs) rotating at low tip speed ratios (TSRs) are studied numerically. The turbulent flows are simulated using the large-eddy simulation (LES) model, and the rotating turbine blades are modeled using the actuator line method. The helical VAWT has identical key parameters as the straight VAWT except for the 135° helical twist of the blades over the 0.3m vertical span. A set of LES runs are performed for two TSRs, 0.6 and 0.4, and the results are reported and analyzed. At these low TSRs, the wake behind the straight-bladed VAWT exhibits two-dimensional dominant flow motions (in the horizontal plane perpendicular to the straight blades) in the near-wake region that cause considerable spanwise expansion of the wake as it extends downstream. In contrast, the helical-bladed VAWT generates highly three-dimensional (3D) wake flow structures and upward/downward mean flow motions within the wake that cause the wake to expand mainly in the vertical direction. Turbulence statistical analyses also show that the 3D wake flow features induced by the helical blades accelerate the wake transition to turbulence and enhance the small-scale turbulent dissipation (as shown by the subgrid-scale turbulent dissipation in the LES), which leads to a more rapid decay of the wake turbulence intensity than that in the straight-bladed VAWT case at the same TSR. Finally, compared with the straight-bladed VAWT, the helical-bladed VAWT also exhibits much smaller temporal variations for the torque and power coefficients during the rotation cycle, which can be beneficial for wind power generation.

17 WIND ENERGY↗

In-cylinder spray evolution in a motored central-injection gasoline engine: Imaging and simulating the effects of flash-boiling and intake crossflow

Accurate predictions of fuel spray behavior and mixture formation in simulations of direct-injection spark-ignition (DISI) engines are fundamental to ensure proper description of all subsequent processes including ignition, combustion, and emissions. In this work, the spray evolution in a single-cylinder optical DISI engine was studied experimentally and numerically with the goal of enabling predictive computational fluid dynamics (CFD) modeling of in-cylinder sprays. The authors explored a wide range of operating conditions characterized by several fuel injection temperatures and engine speeds, using a well-characterized nine-component gasoline surrogate known as PACE-20. The effect of flash boiling and intake crossflow on the spray is discussed, with a focus on evaluating the ability of the spray models to capture highly transient spray behavior. In the experiments, the fuel temperature was varied between 20°C and 80°C, allowing for non-flash- to flash-boiling transition to emerge with enhanced flashing intensity at the highest temperatures. Spray collapse resulted in vapor-rich regions, owing to the locally lower inertia of the fluid. Varying the engine speed from 650 to 1950 rpm promoted increasingly more turbulent in-cylinder crossflow which interacted with the spray during the injection event and resulted in enhanced spray dispersion. The CFD model was able to capture the spray morphology transition at different fuel temperatures and engine speeds adequately. Further, it is shown that the spray breakup model could capture the transitional spray behavior induced by flash boiling atomization and intake flow via proper initialization of the spray cone angle and calibration of the spray models’ constants.

33 ADVANCED PROPULSION SYSTEMS↗

Hypersonic Fluid–Structure Interaction on a Cone–Slice–Ramp Geometry

In this study, fluid–structure interactions were measured between a representative control surface and the hypersonic flow deflected by it. The control surface is simplified as a spanwise finite ramp placed on a longitudinal slice of a cone. The front surface of the ramp contains a thin panel designed to respond to the unsteady fluid loading arising from the shock-wave/boundary-layer interactions. Experiments were conducted at Mach 5 and Mach 8 with ramps of different angles. High-speed schlieren captured the unsteady flow dynamics and accelerometers behind the thin panel measured its structural response. Panel vibrations were dominated by natural modes that were excited by the broadband aerodynamic fluctuations arising in the flowfield. However, increased structural response was observed in two distinct flow regimes: 1) attached or small separation interactions, where the transitional regime induced the strongest panel fluctuations. This was in agreement with the observation of increased convective undulations or bulges in the separation shock generated by the passage of turbulent spots, and 2) large separated interactions, where shear layer flapping in the laminar regime produced strong panel response at the flapping frequency. In addition, panel heating during the experiment caused a downward shift in its natural mode frequencies.

33 ADVANCED PROPULSION SYSTEMS↗

Pore‐Scale Modeling of Carbon Dioxide and Hydrogen Transport During Geologic Gas Storage

Abstract Geologic storage of CO 2 and H 2 are climate‐positive techniques for meeting the energy transition. While similar formations could be considered for both gases, the flow dynamics could differ due to differences in their thermophysical properties. We conduct a rigorous pore‐scale study of water/CO 2 and water/H 2 systems at relevant reservoir conditions in a Bentheimer rock sample using the lattice Boltzmann method to quantify the effects of capillary, viscous, inertial, and wetting forces during gas invasion. At similar conditions, H 2 invasion is weaker compared to CO 2 due to unfavorable viscosity ratios. Increasing flow rate, however, increases the breakthrough saturation for both gas systems in the range of capillary numbers studied. At isolated conditions of flow rate, viscosity ratio, and wettability, local inertial effects are found to be critical and show consistent increase in the invaded gas saturation. The effect of inertial forces persits for both gases across all field conditions tested.

58 GEOSCIENCES↗

E-transit-bench: simulation platform for analyzing electric public transit bus fleet operations

When electrified transit systems make grid aware choices, improved social welfare is achieved by reducing grid stress, reducing system loss, and minimizing power quality issues. Electrifying transit fleet has numerous challenges like non availability of buses during charging, varying charging costs and so on, that are related the electric grid behavior. However, transit systems do not have access to the information about the co-evolution of the grid's power flow and therefore cannot account for the power grid's needs in its day-to-day operation. In this paper we propose a framework of transportation-grid co-simulation, analyzing the spatio-temporal interaction between the transit operations with electric buses and the power distribution grid. Real-world data for a day's traffic from Chattanooga city's transit system is simulated in SUMO and integrated with a realistic distribution grid simulation (using GridLAB-D) to understand the grid impact due to transit electrification. Charging information is obtained from the transportation simulation to feed into grid simulation to assess the impact of charging. We also discuss the impact to the grid with higher degree of transit electrification that further necessitates such an integrated transportation-grid co-simulation to operate the integrated system optimally. Our future work includes extending the platform for optimizing the charging and trip assignment operations.

Sen, Rishav↗

Coarse grained simulation of convectively driven turbulent mixing, transition, and turbulence decay

Accurate predictions with quantifiable uncertainties are essential to many practical turbulent flow applications exhibiting extreme geometrical complexity and broad ranges of length and time scales. Under-resolved computer simulations are typically unavoidable in such applications, and implicit large-eddy simulation (ILES) often becomes the effective strategy. In this work, we focus on assessing ILES initialized with well-characterized 25 6 3 homogeneous isotropic turbulence datasets generated with direct numerical simulation (DNS). ILES solutions based on the LANL xRAGE code are studied as function of resolution for 6 4 3 , 12 8 3 , 25 6 3 , and 51 2 3 grids. ILES performance of new directionally-unsplit high-order numerical hydrodynamics algorithms in xRAGE is examined in this context. Compared to the initial 25 6 3 DNS, we find longer inertial subranges and higher turbulence Reynolds number (Re) for 25 6 3 and 51 2 3 xRAGE — attributed to having linked DNS (Navier–Stokes based) solutions to nominally inviscid (higher Re) Euler based ILES solutions. For fixed 25 6 3 resolution, we find that significantly larger simulated turbulence Re can be achieved with the higher-order unsplit (vs. split) discretizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Behavior of an air-cooled thermosiphon during a non-air gas ingress event

Here, this study investigates the transient response of an air-cooled thermosiphon subjected to a sudden non-air gas ingress event. Results presented are from experimental testing conducted on a large-scale thermal hydraulic facility that examines heat removal performance of the Reactor Cavity Cooling System (RCCS) concept, a passive safety system for advanced nuclear reactors. Relying on buoyancy driven natural circulation flow, these systems feature a series of riser standpipes and coupled chimney network that provide a pathway for heat rejection to the atmosphere. In this work, the facility was modified to allow a transition of the draft intake from normal ambient air, to draft intake from an open-top volume containing high purity argon. Testing began by establishing normal operation and allowing the facility to reach steady-state thermal hydraulic flow conditions with natural air draft. The transitioned sequence was then initiated which simultaneously changed the inlet boundary from open air to a pathway from the quiescent argon volume. This event caused 1,200 cu-feet of the heavy gas, twice the internal volume of total facility flow path, to be drawn into the inlet plenum and ingress into the heated riser standpipes. With the upper chimney network still containing residual air, there was insufficient density difference to maintain buoyancy driven natural circulation. After a period of only 90 s past the ingress event, the facility experienced complete flow stagnation before entering an extended period of severely degraded system flow. Due to the cessation of bulk fluid movement and subsequent failure of its heat removal function, fluid and structural temperatures began to rise sharply. Re-circulation patterns developed within the multiple parallel riser standpipes, where hot gas near the outlets was observed to travel downward and re-enter the inlet of adjacent channels. After approximately 18-minutes, fluid temperatures and their associated density difference rose to a level sufficient to allow reestablishment of buoyancy driven system flow, and ultimately, recovery of facility operation to normal behavior.

42 ENGINEERING↗

Changes in soil water content and lateral flow exert large effects on soil thermal dynamics across Alaskan landscapes

Both lateral surface and subsurface water flow affect soil moisture dynamics, yet most land surface models only solve subsurface water movement vertically. Here, we use a 3D ecosystem model that considers both land surface and subsurface hydrologic processes to simulate soil moisture, which is then used to drive a 1-D vertical soil thermal model to simulate the soil moisture effects on soil thermal dynamics in central Alaska. Our coupled model improves soil temperature (ST) estimates by 43.5% in comparison with observational data. Soil moisture has little effect on ST during the wet season (-1.5%) and a substantial influence during the dry season (60%). Spatially, water lateral flow has significant impacts on both soil moisture and ST, causing model estimates for thawed areas in the transition season to increase by ~10% in the study area. Our results highlight the importance of considering dynamical soil moisture, as well as lateral flow effects, on soil thermal dynamics in permafrost regions.

54 ENVIRONMENTAL SCIENCES↗

The high-density equation of state in heavy-ion collisions: constraints from proton flow

Abstract A set of different equations of state is implemented in the molecular dynamics part of a non-equilibrium transport simulation (UrQMD) of heavy-ion collisions. It is shown how different flow observables are affected by the density dependence of the equation of state. In particular, the effects of a phase transition at high density are explored, including an expected reduction in mean $$m_T$$ m T . We also show that an increase in $$v_2$$ v 2 is characteristic for a strong softening of the equation of state. The phase transitions with a low coexistence density, $$n_{\text {CE}}<4 n_0$$ n CE < 4 n 0 , show a distinct minimum in the slope of the directed flow as a function of the beam energy, which would be a clear experimental signal. By comparing our results with experimental data, we can exclude any strong phase transition at densities below $$4n_0$$ 4 n 0 .

74 ATOMIC AND MOLECULAR PHYSICS↗

Cyber-Physical Simulation of an Innovative Solid Oxide Electrolysis Cell - Gas Turbine (SOEC-GT) Hybrid Energy System

To produce green hydrogen at scale and at low cost, the solid oxide electrolysis cell (SOEC) systems would be tied to renewable power sources (mainly solar and wind), leveraging zero-carbon electricity at low prices and even nearly free during overgeneration scenarios. However, due to the intermittency of renewable power generation, the SOEC system is subjected to rapid load transitions that occur not only in diurnal cycles but also in short timeframes (e.g., sub-minute). This can result in fast degradation and thus greatly reduce the SOEC’s lifetime. Our team at the National Energy Technology Laboratory has demonstrated that the anode air flow can have a crucial role in SOEC thermal management. Typically, a higher air flow rate can help to mitigate the local temperature gradient distortion during rapid load transitions. To move a large amount of air, the gas turbine (GT; i.e., compressor-turbine-generator set) has been regarded as the most mature and efficient technology. In this presentation, we show the configuration and preliminary results of an innovative SOEC-GT hybrid energy system in a cyber-physical simulation (CPS) approach.

Zhang, Biao↗

Development of PFLOTRAN Transport Capability for Use in the Waste Isolation Pilot Plant Performance Assessment - 20545

Waste Isolation Pilot Plant (WIPP) performance assessment (PA) calculations estimate the probability of radionuclide release from the repository to the land surface and across the land withdrawal boundary for a regulatory period of 10,000 years after facility closure. Simulations of flow and transport in the repository and the surrounding Salado Formation are foundational to the PA. Because proposed additional waste emplacement panels would result in an asymmetric repository layout, the US Department of Energy (DOE) is preparing to transition to use of a three-dimensional (3-D) model domain for simulation of flow and transport instead of the two-dimensional (2-D) flared grid domain currently used. DOE has charged Sandia National Laboratories with developing the capability necessary to simulate processes affecting flow and transport in the WIPP in PFLOTRAN, an open-source massively parallel multi-phase flow and reactive transport code. The new flow and transport capabilities developed in PFLOTRAN incorporate WIPP-specific process models and will replace the 2-D simulators (BRAGFLO and NUTS) that are currently utilized for Salado flow and transport calculations in WIPP PA. The focus of this paper is on the development of a new Nuclear Waste Transport (NWT) mode in PFLOTRAN that has all of the capabilities necessary for Salado transport simulations, including the ability to handle complete dry-out (100% gas saturation) of arbitrary cells in the model domain, radionuclide mass conservation at step changes in porosity associated with borehole intrusion, and the ability to calculate fluxes on a flared grid. The new PFLOTRAN transport capability and a suite of verification tests were designed around a list of functional requirements for WIPP PA calculations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Characterizing the Non-linear Interactions Between Tide, Storm Surge, and River Flow in the Delaware Bay Estuary, United States

Low-lying coastal areas in the mid-Atlantic region are prone to compound flooding resulting from the co-occurrence of river floods and coastal storm surges. To better understand the contribution of non-linear tide-surge-river interactions to compound flooding, the unstructured-grid Finite Volume Community Ocean Model was applied to simulate coastal storm surge and flooding in the Delaware Bay Estuary in the United States. The model was validated with tide gauge data in the estuary for selected hurricane events. Non-linear interactions between tide-surge-river were investigated using a non-stationary tidal analysis method, which decomposes the interactions’ components at the frequency domain. Model results indicated that tide-river interactions damped semidiurnal tides, while the tide-surge interactions mainly influenced diurnal tides. Tide-river interactions suppressed the water level upstream while tide-surge interaction increased the water level downstream, which resulted in a transition zone of damping and enhancing effects where the tide-surge-river interaction was prominent. Evident compound flooding was observed as a result of non-linear tide-surge-river interactions. Furthermore, sensitivity analysis was carried out to evaluate the effect of river flooding on the non-linear interactions. The transition zone of damping and enhancing effects shifted downstream as the river flow rate increased.

54 ENVIRONMENTAL SCIENCES↗

A study of velocity, temperature, and density in the plasma generated by laser-induced breakdowns

The paper presents velocity measurements of shock-induced flow field, leading to vortex generation and plasma deformation in air. Femtosecond-laser electronic excitation tagging (FLEET) velocimetry was performed at the times Δt = 3, 20, 50, and 100 µs post laser-induced breakdown. Emissions over Δt = 3–10 µs showed the propagations of the initially elliptic shock, transitioning into a spherical front. The shock emanated along the laser axis causes the flow outward, and then the pressure gradient generated by the rarefaction wave drives the inward flow at later moments, with the velocity magnitude approaching a steady-state value of 40 m s -1 . Temporal velocity evolution was compared with non-self-similar solutions behind the propagating shock, which are sensitive to the size of the energy deposition, and the use of the measured initial plasma diameter reproduced the experiment. There establishes a region of uniform velocity around 35 m s -1 in the air-flow running through the plasma, which triggers the roll-up of the plasma surface by a large-scale vortex, providing the detail of flow field evolution from the shock propagation to the plasma deformation. A collective Thomson scattering and hybrid fs/ps pure rotational coherent anti-Stokes Raman scattering (CARS) were also performed to gain insight into the high-temperature plasma. An effective electron–ion recombination rate of 2 x 10 -12 cm 3 s -1 was measured at Δt = 0.5–10 µs, during a dynamic plasma expansion and compression. When the shock resides close to the plasma at Δt = 0.5–1 µs, the temperature distributions were found to follow the similarity law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Formation of Magnetically Truncated Accretion Disks in 3D Radiation-transport Two-temperature GRMHD Simulations

Multiwavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nucleus transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole (BH). However, the formation, structure, and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at ~35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semiopaque corona reach, respectively, $T_i$ ≳ 10 10 K and $T_e$ ≳ 5 × 10 8 K. The system produces radiation, powerful collimated jets, and broader winds at the total energy efficiency exceeding 90%, the highest ever energy extraction efficiency from a spinning BH by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium may naturally lead to broadened iron line emission observed in the hard state.

79 ASTRONOMY AND ASTROPHYSICS↗