Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Toward Better Modeling of Supercritical Turbulent Mixing

study was done as part of an effort to develop computational models representing turbulent mixing under thermodynamic supercritical (here, high pressure) conditions. The question was whether the large-eddy simulation (LES) approach, developed previously for atmospheric-pressure compressible-perfect-gas and incompressible flows, can be extended to real-gas non-ideal (including supercritical) fluid mixtures. [In LES, the governing equations are approximated such that the flow field is spatially filtered and subgrid-scale (SGS) phenomena are represented by models.] The study included analyses of results from direct numerical simulation (DNS) of several such mixing layers based on the Navier-Stokes, total-energy, and conservation- of-chemical-species governing equations. Comparison of LES and DNS results revealed the need to augment the atmospheric- pressure LES equations with additional SGS momentum and energy terms. These new terms are the direct result of high-density-gradient-magnitude regions found in the DNS and observed experimentally under fully turbulent flow conditions. A model has been derived for the new term in the momentum equation and was found to perform well at small filter size but to deteriorate with increasing filter size. Several alternative models were derived for the new SGS term in the energy equation that would need further investigations to determine if they are too computationally intensive in LES.

Selle, Laurent↗
DOE OSTI2021

On Two-Phase Pressure and Temperature Equilibration with Mie-Grüneisen Equations of State

In the course of reactive flow modeling, the phases of high explosive (HE) reactants and products introduce additional degrees of freedom to the underlying hydrodynamic equations of compressible reactive fluid mechanics. Hydrodynamic simulations track the bulk (mixture) conservation properties of mass, momentum and energy. When multiple phases are introduced, one must decide on appropriate closure rules to determine the relative phase volumes and energy distribution between the phases. One typical closure is to assume that the phases are in pressure and temperature equilibrium. Here, the requisite equations of p-T equilibration are explored for Mie-Grüneisen Equations of State (EOS). Examples with Davis reactants and Davis products EOSs will be given. The Davis forms are used in both the WSD and AWSD reactive flow models.

42 ENGINEERING↗
DOE OSTI2023

Updated Application of Frequency of Detection Methods for the INL Site Ambient Air Monitoring Network

This report presents a quantitative assessment of the current INL Site air monitoring network using frequency of detection (FD) methods. The first assessment of the INL network was performed in 2015 and made recommendations for improving the network. As a result of changes made in response to the recommendations and the addition of new source locations, the network was modified and reassessed in 2017. Since 2017, administration of the air sampling program has been consolidated under one contractor, which resulted in additional changes to the network (sampler numbers and locations) and changes in radionuclide detection levels. As a result of these changes and others, an updated assessment of the INL Site ambient air monitoring network was performed. The same two exposure scenarios used in previous assessments were used for this assessment: a resident scenario and a shepherd/rancher scenario. The resident was assumed to be continuously present at their residence/business/farm operation outside the INL Site boundary while the shepherd/rancher was assumed to be present 24-hours at the nearest INL grazing allotment boundary in each of the 22.5-degree sectors along the sector centerline from each source. Updates to both the resident and shepherd/rancher receptor locations were included. Other changes include updates to flow rates for stack sources, expansion of the list of important radionuclides based on the most recent National Emission Standards for Hazardous Pollutants (NESHAPs) analysis, and updated dose coefficients. The assessment was conducted to determine whether the current INL monitoring network is capable of detecting releases of important radionuclides from INL Site sources that have the potential to exceed a conservative dose threshold for the two exposure scenarios. The assessment revealed that for the resident scenario, the current network meets the desired performance objective (FD = 95%) for all radionuclides and sources except for Cl-36 from the TRA-770 stack (94.4%). For the shepherd/rancher scenario, the FD performance objective is met for all radionuclides and sources except tritium from MFC-774 and TAN 679 (91% for both). An investigation of reported emissions for the past three years revealed that Cl 36 is not emitted from TRA-770, and routine tritium emissions from MFC-774 and TAN-679 are very small and the sources are likely incapable of emitting enough tritium to cause a release that should be detectable by the monitoring network. This assessment is based on a conservative dose threshold. This coupled with fact that the FD for Cl 36 is only slightly less than the performance objective and Cl-36 is not emitted from TRA-770, modifying the network (i.e. adding another sampler, increasing sampler flow rate, moving samplers) to meet the 95% performance objective for this radionuclide/source/receptor scenario is not warranted. Similarly, because tritium emissions from MFC-774 and TAN-679 are very small and these two sources are likely incapable of causing a dose due to tritium release that should be detectable by the network, modifications to increase tritium detection for these sources is also unwarranted at this time. However, if it is required to meet the performance objective for tritium for all sources and receptor scenarios, additional analysis determined the FD could be raised from 91% to > 99% for both sources by adding two tritium samplers to the network.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Sensitivity Analysis of Transient Critical Heat Flux by RIA Under High-Pressure Flow Boiling Conditions in TRTL

A reactivity-initiated accident (RIA) is a design-basis accident under which critical heat flux (CHF) is likely to be exceeded. The operational margin for RIAs is currently determined using steady-state CHF lookup tables, which provide conservative estimates relative to transient CHF phenomena. The Transient Reactor Test Loop (TRTL) facility at Oregon State University is capable of performing out-of-pile rapid heating experiments representative of a RIA at conditions representative of a pressurized water reactor (PWR). Here, to further our understanding of and ability to predict transient CHF under PWR conditions, we performed a sensitivity analysis on a RELAP5-3D model of the TRTL facility coupled to the RAVEN code framework to define a proposed experimental test matrix to be performed at the TRTL facility. We then implemented a flow boiling CHF correlation into RELAP5-3D and performed a secondary sensitivity analysis inspecting the impact of the built-in RELAP5-3D CHF and heat transfer multipliers on both the prediction of CHF and key safety parameters, such as peak cladding temperature and heat flux. The results show that the multiplier with the highest influence toward the prediction of CHF occurrence and the safety parameters is the transient CHF multiplier. Operational performance envelopes have been developed for each of the test matrix cases and will be used for validation once the experiments are performed. The TRTL facility is currently performing shakedown testing to verify system performance prior to proceeding with the experimental campaign. Restart testing results include pump curve restart testing, pressure tests, and heater rod thermocouple transients.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1991

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form were iterated in time until a steady solution was reached. Evidence was given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at the domain boundaries was proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects were accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. The flow in a developing S-duct was then solved in the laminar regime in a Reynolds number (Re) of 790 and in the turbulent regime at Re equals 40,000 by using the Baldwin-Lomax model. The Stanitz elbow was then solved by using an invicid version of the same code at M sub inlet equals 0.4. Grid dependence and convergence rate were investigated, showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re equals 2.5 times 10(exp 6) was solved with the Baldwin-Lomax and the q-omega models. Both approaches show satisfactory agreement with experiments, although the q-omega model was slightly more accurate.

Michelassi, V.↗
NASA NTRSDate not supplied

Evaluating Changes in Freshwater Health in the Tonlé Sap Basin using Remotely Sensed Data

The Tonlé Sap lake and river basin provide freshwater resources to about 4.5 million people in central Cambodia. Accelerating dam construction, intensifying agriculture, deforestation, and changing climate patterns threaten the health of the lake and the food security of nearby communities. NASA DEVELOP collaborated with Conservation International to calculate multiple components of the ecosystem vitality indicator of the Freshwater Health Index for Tonlé Sap lake from 2000 – 2020 using remote sensing. Landcover classifications, lake level altimetry, groundwater storage and nutrient flows were derived from multiple satellites including Landsat 5, Landsat 7, Landsat 8, Sentinel-2, the Global Precipitation Measurement (GPM)mission, and the Gravity Recovery and Climate Experiment (GRACE). Additionally, the Soil and Water Assessment Tool (SWAT) was used to analyze nutrient flow and water quality in the basin. A breakdown in the volume and regularity of annual lake levels from 2000-2020 was observed, reflecting increased pressure on water supply and agricultural productivity. A decrease in landcover naturalness was also observed and at least 8% of forested areas in the basin were lost while rice harvest intensity increased throughout the study period. Additionally, an initial water quality analysis was demonstrated using SWAT. Continued collaboration with Conservation International and local representatives will guide the remaining water quality analyses and allow them to make informed decisions regarding freshwater management and resource allocation in the Tonlé Sap lake and river basin.

Adriana Le Compte↗
NASA NTRS1977

Numerical turbulence modeling

The underlying bases and developments in two techniques of detailed turbulence modeling are described where the flow is treated in the Eulerian sense, and one technique where the Lagrangian motions of vortices are followed. First, a technique is described for solving the single-point statistically averaged conservation equations. The Reynolds stresses that appear in these equations are evaluated by solving supplemental differential equations which contain terms that are modeled. A sequence of increasingly complex, but also increasingly general, modeling equations is described and computations based on these equations are compared with experimental data. The hierarchy of models described terminates with equations for the individual components of the Reynolds stress tensor. The second Eulerian technique approach to turbulence modeling is the direct numerical simulation of turbulent fields. In this approach, all three dimensional eddies between a predetermined range of sizes are computed in time within a specified volume of flow. Present day computers require a tradeoff between the size of the volume that can be considered and the degree of resolution of the turbulent eddies. Techniques of modeling the smallest eddies are described that permit enlarging the volume, or Reynolds number, that can be considered.

Rubesin, M. W.↗
NASA NTRS1981

A preliminary study of relaxation methods for the inviscid conservative gasdynamics equations using flux splitting

Plus minus flux vector split schemes are combined with the multigrid relaxation method to obtain fast rates of iterative (or steady state) convergence. On model nonlifting biconvex airfoil flow, steady state convergence is obtained using multigrid in 40 iteractions for subsonic cases and in 120 iterations for transonic cases. A 65 x 33 stretched grid is used. However, the multigrid scheme used here (and there are numerous variations) sometimes could not be made to work for transonic cases with a large number of supersonic points. Moreover, the multigrid method is only several times faster than the conventional algorithm which itself could perhaps be accelerated by other means. Nethertheless, with proper adjustment the multigrid method will ultimately lead to very fast rates of iterative convergence.

Steger, J. L.↗
NASA NTRS1984

Numerical study of two-dimensional non-plane MHD wave propagation in a supersonic, superalfvenic magnetohydrodynamic flow

The features of a 2.5-dimensional time-dependent MHD numerical code used to simulate the propagation of finite amplitude MHD waves through an inhomogeneous, supersonic superalfvenic medium are described. Basic equations for conservation of mass, momentum, and free energy in a unit volume plasma gas and for magnetic induction are defined. Initial conditions are functions of the radial coordinates and disturbances are introduced at the lower boundary. A set of finite difference equations based on a Lax-Wendroff scheme is used for the simulation. The model is applied to analyzing a solar flare shock wave in steady-state and global transient conditions while propagating at 1 AU heliolongitude.

Han, S. M.↗
DOE OSTI2024

Path-Integrated Stereo X-Ray Digital Image Correlation: Resolving the Violation of Conservation of Intensity

In this study, x-ray imaging addresses many challenges with visible light imaging in extreme environments, where optical diagnostics such as digital image correlation (DIC) and particle image velocimetry (PIV) suffer biases from index of refraction changes and/or cannot penetrate visibly occluded objects. However, conservation of intensity—the fundamental principle of optical image correlation algorithms—may be violated if ancillary components in the X-ray path besides the surface or fluid of interest move during the test. The test series treated in this work sought to characterize the safe use of fiber-epoxy composites in aerospace and aviation industries during fire accident scenarios. Stereo X-ray DIC was employed to measure test unit deformation in an extreme thermal environment involving a visibly occluded test unit, incident surface heating to temperatures above 600°C, and flames and soot from combusting epoxy decomposition gasses. The objective of the current work is to evaluate two solutions to resolve the violation of conservation of intensity that resulted from both the test unit and the thermal chamber deforming during the test. The first solution recovered conservation of intensity by pre-processing the path-integrated X-ray images to isolate the DIC pattern of the test unit from the thermal chamber components. These images were then correlated with standard, optical DIC software. The second solution, called Path-Integrated Digital Image Correlation (PI-DIC), modified the fundamental matching criterion of DIC to account for multiple, independently-moving components contributing to the final image intensity. The PI-DIC algorithm was extended from a 2D framework to a stereo framework and implemented in a custom DIC software. Both solutions accurately measured the cylindrical shape of the undeformed test unit, recovering radii values of R = 76.20±0.12 mm compared to the theoretical radius of R theor = 76.20 mm. During the test, the test unit bulged asymmetrically as decomposition gasses pressurized the interior and eventually burned in a localized jet. Both solutions measured the heterogeneous radius of this bulge, which reached a maximum of approximately R = 91 mm, with a discrepancy of 2–3% between the two solutions. Two solutions that resolve the violation of conservation of intensity for path-integrated X-ray images were developed. Both were successfully employed in a stereo X-ray DIC configuration to measure deformation of an aluminum-skinned, fiber-epoxy composite test unit in a fire accident scenario.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

The European wood pellets for heating market - Price developments, trade and market efficiency

Competitive international markets imply adjustments towards competitive spatial equilibrium in which excess from one market is transferred to another and prices are equilibrated except for remaining differences that can be assigned to transfer costs. The European market for wood pellets used in small-scale heating systems has been expanding significantly over the past decade. Small scale pellet heating is arguably a mature technology, but whether the market is mature is another question. In this paper we analyse recent data on trade flows and price developments between Italy, Austria, Germany and France to understand the developments of wood pellet market efficiency and to draw conclusions about market function. The objective of this study is to establish a framework to test the European residential wood pellet market for competitive spatial equilibrium using modern trade theory. We find mainly inefficiently integrated markets with remaining positive marginal profits and detectable arbitrageurs’ activity. Based on a thorough discussion of these findings and the underlying data we outline possible methodology advancements and list policy recommendations to secure access and affordability of this renewable heating commodity in the long run.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Computational fluid dynamic modeling to determine the indoor environment of an electron-ion collider service building

The design for the Electron-Ion Collider (EIC) calls for several service buildings that house various power supplies and control electronics for the collider ring itself. In order to operate within specified conditions, the ambient air entering the power supplies needs to be within a certain temperature range while dissipating the heat from losses. Proper cooling is therefore a necessity in the service buildings to ensure that every aspect of the EIC works as intended. Since the EIC is in the design stage, we are evaluating the indoor environment of the service building using the current design specifications. We have researched multiple forms of literature and performed the necessary calculations to compile a list of boundary conditions that accurately represent the situation at hand. We are using computational fluid dynamics modeling to solve the conservation equations for mass, momentum, and energy (Navier-Stokes). This allows us to perform a finite element analysis which will give us the flow distribution in the room as well as temperature profiles throughout the building. We have obtained a simulation result giving us the temperature profiles for the building and it shows that the placement of the racks and supply vents are essential to obtaining an even temperature distribution. This model will provide a basis for design decisions which will affect the overall cooling of the service building without extending the schedule and avoiding a costly reworking of the cooling system.

43 PARTICLE ACCELERATORS↗
NASA NTRS1981

Shockless design and analysis of transonic blade shapes

A fast computer program was developed to eliminate the shocks by slightly altering portions of the contour of a given airfoil in the cascade. The program can be used in two basic modes: (1) An analysis for steady, transonic, potential flow through a given planar cascade of airfoils and (2) a design for converting a given cascade into a shockless transonic cascade. The design mode can automatically be followed by the analysis mode, which confirms that the flow field is shock free. The program generates its own multilevel boundary conforming computational grids and solves a full potential equation in a fully conservative form. The shockless design is performed by implementing Sobieczky's fictitious-gas elliptic continuation concept.

Dulikravich, D. S.↗
NASA NTRS1982

Impact mechanics of the Cretaceous-Tertiary extinction bolide

An examination of the mechanics of asteroidal, cometary, and meteor swarm impact on the earth determined if the enrichment of projectile material in the K-T layer is consistent with melts and impact breccias on the earth and moon, the size of the impacters, the distribution of the kinetic energy, and the sequence of impacts that could give rise to observed extinction phenomena. Flows resulting from spherical projectile impacts onto layers of air, water, and silicates were modeled and Eulerian finite difference algorithms were employed to solve conservation equations and equations of state. A range of speeds and impacter densities were considered, along with sizes from 0.17 km, which would be consumed in the atmosphere, to a 10 km object, which would have had a diameter greater than a reference 7.1 km atmosphere depth. It is concluded that an impact of the K-T bolide could result in global biotic extinction and worldwide material deposition.

Okeefe, J. D.↗
NASA NTRS1982

Small particle transport across turbulent nonisothermal boundary layers

The interaction between turbulent diffusion, Brownian diffusion, and particle thermophoresis in the limit of vanishing particle inertial effects is quantitatively modeled for applications in gas turbines. The model is initiated with consideration of the particle phase mass conservation equation for a two-dimensional boundary layer, including the thermophoretic flux term directed toward the cold wall. A formalism of a turbulent flow near a flat plate in a heat transfer problem is adopted, and variable property effects are neglected. Attention is given to the limit of very large Schmidt numbers and the particle concentration depletion outside of the Brownian sublayer. It is concluded that, in the parameter range of interest, thermophoresis augments the high Schmidt number mass-transfer coefficient by a factor equal to the product of the outer sink and the thermophoretic suction.

Rosner, D. E.↗
NASA NTRS1984

Metric-discontinuous zonal grid calculations using the Osher scheme

Computations on zonal grids - in particular, grids with metric discontinuities resulting from the interspersion of highly clustered regions with coarse regions - are possible using a fully conservative form of the Osher upwind scheme. These zonal grids can result from an abrupt clustering of points near solution discontinuities or near other flow features that require improved resolution. The zonal approach is shown to capture shocks with almost 'shock-fitting' quality but with minimal effort. Results for inviscid flow, including quasi-one-dimensional nozzle flow, supersonic flow over a cylinder, and blast-wave diffraction by a ramp, are presented. These calculations demonstrate the powerful capabilities of the Osher scheme used in conjunction with zonal grids in simulating flow fields with complex shock patterns.

Rai, M. M.↗
NASA NTRS1985

Theoretical and experimental flow fields for a supersonic cruise fighter forebody

The capability of two numerical methods to predict the flow field about a representative supersonic cruise aircraft was examined. The two codes were a small-disturbance transonic program of Boppe and a conservative form full-potential transonic program of Shankar and Szema. For the former code, comparisons were made with wind-tunnel data at Mach numbers of 0.6, 0.9, and 1.2 and angles of attack of 0 deg, 5 deg, and 10 deg (7.5 deg instead of 10 deg at a Mach number of 1.2). Predictions from the two codes were compared at a Mach number of 1.2 and an angle of attack of 7.5 deg. The comparison criteria were contours of local angle of attack, local angle of sideslip, and local Mach number. The comparisons indicated that both codes may be considered useful for design applications, depending on the degree of accuracy required by the user's solution. Both solutions show an inaccuracy in their predictions, particularly as Mach number and/or angle of attack increases, because of their lack of viscous effects and any mechanisms to predict vortex development.

Yaros, S. F.↗
NASA NTRS1986

Computational simulation of liquid rocket injector anomalies

A computer model has been developed to analyze the three-dimensional two-phase reactive flows in liquid fueled rocket combustors. The model is designed to study the influence of liquid propellant injection nonuniformities on the flow pattern, combustion and heat transfer within the combustor. The Eulerian-Lagrangian approach for simulating polidisperse spray flow, evaporation and combustion has been used. Full coupling between the phases is accounted for. A nonorthogonal, body fitted coordinate system along with a conservative control volume formulation is employed. The physical models built into the model include a kappa-epsilon turbulence model, a two-step chemical reaction, and the six-flux radiation model. Semiempirical models are used to describe all interphase coupling terms as well as chemical reaction rates. The purpose of this study was to demonstrate an analytical capability to predict the effects of reactant injection nonuniformities (injection anomalies) on combustion and heat transfer within the rocket combustion chamber. The results show promising application of the model to comprehensive modeling of liquid propellant rocket engines.

Przekwas, A. J.↗