SEARCH · Engineering Papers
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.
Papers are presented dealing with interdisciplinary research in the fields of geodesy, solid earth and ocean physics. Topics discussed include: solid earth and ocean tides; the rotation of the earth and polar motion; vertical crustal motions; the geoid and ocean surface; earthquake mechanism; sea level changes; and lunar dynamics.
Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem.
Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem. Previously announced in STAR as N83-30109
A series of time-dependent numerical computations have been performed for the flow of an unsteady boundary layer over an embedded cavity, which is aligned normal to the flow direction. The unsteady flow is the result of the imposition of a Stuart vortex, at the inflow boundary of the computational domain, onto an otherwise unperturbed two-dimensional laminar boundary layer flow. This produces a strongly vortical, but spatially monochromatic, motion at inflow which is allowed to evolve downstream through the time-dependent numerical solution of appropriate conservation equations. As the imposed vortex evolves over the solid boundary, a weaker induced vortical motion is produced and together these perturbing motions interact with the fluid motion in the embedded cavity. An analysis of the pertinent dynamic variables as well as the time-dependent drag characteristics for this type of non-planar geometry is performed; with particular emphasis on the relative contribution of both pressure drag and frictional drag at various times in the evolution process.
A two-blade rotor having a diameter of 4 feet and a solidity of 0.037 was subjected to sharp-edge vertical gusts while being operated at various forward speeds to study the effect of the gusts on the blade periodic bending moments and flapping angles. Variables studied included gust velocity, collective pitch angle, flapping hinge offset, and tip-speed ratio. Dimensionless coefficients are derived for the periodic components of the incremental changes in blade flapping angles and bending moments which arise when a rotor blade penetrates a sharp-edge gust. Mental changes in both the flapping angles and bending moments are essentially proportional to gust velocity, and the coefficients express the ratio of these increments to gust velccity. The results show that the flapping coefficient usually increases with an increase in collective pitch angle, is generally dependent on tip-speed ratio, and is essentially independent of the amount of flapping hinge offset. The bending-moment coefficient is also dependent on collective pitch angle and tip-speed ratio. Expected reductions in bending moments are realized by the use of flapping hinges, and further reductions in bending moments are achieved as the amount of flapping hinge offset is increased. Comparison of the experimental results of this investigation with limited available theoretical results shows substantial agreement but indicates that the assumption that the response of the rotor to a sharp-edge gust is independent of the collective pitch angle prior to gust entry is probably inadequate.
Molten salts, such as lithium fluoride, are attractive candidates for thermal energy storage in solar dynamic space power systems because of their high latent heat of fusion. However, these same salts have poor thermal conductivities which inhibit the transfer of heat into the solid phase and out of the liquid phase. One concept for improving the thermal conductivity of the thermal energy storage system is to add a conductive filler material to the molten salt. High thermal conductivity pitch-based graphite fibers are being considered for this application. Although there is some information available on the thermal conductivity of lithium fluoride solid, there is very little information on lithium fluoride liquid, and no information on molten salt graphite fiber composites. This paper describes a hot filament technique for determining the thermal conductivity of molten salts. The hot filament technique was used to find the thermal conductivity of molten lithium fluoride at 930 C, and the thermal conductivity values ranged from 1.2 to 1.6 W/mK. These values are comparable to the slightly larger value of 5.0 W/mK for lithium fluoride solid. In addition, two molten salt graphite fiber composites were characterized with the hot filament technique and these results are also presented.
A two-year program to develop advanced 3D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades, and vanes is described. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulation models were developed: an eight-noded midsurface shell element; a nine-noded midsurface shell element; and a twenty-noded isoparametric solid element. A separate computer program has been developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.
In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.
In November of 2003 a rollout test was performed to gain a better understanding of the dynamic environment for the Space Shuttle during transportation from the Vehicle Assembly Building to the launch pad. This was part of a study evaluating the methodology for including the rollout dynamic loads in the Space Shuttle fatigue life predictions. The rollout test was conducted with a partial stack consisting of the Crawler Transporter, Mobile Launch Platform, and the Solid Rocket Boosters with an interconnecting crossbeam. Instrumentation included over 100 accelerometers. Data was recorded for steady state speeds, start-ups and stops, and ambient wind excitations with the vehicle at idle. This paper will describe the operating deflection shape analysis performed using the measured acceleration response data. The response data for the steady state speed runs were dominated by harmonics of the forcing frequencies, which were proportional to the vehicle speed. Assuming a broadband excitation for the wind, analyses of the data sets with the vehicle at idle were used to estimate the natural frequencies and corresponding mode shapes. Comparisons of the measured modal properties with numerical predictions are presented.
Experience with the dynamics and data analyses for earth and moon reveals both similarities and differences. Analysis of Lunar Laser Ranging (LLR) data provides information on the lunar orbit, rotation, solid-body tides, and retroreflector locations.
The study of Planet earth is a very complex problem. It has many non-linear and chaotic systems operating in parallel and have interdependencies. In reality, these systems/phenomena s are not well understood or mathematically modeled because of our lack of knowledge of such intricate processes. However, in order to further the subject of Earth as an integrated system, space provides excellent vantage points to look at these processes in multidimensional framework. For example, we can make strives to understand the global water cycle, carbon cycle, atmospheric chemistry, biomass changes, oceans and solid Earth variations by making multitude of global measurements such as soil moisture, precipitation, tropospheric and stratospheric gases, aerosols, tropospheric winds, ocean salinity, ocean color, vegetation cover, crustal dynamics and many more. Such suites of measurements derive the coupled models so we may predict the changes due to natural and anthropogenic forcing. NASA along with other international space agencies have made tremendous investments in recent years in developing and flying remote sensing space borne sensors to enable these measurements. These science measurements and products are further used to address pressing issues such as coastal zone erosion, air quality, severe weather, water availability and quality, public health, fires, earthquakes, land slides and others for societal benefits. This presentation provides a comprehensive overview of NASA s science investigations, related technologies and satellites/sensors and applications.
In order to protect the Solar Dynamics Observatory's instruments from blowing their fuses and being rendered unusable, Goddard Space Flight Center worked with Micropac Industries Inc., based in Garland, Texas, to develop solid-state power controllers, which can depower and then resupply power to an instrument in the event of an electric surge. The company is now selling the technology for use in industrial plants.
The loss of mass from the Greenland ice sheet causes an increasing influx of freshwater to the Greenlandic fjords and the oceans. Freshwater fluxes from marine-terminating glaciers are important to understand fjord circulation and ecosystem dynamics. Here, we present a data set constructed by reformulating existing products into a shared temporal and spatial framework. We combine three publicly available data sets of solid-ice discharge (iceberg), liquid-surface runoff (runoff) and basal melt to present a cohesive overview of the flow of freshwater from marine-terminating glaciers to the Greenlandic fjords. We also calculate glacier drainage basins and compare our findings to previous studies showing that drainage-basin sizes may vary considerably depending on how they were reconstructed. The data set will be a valuable asset to oceanographic, glaciological and marine biological research activities.
In recent years, there has been a resurgence of interest in the development of hybrid rocket engines for advanced launch vehicle applications. Hybrid propulsion systems use a solid fuel such as hydroxyl-terminated polybutadiene (HTPB) along with a gaseous/liquid oxidizer. The performance of hybrid combustors depends on the convective and radiative heat fluxes to the fuel surface, the rate of pyrolysis in the solid phase, and the turbulent combustion processes in the gaseous phases. These processes in combination specify the regression rates of the fuel surface and thereby the utilization efficiency of the fuel. In this paper, we employ computational fluid dynamics (CFD) techniques in order to gain a quantitative understanding of the physical trends in hybrid rocket combustors. The computational modeling is tailored to ongoing experiments at Penn State that employ a two dimensional slab burner configuration. The coordinated computational/experimental effort enables model validation while providing an understanding of the experimental observations. Computations to date have included the full length geometry with and with the aft nozzle section as well as shorter length domains for extensive parametric characterization. HTPB is sed as the fuel with 1,3 butadiene being taken as the gaseous product of the pyrolysis. Pure gaseous oxygen is taken as the oxidizer. The fuel regression rate is specified using an Arrhenius rate reaction, which the fuel surface temperature is given by an energy balance involving gas-phase convection and radiation as well as thermal conduction in the solid-phase. For the gas-phase combustion, a two step global reaction is used. The standard kappa - epsilon model is used for turbulence closure. Radiation is presently treated using a simple diffusion approximation which is valid for large optical path lengths, representative of radiation from soot particles. Computational results are obtained to determine the trends in the fuel burning or regression rates as a function of the head-end oxidizer mass flux, G=rho(e)U(e), and the chamber pressure. Furthermore, computation of the full slab burner configuration has also been obtained for various stages of the burn. Comparisons with available experimental data from small scale tests conducted by General Dynamics-Thiokol-Rocketdyne suggest reasonable agreement in the predicted regression rates. Future work will include: (1) a model for soot generation in the flame for more quantitative radiative transfer modelling, (2) a parametric study of combustion efficiency, and (3) transient calculations to help determine the possible mechanisms responsible for combustion instability in hybrid rocket motors.
Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.
The impact of passage rotation on the gas dynamic wave processes is analyzed through a numerical simulation of ideal shock-tube flow in a closed rotating-channel. Initial conditions are prescribed by assuming homentropic solid-body rotation. Relevant parameters of the problem such as wheel Mach number, hub-to-tip radius ratio, length-to-tip radius ratio, diaphragm temperature ratio, and diaphragm pressure ratio are varied. The results suggest possible criteria for assessing the consequences of passage rotation on the wave processes, and they may therefore be applicable to pressure-exchange wave rotors. It is shown that for a fixed geometry and initial conditions, the contact interface acquires a distorted three-dimensional time-dependent orientation at non-zero wheel Mach numbers. At a fixed wheel Mach number, the level of distortion depends primarily on the density ratio across the interface as well as the hub-to-tip radius ratio. Rarefaction fronts, shocks, and contact interfaces are observed to propagate faster with increasing wheel Mach number.
The solid to liquid phase transformation of LiOH at 744.3 K is considered to be an ideal candidate thermal energy storage (TES) mechanism for a Rankine heat engine based solar dynamic system operating at approximately 682 K. While pure nickel is thought to be a suitable containment material for LiOH, long term containment is of concern because molten hydroxides are usually corrosive. Two commercially pure nickel alloys, Ni-200 and Ni-201, were exposed to molten LiOH, its vapor, and vacuum at 775 K for periods ranging from 50 to 5000 h, and simple mechanical property measurements (77 to 900 K tensile and 750 K creep rupture) of exposed alloys were undertaken. The mechanical property test procedures are described and tabular lists of the test data are presented.
A simple dynamic systems analysis is used to give examples of strong, weak overdriven, and weak underdriven oblique detonations. Steady oblique detonations consisting of a straight lead shock attached to a solid wedge followed by a resolved reaction zone structure are admitted as solutions to the reactive Euler equations. This is demonstrated for a fluid that is taken to be an inviscid, calorically perfect ideal gas that undergoes a two-step irreversible reaction with the first step exothermic and the second step endothermic. This model admits solutions for a continuum of shock wave angles for two classes of solutions identified by a Rankine-Hugoniot analysis: strong and weak overdriven waves. The other class, weak underdriven, is admitted for eigenvalue shock-wave angles. Chapman-Jouguet waves, however, are not admitted. These results contrast those for a corresponding onestep model that, for detonations with a straight lead shock, only admits strong, weak overdriven, and Chapman-Jouguet solutions.