Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable geometry”

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 649 records · Page 36

Diffusion Flame Extinction in a Low Strain Flow

Diffusion flames are of great interest in fire safety and many industrial processes. Many parameters significantly affect the flame structure, shape and stability, of particular importance are the constraints imposed by geometrical boundaries. Physical boundaries determine the characteristics of the flow, affect heat, fuel, and oxidizer transport from and towards the flame and can act as heat sinks or heat sources. As a result, the existence of a flame, its shape and nature are intimately related to the geometrical characteristics of the environment that surrounds it. The counter-flow configuration provides a constant strain flow, therefore, is ideal to study the structure of diffusion flames. Most studies have concentrated on the high velocity, high strain limit, since buoyantly induced instabilities will disintegrate the planar flame as the velocity decreases. Only recently, experimental studies in micro-gravity conditions have begun to explore the low strain regimes. The main objective of these on-going studies is to determine the effect of radiative heat losses and variable strain on the structure and radiation-induced extinction of diffusion flames. For these programs, size, geometry, and experimental conditions have been chosen to keep the flame unaffected by the physical boundaries. Whether is the burning of condensed or gaseous fuels, for most real situations the boundaries impose a significant effect on the nature of the flame. There is, therefore, a need to better understand the effect that geometrical constraints (i.e. flow nonperpendicular to a fuel surface, heat losses to the boundaries, etc.) might have on the final characteristics of a diffusion flame. Preliminary experiments have shown that, in the absence of gravity, and depending on the distance from the flame to the boundary, three characteristically different regimes can be observed. Close to the boundary, the flame is parabolic, very thin and blue, almost soot-less. Diffusion is the main mechanism controlling fuel transport to the reaction zone, conduction towards the inlets is the main source of heat losses. As the distance increases the flame becomes linear and thickens, remaining blue at the oxidizer side and turning yellow at the fuel side. Here, convection brings fuel and oxidizer together and the reaction occurs in the viscous layer formed between the fuel and oxidizer streams. This region corresponds to the characteristic counter-flow flame where conduction and convection become negligible forms of heat losses and radiation becomes dominant. The flame in the third (mixed) region, between the two others, results from the combination of the scenarios presented above.

Sutula, Jason↗

Characterization of Plasma Discharges in a High-Field Magnetic Tandem Mirror

High density magnetized plasma discharges in open-ended geometries, like Tandem Mirrors, have a variety of space applications. Chief among them is the production of variable Specific Impulse (I(sub sp)) and variable thrust in a magnetic nozzle. Our research group is pursuing the experimental characterization of such discharges in our high-field facility located at the Advanced Space Propulsion Laboratory (ASPL). These studies focus on identifying plasma stability criteria as functions of density, temperature and magnetic field strength. Plasma heating is accomplished by both Electron and Ion Cyclotron Resonance (ECR and ICR) at frequencies of 2-3 Ghz and 1-30 Mhz respectively, for both Hydrogen and Helium. Electron density and temperature has measured by movable Langmuir probes. Macroscopic plasma stability is being investigated in ongoing research.

Chang-Diaz, Franklin R.↗

Sources of Sodium in the Lunar Exosphere: Modeling Using Ground-Based Observations of Sodium Emission and Spacecraft Data of the Plasma

Observations of the equatorial lunar sodium emission are examined to quantify the effect of precipitating ions on source rates for the Moon's exospheric volatile species. Using a model of exospheric sodium transport under lunar gravity forces, the measured emission intensity is normalized to a constant lunar phase angle to minimize the effect of different viewing geometries. Daily averages of the solar Lyman alpha flux and ion flux are used as the input variables for photon-stimulated desorption (PSD) and ion sputtering, respectively, while impact vaporization due to the micrometeoritic influx is assumed constant. Additionally, a proxy term proportional to both the Lyman alpha and to the ion flux is introduced to assess the importance of ion-enhanced diffusion and/or chemical sputtering. The combination of particle transport and constrained regression models demonstrates that, assuming sputtering yields that are typical of protons incident on lunar soils, the primary effect of ion impact on the surface of the Moon is not direct sputtering but rather an enhancement of the PSD efficiency. It is inferred that the ion-induced effects must double the PSD efficiency for flux typical of the solar wind at 1 AU. The enhancement in relative efficiency of PSD due to the bombardment of the lunar surface by the plasma sheet ions during passages through the Earth's magnetotail is shown to be approximately two times higher than when it is due to solar wind ions. This leads to the conclusion that the priming of the surface is more efficiently carried out by the energetic plasma sheet ions.

Sarantos, Menelaos↗

Burning and Suppression of Solids–II (BASS-II) Summary Report

The Burning and Suppression of Solids (BASS) experiment hardware is a small flow duct that provides containment for small scale burning of solid samples within the Microgravity Science Glovebox (MSG) aboard the International Space Station (ISS). A video camera with a data overlay and a 35-mm still camera record the combustion events. The controls (ignition, fan speed, etc.) are operated by an astronaut while the principal investigator team monitors the experiment from the ground and communicates directly with the astronaut. For the first time on ISS, BASS–II utilized MSG working volume dilution with gaseous N2. We developed a perfectly stirred reactor model to determine the N2 flow time and flow rate to obtain the desired reduced O2 concentration in the working volume for each test. We calibrated the model with the Compound Specific Analyzer-Combustion Products (CSA-CP) O2 readings offset using the Major Constituents Analyzer reading of the ISS ambient atmosphere data for that day. This worked out extremely well for operations, and added a new vital variable, ambient O2 level, to our test matrices. The main variables tested in BASS–II were ambient O2 concentration, ventilation flow velocity, and fuel type, thickness, and geometry. BASS–II also utilized the onboard CSA-CP for O2 and CO readings, and the Carbon Dioxide Monitor for CO2 readings before and after each test. Readings from these sensors allow us to evaluate the completeness of the combustion. The O2 and CO2 readings before and after each test were analyzed and compared very well to stoichiometric ratios for a one-step gas-phase reaction. The CO versus CO2 followed a linear trend for some datasets, but not for all the different geometries of fuel and flow tested. We calculated the heat release rates during each test from the O2 consumption and burn times, using the constant 13.1 kJ of heat released per gram of O2 consumed. The results showed that most of the tests had heat release rates well below 100 W. Lastly, the global equivalence ratio for the tests is estimated to be fuel rich, 1.3 on average using mass loss and O2 consumption data.

Microgravity↗

Leveraging ­ CPF Spectral Information for Effective Angular Corrections in Sensor Inter-calibration Studies

The Climate Absolute Radiance and Refractivity Observatory Pathfinder (CPF) mission will provide a high-accuracy (0.3% radiometric uncertainty at k=1) SI-traceable on-orbit calibration reference for sensor intercalibration in reflective solar spectral region. CPF inter-calibration algorithms have been developed to address various sampling differences between CPF and target sensors, particularly in spectral, angular, spatial and temporal domains. CPF provides unique hyper-spectral measurements that allow the utilization of rich spectral information for critical inter-calibration procedures including angular correction and spectral-gap filling. The angular correction scheme of CPF mission has been constructed based on the spectral correlation relationship between CPF observed spectra and the spectra to be measured at the viewing/sun geometry angles of target sensors. The angular correction relationship is both angular and scene dependent. The scene variability issue can impose undesired uncertainties that complicates the correction for those errors. We have demonstrated that hyper-spectral radiance information can be well utilized to identify different scene types associated with CPF observations. Using spectral radiances-based scene clustering analysis can greatly reduce the uncertainty associated with the angular correction algorithm. Additionally, a reliable quality control scheme has been established, utilizing spectral similarity and regression-prediction error analyses, to ensure the successful and efficient implementation of the angular correction method.

Wan Wu↗

Cosmological tests of the Hoyle-Narlikar conformal gravity

For the first time the Hoyle-Narlikar theory with creation of matter and a variable gravitational constant G, is subjected to the following cosmological tests: (1) the magnitude versus z relation, (2) the N(m) versus m relation for quasars, (3) the metric angular diameters versus z relation, (4) the isophotal angles versus z relation, (5) the log N-log S radio source count, and finally (6) the 3 K radiation. It is shown that the theory passes all these tests just as well as the standard cosmology, with the additional advantage that the geometry of the universe is uniquely determined, with a curvature parameter equal to zero. It is also interesting to note that the variability of G affects the log N-log S curve in a way similar to the density evolution introduced in standard cosmologies. The agreement with the data is therefore achieved without recourse to an ad hoc density evolution.

Canuto, V. M.↗

Two-and-a-half-dimensional magnetohydrodynamic turbulence

The homogeneous turbulence for which fluctuating magnetic fields and velocity fields are independent of one spatial coordinate but still possess all three components are studied to form a generalized two-and-half dimensional geometry. The integral of the z-vector potential and the magnetic helicity are shown to be ideal invariants, and the basic dynamical variables and the equations of uniform-density incompressible magnetohydrodynamics are defined. The possibility of simultaneous inverse cascades is considered, and Kolmogoroff dimensional analysis is employed to infer omnidirectional inverse cascade spectra. Implications of a selective decay hypothesis, where a tendency exists in the initial value problem for the ideal invariants directly cascaded to higher wavenumbers to be selectively dissipated, and the possibility of a dynamo action in the considered geometry, are examined.

Montgomery, D.↗

Seasonal and Dayurnal Planetary Albedo Variability from Six Years of DSCOVR EPIC Data

Deep Space Climate Observatory (DSCOVR) measurements of Earth’s reflected solar radiation from the Lissajous orbital position near the Lagrangian L1 point provide continuous monitoring of the Earth’s sunlit hemisphere. CERES-based angle models were used to convert the near-hourly reflected radiances of the EPIC images into a climate-style planetary albedo data-point over the sunlit hemisphere. Integration over the sunlit hemispheres averages out the meteorological weather noise, but retains the seasonal and planetary-scale variability. As the Earth rotates, this generates variations in the Earth’s planetary albedo that are precisely aligned in longitude, which constitutes the dayurnal cycle. This dayurnal variability in the planetary albedo arises from planetary-scale changes in cloud radiative properties that can be directly compared to similarly sampled climate GCM output data. Six years of EPIC data have been analyzed, showing characteristic patterns in the seasonal and dayurnal variability of the Earth’s planetary albedo. Much of the seasonal change in planetary albedo is associated with the changing DSCOVR viewing geometry and the change in solar declination. But throughout the year, the highest planetary albedos are observed over the Central Asia (Iraq) longitude, while the lowest planetary albedos occur over the Central Pacific longitude. For these longitudes, the relative seasonal changes in the planetary albedo are slowly varying and anti-correlated. Dayurnal amplitude maxima tend to occur during the July-September timeframe, with April-May and December exhibit distinct minima in the dayurnal amplitude. West Africa and the West Pacific longitudes likewise exhibit anti-correlated seasonal variability, while they also undergo anti-correlated short period oscillations. On the other hand, in nearby longitudes, there are short period spikes in planetary albedo of a few-days duration, as well as longer period oscillations that may range from a week to several months, that generally tend to be correlated.

Planetary Albedo↗

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗

Comparison of Two Algebraic Methods for Curve/curve Intersection

Most geometric modeling systems use either polynomial or rational functions to represent geometry. In such systems most computational problems can be formulated as systems of polynomials in one or more variables. Classical elimination theory can be used to solve such systems. Here Cayley's method of elimination is summarized and it is shown how it can best be used to solve the curve/curve intersection problem. Cayley's method was found to be a more straightforward approach. Furthermore, it is computationally simpler, since the elements of the Cayley matrix are one variable instead of two variable polynomials. Researchers implemented and tested both methods and found Cayley's to be more efficient. Six pairs of curves, representing mixtures of lines, circles, and cubic arcs were used. Several examples had multiple intersection points. For all six cases Cayley's required less CPU time than the other method. The average time ratio of method 1 to method 2 was 3.13:1, the least difference was 2.33:1, and the most dramatic was 6.25:1. Both of the above methods can be extended to solve the surface/surface intersection problem.

Demontaudouin, Y.↗

A multistage mesh generator for solving the average-passage equation system

One means of numerically simulating the 3-D flow field within a multistage turbomachine is through the solution of the average-passage equation system. One requirement of a current algorithm used to solve this system of equations has been the ability to generate multiple blade row meshes which satisfy specific geometrical constraints. In addition to meeting this criterion, one desires a mesh generation code which requires minimal user input, utilizes variable mesh control parameters, generates diagnostics helpful to the user, and possesses the capability to handle widely varying geometries. A mesh generation code with these features was written and has been used in solving the inviscid form of the average-passing equation system for both ducted and unducted multiple blade row geometries. This paper serves as a user reference guide, with a description of the mesh generation algorithm, a sample input file, and examples of typical meshes generated.

Mulac, Richard A.↗

Two-scale model of a geomagnetic field variation

The effect of the vertical scale is investigated by considering a simple kinematic two-scale model of fluid flow inducing a variable magnetic field. Depending on the time constant, the induced magnetic field displays a variety of behaviors and geometries. In the high-frequency case, for example, a strong magnetic field tangential to the core mantle boundary, and hidden in the Delta layer, can be generated. A detailed computation and description of this magnetic field are presented. Some possible features of the secular variation of the actual geomagnetic field are discussed in the light of the model proposed here.

Braginsky, S. I.↗

Multi-grid domain decomposition approach for solution of Navier-Stokes equations in primitive variable form

The new multigrid (or adaptive) pseudospectral element method was carried out for the solution of incompressible flow in terms of primitive variable formulation. The desired features of the proposed method include the following: (1) the ability to treat complex geometry; (2) high resolution adapted in the interesting areas; (3) requires minimal working space; and (4) effective in a multiprocessing environment. The approach for flow problems, complex geometry or not, is to first divide the computational domain into a number of fine-grid and coarse-grid subdomains with the inter-overlapping area. Next, it is necessary to implement the Schwarz alternating procedure (SAP) to exchange the data among subdomains, where the coarse-grid correction is used to remove the high frequency error that occurs when the data interpolation from the fine-grid subdomain to the coarse-grid subdomain is conducted. The strategy behind the coarse-grid correction is to adopt the operator of the divergence of the velocity field, which intrinsically links the pressure equation, into this process. The solution of each subdomain can be efficiently solved by the direct (or iterative) eigenfunction expansion technique with the least storage requirement, i.e. O(N(exp 3)) in 3-D and O(N(exp 2)) in 2-D. Numerical results of both driven cavity and jet flow will be presented in the paper to account for the versatility of the proposed method.

Ku, Hwar-Ching↗

Independent pixel and Monte Carlo estimates of stratocumulus albedo

Monte Carlo radiative transfer methods are employed here to estimate the plane-parallel albedo bias for marine stratocumulus clouds. This is the bias in estimates of the mesoscale-average albedo, which arises from the assumption that cloud liquid water is uniformly distributed. The authors compare such estimates with those based on a more realistic distribution generated from a fractal model of marine stratocumulus clouds belonging to the class of 'bounded cascade' models. In this model the cloud top and base are fixed, so that all variations in cloud shape are ignored. The model generates random variations in liquid water along a single horizontal direction, forming fractal cloud streets while conserving the total liquid water in the cloud field. The model reproduces the mean, variance, and skewness of the vertically integrated cloud liquid water, as well as its observed wavenumber spectrum, which is approximately a power law. The Monte Carlo method keeps track of the three-dimensional paths solar photons take through the cloud field, using a vectorized implementation of a direct technique. The simplifications in the cloud field studied here allow the computations to be accelerated. The Monte Carlo results are compared to those of the independent pixel approximation, which neglects net horizontal photon transport. Differences between the Monte Carlo and independent pixel estimates of the mesoscale-average albedo are on the order of 1% for conservative scattering, while the plane-parallel bias itself is an order of magnitude larger. As cloud absorption increases, the independent pixel approximation agrees even more closely with the Monte Carlo estimates. This result holds for a wide range of sun angles and aspect ratios. Thus, horizontal photon transport can be safely neglected in estimates of the area-average flux for such cloud models. This result relies on the rapid falloff of the wavenumber spectrum of stratocumulus, which ensures that the smaller-scale variability, where the radiative transfer is more three-dimensional, contributes less to the plane-parallel albedo bias than the larger scales, which are more variable. The lack of significant three-dimensional effects also relies on the assumption of a relatively simple geometry. Even with these assumptions, the independent pixel approximation is accurate only for fluxes averaged over large horizontal areas, many photon mean free paths in diameter, and not for local radiance values, which depend strongly on the interaction between neighboring cloud elements.

Cahalan, Robert F.↗

Non-Coalescence Effects in Microgravity

Non-coalescence of two bodies of the same liquid and the suppression of contact between liquid drops and solid surfaces is being studied through a pair of parallel investigations being conducted at the Georgia Institute of Technology and the Microgravity Research and Support (MARS) Center in Naples, Italy. Both non-coalescence and contact suppression are achieved by exploiting the mechanism of thermocapillary convection to drive a lubricating film of surrounding gas (air) into the space between the two liquid free surfaces (non-coalescence) or between the drop free surface and the solid (contact suppression). Earlier experiments performed included flow-visualization experiments in both axisymmetric and (nearly) two-dimensional geometries and quantitative measurements of film thickness in the contact-suppression case in both geometries. Work done in the second year has focused on obtaining quantitative results relating to the effects of variable air pressure, development of analytical and numerical models of non-coalescing droplets and to pursuing potential applications of these self-lubricated systems.

Neitzel, G. Paul↗

The Oxidation Kinetics of Continuous Carbon Fibers in a Cracked Ceramic Matrix Composite

Experimental observations and results suggest two primary regimes as a function of temperature, i.e., diffusion and reaction controlled kinetics. Thermogravimetric analysis of carbon fiber in flowing oxygen gave an activation energy of 64.1 kJ/mol in the temperature range of 500 to 600 C and an apparent activation energy of 7.6 kJ/mol for temperatures from 600 to 1400 C. When C/SiC composite material was unstressed, matrix effects at temperatures from 900 to 1400 C protected the internal fibers. When under stress, self-protection was not observed. Increasing the stress from 10 to 25 ksi caused a 67 to 82 percent reduction in times to failure at temperatures from 750 to 1500 C. Based on experimental results, observation, and theory, a finite difference model was developed, which simulates the diffusion of oxygen into a matrix crack that is bridged by carbon fibers. The model allows the influence of important variables on oxidation kinetics to be studied systematically, i.e., temperature, reaction rate constant, diffusion coefficient, environment, and sample geometry.

Halbig, Michael C.↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called 'structured' (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called "structured" (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗