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,369 records · Page 76

Viscous shock-layer flowfield analysis by an explicit-implicit method

The present work extends the recently reported implicit analogue of MacCormack's earlier widely-used explicit method to external axisymmetric laminar flows with strong entropy gradients. The details of the 'numerics' of the implicit part are provided in a body-oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analogue may appeal to researchers for using it in the analysis of three-dimensional flow problems.

Gupta, R. N.↗

Development of a Regional Groundwater Flow Model for the U.S. Department of Energy Oak Ridge Reservation: Numerical Model using MODFLOW-USG - 17368

A regional scale groundwater flow model (GFM) has been developed for the U.S. Department of Energy (DOE) Oak Ridge Reservation (ORR). Development of the model implements a key recommendation of the ORR Groundwater Strategy. The model will be used as the single calibrated flow model for the ORR and as the framework for future, smaller scale, modeling efforts to support cleanup actions and decisions under the Comprehensive Environmental Response, Compensation, and Liability Act of 1980 (CERCLA). The Oak Ridge site is located in a geologically complex region and encompasses three large government facilities, including the Y-12 National Security Complex (Y-12), Oak Ridge National Laboratory (ORNL), and the East Tennessee Technology Park (ETTP). An interagency approach for addressing legacy groundwater contamination from past operations at these facilities has been developed, resulting in an ORR Groundwater Strategy that was agreed to by DOE, the U.S. Environmental Protection Agency (EPA), and Tennessee Department of Environment and Conservation (TDEC). The ORR Groundwater Strategy provides a comprehensive approach to addressing and prioritizing groundwater issues across the DOE reservation. One component of the strategy is to 'develop and maintain an ORR-wide regional groundwater flow model to ensure a single, regional, calibrated model to support groundwater characterization, decision-making, and remediation.' As part of the cleanup plans for the Oak Ridge site under CERCLA and the Federal Facility Agreement (FFA), there are six final watershed-scale groundwater Records of Decision (RODs) currently planned. In preparation for future projects, the regional GFM has been constructed to serve as the single, calibrated regional flow model to be used as the hydrologic base for the groundwater plume-specific modeling developed for the RODs. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Space-Time Conservation Element and Solution Element Method-A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: Numerical Simulation of Shock Waves and Contact Discontinuities - 2

Without resorting to special treatment for each individual test case, the 1D and 2D CE/SE shock-capturing schemes described previously (in Part I) are used to simulate flows involving phenomena such as shock waves, contact discontinuities, expansion waves and their interactions. Five 1D and six 2D problems are considered to examine the capability and robustness of these schemes. Despite their simple logical structures and low computational cost (for the 2D CE/SE shock-capturing scheme, the CPU time is about 2 micro-secs per mesh point per marching step on a Cray C90 machine), the numerical results, when compared with experimental data, exact solutions or numerical solutions by other methods, indicate that these schemes can accurately resolve shock and contact discontinuities consistently.

Wang, Xiao-Yen↗

Factorizable Schemes for the Equations of Fluid Flow

We present an upwind high-resolution factorizable (UHF) discrete scheme for the compressible Euler equations that allows to distinguish between full-potential and advection factors at the discrete level. The scheme approximates equations in their general conservative form and is related to the family of genuinely multidimensional upwind schemes developed previously and demonstrated to have good shock-capturing capabilities. A unique property of this scheme is that in addition to the aforementioned features it is also factorizable, i.e., it allows to distinguish between full-potential and advection factors at the discrete level. The latter property facilitates the construction of optimally efficient multigrid solvers. This is done through a relaxation procedure that utilizes the factorizability property.

Sidilkover, David↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

Dynamically consistent hydrography and absolute velocity in the eastern North Atlantic Ocean

The problem of mapping a dynamically consistent hydrographic field and associated absolute geostrophic flow in the eastern North Atlantic between 24 deg and 36 deg N is related directly to the solution of the so-called thermocline equations. A nonlinear optimization problem involving Needler's P equation is solved to find the hydrography and resulting flow that minimizes the vertical mixing above about 1500 m in the ocean and is simultaneously consistent with the observations. A sharp minimum (at least in some dimensions) is found, apparently corresponding to a solution nearly conserving potential vorticity and with vertical eddy coefficient less than about 10(exp -5) sq m/s. Estimates of `residual' quantities such as eddy coefficients are extremely sensitive to slight modifications to the observed fields. Boundary conditions, vertical velocities, etc., are a product of the optimization and produce estimates differing quantitatively from prior ones relying directly upon observed hydrography. The results are generally insensitive to particular elements of the solution methodology, but many questions remain concerning the extent to which different synoptic sections can be asserted to represent the same ocean. The method can be regarded as a practical generalization of the beta spiral and geostrophic balance inverses for the estimate of absolute geostrophic flows. Numerous improvements to the methodology used in this preliminary attempt are possible.

Wunsch, Carl↗

An improved far field drag calculation method for nonlinear CFD codes

An improved method to calculate drag based on far field conditions is presented and demonstrated. The method is illustrated using two CFD codes, PARC and CFL3D, and examining their ability to preserve lift and drag coefficients along grid line defined integration paths. The flow fields were generated by solving the Euler equations for a NACA 0012 airfoil at a free stream Mach number of 0.8 and angle of attack of 1.25 degrees. In comparison to force coefficients obtained by surface pressure integration, neither code acceptably preserved both force coefficients throughout the near and far fields. An investigation into the relationship between numerical prediction error and calculated force coefficients revealed a direct connection between solution mass conservation error and force coefficients error. A method to correct the predicted force coefficients based on integrated mass conservation error is described and demonstrated. The corrected force coefficients for both codes are shown to be more accurate than the uncorrected values. The correction method is applicable to both two and three dimensions and is independent of the algorithm used to generate the flow field.

Sirbaugh, James R.↗

A viscous shock-layer flowfield analysis by an explicit-implicit method

An implicit analogue of a widely used explicit method to external axisymmetric laminar flows with strong entropy gradients is extended. The details of the "numerics" of the implicit part are provided in a body oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analog may appeal to researchers for using it in the analysis of 3-D flow problems.

Gupta, R. M.↗

Three-dimensional local ALE-FEM method for fluid flow in domains containing moving boundaries/objects interfaces

A three-dimensional finite element method for the numerical simulations of fluid flow in domains containing moving rigid objects or boundaries is developed. The method falls into the general category of Arbitrary Lagrangian Eulerian methods; it is based on a fixed mesh that is locally adapted in the immediate vicinity of the moving interfaces and reverts to its original shape once the moving interfaces go past the elements. The moving interfaces are defined by separate sets of marker points so that the global mesh is independent of interface movement and the possibility of mesh entanglement is eliminated. The results is a fully robust formulation capable of calculating on domains of complex geometry with moving boundaries or devises that can also have a complex geometry without danger of the mesh becoming unsuitable due to its continuous deformation thus eliminating the need for repeated re-meshing and interpolation. Moreover, the boundary conditions on the interfaces are imposed exactly. This work is intended to support the internal combustion engines simulator KIVA developed at Los Alamos National Laboratories. The model's capabilities are illustrated through application to incompressible flows in different geometrical settings that show the robustness and flexibility of the technique to perform simulations involving moving boundaries in a three-dimensional domain.

16 TIDAL AND WAVE POWER↗

A robust and efficient finite volume method for compressible inviscid and viscous two-phase flows [A robust and efficient finite volume method for compressible viscous two-phase flows]

A robust and efficient density-based finite volume method is developed for solving the six-equation single pressure system of two-phase flows at all speeds on hybrid unstructured grids. Unlike conventional approaches where an expensive exact Riemann solver is normally required for computing numerical fluxes at the two-phase interfaces in addition to AUSM-type fluxes for single-phase interfaces in order to maintain stability and robustness in cases involving interactions of strong pressure and void-fraction discontinuities, a volume-fraction coupling term for the AUSM + -up fluxes is introduced in this work to impart the required robustness without the need of the exact Riemann solver. The resulting method is significantly less expensive in regions where otherwise the Riemann solver would be invoked. A transformation from conservative variables to primitive variables is presented and the primitive variables are then solved in the implicit method in order for the current finite volume method to be able to solve, effectively and efficiently, low Mach number flows in traditional multiphase applications, which otherwise is a great challenge for the standard density-based algorithms. Here, a number of benchmark test cases are presented to assess the performance and robustness of the developed finite volume method for both inviscid and viscous two-phase flow problems. The numerical results indicate that the current density-based method provides an attractive and viable alternative to its pressure-based counterpart for compressible two-phase flows at all speeds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Unified Nusselt- and Sherwood-number correlations in axisymmetric finite-gap stagnation and rotating-disk flows

This paper develops a unified analysis of stagnation flow heat and mass transport, considering both semi-infinite domains and finite gaps, with and without rotation of the stagnation surface. An important objective is to derive Nusselt- and Sherwood-number correlations that represent heat and mass transport at the stagnation surface. The approach is based on computationally solving the governing conservation equations in similarity form as a boundary-value problem. The formulation considers ideal gases and incompressible fluids. The correlated results depend on fluid properties in terms of Prandtl, Schmidt, and Damkohler numbers. Heterogeneous chemistry at the stagnation surface is represented as a single first-order reaction. A composite Reynolds number represents the combination of stagnation flows with and without stagnation-surface rotation.

97 MATHEMATICS AND COMPUTING↗

A novel Variable Refrigerant Flow (VRF) heat recovery system model: Development and validation

As one of the latest emerging HVAC technologies, the Variable Refrigerant Flow (VRF) system with heat recovery (HR) configurations has obtained extensive attention from both the academia and industry. Compared with the conventional VRF systems with heat pump (HP) configurations, VRF-HR is capable of recovering heat from cooling zones to heating zones and providing simultaneous cooling and heating operations. Moreover, this could lead to substantial energy saving potential and more flexible zonal control. In this paper, a novel model is developed to simulate the energy performance of VRF-HR systems. It adheres to a more physics-based development with the ability to simulate the refrigerant loop performance and consider the dynamics of more operational parameters, which is crucial for representing more advanced control logics. Another key feature of the model is the introduction of component-level curves for indoor units and outdoor units instead of overall performance curves for the entire system, and thus it requires much fewer user-specified performance curves as model inputs. The validation study reflects good agreements between the simulated energy use from the new VRF-HR model and the laboratory measurement data across all operational modes at sub-hourly time steps. The model has been adopted in the official release of the EnergyPlus simulation program since Version 8.6, which enables more accurate and robust assessments of VRF-HR systems to support their applications in energy retrofit of existing buildings or design of zero-net-energy buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Smoothed Particle Hydrodynamics Continuous Boundary Force method for Navier-Stokes equations subject to Robin boundary condition

Robin boundary condition for the Navier-Stokes equations is used to model slip conditions at the fluid-solid boundaries. A novel Continuous Boundary Force (CBF) method is proposed for solving the Navier-Stokes equations subject to Robin boundary condition. In the CBF method, the Robin boundary condition at boundary is replaced by the homogeneous Neumann boundary condition at the boundary and a volumetric force term added to the momentum conservation equation. Smoothed Particle Hydrodynamics (SPH) method is used to solve the resulting Navier-Stokes equations. We present solutions for two-dimensional and three-dimensional flows in domains bounded by flat and curved boundaries subject to various forms of the Robin boundary condition. The numerical accuracy and convergence are examined through comparison of the SPH-CBF results with the solutions of finite difference or finite element method. Taken the no-slip boundary condition as a special case of slip boundary condition, we demonstrate that the SPH-CBF method describes accurately both no-slip and slip conditions.

Pan, Wenxiao↗

Lagrangian and Dirac constraints for the ideal incompressible fluid and magnetohydrodynamics

The incompressibility constraint for fluid flow was imposed by Lagrange in the so-called Lagrangian variable description using his method of multipliers in the Lagrangian (variational) formulation. An alternative is the imposition of incompressibility in the Eulerian variable description by a generalization of Dirac’s constraint method using noncanonical Poisson brackets. Here it is shown how to impose the incompressibility constraint using Dirac’s method in terms of both the canonical Poisson brackets in the Lagrangian variable description and the noncanonical Poisson brackets in the Eulerian description, allowing for the advection of density. Both cases give the dynamics of infinite-dimensional geodesic flow on the group of volume preserving diffeomorphisms and explicit expressions for this dynamics in terms of the constraints and original variables is given. Since Lagrangian and Eulerian conservation laws are not identical, comparison of the various methods is made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energetic Potential for Demand Response in Detached Single Family Homes in Austin, TX

This study builds a generalized tool to quantify the maximum peak load reduction achievable with residential demand response and distributed energy resources (DERs). The tool is demonstrated by using empirical energy usage data from detached single family households within the service territory of Austin Energy - the local municipal utility in Austin, TX. The demand response optimization algorithm shifts energy usage of four controllable, high-consumption residential devices to off-peak hours: HVAC (heating, ventilation, and air conditioning) systems, electric water heaters, electric vehicles, and pool pumps. In addition, rooftop solar generation on individual residences is considered to reduce net load to the utility, and an energy storage system optimization model performing energy arbitrage further reduces peak demand. Results from the aggregate demand response and battery optimization study indicate that peak demand for the 2017 summer peak in Austin Energy could have been reduced by 26 MW or 2.9%. Here, by considering the combined effect of controllable consumer loads, solar generation, and energy storage systems, this study highlights the potential of strategic demand response and energy arbitrage to reduce peak demand at the distribution level.

14 SOLAR ENERGY↗

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

Transfer and Flushing Evaluation for Potential Precipitated Solids in the 241-AP-106 to EMF Transfer Pipeline for DFLAW

Waste transfer pipeline flushing is necessary to maintain pipeline performance, but added flush fluid increases the nuclear waste inventory. During the Direct Feed Low-Activity Waste (DFLAW) mission, waste transfer lines from the Interim Low-Activity Waste Storage Tank (i.e., 241-AP-106) to the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Facility must have the capability of being flushed to prevent accumulation of solids and to mitigate corrosion concerns. Operational experience at Hanford has shown that the solids precipitation occurs during processing operations from dilute Hanford liquids. Pacific Northwest National Laboratory was requested to evaluate the potential precipitated solid particle transport and flushing operational capabilities of the AP-106 through the Effluent Management Facility (EMF) low-point drain LAW feed pipeline. The evaluation results do not provide operational requirements, but rather provide a scoping basis for understanding the potential operational significance of solids precipitation in the pipeline. The salt and aluminum phase solids most likely to precipitate from the LAW feeds during cooling, evaporation, or mixing have particle densities that are estimated to range from 1.62 to 2.78 g/mL and a spherical particle size range of 8 to 2100 μm, depending on the solid phase considered. Application of the method to calculate the critical deposition velocity required by TFC-ENG-STD-26 at conservatively bounding estimates for the solids concentrations demonstrates that potentially precipitated solids may deposit on the bottom of the transfer pipe invert at the lower LAW feed pipeline flow rate, but the upper LAW feed pipeline flow rate will likely prevent deposition. The system pressure limit for the AP-106 through the EMF low-point drain LAW feed pipeline likely exceeds any pressure loss at a calculated critical deposition velocity even for the most adverse potential precipitated solid. However, should solids settle as a result of no-flow conditions, the flow capabilities of the AP-106 through the EMF low-point drain LAW feed pipeline are shown, based on the available literature, to be potentially inadequate for effective flushing operations to remove solids, even at solids concentrations below the maximum specified for LAW feed. Thus, stepwise solids accumulations over multiple transfers may be an issue, and eventual line plugging may occur. Further inadequacies may potentially be realized should the cohesive nature of certain precipitates be accounted for.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗