Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUATIONS - DRAG, PROFILE - BICKNELL”

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 271 records · Page 15

A Topside Equatorial Ionospheric Density and Composition Climatology During and After Extreme Solar Minimum

During the recent solar minimum, solar activity reached the lowest levels observed during the space age. This extremely low solar activity has accompanied a number of unexpected observations in the Earth's ionosphere and thermosphere when compared to previous solar minima. Among these are the fact that the ionosphere is significantly contracted beyond expectations based on empirical models. Climatological altitude profiles of ion density and composition measurements near the magnetic dip equator are constructed from the C/NOFS satellite to characterize the shape of the top side ionosphere during the recent solar minimum and into the new solar cycle. The variation of the profiles with respect to local time, season, and solar activity are compared to the IRI-2007 model. Building on initial results reported by Heelis et al. [2009], here we describe the extent of the contracted ionosphere, which is found to persist throughout 2009. The shape of the ionosphere during 2010 is found to be consistent with observations from previous solar minima.

Klenzing, J. H.↗

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗

Relations between VAS physical retrievals and their first guess input

The Smith (1983) physical retrieval algorithm, which is designed to yield mesoscale soundings through a blending of first-guess limited fine mesh-derived profiles with VISSR Atmospheric Sounder (VAS) detected radiance data by means of an inversion of the radiative transfer equation, yields retrievals that may be better related to ground truth data than the guess profiles at some levels. At other levels, the first-guess profiles were slightly superior to the retrievals. Correlation functions related to the window and midtropospheric temperature channels indicate that magnitudes and orientations of retrieval gradients were intermediate to those of the guess data and brightness temperatures.

Beven, John L.↗

The composition and vertical structure of the lower cloud deck on Venus

Excluding the insertion of a probe into the lower Venus atmosphere, the only means of examining the atmospheric properties below the upper cloud layer are by passive or active radio techniques. Particularly, in a radio occultation experiment, the amplitude and Doppler data of a signal which has been refracted by the dense Venus atmosphere provide temperature, pressure, atmospheric density, cloud mass content, and some information about the cloud composition. The Venus orbiting spacecraft Pioneer Venus has provided radio occultation data for more than two years. The present investigation has the primary goal to use the absorption coefficient profiles derived from the radio occultation experiment to determine the mass content and composition of the Venusian lower deck. Attention is given to an experiment description, liquid content calculations, the mass content equation, the subtraction of gaseous absorbers, dielectric constant measurements, the concentration profile, the spherical shell model, and mass content profiles for all orbits.

Cimino, J.↗

Phenomenological and experimental study of the thermal response of low density silica ablators to high enthalpy plasma flows

A complete analysis of the melting-ablator governing equations is carried out for porous ablators with the purpose of establishing an approximate analytical model capable of determining surface recession, temperature, and heat flux as well as in-depth temperature profiles. A detailed dimensional analysis of the melt-flow governing equations is performed for parameters applicable to high-enthalpy flows. Results of arc-jet flow experiments performed on NASA Reusable Surface Insulation materials of different densities and over a range of stagnation pressures and heat fluxes are discussed. A dimensionless correlation for melt run-offs is derived in terms of material and flow parameters; however, it is pointed out that the accuracy of the current data is sufficient to conclude the compliance with the model.

Henline, William D.↗

Mapping and Geological Analysis of Mercury and Venus Radar Ranging Data

Although many radar profiles and images of the areas within 20 deg of Mercury's equator had been obtained from 1971, at both Goldstone and Arecibo radar facilities, surprisingly little geological analysis had been done with these data until recently. Topographic profiles and radar roughness reflectivity images which can be derived from these data will be crucial in completing the geological mapping of Mercury now underway at the U.S. Geological Survey. Earth based radar observations of Mercury and the other terrestrial planets as well are a potentially very valuable tool in the determination of the physical nature of their surfaces. Processing of available radar data must be completed to establish any systematic relationship between raar reflectivities and roughness, density, dielectric constant, and other related geological parameters.

Clark, P. E.↗

Computation of Large-Scale Structure Jet Noise Sources With Weak Nonlinear Effects Using Linear Euler

An approximate technique is presented for the prediction of the large-scale turbulent structure sound source in a supersonic jet. A linearized Euler equations code is used to solve for the flow disturbances within and near a jet with a given mean flow. Assuming a normal mode composition for the wave-like disturbances, the linear radial profiles are used in an integration of the Navier-Stokes equations. This results in a set of ordinary differential equations representing the weakly nonlinear self-interactions of the modes along with their interaction with the mean flow. Solutions are then used to correct the amplitude of the disturbances that represent the source of large-scale turbulent structure sound in the jet.

Dahl, Milo D.↗

Topside Equatorial Ionospheric Density and Composition During and After Extreme Solar Minimum

During the recent solar minimum, solar activity reached the lowest levels observed during the space age. This extremely low solar activity has accompanied a number of unexpected observations in the Earth s ionosphere-thermosphere system when compared to previous solar minima. Among these are the fact that the ionosphere is significantly contracted beyond expectations based on empirical models. Altitude profiles of ion density and composition measurements near the magnetic dip equator are constructed from the Communication/Navigation Outage Forecast System (C/NOFS) satellite to characterize the shape of the topside ionosphere during the recent solar minimum and into the new solar cycle. The variation of the profiles with respect to local time, season, and solar activity are compared to the IRI-2007 model. Building on initial results reported by Heelis et al. (2009), here we describe the extent of the contracted ionosphere, which is found to persist throughout 2009. The shape of the ionosphere during 2010 is found to be consistent with observations from previous solar minima.

Klenzing, J.↗

Mapping and geological analysis of Mercury radar data

Although many radar profiles and images of the area within 20 deg of Mercury's equator had been obtained from 1971 to 1981, at both Goldstone and Arecibo radar facilities, surprisingly little geological analysis had been done with these data until recently. Topographic profiles and radar roughness reflectivity images which can be derived from these data will be crucial in completing the geological mapping of Mercury now underway at the U.S. Geological Survey. Processing of available radar data must be completed to establish any systematic relationship between radar reflectivities and roughness, density, dielectric constant, and other related geological parameters. Specific tasks accomplished for these purposes include the following. Documentation was located and searched to establish the type and quantity of Goldstone 12.5 cm radar observations which were available for Mercury. Data has been collected during approximately 50 observation periods from 1971 to 1981. About half of the data, collected during 1972 and 1973, have been processed, but without adequate documentation. A standardized, well-documented procedure for processing and analysis for all Goldstone Earth-based observations of Mercury was established.

Clark, P. E.↗

The structure of a three-dimensional turbulent boundary layer

The three-dimensional turbulent boundary layer is shown to have a self-consistent two-layer asymptotic structure in the limit of large Reynolds number. In a streamline coordinate system, the streamwise velocity distribution is similar to that in two-dimensional flows, having a defect-function form in the outer layer which is adjusted to zero at the wall through an inner wall layer. An asymptotic expansion accurate to two orders is required for the cross-stream velocity which is shown to exhibit a logarithmic form in the overlap region. The inner wall-layer flow is collateral to leading order but the influence of the pressure gradient, at large but finite Reynolds numbers, is not negligible and can cause substantial skewing of the velocity profile near the wall. Conditions under which the boundary layer achieves self-similarity and the governing set of ordinary differential equations for the outer layer are derived. The calculated solution of these equations is matched asymptotically to an inner wall-layer solution and the composite profiles so formed describe the flow throughout the entire boundary layer. The effects of Reynolds number and cross-stream pressure gradient on the crossstream velocity profile are discussed and it is shown that the location of the maximum cross-stream velocity is within the overlap region.

Degani, A. T.↗

Methods of Stress Calculation in Rotating Disks

The paper describes nethods of computing the stresses in disks of a given profile as well as methods of choosing the disk profiles for a given stress distribution for turhines, turbo blowers, and so forth. A new method of in tegrating the differential equations of Stodola leads to a simplification of the computation for disks of hyperbolic profile.

Tumarkin, S.↗

Sporadic-E layers and unstable wind shears

Electron density profiles of sporadic-E layers have been observed with good height resolution using rocket-borne probes. These generally show a simple shape consistent with the effect of a linear wind shear acting on metallic ions. Occasionally more complex shapes have been recorded, including double peaks and, on one occasion, a nearly rectangular profile. A direct method of obtaining the wind profile from the concentration profile of metallic ions has been developed. The metallic ion concentration profile itself is obtained from the electron density profile. Both procedures derive from the steady-state continuity equation. For linear wind shears it is found that the maximum value of the shear is about 50 m/s/km which corresponds to a Richardson number of 1/4. Layers of complex shape are associated with non-linear wind shears in which the maximum shear considerably exceeds this value. It is concluded that the complex profiles of sporadic-E layers can be interpreted as an effect of unstable wind shears.

Smith, L. G.↗

Numerical Issues for Circulation Control Calculations

Steady-state and time-accurate two-dimensional solutions of the compressible Reynolds-averaged Navier- Stokes equations are obtained for flow over the Lockheed circulation control (CC) airfoil and the General Aviation CC (GACC) airfoil. Numerical issues in computing circulation control flows such as the effects of grid resolution, boundary and initial conditions, and unsteadiness are addressed. For the Lockheed CC airfoil computed solutions are compared with detailed experimental data, which include velocity and Reynolds stress profiles. Three turbulence models, having either one or two transport equations, are considered. Solutions are obtained on a sequence of meshes, with mesh refinement primarily concentrated on the airfoil circular trailing edge. Several effects related to mesh refinement are identified. For example, sometimes sufficient mesh resolution can exclude nonphysical solutions, which can occur in CC airfoil calculations. Also, sensitivities of the turbulence models with mesh refinement are discussed. In the case of the GACC airfoil the focus is on the difference between steady-state and time-accurate solutions. A specific objective is to determine if there is self-excited vortex shedding from the jet slot lip.

Swanson, Roy C., Jr.↗

Plasma property and performance prediction for mercury ion thrusters

The discharge chambers of mercury ion thrusters are modelled so the principal effects and processes which govern discharge plasma properties and thruster performance are described. The conservation relations for mass, charge and energy when applied to the Maxwellian electron population in the ion production region yield equations which may be made one-dimensional by the proper choice of coordinates. Solutions to these equations with the appropriate boundary conditions give electron density and temperature profiles which agree reasonably well with measurements. It is then possible to estimate plasma properties from thruster design data and those operating parameters which are directly controllable. By varying the operating parameter inputs to the computer code written to solve these equations, perfromance curves are obtained which agree quite well with measurements.

Longhurst, G. R.↗

Inference of dust opacities for the 1977 Martian great dust storms from Viking Lander 1 pressure data

The tidal heating components for the dusty Martian atmosphere are computed based on dust optical parameters estimated from Viking Lander imaging data, and used to compute the variation of the tidal surface pressure components at the Viking Lander sites as a function of season and the total vertical extinction optical depth of the atmosphere. An atmospheric tidal model is used which is based on the inviscid, hydrostatic primitive equations linearized about a motionless basic state the temperature of which varies only with height, and the profiles of the tidal forcing components are computed using a delta-Eddington approximation to the radiative transfer equations. Comparison of the model results with the observed variations of surface pressure and overhead dust opacity at the Viking Lander 1 site reveal that the dust opacities and optical parameters derived from imaging data are roughly representative of the global dust haze necessary to reproduce the observed surface pressure amplitudes, with the exception of the model-inferred asymmetry parameter, which is smaller during the onset of a great storm. The observed preferential enhancement of the semidiurnal tide with respect to the diurnal tide during dust storm onset is shown to be due primarily to the elevation of the tidal heating source in a very dusty atmosphere.

Zurek, R. W.↗

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗

The differential rotation of the solar surface

The large-scale flow in the solar convection zone is discussed. The objective is to deduce from observation the principal physical balances in the governing equations. The simplest set of equations that seem potentially realistic are utilized. It is assumed that magnetic fields are negligible, that mixing-length theory gives an accurate representation of the mean structure, and that rigid-body rotation exists at the base. It is deduced that the zonal glow is geostrophic, the meridional flow is controlled by friction, and diffusive heating balances advective cooling due to vertical motion. Next, a detailed calculation of the latitude profile of surface angular velocity is performed, based on approximate equations containing only the principal balances, and agrees well with observation. Finally, it is demonstrated that the amplitude of the flow may be determined by the constraint that the net equatorward angular momentum flux vanish. The conclusions are consistent with the hypothesis that the flow is primarily a single large axisymmetric convection cell in each hemisphere.

Gierasch, P. J.↗