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 Solid Oxide Fuel Cell (SOFC) system having a hot zone with a center cathode air feed tube for improved reactant distribution, a CPOX reactor attached at the anode feed end of the hot zone with a tail gas combustor at the opposing end for more uniform heat distribution, and a counter-flow heat exchanger for efficient heat retention.
Abstract We establish and develop a novel methodology to treat higher-order non-linear effects of gravitational radiation that is scattered from binary inspirals, which employs modern scattering-amplitudes methods on the effective picture of the binary as a composite particle. We spell out our procedure to study such effects: assembling tree amplitudes via generalized-unitarity methods and employing the closed-time-path formalism to derive the causal effective actions, which encompass the full conservative and dissipative dynamics. We push through to a new state of the art for these higher-order effects, up to the third subleading tail effect, at order$$ {G}_N^5 $$ G N 5 and the 5-loop level, which corresponds to the 8.5PN order. We formulate the consequent dissipated energy for these higher-order corrections, and carry out a renormalization analysis, where we uncover new subleading RG flow of the quadrupole coupling. For all higher-order tail effects we find perfect agreement with partial observable results in PN and self-force theories, where available.
As a follow-up on the drift-kinetic study of the non-local bootstrap current in the steep edge pedestal of tokamak plasma by Koh et al. [Phys. Plasmas 19, 072505 (2012)], a gyrokinetic neoclassical study is performed with gyrokinetic ions and drift-kinetic electrons. Besides the gyrokinetic improvement of ion physics from the drift-kinetic treatment, a fully non-linear Fokker-Planck collision operator—that conserves mass, momentum, and energy—is used instead of Koh et al.'s linearized collision operator in consideration of the possibility that the ion distribution function is non-Maxwellian in the steep pedestal. An inaccuracy in Koh et al.'s result is found in the steep edge pedestal that originated from a small error in the collisional momentum conservation. The present study concludes that (1) the bootstrap current in the steep edge pedestal is generally smaller than what has been predicted from the small banana-width (local) approximation [e.g., Sauter et al., Phys. Plasmas 6, 2834 (1999) and Belli et al., Plasma Phys. Controlled Fusion 50, 095010 (2008)], (2) the plasma flow evaluated from the local approximation can significantly deviate from the non-local results, and (3) the bootstrap current in the edge pedestal, where the passing particle region is small, can be dominantly carried by the trapped particles in a broad trapped boundary layer. A new analytic formula based on numerous gyrokinetic simulations using various magnetic equilibria and plasma profiles with self-consistent Grad-Shafranov solutions is constructed.
In this paper, the Navier-Stokes (NS) equations are approximated with a reduced NS system, that represents the lowest-order terms in an asymptotic Re expansion. This system allows for simplified boundary conditions, more generality in the location of the outflow boundary, and insures mass conservation in all subdomain grid interfaces, and also at the outflow boundary. The higher-order NS diffusion terms are included through a deferred corrector, in chosen subdomains, when required.
Centered finite-difference schemes are commonly used for high-fidelity turbulent flow simulations in canonical configurations because of their non-dissipative property and computational efficiency. However, their use in flow simulations over complex geometries is limited by the requirements of a structured grid and a stable boundary treatment in the absence of artificial (numerical) dissipation. Cartesian embedded boundary (EB) approaches provide an efficient structured-grid framework to apply difference schemes over complex domains. However, they are often restricted to low orders of accuracy because of numerical instabilities at the embedded boundaries and the issues of small-cell problem that are difficult to address with high-order accuracy. The present work discusses a systematic approach to obtain high-order EB methods with non-dissipative centered schemes in the interior. This approach, based on satisfying the primary and secondary conservation conditions, is employed to derive EB schemes that are up to sixth-order accurate in the interior and fourth-order accurate globally for hyperbolic, parabolic as well as incompletely parabolic problems. The proposed finite-difference discretization is, by construction, dimensionally split and addresses the small-cell problem without any cell/geometry transformations, thus, highly simplifying implementation in a flow solver. Various linear and non-linear numerical tests are performed to evaluate the stability and the accuracy of the proposed EB schemes.
This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.
Numerical simulation is essential for the prediction and evaluation of hydrocarbon reservoir performance. Numerical simulators developed for the description of the behavior of hydrates under production and the corresponding flow of fluids and heat accounting for all known processes are powerful, but they need validation through comparison to field or experimental data in order to instill confidence in their predictions. In this study, we analyze by means of numerical simulation the results of an experiment of methane hydrate dissociation by thermal stimulation in unconsolidated porous media heated through the vessel walls. The physics captured by the model include multicomponent heat and mass transfer, multiphase flow through porous media, and the phase behavior of the CH 4 + H 2 O system involved in methane hydrate formation and dissociation. The set of governing equations consists of the mass and energy conservation equations coupled with constitutive relationships, i.e., the dissolution of gas in H 2 O, relative permeability and capillary pressure models, composite thermal conductivity models, and methane hydrate phase equilibria. The model geometry describes accurately the hydrate reactor used in a recent experimental study investigating methane hydrate dissociation behavior [Chong et. al. Appl. Energy 2016, 177, 409-421]. The cumulative gas production is estimated and validated against three tests of experimental data involving different boundary temperatures, showing a good agreement between observations and numerical predictions. The predicted evolution of the spatial distributions of different phases over time shows that hydrate dissociation progresses inward from the reactor boundary to the center, methane gas accumulates to the top of the reactor because of buoyancy, and water migrates down to the bottom of the reactor because of gravity. A sharp hydrate dissociation front is predicted, and the estimated location of hydrate dissociation front suggests a linear relationship with the square root of time. A sensitivity analysis on the thermal conductivity of sand under fully saturated conditions is conducted to elucidate its effect on the gas production behavior. In addition, the energy efficiency ratio computed from the simulation of this boundary-wall heating technique varies from 14.0 to 16.2. Deviations between observations and predictions of the evolution of the temperature profile are attributed to initial heterogeneous distribution of the hydrate phase in the hydrate reactor.
We propose an algorithm for low Mach number reacting flows subjected to electric field that includes the chemical production and transport of charged species. This work is an extension of a multi-implicit spectral deferred correction (MISDC) algorithm designed to advance the conservation equations in time at scales associated with advective transport. The fast and nontrivial interactions of electrons with the electric field are treated implicitly using a Jacobian-Free Newton Krylov approach for which a preconditioning strategy is developed. Within the MISDC framework, this enables a close and stable coupling of diffusion, reactions and dielectric relaxation terms with advective transport and is shown to exhibit second-order convergence in space and time. The algorithm is then applied to a series of steady and unsteady problems to demonstrate its capability and stability. Although developed in a one-dimensional case, the algorithmic ingredients are carefully designed to be amenable to multidimensional applications.
DiscoFluxM is software for modeling the deformation behavior of crystalline materials within the third-party finite-element based software, MOOSE. The software enables full-field crystal plasticity finite element simulations of deformation processes of metals and other crystalline materials. DiscoFluxM is an implementation of the DiscoFlux theoretical and algorithmic framework for coupling single-crystal plasticity with dislocation transport. The continuum mesoscale modeling framework is founded upon a nonlocal dislocation-density based crystal plasticity theory. The nonlocal theory couples continuum dislocation transport with an otherwise spatially-local single-crystal model employing nonlinear thermoelasticity and crystallographic plasticity. Dislocation transport is modeled by enforcing dislocation conservation at a slip-system level through the solution of advection-diffusion equations. The configuration of geometrically necessary dislocation density gives rise to a back-stress that inhibits or accentuates the flow of dislocations.
The results of a theoretical investigation of thermodynamic effects on developed cavitation are presented. An approximate solution to the conservation equations for a two-phase laminar boundary layer is obtained. This analysis produces an expression for the temperature difference between the liquid and vapor phases which can be applied to developed cavity flows. Experimental data of cavity temperature depressions are correlated using this result. In addition, a theoretical estimate of the Nusselt number for the cavity is made using a turbulent boundary layer cavity model proposed by Brennen. The result agrees in part with empirically determined expressions for the cavity Nusselt number.
A revised version of a split-velocity method for numerical calculation of compressible duct flow was developed. The revision incorporates balancing of mass flow rates on each marching step in order to maintain front-to-back continuity during the calculation. The (checkerboard) zebra algorithm is applied to solution of the three-dimensional continuity equation in conservative form. A second-order A-stable linear multistep method is employed in effecting a marching solution of the parabolized momentum equations. A checkerboard successive overrelaxation iteration is used to solve the resulting implicit nonlinear systems of finite-difference equations which govern stepwise transition.
Numerical studies of time dependent cosmic ray shock structures in planar geometry are interesting because analytical time-independent solutions are available which include the non-linear reactions on the plasma flow. A feature of these time asymptotic solutions is that for higher Mach numbers (M approximately 5) and for a low cosmic ray upstream pressure the solution is not uniquely determined by the usual conservation laws of mass, momentum and energy. These numerical solutions clearly indicate that much work needs to be done before we understand shock acceleration as a time dependent process. The slowness of the process is possibly due to the fact that there is a diffusive flux into the downstream region in addition to the usual advective losses. Analytic investigations of thie phenomenon are required.
Advances in the mechanics of boundary layer flow are reported. The physical problems of large scale coherent structures in real, developing free turbulent shear flows, from the nonlinear aspects of hydrodynamic stability are addressed. The presence of fine grained turbulence in the problem, and its absence, lacks a small parameter. The problem is presented on the basis of conservation principles, which are the dynamics of the problem directed towards extracting the most physical information, however, it is emphasized that it must also involve approximations.
A theory is presented which predicts tides in turbulent, self-gravitating, and loading oceans possessing linearized bottom friction, realistic bathymetry, and continents (at coastal boundaries no-flow conditions are imposed). The theory is phrased in terms of spherical harmonics, which allows the tide equations to be reduced to linear matrix equations. This approach also allows an ocean-wide mass conservation constraint to be applied. Solutions were obtained for 32 long and short period luni-solar tidal constituents (and the pole tide), including the tidal velocities in addition to the tide height. Calibrating the intensity of bottom friction produces reasonable phase lags for all constituents; however, tidal amplitudes compare well with those from observation and other theories only for long-period constituents. In the most recent stage of grant research, traditional theory (Liouville equations) for determining the effects of angular momentum exchange on Earth's rotation were extended to encompass high-frequency excitations (such as short-period tides).
A Gauss-Seidel implicit aerothermodynamic Navier-Stokes computational CFD code with thermochemical surface conditions (GIANTS) and with appropriate boundary conditions was developed for the analysis of hypersonic arcjet flows over large angle blunt bodies. Two-dimensional axisymmetric multicomponent full Navier-Stokes equations were solved, together with boundary conditions, to account for the surface slip, surface catalysis, and energy and mass conservation at the body surface. It is shown that the predictions for isothermal, noncatalytic, and no-slip surface conditions which are made using the GIANTS code with simple boundary conditions agree with those of the LAURA code, and the predicted normalized heating distributions over the forebody surface with boundary conditions for various cone angles are consistent with the trends of data from arcjet experiments.
PHASTRAN analyzes heat-transfer and flow behaviors of materials undergoing phase changes. Many phase changes operate over range of accelerations or effective gravitational fields. To analyze such thermal systems, it is necessary to obtain simultaneous solutions for equations of conservation of energy, momentum, and mass, and for equation of state. Written in APL2.
In this study, the NCAR CLM version 2.0 land-surface model was integrated into the NASA/NCAR fvGCM. The CLM was developed collaboratively by an open interagency/university group of scientists and based on well-proven physical parameterizations and numerical schemes that combine the best features of BATS, NCAR-LSM, and IAP94. The CLM design is a one-dimensional point model with 1 vegetation layer, along with sub-grid scale tiles. The features of the CLM include 10-uneven soil layers with water, ice, and temperature states in each soil layer, and five snow layers, with water flow, refreezing, compaction, and aging allowed. In addition, the CLM utilizes two-stream canopy radiative transfer, the Bonan lake model and topographic enhanced streamflow based on TOPMODEL. The DAO fvGCM uses a genuinely conservative Flux-Form Semi-Lagrangian transport algorithm along with terrain- following Lagrangian control-volume vertical coordinates. The physical parameterizations are based on the NCAR Community Atmosphere Model (CAM-2). For our purposes, the fvGCM was run at 2 deg x 2.5 deg horizontal resolution with 55 vertical levels. The 10-year climate from the fvGCM with CLM2 was intercompared with the climate from fvGCM with LSM, ECMWF and NCEP. We concluded that the incorporation of CLM2 did not significantly impact the fvGCM climate from that of LSM. The most striking difference was the warm bias in the CLM2 surface skin temperature over desert regions. We determined that the warm bias can be partially attributed to the value of the drag coefficient for the soil under the canopy, which was too small resulting in a decoupling between the ground surface and the canopy. We also discovered that the canopy interception was high compared to observations in the Amazon region. A number of experiments were then performed focused on implementing model improvements. In order to correct the warm bias, the drag coefficient for the soil under the canopy was considered a function of LAI (Leaf Area Index). Analysis of the results revealed that there was a substantial impact, and the warm and dry bias in the CLM2 was significantly reduced. For the interception scheme, the canopy throughfall was increased to allow for more infiltration of precipitation into the soil, resulting in increased low-level moisture and a decrease in the interception loss ratio (canopy evaporation to precipitation).