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 991 records · Page 55

Seasonal mixing from intermittent flow drives concentration-discharge behaviour in a stream affected by coal mine drainage

Abandoned mining operations continue to severely degrade many ecosystems worldwide by releasing acidic water and/or heavy metals into surface and groundwater. Contaminant concentrations in affected streams vary with discharge in patterns that reflect both geochemical reactions and variable mixing of contaminated and non-contaminated waters. However, controls on concentration-discharge (C-Q) patterns remain unclear, particularly for constituents that experience changing solubility across redox and pH gradients. Understanding the C-Q behaviour of contaminants aids in predicting both downstream transport and effects on aquatic life under variable flow. Here, we examined the C-Q behaviours of non-reactive (Na, K, Ca, Mg, Cl - ) and reactive (Fe, Mn, Al, H + , SO 4 2- ) solutes in a stream contaminated with acid mine drainage in northeastern Ohio, USA. Concentration-discharge patterns at the watershed outlet primarily reflected mixing of contaminated baseflow with intermittent inputs of high pH water draining from a passive limestone treatment system into the stream. The treatment system acted as an ephemeral tributary that mitigated contamination in the stream by diluting solutes, raising pH, and driving metal precipitation, but only when flow was present during wet seasons. Consequently, AMD-derived reactive solutes (H + , Fe, Mn, Al) decreased with increasing stream discharge while relatively conservative solutes (e.g., Ca, Mg, K, Na) decreased only slightly or were chemostatic. This study highlights both the unique C-Q patterns of reactive solutes when compared to those of non-reactive solutes and the potential for intermittent streams to control C-Q behaviour in headwater catchments.

58 GEOSCIENCES↗

Krypton-81 dating of the deep Continental Intercalaire aquifer with implications for chlorine-36 dating

Deep groundwater samples from the Continental Intercalaire (CI) aquifer in the Northern Tunisian Sahara have been analyzed for noble gases 3 He, 4 He, Ne and 81 Kr, and for 14 C to better constrain the groundwater residence time of this large transboundary aquifer. Its significant radiogenic 4 He content and background-level 14 C both indicate water older than a few tens of thousands of years. Distinct helium concentrations and Ne/He ratios suggest different groundwater flow paths through the Tozeur and Kebili regions. 81 Kr is applied for the first time on the CI aquifer, providing more direct evidence for the presence of old groundwater with 81 Kr/Kr ratios at 63% – 15% of the atmospheric value, corresponding to residence times of 150 – 630 ka. These newly obtained 81 Kr ages enable a re-evaluation of the previously reported data from the Tozeur region in order to better constrain the input values of 36 Cl/Cl and Cl. The initial 36 Cl/Cl ratio is found to be 4 – 6 times larger than previously assumed, and the initial Cl concentration 6 – 15 times smaller. The results of 81 Kr dating and the recalibrated 36 Cl dating define an eastward age progression across > 500 km of flow path from the Algerian to Tunisian sections. Here, this study exemplifies the utility of 81 Kr not only as an independent and conservative chronometer, but also as a supplemental measurement to better calibrate the parameters required for other age tracers, including cosmogenic 36 Cl and radiogenic 4 He.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Neoclassical transport in strong gradient regions of large aspect ratio tokamaks

We present a new neoclassical transport model for large aspect ratio tokamaks where the gradient scale lengths are of the size of the ion poloidal gyroradius. Previous work on neoclassical transport across transport barriers assumed large density and potential gradients but a small temperature gradient, or neglected the gradient of the mean parallel flow. Using large aspect ratio and low collisionality expansions, we relax these restrictive assumptions. We define a new set of variables based on conserved quantities, which simplifies the drift kinetic equation whilst keeping strong gradients, and derive equations describing the transport of particles, parallel momentum and energy by ions in the banana regime. The poloidally varying parts of density and electric potential are included. Studying contributions from both passing and trapped particles, we show that the resulting transport is dominated by trapped particles. We find that a non-zero neoclassical particle flux requires parallel momentum input which could be provided through interaction with turbulence or impurities. We derive upper and lower bounds for the energy flux across a transport barrier in both temperature and density and present example profiles and fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗

Development of fracture diagnostic methods for fluid distribution based on quantitative interpretation of distributed acoustic sensing and distributed temperature sensing

Multistage hydraulic fracturing design on horizontal wells has significantly evolved with larger fluid volume, more fracturing stages, and tighter perforation cluster spacing to efficiently stimulate unconventional reservoirs. From the published field observations, the recent fracturing design results in complex fracture networks or swarm of fractures. Fracture treatment evaluation is extremely challenging in such a case because of the large amount of variables in well completion and stimulation design. Combined measurements from different technologies can help in fracture diagnosis. Fluid distribution, either during fracture injection or during production, directly relates to the stimulation efficiency at the cluster level and at the stage level. Because it is unlikely in the real world to distribute the injected fluid uniformly among all the clusters, we need diagnostic techniques to generate the flow profile along a lateral. Fiber-optic measurements, such as distributed acoustic sensing (DAS) and distributed temperature sensing (DTS), are currently used to diagnose downhole flow conditions. This technology allows us to qualitatively confirm the fluid flow profile and other issues occurring downhole during fracturing such as leakage through plugs. For optimizing a fracturing design, we also need to understand how the design parameters are correlated with the stimulation efficiency. In this study, we combine two sets of models of DAS and DTS data interpretation for injected fluid volume distribution. The DAS is interpreted based on an empirical correlation between fluid flow rates and frequency band energy from the acoustic signals. The DTS is interpreted by performing temperature history match-based thermal energy conservation. Because of the completely different physics behind the interpretations, the confirmation of two interpretations provides confidence in fluid distribution.

Geochemistry & Geophysics↗

Pulse Detonation Engine for Advanced Oxy-Combustion of Coal-Based Fuel for Direct Power Extraction Applications

Pressure gain combustion and magnetohydrodynamic (MHD) systems have the potential to provide a step increase in the efficiency of combined-cycle power plants. Specific advantages include a net pressure increase to the system instead of a pressure drop, the high temperature of the detonation waves can increase the efficiency of power extraction compared to other processes, significant thermal energy can be released in a compact region, and the high velocities of the flow increase extraction of electrical power. In summary, a pressure gain combustor coupled with a MHD has the potential to be transformative. Despite the potential advantages, relatively little research has been conducted considering coupled pressure gain combustion systems with MHD systems. With this background and motivation, the overall goal of this effort was to advance the knowledge, technology, and computational tools associated with coupled detonation and MHD systems. A joint experimental and computational approach was used while seeking to accomplish the goals of this work. Specifically, two pulse detonation engines were developed and used for the experiments to produce detonations. Detonation speeds were measured for a variety of flow and fuel conditions (e.g., methane, propane, with coal particles). Preliminary electrical conductivity measurements were collected. An extensive amount of research was performed to identify sensitivities of detonation behaviors to the presence of combustion products. Computationally, a twofold approach was used in this work. First, a solver was developed for solving the governing equations for a reactive flow with coupled detonation and ionization chemistry. The solver was applied to study the impacts of seed material ionization on detonation. Second, a conservation element-solution element (CE-SE) based numerical solver for detonation studies with a reduced reaction mechanism for oxy-methane combustion was developed and verified on standard test cases. Key findings and contributions from this work are as follows. A system was developed for injecting powderized coal, or other seeding material, into a pulse-detonation system. The influence of a combustion product (i.e., CO 2 ) on detonation behavior was identified. Knowledge gained from this work is applicable to devices such as rotating detonation engines, where combustion products mix with fresh reactants. A system for measuring the electrical conductivity of the exhaust from a pulse-detonation engine was developed. The open-source solver, Clawpack, was extended to solve the reactive Euler equations for simulating detonations. A coupled combustion and ionization chemistry was developed in a single chemical kinetic model for methane oxidation. This model can be used to solve coupled MHD and detonation simulations. It was found that parasitic interactions from ionization chemistry with the magnetic field can reduce the detonation velocity by up to 8%, with a potential impact on power extraction of 15%. It is recommended that interactions between the detonation front and MHD field be considered. Higher gas temperatures and velocities were achieved owing to oxy-fuel detonations. Use of radical dissociation reactions in the reduced reaction mechanism, was found to be critical in predicting detonation temperature and velocity accurately.

01 COAL, LIGNITE, AND PEAT↗

Full potential solution of a transonic quasi-3-D flow through a cascade using artificial compressibility

A reliable method is presented for calculating the flowfield about a cascade of arbitrary 2-D airfoils. The method approximates the three-dimensional flow in a turbomachinery blade row by correcting for streamtube convergence and radius change in the throughflow direction. The method is a fully conservative solution of the full potential equation incorporating the finite volume technique on a body-fitted periodic mesh, with an artificial density imposed in the transonic region to ensure stability and the capture of shock waves. Comparison of results for several supercritical blades shows good agreement with their hodograph solutions. Other calculations for these profiles as well as standard NACA blade sections indicate that this is a useful scheme for analyzing both the design and off-design performance of turbomachinery blading.

Farrell, C.↗

The evolution of Ertel's potential vorticity during stratospheric sudden warmings

In the winter stratosphere of the Northern Hemisphere, the disruption of the westerly vortex and associated warming of polar latitudes is a well known phenomenon. It has become apparent that some important dynamical processes in the stratosphere are highly nonlinear and are best thought of locally rather than in terms of the interaction between the zonal-mean flow and harmonic waves around latitude circles. The importance of nonlinear processes was suggested by McIntyre and Palmer (1983, 1984) who used isentropic maps of Ertel's Potential Vorticity to show that during disturbed episodes material lines may become strongly and irreversibly deformed in certain places. They adopted the term planetary wave breaking to describe this process. Isentropic maps of Q are used to follow the evolution of a Canadian warming in November - December 1981 and a particularly strong warming in January 1982. The advection of Q over large distances on isentropic surfaces was a striking feature of the flow during each event. This could be identified because of our ability to follow the movement of material lines due to the approximate conservation of Q over several days. The advection of Q was a nonlinear process because its changing distribution affected the advecting wind field. The Canadian warming did not lead to a permanent change in the structure of the westerly vortex, as defined by the coarse-grain field of Q, whereas the January event was accompanied by a substantial loss of resolved Q which was never fully recovered.

Fairlie, T. D. A.↗

Finite difference simulation of unsteady interactive flows

The development of time-dependent numerical simulations of unsteady interactive flows of an aerodynamic nature is reviewed with emphasis on compressible flows at flight Reynolds numbers and noniterative schemes based on Navier-Stokes equations. The importance of writing the equations in strong conservation-law form for a generalized body-oriented coordinate system is pointed out. The discussion covers time and length scales and numerical methods currently in use. Some computed results are presented and compared with experimental data.

Deiwert, G. S.↗

Towards a rational theory for CFD global stability

The fundamental notion of the consistent stability of semidiscrete analogues of evolution PDEs is explored. Lyapunov's direct method is used to develop CFD semidiscrete algorithms which yield the TVD constraint as a special case. A general formula for supplying dissipation parameters for arbitrary multidimensional conservation law systems is proposed. The reliability of the method is demonstrated by the results of two numerical tests for representative Euler shocked flows.

Baker, A. J.↗

Numerical Modeling of Pressurization of a Propellant Tank

An unsteady finite volume procedure has been developed to predict the history of pressure, temperature and mass flow rate of the pressurant and propellant during the expulsion of the propellant from a tank. The time dependent mass, momentum and energy conservation equations are solved at the ullage space. The model accounts for the change in the ullage volume due to expulsion of the propellant. It also accounts for the heat transfer from the tank wall and propellant to the ullage gas. The procedure was incorporated in the Generalized Fluid System Simulation Program (GFSSP). The results of several test cases were then compared with a published correlation of pressurant requirements for a given displacement of propellant. The agreement between the predictions and the correlation was found to be satisfactory.

Majumdar, Alok↗

Numerical Modeling of Pressurization of a Propellant Tank

An unsteady finite volume procedure has been developed to predict the history o pressure, temperature and mass flow rate of the pressurant and propellant during the expulsion of the propellant from a tan. The time dependent mass, momentum and energy conservation equations are solved at the ullage space. The model accounts for the change in the ullage volume due to expulsion of the propellant. It also accounts for the heat transfer from the tank wall and propellant to the ullage gas. The procedure was incorporated in the Generalized Fluid System Simulation Program (GFSSP). The results of several test cases were then compared with a published correlation of pressurant requirements for a given displacement of propellant. The agreement between the predictions and the correlation was found to be satisfactory.

Majumdar, Alok↗

Investigation of Transitional Shock-Wave/Boundary Layer Interactions Using Direct Numerical Simulations

Interaction of a transitional boundary layer with a normal-shock is investigated using unstructured tetrahedral meshes under the numerical framework of the space-time conservation element, solution element (CESE) method. The computations mimic recent experimental efforts at the University of Tennessee Space Institute, where a Mach 2.0 flow interacts with a tall cylinder attached to a flat plate. The location of the cylinder with respect to the flat plate leading edge determines if the incoming boundary layer is laminar, transitional or fully turbulent. Four representative flow conditions exemplifying laminar and transitional boundary layers are analyzed by direct numerical simulations. Similar to what was observed in the experiments for the case of transitional interaction, the computations reveal an intermittent upstream influence (UI) shock that repeatedly travels upstream from the lambda-foot toward the leading edge before vanishing. Through detailed unsteady flow analysis obtained using Fourier analysis and dynamic mode decomposition techniques, the presence of disturbances with similar frequencies as those measured in experiments were identified in the flow along with locations that appear to influence the dynamics of the flow.

Venkatachari, Balaji Shankar↗

Reduced Navier Stokes Relaxation Procedures for Internal Flows

In spite of significant advancement in the field of high speed computing, flow calculations involving complex geometries and/or flow behavior still require large amounts of CPU time and memory. In order to predict such flows without sacrificing grid convergence and accuracy, adaptive gridding techniques that provide optimal resolution are highly desirable. The present work combines multigrid techniques and domain decomposition concepts to provide local, solution adaptive, grid refinement. Several viscous compressible and incompressible, two and three-dimensional, flows with strong inviscid interaction and/or axial flow reversal, are considered with a segmented multigrid domain decomposition (SMGDD) procedure for which uniform meshes result in each domain. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable non-staggered grid computation. The SMGDD technique presented herein has previously been applied for incompressible two dimensional flows. The present work offers improvement in the gridding strategy, by allowing for disjoint subdomains that provide optimal resolution of disparate flow features. It also extends the SMGDD technique to three dimensional compressible flows. Laminar and turbulent flow in a backward facing step channel is considered; although the procedure is applicable to more severe geometries. The standard K-epsilon model is applied for turbulence closure. For Re greater than 400, differences between two-dimensional theory and experiment are resolved through a three dimensional simulation, which confirms the experimentally observed three dimensionality of the recirculation patterns on the upper and lower surfaces.

Rubin, Stanley G.↗

A global marching technique for the prediction of separated flows over arbitrary airfoils

The flow over arbitrary airfoils is analyzed using a boundary-layer-type solver developed for studying separated flows. The scheme includes extension to a general coordinate system and use of a more general zonal technique for solving the coupled equations. In order to be able to consider arbitrary geometries, second order accurate (in space) conservative differences are generated by considering the integral formulation of the governing equations in a general coordinate system. The general coordinate system is handled in as general a manner as possible to allow for the use of either analytically or numerically generated coordinate system. Different C type grids (orthogonal and nonorthogonal) were generated using different grid generation techniques to analyze the flow over arbitrary airfoils (J012, J025, NACA0012). A comparison with the experiment is demonstrated.

Halim, Ahmad A. M.↗

The efficient solution of transonic wing flow fields

An evaluation of the transonic-wing-analysis computer code TWING is presented. TWING utilizes a fully implicit, approximate-factorization iteration scheme to solve the full-potential equation in conservative form. A numerical elliptic-solver grid-generation scheme is used to generate the required finite-difference mesh. Several wing configurations have been analyzed, and comparisons of computed results have been made with available experimental data. Results indicate that the code is robust, accurate (when significant viscous effects are not present), and efficient. TWING generally produces solutions an order of magnitude faster than other conservative, full-potential codes using successive-line overrelaxation. The present method is applicable to a wide range of isolated wing configurations, including high-aspect-ratio transport wings and low-aspect-ratio, high-sweep, fighter configurations.

Holst, T. L.↗

A multidomain global modeling technique for analysis of Space Shuttle main engine

A multidomain global modeling technique is presented. The purpose of this multidomain methodology is to simultaneously analyze the three-dimensional compressible flow field in Space Shuttle main engine components. This technique provides flexibility and modularity for individual domain and also simplifies the grid generation task for each domain. A two-way coupling procedure is adopted for exchanging interface boundary conditions. This algorithm which requires that the adjacent domains share an overlapping grid slab, is fully conservative and implicit. For non-matching grid systems and for non-overlapping interfaces, an interpolation and/or extrapolation procedure is used when transferring the interface boundary conditions at the adjacent domains. To maintain overall conservation, a flux correction method is subsequently used. The coupling procedure was verified by modeling several geometries with features similar to SSME components. Satisfactory results for all test cases were obtained. Finally, this algorithm has been applied to the SSME flow simulation.

Tam, L. T.↗

An incompressible Navier-Stokes flow solver in three-dimensional curvilinear coordinate systems using primitive variables

An implicit, finite-difference computer code has been developed to solve the incompressible Navier-Stokes equations in a three-dimensional, curvilinear coordinate system. The pressure-field solution is based on the pseudo compressibility approach in which the time derivative pressure term is introduced into the mass conservation equation to form a set of hyperbolic equations. The solution procedure employs an implicit, approximate factorization scheme. The Reynolds stresses, that are uncoupled from the implicit scheme, are lagged by one time-step to facilitate implementing various levels of the turbulence model. Test problems for external and internal flows are computed, and the results are compared with existing experimental data. The application of this technique for general three-dimensional problems is then demonstrated.

Kwak, D.↗