SEARCH · Search NASA

Results for “Particle 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.

446 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Emergent viscous hydrodynamics from a single quantum particle

We investigate an explicit example of how spatial decoherence can lead to hydrodynamic behavior in the late-time, long-wavelength regime of open quantum systems. We focus on the case of a single nonrelativistic quantum particle linearly coupled to a thermal bath of noninteracting harmonic oscillators at temperature T , a la Caldeira and Leggett. Taking advantage of decoherence in the position representation, we expand the reduced density matrix in powers of the off-diagonal spatial components, so that high-order terms are suppressed at late times. Truncating the resulting power series at second order leads to a set of dissipative transient hydrodynamic equations similar to the nonrelativistic limit of equations widely used in simulations of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Transport coefficients are directly determined by the damping constant γ , which quantifies the influence of the environment. The asymptotic limit of our hydrodynamic equations reduces to the celebrated Navier-Stokes equations for a compressible fluid in the presence of a drag force. Furthermore, our results shed new light on the onset of hydrodynamic behavior in open quantum systems where a system with few degrees of freedom is coupled to a large thermal environment.

Hydrodynamics↗

NASA NTRS2005

Forecasting the Impact of an 1859-calibre Superstorm on Satellite Resources

We have assembled a database of operational satellites in orbit as of 2004, and have developed a series of simple models to assess the economic impacts to this resource caused by various scenarios of superstorm events possible during the next sunspot cycle between 2010 and 2014. Despite the apparent robustness of our satellite assets against the kinds of storms we have encountered during the satellite era, our models suggest a potential economic loss exceeding $10(exp 11) for satellite replacement and lost profitability caused by a once a century single storm similar to the 1859 superstorm. From a combination of power system and attitude control system (the most vulnerable) failures, we estimate that 80 satellites (LEO, MEO, GEO) may be disabled as a consequence of a superstorm event. Additional consequences may include the failure of many of the GPS, GLONASS and Galileo satellite systems in MEO. Approximately 98 LEO satellites that normally would not have re-entered for many decades, may prematurely de-orbit in ca 2021 as a result of the temporarily increased atmospheric drag caused by the superstorm event occurring in 2012. The $10(exp 11) International Space Station may lose at least 15 kilometers of altitude, placing it in critical need for re-boosting by an amount that is potentially outside the range of typical Space Shuttle operations during the previous solar maximum in ca 2000, and at a time when NASA plans to decommission the Space Shuttle. Several LEO satellites will unexpectedly be placed on orbits that enter the ISS zone of avoidance, requiring some action by ground personnel and ISS astronauts to avoid close encounters. Radiation effects on astronauts have also been considered and could include a range of possibilities from acute radiation sickness for astronauts inside spacecraft, to near-lethal doses during EVAs. The specifics depends very sensitively on the spectral hardness of the accompanying SPE event. Currently, the ability to forecast extreme particle events and coronal mass ejections, or predict their fluences and geo-severity in the 24-hrs prior to the event, appears to be no better than 50/50. If the events of the 1859 superstorm serve as a guide, the scope of a contemporary superstorm will most certainly be an awesome event, but one that the vast majority of our other satellite resources may reasonably be expected to survive.

Odenwald, Sten↗

NASA NTRS2011

Heliocentric Distance of Coronal Mass Ejections at the Time of Energetic Particle Release: Revisiting the Ground Level Enhancement Events of Solar Cycle 23

Using the kinematics of coronal mass ejections (CMEs), onset time of soft X-ray flares, and the finite size of the pre-eruption CME structure, we derive the heliocentric distan~e at which the energetic particles during the ground level enhancement (GLE) events of Solar Cycle 23. We find that the GLE particles are released when the CMEs reach an average heliocentric distance of approx.3.25 solar radii (Rs). From this we infer that the shocks accelerating the particles are located at similar heights. Type II radio burst observations indicate that the CMEs are at much lower distances (average approx.1.4 Rs) when the CME-driven shock first forms. The shock seems to travel approx.1.8 Rs over a period of approox.30 min on the average before releasing the GLE particles. In deriving these results, we made three assumptions that have observational support: (i) the CME lift off occurs from an initial distance of about 1.25 Rs; (ii) the flare onset and CME onset are one and the same because these are two different manifestations of the same eruption; and (iii) the CME has positive acceleration from the onset to the first appearance in the coronagraphic field of view (2.5 to 6 Rs). Observations of coronal cavities in eclipse pictures and in coronagraphic images justify the assumption (i). The close relationship between the flare reconnection magnetic flux and the azimuthal flux of interplanetary magnetic clouds justify assumption (ii) consistent with the standard model (CSHKP) of solar eruption. Coronagraphic observations made close to the solar surface indicate a large positive acceleration of CMEs to a heliocentric distance of approx.3 Rs before they start slowing down due to the drag force. The inferred acceleration (approx.1.5 km/s/s) is consistent with reported values in the literature.

Gopalswamy, Natchimuthuk↗

NASA NTRS2017

Wake Structure Induced by Seal Whisker: Effects of Angle of Incidence

Seal whiskers have been found to produce unique wake flow structures that minimize self-induced vibration and reduce drag. The cause of these wake features are due to the peculiar three-dimensional morphology of the whisker surface. The whisker can be defined as an elliptical cross section with variation in the major and minor axis of the ellipse along the length of the whisker as well as rotation of the elliptical plane with respect to the whisker axis, angle of incidence. While the effects of several dominant parameters of the whisker morphology have been studied, the effect of the rotation of the elliptical plane has not been well understood. This paper will examine the influence of the angle of incidence on the wake flow structure through series of water channel studies. Several models of whisker-like geometries will be tested which isolate the rotation angle as the only variation between models. The wake structure behind each seal whisker model will be measured through particle image velocimetry (PIV). The different models’ wake structures will be compared identifying the effect of angle of incidence on the wake structure. Angle of incidence was found to influence the wake structure through reorganization of velocity patterns, reduction of recovery length and modification of magnitude of Tu. These results help provide a more complete understanding of the seal whisker geometry relationship to wake structure and can provide insight into design practices for application of whisker geometery to various engineering problems.

vortex induced vibration↗

DOE OSTI2020

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

"Hall" transport of liquid crystal solitons in Couette flow

Topology establishes a unifying framework for a diverse range of scientific areas including particle physics, cosmology, and condensed matter physics. One of the most fascinating manifestations of topology in the context of condensed matter is the topological Hall effect, and its relative: the Skyrmion Hall effect. Skyrmions are stable vortex-like spin configurations in certain chiral magnets, and when subject to external electric currents can drift in the transverse direction to the current. These quasi-particles are characterised by a conserved topological charge which in the Skyrmion Hall effect plays the role of electric charges in the ordinary Hall effect. Recently, it has been shown that liquid crystals endowed with chiral properties serve as an ideal testbed for the fundamental investigation of topological solitons, including their two- and three-dimensional realisations. Here, we show experimentally and numerically that three-dimensional solitons aka "torons" exhibit a Hall-like effect when driven by shear flows: the torons are deflected in the direction perpendicular to the shear plane. The experimental results are rationalised in terms of the dynamic Ericksen-Leslie equations, which predict the emergence of the transverse component of the net mass flow, the magnitude of which scales as the 3rd power of the shear rate. The perturbation analysis highlights an interplay of the viscous and chiral elastic torques as the mechanism for the emergence of net transverse currents. Numerical simulations demonstrate, however, that torons are not merely dragged by the flow but move with their own transverse speed, much larger than the average flow velocity in the transverse direction. Our findings may enable responsive microfluidic applications relying on soft topological solitons.

FOS: Physical sciences↗

NASA NTRS1999

Computer Modeling of Sand Transport on Mars Using a Compart-Mentalized Fluids Algorithm (CFA)

It has been postulated that aeolian transport on Mars may be significantly different from that on Earth. From laboratory experiments simulating martian grain transport [2], it has been observed that (saltating) grains striking the bed can cause hundreds of secondary reptation trajectories when impact occurs at speeds postulated for Mars. Some of the ballistically induced trajectories "die ouf' and effectively join the ranks on the creep population that is merely nudged along by impact. Many of the induced reptation trajectories, however, are sufficiently high for the grains to become part of the saltation load (it is irrelevant to the boundary layer how a grain attained its initial lift force). When these grains, in turn, strike the surface, they too are capable of inducing more reptating grains. This cascading effect has been discussed in connection with terrestrial aeolian transport in an attempt to dispel the notion that sand motion is divisible only into creep and saltation loads. On Earth, only a few grains are splashed by impact. On Mars, it may be hundreds. We developed a computer model to address this phenomenon because there are some important ramifications: First, this ratio may mean that martian aeolian transport is dominated by reptation flux rather than saltation. On Earth, the flux would be a roughly balanced mixture between reptation/creep and saltation. On Venus, there would be no transport other than by saltation. In other words, an understanding of planetary aeolian processes may not be necessarily understood by extrapolating from the "Earth case", with only gravity and atmospheric density/viscosity being considered as variables. Second, the reptation flux on Mars may be self sustaining, so that little input is required by the wind once transport has been initiated. The number of grains saturating the boundary layer near the bed may mean that average grain speed on Mars might conceivably be less than that on Earth. This would say much for models of sand comminution on Mars. A multiple-grain transport model using just the equations of grain motion describing lift and drag is impossible to develop owing to stochastic effects --the very effects we wish to model. Also, unless we were to employ supercomputing techniques and extremely complex computer codes that could deal with millions of grains simultaneously, it would also be difficult to model grain transport if we attempted to consider every grain in motion. No existing computer models were found that satisfactorily used the equations of motion to arrive at transport flux numbers for the different populations of saltation and reptation. Modeling all the grains in a transport system was an intractable problem within our resources, and thus we developed what we believe to be a new modeling approach to simulating grain transport. The CFA deals with grain populations, but considers them to belong to various compartmentalized fluid units in the boundary layer. In this way, the model circumvents the multigrain problem by dealing primarily with the consequences of grain transport --momentum transfer between air and grains, which is the physical essence of a dynamic grain-fluid mixture. We thus chose to model the aeolian transport process as a superposition of fluids. These fluids include the air as well as particle populations of various properties. The prime property distinguishing these fluids is upward and downward grain motion. In a normal saltation trajectory, a grain's downwind velocity increases with time, so a rising grain will have a smaller downwind velocity than a failing grain. Because of this disparity in rising and falling grain proper-ties, it seemed appropriate to track these as two separate grain populations within the same physical space. The air itself can be considered a separate fluid superimposed within and interacting with the various grain-cloud "fluids". Additional informaiton is contained in the original.

Marshall, J.↗

DOE OSTI2021

Investigation of theoretical scaling laws using large eddy simulations for airborne spreading of viral contagion from sneezing and coughing

Using a set of large eddy point-particle simulations, we explore the fluid dynamics of an ejected puff resulting from a cough/sneeze. The ejection contains over 61,000 potentially virus-laden droplets at an injection Reynolds number of about 46,000, comparable to an actual cough/ sneeze. We observe that global puff properties, such as centroid, puff volume, momentum, and buoyancy vary little across realizations. Other properties, such as maximum extent, shape, and edge velocity of the puff, may exhibit substantial variation. In many realizations, a portion of the puff splits off and advances along a random direction, while keeping airborne droplet nuclei afloat. This peeled-off portion provides a mechanism for virus-laden droplets to travel over large distances in a short amount of time. We also observe that the vast majority of droplets remain suspended within the puff after all liquid has evaporated. The main objectives of the study are to (i) evaluate assumptions of Balachandar’s et al. theory [Int. J. Multiphase Flow 132, 103439 (2020)], which include buoyancy effects, shape of the puff, and droplet evaporation rate, (ii) obtain values of closure parameters, which include location and time of the virtual origin, and puff entrainment and drag coefficients, and (iii) evaluate the accuracy of the theory in predicting the shape, size, and location of the puff, as well as droplet number density long after ejection. The theory adequately predicts global puff properties including size, velocity, and distance traveled, the largest size of droplets that exit the puff due to settling, and the droplet size distribution within the puff long after ejection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗

DOE OSTI2024

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

SSAPy - Space Situational Awareness for Python

SSAPy is a fast and flexible orbit modeling and analysis tool for orbits spanning from low-Earth into the cislunar regime. Orbits can be flexibly specified from common input formats such as Keplerian elements or two-line element (TLE) data files. SSAPy allows users to model satellites and specify parameters such as satellite area, mass, and drag coefficients. SSAPy includes a customizable force-propagation with a range of Earth, Lunar, radiation, atmospheric, and maneuvering models. SSAPy makes use of various community integration methods and can calculate time-evolved orbital quantities, including satellite magnitudes and state vectors. Users can specify various space- and ground-based observation models with support for multiple coordinate and reference frames. SSAPy also supports orbit analysis and propagation methods such as multiple hypothesis tracking and has built-in uncertainty quantification. The majority of SSAPy’s methods are vectorized and parallelizable, allowing for effective use of high-performance computer (HPC) systems. Finally, SSAPy has plotting functionality, allowing users to visualize orbits and trajectories. Examples are shown in Figure 1 and Figure 2.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Cosmic recombination in the presence of primordial magnetic fields

Primordial magnetic fields (PMFs) may explain observations of magnetic fields on extragalactic scales. They are most cleanly constrained by measurements of cosmic microwave background radiation (CMB) anisotropies. Their effects on cosmic recombination may even be at the heart of the resolution of the Hubble tension. We present the most detailed analysis of the effects of PMFs on cosmic recombination to date. To this end we extend the public magneto-hydrodynamic code ENZO with a new cosmic recombination routine, Monte-Carlo simulations of Lyman-α photon transport, and a Compton drag term in the baryon momentum equation. The resulting code allows us, for the first time, to realistically predict the impact of PMFs on the cosmic ionization history and the clumping of baryons during cosmic recombination. Our results identify the importance of mixing of Lyman-α photons between overdense- and underdense- regions for small PMF strength. This mixing speeds up recombination beyond the speed-up due to clumping. We also investigate the effects of pecuilar flows on the recombination rate and find it to be small for small PMF strengths. For non-helical PMFs with a Batchelor spectrum we find a surprising dependency of results on ultra-violet magnetic modes. We further show that the increase in the ionization fraction at low redshift by hydrodynamic baryon heating due to PMF dissipation is completely compensated by the faster recombination from baryon clumping. In conclusion, the present study shall serve as a theoretical foundation for a future precise comparison of recombination with PMFs to CMB data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Nonlinear dynamics and dissipation of a curvilinear vortex driven by a strong time-dependent Meissner current

Here, we report numerical simulations of large-amplitude oscillations of a trapped vortex line under a strong ac magnetic field $\textit{H(t) = H}$ sin$\textit{ωt}$ parallel to the surface. The power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was calculated by taking into account the nonlinear vortex line tension, vortex mass, and a nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient $\textit{η(v)}$. We show that the LO decrease of $\textit{η(v)}$ with the vortex velocity $\textit{v}$ can radically change the field dependence of the surface resistance R i ( H ) caused by trapped vortices. At low frequencies, $R_i(H)$ exhibits a conventional increase with $\textit{H}$, but as $\textit{ω}$ increases, the surface resistance becomes a nonmonotonic function of $\textit{H}$ which decreases with $\textit{H}$ at higher fields. The effects of frequency, pin spacing, and the mean-free path $l_i$ on the field dependence of $R_i(H)$ were calculated. It is shown that, as the surface gets dirtier and $l_i$ decreases, the anomalous drop of $R_i(H)$ with $\textit{H}$ shifts to lower fields which can be much smaller than the lower critical magnetic field. Our numerical simulations also show that the LO decrease of $\textit{η(v)}$ with $\textit{v}$ can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along the vortex. Measurements of $R_i(H)$ caused by sparse vortices trapped perpendicular to the surface can offer opportunities to investigate an extreme nonlinear dynamics of vortices driven by strong current densities up to the depairing limit at low temperatures. The behavior of $R_i(H)$, which can be tuned by varying the rf frequency or concentration of nonmagnetic impurities, is not masked by strong heating effects characteristic of dc or pulse transport measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Anomalous resistivity and electron heating by lower hybrid drift waves inside reconnecting current sheets

Inside an electron diffusion region of laboratory reconnection experiments, the quasi-electrostatic lower hybrid drift wave (ES-LHDW) is observed when a significant guide field component is present. Through direct measurement of the anomalous drag term and quasilinear analysis, it is shown that ES-LHDW can account for approximately 20% of the mean reconnection electric field in a case with moderate guide field. This value exceeds the contribution from classical resistivity, which is around 10%. The effects of the Lorentz force term, often neglected for electrostatic waves, are crucial for the observed correlation between electric field and density fluctuations. Anomalous electron heating by the perturbed current and resistivity (2.6 MW/m 3 ) also surpasses the classical Ohmic heating, which is about 2.0 MW/m 3 . For the case with a high guide field, significantly higher local electron temperatures were observed during periods of strong ES-LHDW activity. A statistical analysis further supports electron heating by LHDW, showing a larger increase in electron temperature with a high guide field. Finally, data from the Magnetospheric Multiscale mission provide evidence of Landau damping of ES-LHDW, suggesting that ES-LHDW may contribute to the generation of nonthermal electrons along the direction parallel to the magnetic field.

Magnetic reconnection↗