Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “drag estimation”

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 451 records · Page 25

Coupled Ablation, Heat Conduction, Pyrolysis, Shape Change and Spallation of the Galileo Probe

The Galileo probe enters the atmosphere of Jupiter in December 1995. This paper presents numerical methodology and detailed results of our final pre-impact calculations for the heat shield response. The calculations are performed using a highly modified version of a viscous shock layer code with massive radiation coupled with a surface thermochemical ablation and spallation model and with the transient in-depth thermal response of the charring and ablating heat shield. The flowfield is quasi-steady along the trajectory, but the heat shield thermal response is dynamic. Each surface node of the VSL grid is coupled with a one-dimensional thermal response calculation. The thermal solver includes heat conduction, pyrolysis, and grid movement owing to surface recession. Initial conditions for the heat shield temperature and density were obtained from the high altitude rarefied-flow calculations of Haas and Milos. Galileo probe surface temperature, shape, mass flux, and element flux are all determined as functions of time along the trajectory with spallation varied parametrically. The calculations also estimate the in-depth density and temperature profiles for the heat shield. All this information is required to determine the time-dependent vehicle mass and drag coefficient which are necessary inputs for the atmospheric reconstruction experiment on board the probe.

Milos, Frank S.↗

Interplanetary Propagation of Coronal Mass Ejections

Although more than ten thousand coronal mass ejections (CMEs) are produced during each solar cycle at the Sun, only a small fraction hits the Earth. Only a small fraction of the Earth-directed CMEs ultimately arrive at Earth depending on their interaction with the solar wind and other large-scale structures such as coronal holes and CMEs. The interplanetary propagation is essentially controlled by the drag force because the propelling force and the solar gravity are significant only near the Sun. Combined remote-sensing and in situ observations have helped us estimate the influence of the solar wind on the propagation of CMEs. However, these measurements have severe limitations because the remote-sensed and in-situ observations correspond to different portions of the CME. Attempts to overcome this problem are made in two ways: the first is to model the CME and get the space speed of the CME, which can be compared with the in situ speed. The second method is to use stereoscopic observation so that the remote-sensed and in-situ observations make measurements on the Earth-arriving part of CMEs. The Solar Terrestrial Relations Observatory (STEREO) mission observed several such CMEs, which helped understand the interplanetary evolution of these CMEs and to test earlier model results. This paper discusses some of these issues and updates the CME/shock travel time estimates for a number of CMEs.

Gopalswamy, Nat↗

Simulation Study of a Follow-on Gravity Mission to GRACE

The gravity recovery and climate experiment (GRACE) has been providing monthly estimates of the Earth's time-variable gravity field since its launch in March 2002. The GRACE gravity estimates are used to study temporal mass variations on global and regional scales, which are largely caused by a redistribution of water mass in the Earth system. The accuracy of the GRACE gravity fields are primarily limited by the satellite-to-satellite range-rate measurement noise, accelerometer errors, attitude errors, orbit errors, and temporal aliasing caused by unmodeled high-frequency variations in the gravity signal. Recent work by Ball Aerospace and Technologies Corp., Boulder, CO has resulted in the successful development of an interferometric laser ranging system to specifically address the limitations of the K-band microwave ranging system that provides the satellite-to-satellite measurements for the GRACE mission. Full numerical simulations are performed for several possible configurations of a GRACE Follow-On (GFO) mission to determine if a future satellite gravity recovery mission equipped with a laser ranging system will provide better estimates of time-variable gravity, thus benefiting many areas of Earth systems research. The laser ranging system improves the range-rate measurement precision to approximately 0.6 nm/s as compared to approx. 0.2 micro-seconds for the GRACE K-band microwave ranging instrument. Four different mission scenarios are simulated to investigate the effect of the better instrument at two different altitudes. The first pair of simulated missions is flown at GRACE altitude (approx. 480 km) assuming on-board accelerometers with the same noise characteristics as those currently used for GRACE. The second pair of missions is flown at an altitude of approx. 250 km which requires a drag-free system to prevent satellite re-entry. In addition to allowing a lower satellite altitude, the drag-free system also reduces the errors associated with the accelerometer. All simulated mission scenarios assume a two satellite co-orbiting pair similar to GRACE in a near-polar, near-circular orbit. A method for local time variable gravity recovery through mass concentration blocks (mascons) is used to form simulated gravity estimates for Greenland and the Amazon region for three GFO configurations and GRACE. Simulation results show that the increased precision of the laser does not improve gravity estimation when flown with on-board accelerometers at the same altitude and spacecraft separation as GRACE, even when time-varying background models are not included. This study also shows that only modest improvement is realized for the best-case scenario (laser, low-altitude, drag-free) as compared to GRACE due to temporal aliasing errors. These errors are caused by high-frequency variations in the hydrology signal and imperfections in the atmospheric, oceanographic, and tidal models which are used to remove unwanted signal. This work concludes that applying the updated technologies alone will not immediately advance the accuracy of the gravity estimates. If the scientific objectives of a GFO mission require more accurate gravity estimates, then future work should focus on improvements in the geophysical models, and ways in which the mission design or data processing could reduce the effects of temporal aliasing.

Loomis, Bryant D.↗

Gravitational clock: A proposed experiment for the measurement of the gravitational constant G

The increased importance and the fundamental significance of accurately measuring the gravitational constant G are discussed along with recent or proposed experimental measurements of G. The method of using mutually gravitating bodies in the clock mode in a drag-free satellite is described. A satellite experiment consisting of the flat-plate spherical mass oscillator proposed combines the mathematical and experimental conveniences most simply. It is estimated that accuracies of 1 part in 1,000,000 are easily obtainable by careful fabrication of parts. The use of cryogenic techniques, thin films, and superconductors allows increased accuracies of two or three orders of magnitude or better. These measurements can be increased to the level of 1 part in 10 to the 11th power at which time-variations, and other variations, in G can be observed.

Smalley, L. L.↗

On the electrostatic charging of the cometary nucleus

The paper demonstrates that a water-dominated cometary nucleus is not shielded from solar UV radiation or solar wind by protective atmospheres beyond a heliocentric distance of 5 AU. It is also shown that submicron dust particles can be electrostatically levitated and blown off the cometary surface. An upper limit on the total mass that can be lost from the sunlit hemisphere is derived using two solar wind models, and the dust mantle erosion is estimated at 5-12 cm. It is noted that the erosion is comparable to the mass of dust blown off a typical medium bright comet by the viscous drag of the sublimating gases during perihelium passage between 3.0 and 0.5 AU, and it is therefore proposed, that the surface dust blow-off be included in the current balance equation.

Mendis, D. A.↗

NDARC-NASA Design and Analysis of Rotorcraft Theoretical Basis and Architecture

The theoretical basis and architecture of the conceptual design tool NDARC (NASA Design and Analysis of Rotorcraft) are described. The principal tasks of NDARC are to design (or size) a rotorcraft to satisfy specified design conditions and missions, and then analyze the performance of the aircraft for a set of off-design missions and point operating conditions. The aircraft consists of a set of components, including fuselage, rotors, wings, tails, and propulsion. For each component, attributes such as performance, drag, and weight can be calculated. The aircraft attributes are obtained from the sum of the component attributes. NDARC provides a capability to model general rotorcraft configurations, and estimate the performance and attributes of advanced rotor concepts. The software has been implemented with low-fidelity models, typical of the conceptual design environment. Incorporation of higher-fidelity models will be possible, as the architecture of the code accommodates configuration flexibility, a hierarchy of models, and ultimately multidisciplinary design, analysis and optimization.

Johnson, Wayne↗

Asteroid Fragmentation Approaches for Modeling Atmospheric Energy Deposition

During asteroid entry, energy is deposited in the atmosphere through thermal ablation and momentum-loss due to aerodynamic drag. Analytic models of asteroid entry and breakup physics are used to compute the energy deposition, which can then be compared against measured light curves and used to estimate ground damage due to airburst events. This work assesses and compares energy deposition results from four existing approaches to asteroid breakup modeling, and presents a new model that combines key elements of those approaches. The existing approaches considered include a liquid drop or "pancake" model where the object is treated as a single deforming body, and a set of discrete fragment models where the object breaks progressively into individual fragments. The new model incorporates both independent fragments and aggregate debris clouds to represent a broader range of fragmentation behaviors and reproduce more detailed light curve features. All five models are used to estimate the energy deposition rate versus altitude for the Chelyabinsk meteor impact, and results are compared with an observationally derived energy deposition curve. Comparisons show that four of the five approaches are able to match the overall observed energy deposition profile, but the features of the combined model are needed to better replicate both the primary and secondary peaks of the Chelyabinsk curve.

Modeling↗

Atomistic Investigation of Plastic Deformation and Dislocation Motion in Uranium Mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. However, many aspects of its mechanical behavior, particularly at reactor-relevant conditions, remain unclear. In this study, molecular dynamics (MD) simulations were employed to investigate the deformation behavior and dislocation motion in UN. We found that the Kocevski potential predicts the principal slip system as $\frac{1}{2}$ $\langle110\rangle${110}, aligning with experimental data. On the other hand, the Tseplyaev potential predicts slip to primarily occur on $\frac{1}{2}$ $\langle110\rangle${111}. MD simulations of stress–strain behavior were used to estimate the nanoindentation hardness, revealing that the Kocevski potential accurately predicts hardness even though it fails to model dynamic plasticity. Complete dislocation mobility functions have been fitted for the edge and screw dislocations in both the thermally activated and phonon-drag regimes. The 300 K linear mobility of the edge dislocation using the Tseplyaev potential was found to be 817 Pa -1 ·s -1 , whereas that of the screw dislocation using the Kocevski potential was found to be 4546 Pa -1 ·s -1 . At intermediate stresses, we observed that the subsonic steady-state motion of the edge dislocation in UN is intermittently interrupted by velocity jumps, reaching the average sound velocity. Finally, the threshold Schmid stress is calculated as 179–197 MPa, which gives an upper-limit estimate of the uniaxial yield stress of polycrystalline UN of 548–603 MPa. These findings, including the fitted dislocation mobility function, provide essential input for future plasticity and dislocation dynamics models of nuclear fuels.

36 MATERIALS SCIENCE↗

Enhanced viscous flow drag reduction using acoustic excitation

Large eddy break up devices (LEBUs) constitute a promising method of obtaining drag reduction in a turbulent boundary layer. Enhancement of the LEBU effectiveness by exciting its trailing edge with acoustic waves phase locked to the large scale structure influencing the momentum transfer to the wall is sought. An initial estimate of the required sound pressure level for an effective pulse was obtained by considering the magnitude of the pressure perturbations at the near wake of a thin plate in inviscid flow. Detailed skin friction measurments were obtained in the flow region downstream of a LEBU excited with acoustic waves. The data are compared with skin friction measurements of a simply manipulated flow, without acoustic excitation and with a plain flow configuration. The properties and the scales of motion in the flow regime downstream of the acoustically excited LEBU are studied. A parametric study based upon the characteristics of the acoustic input was pursued in addition to the careful mapping of the drag reduction phenomenon within the acoustically manipulated boundary layer. This study of boundary layer manipulation has lead to improved skin friction drag reduction and further understanding of the turbulent boundary layer.

Nagel, R. T.↗

Orbital, Rotational and Climatic Interactions: Energy Dissipation and Angular Momentum Exchange in the Earth-Moon System

A numerical ocean tide model has been developed and tested using highly accurate TOPEX/Poseidon (T/P) tidal solutions. The hydrodynamic model is based on time stepping a finite difference approximation to the non-linear shallow water equations. Two novel features of our implementation are a rigorous treatment of self attraction and loading (SAL), and a physically based parameterization for internal tide (IT) radiation drag. The model was run for a range of grid resolutions, and with variations in model parameters and bathymetry. For a rational treatment of SAL and IT drag, the model run at high resolution (1/12 degree) fits the T/P solutions to within 5 cm RMS in the open ocean. Both the rigorous SAL treatment and the IT drag parameterization are required to obtain solutions of this quality. The sensitivity of the solution to perturbations in bathymetry suggest that the fit to T/P is probably now limited by errors in this critical input. Since the model is not constrained by any data, we can test the effect of dropping sea-level to match estimated bathymetry from the last glacial maximum (LGM). Our results suggest that the 100 m drop in sea-level in the LGM would have significantly increased tidal amplitudes in the North Atlantic, and increased overall tidal dissipation by about 40%. However, details in tidal solutions for the past 20 ka are sensitive to the assumed stratification. IT drag accounts for a significant fraction of dissipation, especially in the LGM when large areas of present day shallow sea were exposed, and this parameter is poorly constrained at present.

Egbert, Gary D.↗

Gamma-ray bursts from remnant neutron star disks

The consequences of a disk of matter orbiting an old neutron star are examined. When the inner edge of the disk approaches close to the star, due to internal viscous drag, runaway ionization of the disk occurs and the resulting plasma is precipitated to the surface of the neutron star, thereby producing a gamma-ray burst. Rough numerical estimates of the occurrence rate are given and found to be consistent with gamma-ray burst observations. The estimates indicate that energies of 10 to the 39th ergs or more could be released with rise times as fast as 0.3 ms. Consideration is given to explanations of the March 5, 1979 event (Cline et al., 1980). Some possible observational searches for optical or IR emission from such disks at the locations of known burst sources and pulsars are discussed.

Michel, F. C.↗

Configuration trade and code validation study on a conical hypersonic vehicle

A study has been conducted on a generic wing-cone transatmospheric vehicle at Mach numbers form 2.5 to 4.5. The objectives of the study were to experimentally define the aerodynamic characteristics of the vehicle and evaluate several computational aerodynamic prediction methods through comparison with the experimental results. The baseline wing-cone configuration fuselage consisted of a 5 deg half-angle cone forebody, cylindrical midbody, and 9 deg truncated cone afterbody. The 4-percent-thick diamond airfoil wing had an aspect ratio of 1. Several configuration variables were investigated to provide trade information on canard, wing-position and incidence, vertical tail, and nose bluntness effects. Results of the study show that wing-position and wing-incidence effects on the longitudinal aerodynamic characteristics can be significantly influenced by wing-body interference. The use of positive wing incidence to provide favorable forebody orientation for possible inlet performance improvement is accompanied by trim drag and lift-drag ratio penalties. The lateral-directional stability characteristics were strongly influenced by the location of the vertical tails. The higher-order full-potential method provided better estimates of the aerodynamic characteristics than either the linearized supersonic potential method or the tangent-cone/tangent-wedge/shock-expansion on method.

Covell, Peter F.↗

Longitudinal aerodynamic characteristics of a subsonic, energy-efficient transport configuration in the National Transonic Facility

The Reynolds number, aeroelasticity, boundary layer transition, and nonadiabatic wall temperature effects, and data repeatability was determined in the National Transonic Facility (NTF) for a subsonic, energy efficient transport model. The model was tested over a Mach number range of 0.50 to 0.86 and a Reynolds number range of 1.9 million to approximately 23.0 million (based on mean geometric chord). The majority of the data was taken using cryogenic nitrogen (data at 1.9 million Reynolds number was taken in air). Force and moment, wing pressure, and wing thermocouple data are presented. The data indicate that increasing Reynolds number resulted in greater effective camber of the supercritical wing and horizontal tail, resulting in greater lift and pitching moment coefficients at nearly all angles of attack for M = 0.82. As Reynolds number was increased, untrimmed L/D increased, the angle of attack for maximum L/D decreased, drag creep was reduced significantly, and drag divergence Mach number increased slightly. Data repeatability for both modes of operation of the NTF (air and cryogenic nitrogen) was generally very good, and nonadiabatic wall effects were estimated to be small. Transition-free and transition-fixed configurations had significantly different force and moment data at M = 0.82 for low Reynolds number, and very small differences were noted at high Reynolds numbers.

Jacobs, Peter F.↗

Hydrodynamic and Aerodynamic Tests of Four Models of Outboard Floats : (N.A.C.A. Models 51-A, 51-B, 51-C, and 51-D)

Four models of outboard floats (N.A.C.A. models 51-A, 51-B, 51-C, and 51-D) were tested in the N.A.C.A. tank to determine their hydrodynamic characteristics and in the 20-foot wind tunnel to determine their aerodynamic drag. The results of the tests, together with comparisons of them, are presented in the form of charts. From the comparisons, the order of merit of the models is estimated for each factor considered. The best compromise between the various factors seems to be given by model 51-D. This model is the only one in the series with a transverse step.

Dawson, John R↗

The Use of the Satellite Breakup Risk Assessment Model (SBRAM) to Characterize Collision Risk to Manned Spacecraft

NASA uses environment models such as ORDEM96 to characterize the long-term orbital debris collision hazard for spacecraft in LEO. Occasionally, however, there are breakups of satellites or rocket bodies that create enhanced collision hazard for a period of time. This enhanced collision hazard can pose increased risks to space operations - especially those involving manned missions where the tolerance for risk is very low. NASA has developed SBRAM to simulate the enhanced debris environment in the days and weeks that follow such a breakup. This simulation provides the kind of risk probabilities that can be used by mission planners to consider if changes are warranted for the mission. Announcements of breakups come to NASA from US Space Command as soon as they are identified. The pre-breakup orbit and time of breakup are used to determine the initial conditions of the explosion. SBRAM uses the latest explosion models developed at NASA to simulate a debris cloud for the breakup. The model uses a Monte Carlo technique to create a random debris cloud from the probability distributions in the breakup model. Each piece of debris randomly created in the cloud is propagated in a deterministic manner to include the effects of drag and other orbital perturbations. The detailed geometry of each simulated close approach to the target spacecraft is noted and logged and the collision probability is computed using an estimated probability density in down-range and cross-range positions of both the target spacecraft and debris object. The collision probability is computed from the overlap of these probability densities for each close-approach geometry and summed over all computed conjunctions. Cloud propagation runs over the desired time interval are then repeated until the scale of the collision risk can be estimated to a desired precision. This paper presents an overview of the SBRAM model and a number of examples, both real and hypothetical, to demonstrate its use. In addition, a number of different examples are shown how the data can be used by decision makers on issues such as spacecraft orientation and timing of EVAs.

Matney, Mark J.↗

Nonlinear and adaptive estimation in reentry.

The problem of real-time estimation of a lifting reentry vehicle trajectory of the shuttle orbiter type is considered. Simulations feature large position and velocity uncertainties at radar acquisition and realistic model errors in lift, drag and other model parameters. Radar tracking and accelerometer data are simulated. Significant nonlinearities are found to exist on spacecraft acquisition. An iterated nonlinear filter is shown to perform optimally during the radar acquisition phase. An adaptive filter is shown to track time-varying model errors, such as errors in the lift and drag coefficients, down to the noise level. Such real-time model tracking (identification) is frequently required for guidance and control implementation.

Jazwinski, A. H.↗

Computational flow predictions for hypersonic drag devices

The effectiveness of two types of hypersonic decelerators are computationally examined: mechanically deployable flares and inflatable ballutes. CFD is used to predict the flowfield around a solid rocket motor (SRM) with a deployed decelerator. The computations are performed with an ideal gas solver using an effective specific heat ratio of 1.15. The surface pressure coefficients, the drag, and the extent of the compression corner separation zone predicted by the ideal gas solver compare well with those predicted by the nonequilibrium solver. The ideal gas solver is computationally inexpensive and is shown to be well suited for preliminary design studies. The computed solutions are used to determine the size and shape of the decelerator that are required to achieve a drag coefficient of 5 in order to assure that the SRM will splash down in the Pacific Ocean. Heat transfer rates to the SRM and the decelerators are predicted to estimate the amount of thermal protection required.

Tokarcik, Susan↗

Feasibility Study for Low Drag Acoustic Liners Final Report

This report documents the design and structural analysis as a final deliverable for the Phase 1 contract activity. Also included is a community noise test plan, which is a key deliverable for Phase 2. Finally, a high-level estimate (Phase 3 deliverable) is provided for the work statement of Phases 2-4, which covers the build of two inlet test articles, planning and execution of a flight test with the test inlets, as well as data analysis and final documentation. The two test inlets will be compared to the production baseline inlet configuration. There is also a plan to test one of the inlets "hardwalled" using speed tape or some other similar tape to block the acoustic perforations.

Riedel, Brian↗