SEARCH · Engineering Papers

Results for “drag”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2023

A gas dynamic perspective on particle lifting in electrostatic discharge-like devices using multiphase particle-in-cell approach

Although particle–laden electrostatic discharges are widely used in laboratories as well as in industrial applications, the mechanism of particle lifting for particles initially at rest in such highly unsteady systems is not well understood. A multiphase gas–particle solver is developed using the multiphase particle-in-cell (MP-PIC) approach to emulate the interaction of a compressible shock-dominated gas phase with the dense particle phase. First, the two-way coupled solver is initially used to study the interaction of a planar traveling shock with a vertical curtain of particulates. The gas and particle phase evolution was found to be in good agreement with a similar experimental study in Ling et al. [Phys. Fluids 24, 113301 (2012)]. Second, the MP-PIC code is used to study the interaction of an expanding blast wave with a thick bed of particles. The simulation considered forces such as quasi-steady drag, pressure-gradient, added-mass, Saffman, and Magnus forces. We observe that the vertical liftoff particles close to the shock impingement point in this configuration are associated with the quasi-steady drag, pressure gradient, and added-mass forces. Also, the Saffman lift and Magnus forces contribute to lifting particles located radially farther away from the shock impingement point. In addition, the study finds a decrease in particle lifting efficiency with decreasing plasma kernel length and shock strength.

42 ENGINEERING↗

DOE OSTI2024

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Driven probe particle dynamics in a bubble and pattern forming system

We numerically examine the dynamics of a probe particle driven at a constant force through an assembly of particles with competing long-range repulsion and short-range attraction that forms a bubble or stripe state. In the bubble regime, we identify several distinct types of motion, including an elastic or pinned regime where the probe particle remains inside a bubble and drags all other bubbles with it. There is also a plastic bubble phase where the bubble in which the probe particle is trapped is able to move past the adjacent bubbles. At larger drives, there is a breakthrough regime where the probe particle jumps from bubble to bubble and, in some cases, can induce correlated rotations or plastic rearrangements of the particles within the bubbles. At the highest drives, the probe particle moves sufficiently rapidly that the background particles undergo only small distortions. The distinctive dynamic flow states and the transitions between them are accompanied by signatures in the effective drag on the driven particle, jumps in the velocity–force curves, and changes in the time-dependent velocity fluctuations. We map the dynamic phase diagram for this system for varied interaction lengths, bubble sizes, and densities.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Gravitomagnetic vorticity generation in black hole accretion discs: a potential spatial constraint on plasma flow stability

We calculate the vorticity generation rate in the accretion disc near a slowly rotating black hole in the low-velocity, weak-field limit of general relativity. Specifically, we find that the frame-dragging effect due to the black hole’s rotation – manifested through the gravitomagnetic field – can generate vorticity in a moving plasma in the accretion disc. The mechanism remains operational as long as the accretion disc has non-negligible vertical height and is independent of the exact thermodynamical profile of the disc. The enstrophy density generation rate, as a measure of turbulence and dissipation, is presented, which indicates that the frame-dragging effect can disrupt the stability of the disc away from the z = 0 plane.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Interlayer Exciton Condensates between Second Landau Level Orbitals in Double Bilayer Graphene

Here, we present Coulomb-drag measurements on a heterostructure comprising two Bernal-stacked bilayer graphene (BLG) sheets separated by a 2.5 nm hexagonal boron nitride (hBN) spacer in the quantum Hall (QH) regime. Using top and bottom gate control, together with an interlayer bias, we independently tune the two BLG layers into either the lowest (N = 0) or second (N = 1) Landau level (LL) orbital and probe their interlayer QH states. When both layers occupy the N = 0 orbital, we observe both interlayer exciton condensates (ECs) at integer total filling and interlayer fractional QH states, echoing the results in double monolayer graphene. In contrast to previous studies, however, when both BLG layers occupy the N = 1 orbital, we also observe quantized drag signals, signifying an interlayer exciton condensate formed between the second LLs. By tuning the layer degree of freedom, we find that this N = 1 EC state arises only when the N = 1 wave function in each BLG is polarized toward the hBN interface to maximize the interlayer Coulomb interaction.

Electrical conductivity↗

DOE OSTI2019

Magnetic domain dynamics in an insulating quantum ferromagnet

The statistics and form of avalanches in a driven system reveal the nature of the underlying energy landscape and dynamics. In conventional metallic ferromagnets, eddy-current back action can dominate the dynamics. Here, we study Barkhausen noise in Li(Ho,Y)F 4 , an insulating Ising ferromagnet that cannot sustain eddy currents. In addition, for large avalanches at temperatures approaching the Curie point, we find a symmetric response free of drag effects. In the low-temperature limit, drag effects contribute to the dynamics, which we link to enhanced pinning from local random fields that are enabled by the microscopic dipole-coupled Hamiltonian (the Ising model in transverse field).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

How the Universe postpones the evaporation and curtails the quantum spreading of black holes

Black holes are expected to evaporate through the process of Hawking radiation. This process is expected to cause the uncertainty in a black hole’s position to grow to ~ $M^4/M^3_{pl}$ over the course of it’s lifetime, even as its momentum spreads only by ~ $M_{Pl}$ For the black holes that have been observed, which have M ≥ M ⊙ this greatly exceeds the Hubble volume. However, the decay of black holes and their quantum spreading, are delayed in the Universe while the influx of energy into the black holes exceeds their Hawking luminosity. We show that for these M ≥ M ⊙ black holes, their decay outside galaxies and clusters is prevented far longer than it takes the black holes to be dragged well beyond the Hubble horizon, where their eventual decay occurs away from the prying eyes of any observer who has not hitched a ride with them. Meanwhile, black holes in an observer’s galaxy or cluster are themselves prevented from decaying long past the extinction of the last stars, and at least until their galaxy/cluster is swept clean of dark matter, in $\gg$ 10 25 y. Even then, if the black holes become unbound, then they are dragged beyond the Hubble radius before undergoing significant decay; if not, they remain in bound orbits, spreading at most over a subvolume of the galaxy/cluster, and long localized by scatterings to much smaller volumes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Theoretical study on the massively augmented electro-osmotic water transport in polyelectrolyte brush functionalized nanoslits

In this paper, we demonstrate that functionalizing nanoslits with pH-responsive polyelectrolyte (PE) brushes can lead to extremely fast electroosmotic (EOS) water transport, where the maximum centreline velocity and the volume ow rate can be an order of magnitude larger than these quantities in identically charged brush-free nanochannels for a wide range of system parameters. Such an enhancement is most remarkable given that the brushes have been known to retard the transport by imparting additional drag on the fluid flow. We argue that this enhancement stems from the localization of the charge density of the brush-induced electric double layer (and hence the EOS body force) away from the nanochannel wall (or the location of the wall-induced drag force). This ensures a much larger impact of the EOS body force triggering such fast water transport. Lastly, the calculated flux values for the present brush-grafted nanochannels are found to be signi ficantly larger than those for a wide range of nanofluidic membranes and channels, suggesting that the brush functionalization can be considered as a mechanism for enabling such superfast nanofluidic transport.

42 ENGINEERING↗

DOE OSTI2020

Enhanced tuneable rotatory power in a rotating plasma

We report that the gyrotropic properties of a rotating magnetized plasma are derived analytically. Mechanical rotation leads to a new cutoff for wave propagation along the magnetic field, and polarization rotation above this cutoff is the sum of the classical magneto-optical Faraday effect and the mechanico-optical polarization drag. Exploiting the very large effective group index near the cutoff, we expose here that polarization drag can be 10 4 larger than Faraday rotation at GHz frequency. The rotation leads to weak absorption while allowing direct frequency control, demonstrating the unique potential of rotating plasmas for nonreciprocal elements. The very large rotation frequency of a dense non-neutral plasma could enable unprecedented gyrotropy in the THz regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Characterizing Minimum Admissible Separation Distances in Heavy Duty Vehicle Platoons

Heavy duty (HD) vehicles platooning under highway operating conditions, have been projected to provide significant fuel economy (FE) gains based on aerodynamic drag reduction of the platooning vehicles. Realizing these benefits under real-world operating conditions have several challenges. This paper (the second as part of series) quantifies the minimum admissible separation distance between platooning vehicles in HD commercial vehicle (CV) applications. This distance must be set up such that no collision occurs with leading vehicles under all deceleration conditions. In this paper key interactions to characterize this distance are demonstrated. These include, road surface conditions, grade, inter-vehicle communication delays, air-brake lag time, vehicle speed, weight, and aerodynamic drag. It is seen that characterizing the vehicle and road dynamics will be critical in establishing required deceleration or stopping distances for each vehicle in a platoon. The results of the study also necessitate vehicle control systems that can dynamically adjust not only the separation distances but also the braking reference points based on road conditions. The narrative provides application decision personnel with a well-defined set of admissible separation distance impact factors to consider as part of their architecture selection process.

Sujan, Vivek↗

DOE OSTI2022

Modeling of Cube Array Roughness: RANS, Large Eddy Simulation, and Direct Numerical Simulation

Abstract Flow over arrays of cubes is an extensively studied model problem for rough wall turbulent boundary layers. While considerable research has been performed in computationally investigating these topologies using direct numerical simulation (DNS) and large eddy simulation (LES), the ability of sublayer-resolved Reynolds-averaged Navier–Stokes (RANS) to predict the bulk flow phenomena of these systems is relatively unexplored, especially at low and high packing densities. Here, RANS simulations are conducted on six different packing densities of cubes in aligned and staggered configurations. The packing densities investigated span from what would classically be defined as isolated, up to those in the d-type roughness regime, filling in the gap in the present literature. Three different sublayer-resolved turbulence closure models were tested for each case: a low Reynolds number k–ϵ model, the Menter k–ω SST model, and a full Reynolds stress model. Comparisons of the velocity fields, secondary flow features, and drag coefficients are made between the RANS results and existing LES and DNS results. There is a significant degree of variability in the performance of the various RANS models across all comparison metrics. However, the Reynolds stress model demonstrated the best accuracy in terms of the mean velocity profile as well as drag partition across the range of packing densities.

Engineering↗

DOE OSTI2022

Effect of the J-Shaped Wind Turbine Airfoil Opening Ratio and Thickness on the Performance of Symmetrical Airfoils

Abstract In this work, the effect of the inner opening ratio on the J-shaped airfoils aerodynamic performance was studied and documented for symmetrical airfoils. Three different airfoil thicknesses were investigated: small (NACA0008), medium (NACA0015), and large (NACA0024). For each airfoil thickness, effects of three inner opening ratios were analyzed: one-third, one-half, and two-thirds. The performance of each opening ratio was compared with the performance of the solid airfoil “zero opening ratio” for different angles of attack between 5 deg and 20 deg. All designs were simulated using the computational fluid dynamics (CFD) technology against experimental results for solid NACA4412 airfoil in the University of Wisconsin-Milwaukee (UWM) wind tunnel facility and other published experimental data. It was found that large eddy simulation yields accurate solutions with a smaller number of mesh cells compared to the k–ω turbulence model but with much longer computational time. The lift-to-drag ratio for all studied airfoils has a maximum value for solid airfoils compared to those equipped with openings. For airfoils equipped with 00.00% opening ratio “solid,” NACA0015 airfoil has the maximum lift-to-drag ratio. Furthermore, it was found that NACA0008 equipped with a 33.33% opening ratio has the best performance of all studied J-shaped airfoils.

Energy & Fuels↗

DOE OSTI2022

Novel Tube Design for Superheater Heat Exchanger Enabled Via Additive Manufacturing

Superheater tubes are critical boiler components that operate at relatively higher temperatures and pressure. Amongst the primary concerns for these tubes is the deposition of ash particles on the tube surface, leading to the reduced thickness of the tube due to material corrosion, consequently causing early creep failure of the component. In this research, a novel tube design has been proposed which resembles a teardrop or ogive shape to reduce the drag and concurrently improve the creep life of the superheater tubes. To administer the practicality of novel tubes, metal additive manufacturing (AM), for instance, laser-powder bed fusion (L-PBF), has been proposed. These unconventional designs were assessed and compared with the baseline circular tube design for mechanical design requirements (hoop stress and creep life) and the particle and flue gas flow characteristics around the differently shaped tubes. A thermomechanical finite element (FE) analysis was performed for hoop stress calculations. This study also emphasizes on effect of circumferential thermal variation on hoop stress distribution in tubes. Therefore, a detailed two-dimensional (2D) thermal simulation has been performed to report the circumferential thermal variation on the tube. A computational fluid dynamics (CFD) analysis coupled with particle tracing was performed for gas flow visualization and particle tracing around the proposed shapes and baseline circular-shaped tube design. The Schlieren optic setup was built and leveraged for qualitative validation of the proposed design. The complete design methodology established in the paper shows teardrop-shaped tubes better in terms of drag and creep life in contrast to the circular-shaped tube.

42 ENGINEERING↗

DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗

DOE OSTI2023

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗

DOE OSTI2025

An infrasound source analysis of the OSIRIS-REx sample return capsule hypersonic re-entry

The OSIRIS-REx sample return capsule's hypersonic re-entry into the atmosphere is a rare opportunity to test a variety of sonic boom source models since the projectile dimensions are well characterized. While the as-flown flight path is unknown, the predicted flight path enables a rough approximation of the source Mach number and location. Six infrasound microphones deployed in the boom carpet along the predicted flight path recorded impulsive signals from the OSIRIS-REx re-entry. Using a suite of atmosphere profiles and the geometric acoustics approximation, we estimate locations with uncertainty estimates along the flight path from which the signals were emitted. Acoustic overpressure and signal duration predictions from Whitham's far field theory, Carlson's simplified sonic boom prediction method, and a drag-dominated hypersonic model are analyzed with uncertainty estimates from the location estimate. While the Carlson simplified sonic boom prediction method could be accurate, our preference is for the drag-dominated source model. Using this source model with an inviscid Burgers's equation solver for propagation, we obtained an excellent match to the recorded data. In conclusion, these results will help better inform future sample return capsule re-entry observation campaigns as well as contribute to a better understanding of high altitude infrasonic sources.

58 GEOSCIENCES↗

DOE OSTI2021

Hyperdiffusion of dust particles in a turbulent tokamak plasma

The effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad-hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number $St\propto a_d/n_0$, with $a_d$ the dust radius and $n_0$ the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories $\Delta r^2$ is shown to depend on the ratio $Ku/St$, with $Ku\propto u_{rms}$ the Kubo number and $u_{rms}$ the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting $\Delta r^2\propto (Ku/St)^2t^3$, or $\Delta r^2\propto (u_{rms}n_0/a_d)^2t^3$. As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible, and drives the transition from a hyperdiffusive regime towards a ballistic one.

dust↗