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,747 records · Page 97

High-frequency monitoring of water fluxes and nutrient loads to assess the effects of controlled drainage on water storage and nutrient transport

High nitrogen (N) and phosphorus (P) fluxes from upstream agriculture threaten aquatic ecosystems in surface waters and estuaries, especially in areas characterized by high agricultural N and P inputs and densely drained catchments like the Netherlands. Controlled drainage has been recognized as an effective option to optimize soil moisture conditions for agriculture and to reduce unnecessary losses of fresh water and nutrients. This is achieved by introducing control structures with adjustable overflow levels into subsurface tube drain systems. A small-scale (1 ha) field experiment was designed to investigate the hydrological and chemical changes after introducing controlled drainage. Precipitation rates and the response of water tables and drain fluxes were measured in the periods before the introduction of controlled drainage (2007–2008) and after (2009–2011). For the N and P concentration measurements, auto-analyzers for continuous records were combined with passive samplers for time-averaged concentrations at individual drain outlets. The experimental setup enabled the quantification of changes in the water and solute balance after introducing controlled drainage. The results showed that introducing controlled drainage reduced the drain discharge and increased the groundwater storage in the field. To achieve this, the overflow levels have to be elevated in early spring, before the drain discharge stops due to dryer conditions and falling groundwater levels. The groundwater storage in the field would have been larger if the water levels in the adjacent ditch were controlled as well by an adjustable weir. The N concentrations and loads increased, which was largely related to elevated concentrations in one of the three monitored tube drains. The P loads via the tube drains reduced due to the reduction in discharge after introducing controlled drainage. Furthermore, this may be counteracted by the higher groundwater levels and the larger contribution of N- and P-rich shallow groundwater and overland flow to the surface water.

54 ENVIRONMENTAL SCIENCES↗

Quasar Feedback in the Ultraluminous Infrared Galaxy F11119+3257: Connecting the Accretion Disk Wind with the Large-scale Molecular Outflow

In Tombesi et al., we reported the first direct evidence for a quasar accretion disk wind driving a massive (>100 M {sub ⊙} yr{sup −1}) molecular outflow. The target was F11119+3257, an ultraluminous infrared galaxy (ULIRG) with unambiguous type 1 quasar optical broad emission lines. The energetics of the accretion disk wind and molecular outflow were found to be consistent with the predictions of quasar feedback models where the molecular outflow is driven by a hot energy-conserving bubble inflated by the inner quasar accretion disk wind. However, this conclusion was uncertain because the mass outflow rate, momentum flux, and mechanical power of the outflowing molecular gas were estimated from the optically thick OH 119 μm transition profile observed with Herschel. Here, we independently confirm the presence of the molecular outflow in F11119+3257, based on the detection of ∼±1000 km s{sup −1} blue- and redshifted wings in the CO(1−0) emission line profile derived from deep ALMA observations obtained in the compact array configuration (∼2.″8 resolution). The broad CO(1−0) line emission appears to be spatially extended on a scale of at least ∼7 kpc from the center. Mass outflow rate, momentum flux, and mechanical power of (80–200) R{sub 7}{sup −1} M {sub ⊙} yr{sup −1}, (1.5–3.0) R{sub 7}{sup −1} L {sub AGN}/c, and (0.15–0.40)% R{sub 7}{sup −1} L{sub AGN}, respectively, are inferred from these data, assuming a CO-to-H{sub 2} conversion factor appropriate for a ULIRG (R {sub 7} is the radius of the outflow normalized to 7 kpc, and L {sub AGN} is the AGN luminosity). These rates are time-averaged over a flow timescale of 7 × 10{sup 6} yr. They are similar to the OH-based rates time-averaged over a flow timescale of 4 × 10{sup 5} yr, but about a factor of 4 smaller than the local (“instantaneous”; ≲10{sup 5} yr) OH-based estimates cited in Tombesi et al. The implications of these new results are discussed in the context of time-variable quasar-mode feedback and galaxy evolution. The need for an energy-conserving bubble to explain the molecular outflow is also reexamined.

79 ASTRONOMY AND ASTROPHYSICS↗

Integrability technique for fluid flow induced deformation of a boundary hair

The deformation of a dense carpet of hair due to Stokes flow in a channel can be described by a nonlinear integrodifferential equation for the shape of a single hair, which possesses several solutions for a given choice of parameters. Although it was posed in a previous study and it bears a resemblance to the pendulum problem from mechanics, this equation has not been analytically solved until now. Despite the presence on an integral with a nonlinear functional dependence on the dependent variable, the system is integrable. We compare the analytically obtained solution to a finite-difference numerical approach, identify the physically realizable solution branch, and briefly study the solution structure through a conserved energylike quantity. Time-dependent fluid-structure interactions are a rich and complex subject to investigate, and we argue that the solution discussed herein can be used as a basis for understanding these systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗

Conservation Voltage Reduction (CVR) via Two-Timescale Control in Unbalanced Power Distribution Systems

Voltage control devices are employed in power distribution systems to reduce the power consumption by operating the system closer to the lower acceptable voltage limits; this technique is called conservation voltage reduction (CVR). The different modes of operation for system’s legacy devices (with binary control) and new devices (e.g. smart inverters with continuous control) coupled with variable photovoltaic (PV) generation results in voltage fluctuations which makes it challenging to achieve CVR objective. This paper presents a two-timescale control of feeder’s voltage control devices to achieve CVR that includes (1)a centralized controller operating in a slower time-scale to coordinate voltage control devices across the feeder and (2) local controllers operating in a faster timescale to mitigate voltage fluctuations due to PV variability. The centralized controller utilizes a three-phase optimal power flow model to obtain the decision variables for both legacy devices and smart inverters. The local controllers operate smart inverters to minimize voltage fluctuations and restore nodal voltages to their reference values by adjusting the reactive power support. The proposed approach is validated using the IEEE-123 bus (medium-size) and R3-12.47-2 (large-size) feeders. It is demonstrated that the proposed approach is effective in achieving the CVR objective for unbalanced distribution systems.

optimal control, distributed generation, fluctuati↗

Efficient Monte Carlo With Graph-Based Subsurface Flow and Transport Models

Simulating flow and transport in fractured porous media frequently involves solving numerical discretizations of partial differential equations with a large number of degrees of freedom using discrete fracture network (DFN) models. Uncertainty in the properties of the fracture network that controls flow and transport requires a large number of DFN simulations to statistically describe quantities of interest. However, the computational cost of solving more than a few realizations of a large DFN can be intractable. As a means of circumventing this problem, we utilize both a high-fidelity DFN model and a graphbased model of flow and transport in combination with a multifidelity Monte Carlo (MC) method to reduce the number of high-fidelity simulations that are needed to obtain an accurate estimate of the quantity of interest. We demonstrate the approach by estimating quantiles of the breakthrough time for a conservative tracer in an ensemble of fractured porous media. Our results demonstrate that a multifidelity MC estimate, whose computational cost is equal to the cost of 10 DFN simulations, can be as accurate as a standard MC estimate that utilizes 1,000 DFN simulations. Thus the combination of our graphbased model with multifidelity MC estimates effectively reduces the computational cost of the problem by a factor of approximately 100.

58 GEOSCIENCES↗

Numerical methods for the calculation of three-dimensional nozzle exhaust flow fields

Numerical codes developed for the calculation of three-dimensional nozzle exhaust flow fields associated with hypersonic airbreathing aircraft are described. Both codes employ reference plane grid networks with respect to three coordinate systems. Program CHAR3D is a characteristic code utilizing a new wave preserving network within the reference planes, while program BIGMAC is a finite difference code utilizing conservation variables and a one-sided difference algorithm. Secondary waves are numerically captured by both codes, while the underexpansion shock and plume boundary are treated discretely. The exhaust gas properties consist of hydrogen-air combustion product mixtures in local chemical equilibrium. Nozzle contours are treated by a newly developed geometry package based on dual cubic splines. Results are presented for simple configurations demonstrating two- and three-dimensional multiple wave interactions.

Dash, S. M.↗

Hydrography synthesis using LANDSAT remote sensing and the SCS models

The land cover requirements of the Soil Conservation Service (SCS) Model used for hydrograph synthesis in urban areas were modified to be LANDSAT compatible. The Curve Numbers obtained with these alternate land cover categories compare well with those obtained in published example problems using the conventional categories. Emergency spillway hydrographs and synthetic flood frequency flows computed for a 21.1 sq. mi. test area showed excellent agreement between the conventional aerial photo-based and the Landsat-based SCS approaches.

Ragan, R. M.↗

Investigation of Perforated Convergent-divergent Diffusers with Initial Boundary Layer

An experimental investigation was made at Mach number 1.90 of the performance of a series of perforated convergent-divergent supersonic diffusers operating with initial boundary layer, which was induced and controlled by lengths of cylindrical inlets affixed to the diffusers. Supercritical mass-flow and peak total-pressure recoveries were decreased slightly by use of the longest inlets (4 inlet diameters in length). Combinations of cylindrical inlets, perforated diffusers, and subsonic diffuser were evaluated as simulated wind tunnels having second throats. Comparisons with noncontracted configurations of similar scale indicated conservatively computed power reductions of 25 percent.

DIFFUSERS, SUPERSONIC↗

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗

Compact Water Vapor Exchanger for Regenerative Life Support Systems

Thermal and environmental control systems for future exploration spacecraft must meet challenging requirements for efficient operation and conservation of resources. Regenerative CO2 removal systems are attractive for these missions because they do not use consumable CO2 absorbers. However, these systems also absorb and vent water to space along with carbon dioxide. This paper describes an innovative device designed to minimize water lost from regenerative CO2 control systems. Design studies and proof-of-concept testing have shown the feasibility of a compact, efficient membrane water vapor exchanger (WVX) that will conserve water while meeting challenging requirements for operation on future spacecraft. Compared to conventional WVX designs, the innovative membrane WVX described here has the potential for high water recovery efficiency, compact size, and very low pressure losses. The key innovation is a method for maintaining highly uniform flow channels in a WVX core built from water-permeable membranes. The proof-of-concept WVX incorporates all the key design features of a prototypical unit, except that it is relatively small scale (1/23 relative to a unit sized for a crew of six) and some components were fabricated using non-prototypical methods. The proof-of-concept WVX achieved over 90% water recovery efficiency in a compact core in good agreement with analysis models. Furthermore the overall pressure drop is very small (less than 0.5 in. H2O, total for both flow streams) and meets requirements for service in environmental control and life support systems on future spacecraft. These results show that the WVX provides very uniform flow through flow channels for both the humid and dry streams. Measurements also show that CO2 diffusion through the water-permeable membranes will have negligible effect on the CO2 partial pressure in the spacecraft atmosphere.

Izenson, Michael G.↗

Recent Development in the CESE Method for the Solution of the Navier-Stokes Equations Using Unstructured Triangular or Tetrahedral Meshes With High Aspect Ratio

In the multidimensional CESE development, triangles and tetrahedra turn out to be the most natural building blocks for 2D and 3D spatial meshes. As such the CESE method is compatible with the simplest unstructured meshes and thus can be easily applied to solve problems with complex geometries. However, because the method uses space-time staggered stencils, solution decoupling may become a real nuisance in applications involving unstructured meshes. In this paper we will describe a simple and general remedy which, according to numerical experiments, has removed any possibility of solution decoupling. Moreover, in a real-world viscous flow simulation near a solid wall, one often encounters a case where a boundary with high curvature or sharp corner is surrounded by triangular/tetrahedral meshes of extremely high aspect ratio (up to 106). For such an extreme case, the spatial projection of a space-time compounded conservation element constructed using the original CESE design may become highly concave and thus its centroid (referred to as a spatial solution point) may lie far outside of the spatial projection. It could even be embedded beyond a solid wall boundary and causes serious numerical difficulties. In this paper we will also present a new procedure for constructing conservation elements and solution elements which effectively overcomes the difficulties associated with the original design. Another difficulty issue which was addressed more recently is the wellknown fact that accuracy of gradient computations involving triangular/tetrahedral grids deteriorates rapidly as the aspect ratio of grid cells increases. The root cause of this difficulty was clearly identified and several remedies to overcome it were found through a rigorous mathematical analysis. However, because of the length of the current paper and the complexity of mathematics involved, this new work will be presented in another paper.

Chang, Sin-Chung↗

Noether currents for Eulerian variational principles in non-barotropic magnetohydrodynamics and topological conservations laws

We derive a Noether current for the Eulerian variational principle of ideal non-barotropic magnetohydrodynamics (MHD). It was shown previously that ideal non-barotropic MHD is mathematically equivalent to a five function field theory with an induced geometrical structure in the case that field lines cover surfaces and this theory can be described using a variational principle. Here we use various symmetries of the flow to derive topological constants of motion through the derived Noether current and discuss their implication for non-barotropic MHD.

42 ENGINEERING↗

Discharge regime of non-ambipolarity with a self-induced steady-state magnetic field in plasma sources with localized radio-frequency power deposition

Involving the idea for the Biermann effect known from space physics as well as recent discussions on non-ambipolarity of the electron and ion fluxes in low-pressure discharges, the study builds the discharge pattern in a source with localized RF power deposition outside the region of high electron density. A vortex dc current flowing in an RF discharge and a steady-state magnetic field induced by this current govern the discharge behavior. Owing to a shift in the positions of the electron-density and plasma-potential maxima, the dc current is driven with the purpose of keeping the conservativity of the dc field in the discharge. The results present the spatial structure of a discharge in a regime of non-ambipolarity of the electron and ion fluxes, including its modifications by the magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CHROTRAN: a mathematical and computational model for in situ heavy metal remediation in heterogeneous aquifers

Groundwater contamination by heavy metals is a critical environmental problem for which in situ remediation is frequently the only viable treatment option. For such interventions, a three-dimensional reactive transport model of relevant biogeochemical processes is invaluable. To this end, we developed a model, CHROTRAN, for in situ treatment, which includes full dynamics for five species: a heavy metal to be remediated, an electron donor, biomass, a nontoxic conservative bio-inhibitor, and a biocide. Direct abiotic reduction by donor-metal interaction as well as donor-driven biomass growth and bio-reduction are modeled, along with crucial processes such as donor sorption, bio-fouling and biomass death. Our software implementation handles heterogeneous flow fields, arbitrarily many chemical species and amendment injection points, and features full coupling between flow and reactive transport. We describe installation and usage and present two example simulations demonstrating its unique capabilities. One simulation suggests an unorthodox approach to remediation of Cr(VI) contamination.

54 ENVIRONMENTAL SCIENCES↗

Supersonic turbulent boundary-layer flows with tangential slot injection.

An implicit finite-difference method has been developed for the solution of the compressible boundary-layer equations. This method is applied to tangential slot injection into supersonic turbulent boundary-layer flows. In addition, the effects induced by the interaction between the boundary-layer displacement thickness and the external pressure field are considered. Three different eddy viscosity models have been used to specify the turbulent momentum exchange. One model depends on the species concentration profile, and the species conservation equation has been included in the system of governing partial differential equations. For air injected into air, the freestream and injected gases are treated as separate species which have common fluid properties. Calculations were made and results were compared with experimental data at freestream Mach numbers of 2.4 and 6.0 and with results of another finite-difference method. Good agreement was obtained for the reduction of wall skin friction with slot injection.

Miner, E. W.↗

Carbon wastewater treatment process

A new powdered-carbon treatment process is being developed for the elimination of the present problems, associated with the disposal of biologically active sewage waste solids, and with water reuse. This counter-current flow process produces an activated carbon, which is obtained from the pyrolysis of the sewage solids, and utilizes this material to remove the adulterating materials from the water. Additional advantages of the process are the elimination of odors, the removal of heavy metals, and the potential for energy conservation.

Humphrey, M. F.↗

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Baker, A. J.↗