Search NASA⌕ Search

SEARCH · Search NASA

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

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

At least 91 records · Page 5

An artificial viscosity method for the design of supercritical airfoils

A numerical technique is presented for the design of two-dimensional supercritical wing sections with low wave drag. The method is a design mode of the analysis code H which gives excellent agreement with experimental results and is widely used in the aircraft industry. Topics covered include the partial differential equations of transonic flow, the computational procedure and results; the design procedure; a convergence theorem; and description of the code.

Mcfadden, G. B.↗

Modeling the dust size distribution in comets with dust fragmentation

A hydrodynamic model was developed of a spherically symmetric dusty gas flow in a cometary atmosphere assuming a single fluid, inviscid, perfect gas. The hydrodynamics for gas and dust, which involves the gas drag force (momentum transfer), heat exchange between gas and dust, photodissociation energy for H2O gas, and radiative heating and cooling terms for dust particles, are solved using the Gear method for stiff, coupled differential equations. Calculations were done with a dust size distribution for radii alpha = 0.01 micron to 10 cm with densities variable with the size. A nucleus size of 4.0 km radius with a density of 0.5 g/cu cm and a total dust-to-gas mass ratio chi = 1 was adopted. There are indications from in situ observations that dust particle fragment into smaller ones. Fragmentation of dust particles was incorporated into the model. This is done by adding source and sink terms in the continuity equations for the dust. Lifetimes for the decay of dust particles were assumed as a function of particle size. It is also assumed that dust particles always fragment only into the next smaller size.

Konno, Ichishiro↗

Self-gravity and dissipation in polar rings

Studies of inclined rings inside galaxy potentials have mostly considered the influence of self-gravity and viscous dissipation separately. In this study, we construct models of highly inclined ('polar') rings in an external potential including both self-gravity and dissipation due to a drag force. We do not include pressure forces and thus ignore shock heating that dominates the evolution of gaseous rings inside strongly nonspherical potentials. We adopt an oblate spheroidal scale-free logarithmic potential with axis ratio q = 0.85 and an initial inclination of 80 deg for the self-gravitating rings. We find that stellar (dissipationless) rings suffer from mass loss during their evolution. Mass loss also drives a secular change of the mean inclination toward the poles of the potential. As much as half of the ring mass escapes in the process and forms an inner and an outer shell of precessing orbits. If the remaining mass is more than approximately 0.02 of the enclosed galaxy mass, rings remain bound and do not fall apart from differential precession. The rings precess at a constant rate for more than a precession period tau(sub p) finding the configuration predicted by Sparke in 1986 which warps at larger radii toward the poles of the potential. We model shear viscosity with a velocity-dependent drag force and find that nuclear inflow dominates over self-gravity if the characteristic viscous inflow time scale tau(sub vi) is shorter than approximately 25(tau(sub p)). Rings with (tau(sub vi))/(tau(sub p)) less than or approximately equal to 25 collapse toward the nucleus of the potential within one precession period independent of the amount of self-gravity. Our results imply that stars and gas in real polar rings exhibit markedly different dynamical evolutions.

Dubinski, John↗

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗

Numerical calculation of the three dimensional transonic flow over a yawed wing.

Results are presented of calculations of the three dimensional steady transonic flow over a finite yawed wing. The full potential flow equation is solved in a transformed coordinate system which permits the boundary conditions to be satisfied exactly. The correct differential properties are enforced by rotating the difference scheme to conform with the flow direction, and fast convergence is assured by simulating a time dependent equation designed to settle quickly to a steady state. Computed lift drag ratios are consistent with the results of wind tunnel tests of a yawed wing conducted by R. T. Jones (1972).-

Jameson, A.↗

Gradiometry coexperiments to the gravity probe B and step missions

The Gravity Probe-B (GP-B) spacecraft, designed to test predictions of general relativity, will fly in the mid 1990s. It will carry four electrostatically suspended gyroscopes in a cryogenic environment and will have a drag-free control system to minimize disturbances on the gyroscopes. The Stanford Test of Equivalence Principle (STEP) spacecraft, to fly later, will carry a set of test masses under very similar conditions. The possibility of using differential measurements of the GP-B gyroscopes suspension forces and the STEP tests mass displacement readout to form single-axis gravity gradiometers is explored. It is shown that the noise in the suspension systems is sufficiently small in the relevant frequency range, and that enough information is collected to compensate for the spacecrafts' attitude motion. Finally, using Breakwell's flat-earth approximation, these experiments are compared to other geodesy experiments and predict the contribution they can make to the knowledge of the Earth's geopotential.

Tapley, M.↗

Parameter estimation with maximal updated densities

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. Further, the solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

42 ENGINEERING↗

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗

Computer-automated opponent for manned air-to-air combat simulations

Two versions of a real-time digital-computer program that operates a fighter airplane interactively against a human pilot in simulated air combat were evaluated. They function by replacing one of two pilots in the Langley differential maneuvering simulator. Both versions make maneuvering decisions from identical information and logic; they differ essentially in the aerodynamic models that they control. One is very complete, but the other is much simpler, primarily characterizing the airplane's performance (lift, drag, and thrust). Both models competed extremely well against highly trained U.S. fighter pilots.

Hankins, W. W., III↗

Theory of idealized two-dimensional ballute in Newtonian hypersonic flow

A differential equation governing the geometry of a two-dimensional ballute in hypersonic flow and its constraining boundary conditions are derived under idealized assumptions. By solving these equations, the shape of the ballute is determined over a range of conditions. Lift, drag, pitching moment, and the allowed limit of center-of-gravity location for stability (meta-center) are then calculated using Newtonian hypersonic approximation. It is shown that the meta-center occurs near the forward end because of compliance of the ballute membrane to the shock layer pressures, especially at low free-stream densities. In order for the vehicle employing the ballute to be stable at all densities, the center of gravity must be within approximately the forward 20 percent of overall length of the vehicle. However, typical flight trajectories of an aeroassisted orbital transfer vehicle employing the ballute for aerobraking show that the vehicle may be able to complete its atmospheric flight without tumbling provided that the center of gravity is located within the forward 43 percent of the vehicle length because of the relatively short duration of flight through the destabilizing low-density regime.

Park, C.↗

Radiation forces on small particles in the solar system

Solar radiation forces on small particles in the solar system are examined, and the resulting orbital evolution of interplanetary and circumplanetary dust is considered. An expression is derived for the effects of radiation pressure and Poynting-Robertson drag on small, spherical particles using the energy and momentum transformation laws of special relativity, and numerical examples are presented to illustrate that radiation pressure and Poynting-Robertson drag are only important for particles within a narrow size range. The orbital consequences of these radiation forces are considered both for heliocentric and planetocentric orbiting particles, and the coupling between particle sizes and dynamics is discussed. A qualitative derivation is presented for the differential Doppler effect, which is due to the differential Doppler shifting of radiation from approaching and receding solar hemispheres, and the Yarkovsky effect, which is important for rotating meter-to kilometer-sized particles, is briefly described.

Burns, J. A.↗

GRACE Mission Design: Impact of Uncertainties in Disturbance Environment and Satellite Force Models

The Gravity Recovery and Climate Experiment (GRACE) primary mission will be performed by making measurements of the inter-satellite range change between two co-planar, low altitude, near-polar orbiting satellites. Understanding the uncertainties in the disturbance environment, particularly the aerodynamic drag and torques, is critical in several mission areas. These include an accurate estimate of the spacecraft orbital lifetime, evaluation of spacecraft attitude control requirements, and estimation of the orbital maintenance maneuver frequency necessitated by differences in the drag forces acting on both satellites. The FREEMOL simulation software has been developed and utilized to analyze and suggest design modifications to the GRACE spacecraft. Aerodynamic accommodation bounding analyses were performed and worst-case envelopes were obtained for the aerodynamic torques and the differential ballistic coefficients between the leading and trailing GRACE spacecraft. These analyses demonstrate how spacecraft aerodynamic design and analysis can benefit from a better understanding of spacecraft surface accommodation properties, and the implications for mission design constraints such as formation spacing control.

Mazanek, Daniel D.↗

Aerothermal Analysis and Design of the Gravity Recovery and Climate Experiment (GRACE) Spacecraft

The Gravity Recovery and Climate Experiment (GRACE) primary mission will be performed by making measurements of the inter-satellite range change between two co-planar, low altitude near-polar orbiting satellites. Understanding the uncertainties in the disturbance environment, particularly the aerodynamic drag and torques, is critical in several mission areas. These include an accurate estimate of the spacecraft orbital lifetime, evaluation of spacecraft attitude control requirements, and estimation of the orbital maintenance maneuver frequency necessitated by differences in the drag forces acting on both satellites. The FREEMOL simulation software has been developed and utilized to analyze and suggest design modifications to the GRACE spacecraft. Aerodynamic accommodation bounding analyses were performed and worst-case envelopes were obtained for the aerodynamic torques and the differential ballistic coefficients between the leading and trailing GRACE spacecraft. These analyses demonstrate how spacecraft aerodynamic design and analysis can benefit from a better understanding of spacecraft surface accommodation properties, and the implications for mission design constraints such as formation spacing control.

Mazanek, Daniel D.↗

Analysis of Lageos' altitude decrease

The paper treats the inverse problem of celestial mechanics which consists of determining the force field or potential from given or observed orbit(s). From the observational information, according to which Lageos loses approximately 1 mm altitude per day, a linear partial differential equation is formulated. The solution of this equation gives the field responsible for the above-mentioned, as yet unexplained, small but well established secular decrease in the semi-major axis. Note that the altitude-loss is not due to air-drag because of the very high altitude of this satellite.

Szebehely, V.↗

Stability criteria for accelerating and decelerating aircraft.

Derivation of the linearized differential equation for the longitudinal oscillations of an aircraft for the general case when the atmospheric density, thrust force, and velocity are varying. Exact solutions and stability criteria are presented for several different conditions. For the zero thrust coasting aircraft that is being decelerated by its drag force, new stability criteria are derived for both ascending and descending flight. A critical altitude is found for ascending aircraft that are attempting to coast out of the atmosphere. Above this critical altitude the oscillations cease, and there is a monotonic increase in the angle of attack. Exact solutions for the oscillations of coasting aircraft are presented in terms of the confluent hypergeometric functions. It is shown that previous attempts to predict a critical altitude for ascending vehicles were in error because they did not include the second type, or logarithmic solution, of the confluent hypergeometric function.

Laitone, E. V.↗

Viscosity Measurements at High Pressures: A Critical Appraisal of Corrections to Stokes' Law

Abstract Fluids and melts in planetary interiors significantly influence geodynamic processes from volcanism to global‐scale differentiation. The roles of these geofluids depend on their viscosities ( η ). Constraining geofluid η at relevant pressures and temperatures relies on laboratory‐based measurements and is most widely done using Stokes' Law viscometry with falling spheres. Yet small sample chambers required by high‐pressure experiments introduce significant drag on the spheres. Several correction schemes are available for Stokes' Law but there is no consensus on the best scheme(s) for high‐pressure experiments. We completed high‐pressure experiments to test the effects of (a) the relative size of the sphere diameter to the chamber diameter and (b) the top and bottom of the chamber, that is, the ends, on the sphere velocities. We examined the influence of current correction schemes on the estimated viscosity using Monte Carlo simulations. We also compared previous viscometry work on various geofluids in different experimental setups/geometries. We find the common schemes for Stokes' Law produce statistically distinct values of η . When inertia of the sphere is negligible, the most appropriate scheme may be the Faxén correction for the chamber walls. Correction for drag due to the chamber ends depends on the precision in the sinking distance and may be ineffective with decreasing sphere size. Combining the wall and end corrections may overcorrect η . We also suggest the uncertainty in η is best captured by the correction rather than propagated errors from experimental parameters. We develop an overlying view of Stokes' Law viscometry at high pressures.

Geochemistry & Geophysics↗

Permeability of a New Parachute Fabric – Measurements, Modeling, and Application

In response to a NASA request, Heathcoat Fabrics Limited has woven a new parachute fabric (Custom Design 1 G-60315-1800-Q01). This fabric was tested to obtain its permeability in air (i.e., flow-through volume of air per area per time) over a range of differential pressures from 0.146 to 25 psf (7 to 1197 Pa). The fabric met its specification permeability of 60 to 100 ft3/ft2/min (30.5 to 50.8 cm3/cm2/s) at the U.S. standard differential pressure of 0.5 inch of water (2.60 psf, 124 Pa). The permeability results were transformed into an effective porosity model for use in calculations related to the total porosity of parachutes. The tested fabric is being considered for use in parachutes for future missions to Mars. Calculations of drag coefficient were performed for two geometrically identical parachutes using either the new fabric or fabric woven to Parachute Industry Specification PIA-C-7020D Type I (Mars Science Laboratory Disk-Gap-Band parachute operating on Mars at a Mach number of 0.41). These calculations indicate essentially no difference in the drag coefficient between the two parachutes.

Cruz, Juan R.↗

The singly averaged differential equations of satellite motion for e greater than or equal to 0 and less than 1

The singly-averaged differential equations of motion of a satellite are developed in terms of parameters valid for all eccentricities less than one. The perturbations included in the acceleration model are due to an aspherical central planet (zonal harmonics up to degree 20 and resonant harmonics up to degree and order 20), atmospheric drag for a time-varying atmosphere, third-body gravity (the sun and moon for an earth satellite), solar radiation pressure with shadowing, and impulsive maneuvers. Analytic averaging is used to remove short-period terms due to the aspherical central planet and third-body gravity. Numerical averaging is used to remove short-period terms due to atmospheric drag and solar radiation pressure.

Dallas, S. S.↗