Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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 757 records · Page 42

Does the Hall Effect Solve the Flux Pileup Saturation Problem?

It is well known that magnetic flux pileup can significantly speed up the rate of magnetic reconnection in high Lundquist number resistive MHD,allowing reconnection to proceed at a rate which is insensitive to the plasma resistivity over a wide range of Lundquist number. Hence, pileup is a possible solution to the Sweet-Parker time scale problem. Unfortunately, pileup tends to saturate above a critical value of the Lundquist number, S_c, where the value ofS_c depends on initial and boundary conditions, with Sweet-Parker scaling returning above S_c. It has been argued (see Dorelli and Bim [2003] and Dorelli [2003]) that the Hall effect can allow flux pileup to saturate (when the scale of the current sheet approaches ion inertial scale, di) before the reconnection rate begins to stall. However, the resulting saturated reconnection rate, while insensitive to the plasma resistivity, was found to depend strongly on the di. In this presentation, we revisit the problem of magnetic island coalescence (which is a well known example of flux pileup reconnection), addressing the dependence of the maximum coalescence rate on the ratio of di in the "large island" limit in which the following inequality is always satisfied: l_eta di lambda, where I_eta is the resistive diffusion length and lambda is the island wavelength.

Dorelli, John C.↗

Design, development, and test of a laser velocimeter for a small 8:1 pressure ratio centrifugal compressor

The instrument was designed as a diagnostic tool for the basic fluid dynamics of the inducer, impeller, and diffuser regions of this type compressor. The LV instrumentation was optimized to measure instantaneous velocities up to approximately 500 m/s, measured in absolute coordinates, within the rotating compressor impeller and in the two dimensional radial plane of the diffuser. Some measurements were made within the diffuser and the impeller inlet flows; however, attempts to make detailed measurements of the velocity field were not successful. Difficulties in maintaining high seed particle rates within the probe volume and the improper operation of the blade gating optics may explain the lack of success. Recommendations are made to further pursue these problems. At 100% speed the stage attained a total static pressure ratio of 7.5:1 at 75% total-static efficiency. Flow range from choke-to-surge was 6.8% of choking mass flow rate. Performance was lower than the design intent of 8:1 pressure ratio at 77% efficiency and 12% flow range. Detailed measurements of the stage components are presented which show the reasons for the stage performance deficiencies.

Dolan, F. X.↗

High accuracy solutions of incompressible Navier-Stokes equations

In recent years, high accuracy finite difference approximations were developed for partial differential equations of elliptic type, with particular emphasis on the convection-diffusion equation. These approximations are of compact type, have a local truncation error of fourth order, and allow the use of standard iterative schemes to solve the resulting systems of algebraic equations. These high accuracy approximations are extended to the solution of Navier-Stokes equations. Solutions are obtained for the model problem of driven cavity and are compared with solutions obtained using other approximations and those obtained by other authors. It is discovered that the high order approximations do indeed produce high accuracy solutions and have a potential for use in solving important problems of viscous fluid flows.

Gupta, Murli M.↗

Convergence of infinite dimensional sampled LQR problems - Theory and numerical results

A theory is developed for the convergence of the closed-loop solution to infinite-dimensional discrete-time linear-quadratic regulator (LQR) problems on the infinite time interval to the solution of a corresponding continuous-time LQR problem as the length of the sampling interval tends toward zero. Convergence of solutions to the operator algebraic Riccati equation and corresponding optimal feedback control gains is guaranteed under appropriate uniform stabilizability and detectability conditions and consistent sampling. Also presented are numerical results involving the optimal LQ control of a heat or diffusion equation, a hereditary or delay differential equation, and a hybrid system of ordinary and partial differential equations describing the transverse vibration of a cantilevered Voigt-Kelvin viscoelastic beam with tip mass.

Rosen, I. G.↗

Magnetic Flux Compression Concept for Nuclear Pulse Propulsion and Power

The desire for fast, efficient interplanetary transport requires propulsion systems having short acceleration times and very high specific impulse attributes. Unfortunately, most highly efficient propulsion systems which are within the capabilities of present day technologies are either very heavy or yield very low impulse such that the acceleration time to final velocity is too long to be of lasting interest, One exception, the nuclear thermal thruster, could achieve the desired acceleration but it would require inordinately large mass ratios to reach the range of desired final velocities. An alternative approach, among several competing concepts that are beyond our modern technical capabilities, is a pulsed thermonuclear device utilizing microfusion detonations. In this paper, we examine the feasibility of an innovative magnetic flux compression concept for utilizing microfusion detonations, assuming that such low yield nuclear bursts can be realized in practice. In this concept, a magnetic field is compressed between an expanding detonation driven diamagnetic plasma and a stationary structure formed from a high temperature superconductor (HTSC). In general, we are interested in accomplishing two important functions: (1) collimation of a hot diamagnetic plasma for direct thrust production; and (2) pulse power generation for dense plasma ignition. For the purposes of this research, it is assumed that rnicrofusion detonation technology may become available within a few decades, and that this approach could capitalize on recent advances in inertial confinement fusion ICF) technologies including magnetized target concepts and antimatter initiated nuclear detonations. The charged particle expansion velocity in these detonations can be on the order of 10 (exp 6)- 10 (exp 7) meters per second, and, if effectively collimated by a magnetic nozzle, can yield the Isp and the acceleration levels needed for practical interplanetary spaceflight. The ability to ignite pure fusion micro-bursts with reasonable levels of input energy is an equally challenging scientific problem. It remains to be seen, however, whether an effective ignition driver can be developed which meets the requirements for practical spaceflight application (namely high power density, compactness, low weight, and low cost). In this paper, system level performance and design issues are examined including generator performance, magnetic flux compression processes, magnetic diffusion processes, high temperature superconductor (HTSC) material properties, plasmadynamic processes, detonation plasma expansion processes, magnetohydrodynamic instabilities, magnetic nozzle performance, and thrust production performance. Representative generator performance calculations based on a simplified skin layer formulation are presented as well as the results of exploratory small-scale laboratory experiments on magnetic flux diffusion in HTSC materials. In addition, planned follow-on scientific feasibility experiments are described which utilize high explosive detonations and high energy gas discharges to simulate the plasma conditions associated with thermonuclear micro-detonations.

Litchford, Ronald J.↗

The Role of Kinetic Processes in Magnetospheric Dynamics

The question of whether the micro scale controls the macro scale or vice-versa remains one of the most challenging problems in plasmas. A particular topic of interest within this context is collisionless magnetic reconnection, where both points of views are espoused by different groups of researchers. This presentation will focus on this topic. We will begin by analyzing the global properties of magnetic reconnection. Following that, we discuss properties of electron diffusion region dynamics both for guide field and antiparallel reconnection, and how they can be scaled to different inflow conditions. As a next step, we will study typical temporal variations of the microscopic dynamics with the objective of understanding the potential for secular changes to the macroscopic system. The research will be based on a combination of analytical theory and numerical modeling.

Hesse, Michael↗

Windblown dust on Venus

The threshold frictional velocity necessary to initiate grain movement on the Venus surface is 1 to 2 cm/s. Particles smaller than 30 or 40 microns in effective diameter will also be so moved and suspended at the threshold of movement. Theoretical arguments and Venera 8 Doppler measurements suggest marginally that dust should not be raised at the Venera 8 landing site (10 deg S), but should be raised at higher latitudes. Dust carried to tens of kilometers altitude will be transported over the entire planet and may make important contributions to the solar energy deposition, general circulation, and cloud chromosphere problems. However, the Venera 8 photometer measurements and the low albedo of reasonable surface materials imply a clear lower atmosphere at 10 deg S, despite the fact that dust raised at high latitudes should contribute to the aerosol content at this latitude. Dust raising on Venus may be inhibited by limited vertical turbulent diffusion or by thermal sintering of particles on the planetary surface.

Sagan, C.↗

Review of Problems in Application of Supersonic Combustion

The problem of air-breathing engines capable of flying at very high Mach numbers is described briefly. Possible performance of supersonic combustion ramjets is outlined briefly and the supersonic combustion process is described. Two mechanisms of combustion are outlined: one is supersonic combustion controlled by convection process, and the second is controlled by diffusion. The parameters related to the combustion process are discussed in detail. Data and analyses of reaction rates and mixing phenomena are represented; the flame mechanism is discussed, and experimental results are presented.

Antonio Ferri↗

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗

User's manual for the REEDM (Rocket Exhaust Effluent Diffusion Model) computer program

The REEDM computer program predicts concentrations, dosages, and depositions downwind from normal and abnormal launches of rocket vehicles at NASA's Kennedy Space Center. The atmospheric dispersion models, cloud-rise models, and other formulas used in the REEDM model are described mathematically Vehicle and source parameters, other pertinent physical properties of the rocket exhaust cloud, and meteorological layering techniques are presented as well as user's instructions for REEDM. Worked example problems are included.

Bjorklund, J. R.↗

A space-time discretization procedure for wave propagation problems

Higher order compact algorithms are developed for the numerical simulation of wave propagation by using the concept of a discrete dispersion relation. The dispersion relation is the imprint of any linear operator in space-time. The discrete dispersion relation is derived from the continuous dispersion relation by examining the process by which locally plane waves propagate through a chosen grid. The exponential structure of the discrete dispersion relation suggests an efficient splitting of convective and diffusive terms for dissipative waves. Fourth- and eighth-order convection schemes are examined that involve only three or five spatial grid points. These algorithms are subject to the same restrictions that govern the use of dispersion relations in the constructions of asymptotic expansions to nonlinear evolution equations. A new eighth-order scheme is developed that is exact for Courant numbers of 1, 2, 3, and 4. Examples are given of a pulse and step wave with a small amount of physical diffusion.

Davis, Sanford↗

The jet engine design that can drastically reduce oxides of nitrogen

The problem is analyzed for the case of hydrogen fuel, taking into account supersonic and hypersonic vehicles using scramjet engines. The combustion in scramjets occurs at very high velocity and in a short time. In scramjet combustor designs, two different criteria can be used to design the engine. The amount of NO formed in the diffusion flame depends substantially on the maximum temperature reached. Effects of changing the mode of combustion from a diffusion flame to a heat conduction flame are considered, giving attention to the amount of NO produced in an engine of a given design.

Ferri, A.↗

Shield topology in lightning transient control

A formalism in which the interaction of a system with an electromagnetic source is described by sets of scattering, penetrations, and propagation functions is reviewed and interpreted in the context of the lightning interaction problem. In this formalism, the system is decomposed into simple volumes separated by closed shield surfaces. These surfaces are nested and interconnected to produce higher levels of shielding and subvolumes within a given level. The interaction analysis uses scattering theory to define current and charge densities on the shield surface in conjunction with the diffusion, apertures, and transmission line analysis to define penetration through shield imperfections and propagation within the protected volume.

Vance, E. F.↗

Filtering analysis of a direct numerical simulation of the turbulent Rayleigh-Benard problem

A filtering analysis of a turbulent flow was developed which provides details of the path of the kinetic energy of the flow from its creation via thermal production to its dissipation. A low-pass spatial filter is used to split the velocity and the temperature field into a filtered component (composed mainly of scales larger than a specific size, nominally the filter width) and a fluctuation component (scales smaller than a specific size). Variables derived from these fields can fall into one of the above two ranges or be composed of a mixture of scales dominated by scales near the specific size. The filter is used to split the kinetic energy equation into three equations corresponding to the three scale ranges described above. The data from a direct simulation of the Rayleigh-Benard problem for conditions where the flow is turbulent are used to calculate the individual terms in the three kinetic energy equations. This is done for a range of filter widths. These results are used to study the spatial location and the scale range of the thermal energy production, the cascading of kinetic energy, the diffusion of kinetic energy, and the energy dissipation. These results are used to evaluate two subgrid models typically used in large-eddy simulations of turbulence. Subgrid models attempt to model the energy below the filter width that is removed by a low-pass filter.

Eidson, T. M.↗

Non-Magnetic, Tough, Corrosion- and Wear-Resistant Knives From Bulk Metallic Glasses and Composites

Quality knives are typically fabricated from high-strength steel alloys. Depending on the application, there are different requirements for mechanical and physical properties that cause problems for steel alloys. For example, diver's knives are generally used in salt water, which causes rust in steel knives. Titanium diver's knives are a popular alternative due to their salt water corrosion resistance, but are too soft to maintain a sharp cutting edge. Steel knives are also magnetic, which is undesirable for military applications where the knives are used as a tactical tool for diffusing magnetic mines. Steel is also significantly denser than titanium (8 g/cu cm vs. 4.5 g/cu cm), which results in heavier knives for the same size. Steel is hard and wear-resistant, compared with titanium, and can keep a sharp edge during service. A major drawback of both steel and titanium knives is that they must be ground or machined into the final knife shape from a billet. Since most knives have a mirrored surface and a complex shape, manufacturing them is complex. It would be more desirable if the knife could be cast into a net or near-net shape in a single step. The solution to the deficiencies of titanium, steel, and ceramic knives is to fabricate them using bulk metallic glasses (or composites). These alloys can be cast into net or near-net shaped knives with a combination of properties that exceed both titanium and steel. A commercially viable BMG (bulk metallic glass) or composite knife is one that exhibits one or all of the following properties: It is based on titanium, has a self-sharpening edge, can retain an edge during service, is hard, is non-magnetic, is corrosion-resistant against a variety of corrosive environments, is tough (to allow for prying), can be cast into a net-shape with a mirror finish and a complex shape, has excellent wear resistance, and is low-density. These properties can be achieved in BMG and composites through alloy chemistry and processing. For each desired property for knife fabrication and performance, there is an alloy development strategy that optimizes behavior. Although BMG knives have been demonstrated as far back as 1995, they never found commercial success because they had to be ground (which presented problems because the alloys contained beryllium), they weren't low cost (because they weren't cast to a net-shape), they were brittle (because they were made with a low-quality commercial material), and they had extremely poor corrosion resistance (because corrosion was not well-understood in these materials). Ultimately, these shortcomings prevented the widespread commercialization. In the current work, the inventors have applied more than a decade of research on BMGs from Caltech and JPL to develop a better understanding of how to make BMG knives that exhibit an optimal combination of properties, processing and cost. Alloys have been developed based in titanium (and other metals), that exhibit high toughness, high hardness, excellent corrosion resistance, no ferromagnetism, edge-retaining selfsharpening, and the ability to be cast like a plastic using commercially available casting techniques (currently used by commercial companies such as Liquidmetal Technologies and Visser Precision Casting). The inventors argue that depending on the application (diving, military, tactical, utility, etc.) there is an optimal combination of design and alloy composition. Moreover, with new casting technologies not available at the inception of these materials, net-shaped knives can be cast into complex shapes that require no aftermarket forming, except for sharpening using water-cooled polishing wheel. These combinations of discoveries seek to make low-cost BMG knives commercially viable products that have no equal among metal or ceramic knives. Current work at JPL focuses on net-shape casting of these alloys and testing their mechanical properties versus commercially available knives to demonstrate their benefits.

Hoffman, Douglas C.↗

A Numerical and Experimental Study of Coflow Laminar Diffusion Flames: Effects of Gravity and Inlet Velocity

In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both computationally and experimentally. A series of flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity conditions with the fuel stream CH4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth vorticity-velocity formulation of the governing equations was employed to describe the reactive gaseous mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional computational domain was employed, where the grid on the axisymmetric domain was a nonuniform tensor product mesh. The governing equations and boundary conditions were discretized on the mesh by a nine-point finite difference stencil, with the convective terms approximated by a monotonic upwind scheme and all other derivatives approximated by centered differences. The resulting set of fully coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton's method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off height, and soot temperature were determined by flame emission images recorded by a digital camera, and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first compared to assess the accuracy of the numerical model. After its validity was established, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of laminar coflow methane-air diffusion flames was explored further by examining quantities derived from the computational results.

microgravity↗

Experimental measurements of the ground cloud growth during the 11 February 1974, Titan-Centaur launch at Kennedy Space Center

The Titan-Centaur was launched from Kennedy Space Center on February 11, 1974 at 0948 eastern daylight time. Ground level effluent measurements were obtained from the solid rocket motors for comparison with NASA diffusion models for predicting effluent ground level concentrations and cloud behavior. The results obtained provide a basis for an evaluation of such key model inputs such as cloud rise rate, stabilization altitude, crosswind growth, volume expansion, and cloud trajectory. Ground level effluent measurements were limited because of changing meteorological conditions, incorrect instrument location, and operational problems. Based on the measurement results, operational changes are defined. Photographs of the ground exhaust clouds are shown. The chemical composition of the exhaust gases was analyzed and is given.

Stewart, R. B.↗

On the structure of contact binaries. I - The contact discontinuity

The problem of the interior structure of contact binaries is reviewed, and a simple resolution of the difficulties which plague the theory is suggested. It is proposed that contact binaries contain a contact discontinuity between the lower surface of the common envelope and the Roche lobe of the cooler star. This discontinuity is maintained against thermal diffusion by fluid flow, and the transition layer is thin to the extent that the dynamical time scale is short in comparison with the thermal time scale. The idealization that the transition layer has infinitesimal thickness allows a simple formulation of the structure equations which are closed by appropriate jump conditions across the discontinuity. The further imposition of the standard boundary conditions suffices to define a unique model for the system once the chemical composition, the masses of the two stars, and the orbital separation are specified.

Shu, F. H.↗