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.
A nonlinear mode-coupling theory of the lower-hybrid-drift instability is presented. It is found that the instability saturates by transferring energy from the growing, long wavelength modes to the damped, short wavelength modes. The saturation energy, mean square of the potential fluctuations, and diffusion coefficient are calculated self-consistently.
For a system of N diffuse, gray body radiation shields which view only adjacent surfaces and space, the net radiation method for enclosures has been used to formulate a system of linear, nonhomogeneous equations in terms of the temperatures to the fourth power of each surface in the coupled system of enclosures. The coefficients of the unknown temperatures in the system of equations are expressed in terms of configuration factors between adjacent surfaces and the emissivities. As an application, a system of four conical radiation shields for a spin stabilized STARPROBE spacecraft has been designed and analyzed with respect to variations of the cone half angles, the intershield spacings, and emissivities.
A first direct measurement of the minority-carrier transit time in a transparent heavily doped emitter layer is reported. The value was obtained by a high-frequency conductance method recently developed and used for low-doped Si. The transit time coupled with the steady-state current enables the determination of the quasi-static charge stored in the emitter and the quasi-static emitter capacitance. Using a transport model, from the measured transit time, the value for the minority-carrier diffusion coefficient and mobility is estimated. The measurements were done using a heavily doped emitter of the Si p(+)-n-p bipolar transistor. The new result indicates that the position-averaged minority-carrier diffusion coefficients may be much smaller than the corresponding majority-carrier values for emitters having a concentration ranging from about 3 x 10 to the 19th per cu cm to 10 to the 20th per cu cm.
This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two-phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection-dominated to diffusion-dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.
Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.
The focus of this research project involved the dry-cast process for polymeric membrane formation, whereby evaporation of solvent from an initially homogeneous polymer/solvent/ nonsolvent solution results in phase separation and the formation of polymer-rich and polymer-lean phases. Under certain conditions the polymer-lean phase gives rise to very large and usually undesirable, tear-drop-shaped pores (size approx. 10 - 50 microns) termed macrovoids (MVs). Although in many cases the presence of MV pores has deleterious effects on membrane performance, there are a number of innovative applications where the presence of such pores is highly desirable. Although researchers have proposed a variety of mechanisms for MV formation over the past three decades, two main hypotheses are currently favored: one asserts that MV growth can be attributed solely to diffusion (the diffusive growth hypothesis), whereas the other states that solutocapillary convection (the SC hypothesis) at the MV interface contributes to growth. The overall goal of this research was to obtain a more comprehensive understanding of the fundamental mechanism of MV growth. This research incorporates a coupled modeling and experimental approach to test a solutocapillary convection hypothesis for the growth of macrovoid pores in polymeric membranes. Specifically, we utilized a modification of the first principles model developed by two of the PIs (ARG and WBK) for dry-cast CA membranes. For the experimental component, two separate and mutually complementary approaches were used to study MV growth. In the first, membranes cast in a zero-g environment aboard the NASA KC-135 aircraft were compared with those cast on the ground to assess the effect of the buoyancy force on membrane morphology and MV size and shape. In the second approach, videomicroscopy flow visualization (VMFV) was utilized to observe MV formation and growth in real time and to assess the effect of surface tension on the MV growth dynamics. As a result of these fundamental studies, our research group advanced a new hypothesis for MV pore development in polymeric membranes.
The propensity for nano-tendril bundle (NTB) vs widespread nano-tendril growth (i.e., fuzz) on polycrystalline W under varying ion energy modulation conditions, from DC to peak-to-peak energy modulation of 42 eV at 13.56 MHz, is correlated with the crystal orientation of the underlying grains. Grains that are vicinal to crystal orientations with high surface diffusivity (e.g., {101} for a body centered cubic crystal structure) exhibit NTB growth at lower ion energy modulation amplitude than grains that are vicinal to low surface diffusivity orientations, such as {100}. Adatom mobility considerations are presented to describe the experimental observations. These results support that surface diffusion or W adatom mobility enhanced by ion bombardment plays a key role in the surface morphology evolution of W under He irradiation.
A formulation of the stability behavior of a finite-length pinch is presented. A general initial perturbation is expressed as a uniformly convergent sum over a complete discrete k set. A variational calculation is then performed, based on the energy principle, in which the end-boundary conditions appear as constraints. The requisite Lagrange multipliers mutually couple the elemental periodic excitations. The resulting extended form of delta-W still admits a proper second-variation treatment so that the minimization and stability considerations of Newcomb remain applicable. Comparison theorems are discussed as is the relevance of this end-effect model to the stability of solar coronal loops.
In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.
This report presents the results from code-to-code comparisons for a prismatic modular high-temperature gas-cooled reactor 350-MW design for steady-state neutron transport calculations. Integral parameters as well as power and flux distributions are provided. Differences between solutions from diffusion and transport operators are analysed. Effects of control rod homogenization between full block and one sixth block are discussed.
Recently, U-shape flammability maps have been constructed showing minimal oxygen vs. flame strain for opposed flame spread in micro-gravity by Olson and Ferkul. The U-shape defines the limiting flammability bounds from radiative extinction and flame blow-off. Here, the minimum of the U corresponds to the minimum possible oxidizer concentration where burning can occur, and is an important quantity of interest for fire safety. While high strain extinction bounds have been well analyzed, low strain radiative extinction has not. To estimate low strain extinction, in this study an analytical theory is developed based on thin flame theory coupled with a heat and mass transfer model for solid fuels. A reaction progress variable based on the Damköhler number is adapted in the theory to account for incomplete combustion at high strain rates and enable the capturing of the full flammability map. The analytical model is compared to a one dimensional numerical model w/ detailed chemical kinetics and coupled radiation heat transfer in planar and spherical geometries. The flammability maps are then qualitatively compared to experimental extinguishment data compiled by Olson and Ferkul for cylindrical rods of PMMA showing similar trends. The results show the newly developed analytics capture the radiative extinction bound compared to the numerical model and qualitatively agrees with microgravity data.
Over the past decade, pressure gain combustion research has promised over 10 percentage-points of increase in power plant thermal efficiency. Alas, to realize such potential gain, one must effectively couple the turbine with the detonation combustor, whose exhaust conditions differ substantially from current state of the art gas turbines. This paper presents a modeling approach that enables a superior thermodynamic cycle with a rotating detonation combustor and a retrofitted gas turbine by including a diffuser downstream of the combustor and by contouring the turbine endwall, while preserving the airfoil geometry. We propose a multi-step optimization strategy, parametrizing the endwall geometry with a few control points, without altering the airfoil geometry, and with the stage turbine efficiency as objective function. In a first step the turbine performance is assessed with steady inlet conditions by solving the steady three-dimensional Reynolds-Averaged Navier-Stokes equations. In a second step, the inlet conditions are unsteady, as predicted from a detonation combustor and a diffuser, and three-dimensional full unsteady simulations are performed with an unsteady Reynolds-Averaged Navier-Stokes solver. By altering the vane endwall, the steady optimization yielded an efficiency increase of 12% relative to the baseline, while the unsteady optimization resulted in 21% increase compared to the datum turbine. Finally, a full engine analysis demonstrated the superiority of pressure gain combustion which included a realistic thermodynamic cycle of the combustor, the diffuser and the optimized turbine components.
Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.
This paper presents results that describe the flow field and surface response conditions for a planetary probe entering the atmosphere of Jupiter. The results are obtained with a viscous-shock-layer analysis which includes mass injection, radiative transfer, diffusion, and viscous effects. Radiative transfer is calculated with an existing nongray radiation model that accounts for molecular band, atomic line, and continuum transitions. Particular emphasis is given to the effect of mass injection on the stagnation and downstream flow. Results for both specified and coupled mass injection are obtained for a carbon heat shield.
While a torsion balance in vacuum can easily measure the momentum transfer from a gas beam impinging on a surface attached to the balance, this measurement depends on the accommodation coefficients of the atoms with the surface and the distribution of the recoil. A torsion balance is described for making absolute flux measurements independent of recoil effects. The torsion balance is a conventional taut suspension wire design and the Young modulus of the wire determines the relationship between the displacement and the applied torque. A compensating magnetic field is applied to maintain zero displacement and provide critical damping. The unique feature is to couple the impinging gas beam to the torsion balance via a Wood's horn, i.e., a thin wall tube with a gradual 90 deg bend. Just as light is trapped in a Wood's horn by specular reflection from the curved surfaces, the gas beam diffuses through the tube. Instead of trapping the beam, the end of the tube is open so that the atoms exit the tube at 90 deg to their original direction. Therefore, all of the forward momentum of the gas beam is transferred to the torsion balance independent of the angle of reflection from the surfaces inside the tube.
Results are reported for a preliminary theoretical study of the coupled mass-, momentum-, and heat-transfer conditions expected within small ampoules used to grow oriented organic solid (OS-) films, by physical vapor transport (PVT) in microgravity environments. It is show that previous studies made restrictive assumptions (e.g., smallness of delta T/T, equality of molecular diffusivities) not valid under PVTOS conditions, whereas the important phenomena of sidewall gas creep, Soret transport of the organic vapor, and large vapor phase supersaturations associated with the large prevailing temperature gradients were not previously considered. Rational estimates are made of the molecular transport properties relevant to copper-phthalocyanine monomeric vapor in a gas mixture containing H2(g) and Xe(g). Efficient numerical methods have been developed and are outlined/illustrated here to making steady axisymmetric gas flow calculations within such ampoules, allowing for realistic realistic delta T/T(sub)w-values, and even corrections to Navier-Stokes-Fourier 'closure' for the governing continuum differential equations. High priority follow-on studies are outlined based on these new results.
A series of numerical calculations are made of a jet in crossflow using progressive mesh refinement up to 2.4 million grid points. A solution algorithm that combined the features of a point-coupled solution of the momentum and continuity equations with multigrid acceleration demonstrated efficient convergence of the governing equations. In comparisons with experimental data, some features of the jet flow field were found to be predicted reasonably well, but the results still remain affected by numerical diffusion.
A slight deviation from Newtonian dynamics can lead to new effects associated with the concept of physical intelligence. Non-Newtonian effects such as deviation from classical thermodynamic as well as quantum-like properties have been analyzed. A self-supervised (intelligent) particle that can escape from Brownian motion autonomously is introduced. Such a capability is due to a coupling of the particle governing equation with its own Liouville equation via an appropriate feedback. As a result, the governing equation is self-stabilized, and random oscillations are suppressed, while the Liouville equation takes the form of the Fokker-Planck equation with negative diffusion. Non- Newtonian properties of such a dynamical system as well as thermodynamical implications have been evaluated.