Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conservative schemes”

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 865 records · Page 48

Validation of a Co-Optimized, Smart Hybrid Heat Pump Control

This study examines a residential dual fuel (DF) heating system that integrates heat pump and gas furnace technology with a home’s forced-air distribution ducts. Specifically, it explores a retrofit where an “add-on” heat pump replaces the air conditioning (A/C) system and is combined with the home’s existing furnace and air handler. The research team developed and validated an innovative control scheme for this integrated DF system that selected the more economical heating source (furnace or heat pump) while maintaining comfort in the home. The project’s field study and subsequent modeling analysis show that this control scheme can decrease energy use, cost, and emissions compared to a furnace-only system. As homes in mixed and cold climates transition to more electric space heating, furnace-to-dual fuel system retrofits offer a solution that can alleviate the retrofit costs, energy cost concerns, and peak load impacts of full electrification.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A direct current rectification scheme for microwave space power conversion using traveling wave electron acceleration

The formation of the Vision-21 conference held three years ago allowed the present author to reflect and speculate on the problem of converting electromagnetic energy to a direct current by essentially reversing the process used in traveling wave tubes that converts energy in the form of a direct current to electromagnetic energy. The idea was to use the electric field of the electromagnetic wave to produce electrons through the field emission process and accelerate these electrons by the same field to produce an electric current across a large potential difference. The acceleration process was that of cyclotron auto-resonance. Since that time, this rather speculative ideas has been developed into a method that shows great promise and for which a patent is pending and a prototype design will be demonstrated in a potential laser power beaming application. From the point of view of the author, a forum such as Vision-21 is becoming an essential component in the rather conservative climate in which our initiatives for space exploration are presently formed. Exchanges such as Vision-21 not only allows us to deviate from the 'by-the-book' approach and rediscover the ability and power in imagination, but provides for the discussion of ideas hitherto considered 'crazy' so that they may be given the change to transcend from the level of eccentricity to applicability.

Manning, Robert M.↗

An Efficient and Accurate Algorithm for Computing Grid-Averaged Solar Fluxes for Horizontally Inhomogeneous Clouds

A computationally efficient method is presented to account for the horizontal cloud inhomogeneity by using a radiatively equivalent plane parallel homogeneous (PPH) cloud. The algorithm can accurately match the calculations of the reference (rPPH) independent column approximation (ICA) results, but use only the same computational time required for a single plane parallel computation. The effective optical depth of this synthetic sPPH cloud is derived by exactly matching the direct transmission to that of the inhomogeneous ICA cloud. The ffective9 scattering asymmetry factor is found from a pre-calculated albedo inverse look-up-table that is allowed to vary over the range from -1.0 to 1.0. In the special cases of conservative scattering and total absorption, the synthetic method is exactly equivalent to the ICA, with only a small bias (about 0.2% in flux) relative to ICA due to imperfect interpolation in using the look-up tables. In principle, the ICA albedo can be approximated accurately regardless of cloud inhomogeneity. For a more complete comparison, the broadband shortwave albedo and transmission calculated from the synthetic sPPH cloud and averaged over all incident directions, have the RMS biases of 0.26% and 0.76%, respectively, for inhomogeneous clouds over a wide variation of particle size. The advantages of the synthetic PPH method are that (1) it is not required that all the cloud subcolumns have uniform microphysical characteristic, (2) it is applicable to any 1D radiative transfer scheme, and (3) it can handle arbitrary cloud optical depth distributions and an arbitrary number of cloud subcolumns with uniform computational efficiency.

cloud inhomogeneity↗

A Runge-Kutta discontinuous finite element method for high speed flows

A Runge-Kutta discontinuous finite element method is developed for hyperbolic systems of conservation laws in two space variables. The discontinuous Galerkin spatial approximation to the conservation laws results in a system of ordinary differential equations which are marched in time using Runge-Kutta methods. Numerical results for the two-dimensional Burger's equation show that the method is (p+1)-order accurate in time and space, where p is the degree of the polynomial approximation of the solution within an element and is capable of capturing shocks over a single element without oscillations. Results for this problem also show that the accuracy of the solution in smooth regions is unaffected by the local projection and that the accuracy in smooth regions increases as p increases. Numerical results for the Euler equations show that the method captures shocks without oscillations and with higher resolution than a first-order scheme.

Bey, Kim S.↗

Thermoelastic response of metal matrix composites with large-diameter fibers subjected to thermal gradients

A new micromechanical theory is presented for the response of heterogeneous metal matrix composites subjected to thermal gradients. In contrast to existing micromechanical theories that utilize classical homogenization schemes in the course of calculating microscopic and macroscopic field quantities, in the present approach the actual microstructural details are explicitly coupled with the macrostructure of the composite. Examples are offered that illustrate limitations of the classical homogenization approach in predicting the response of thin-walled metal matrix composites with large-diameter fibers when subjected to thermal gradients. These examples include composites with a finite number of fibers in the thickness direction that may be uniformly or nonuniformly spaced, thus admitting so-called functionally gradient composites. The results illustrate that the classical approach of decoupling micromechanical and macromechanical analyses in the presence of a finite number of large-diameter fibers, finite dimensions of the composite, and temperature gradient may produce excessively conservative estimates for macroscopic field quantities, while both underestimating and overestimating the local fluctuations of the microscopic quantities in different regions of the composite. Also demonstrated is the usefulness of the present approach in generating favorable stress distributions in the presence of thermal gradients by appropriately tailoring the internal microstructure details of the composite.

Aboudi, Jacob↗

An implicit algorithm for the transonic full-potential equation in conservative form

A fast, implicit approximate factorization algorithm for the solution of the conservative full-potential equation for transonic flow in two and three dimensions is presented. Stability in supersonic regions is maintained by the use of an upwind evaluation of the density coefficient along all coordinate directions, providing an effective upwind difference of the streamwise terms for any orientation of the velocity vector and thereby enhancing the reliability of the algorithm. The algorithm is shown to provide rapid convergence for the computation of certain difficult two-dimensional test cases, including cases with fishtail shock patterns, demonstrating the reliability and efficiency of the procedure. Surface pressure coefficient distributions obtained by the present method are also found to be in good agreement with those computed by successive-line overrelaxation and a hybrid direct-solver/successive-line overrelaxation scheme, with significant reductions in CPU time required. A three-dimensional solution for a swept wing mounted between parallel walls is also presented which demonstrates the high convergence rate of the algorithm in three dimensions as well as two.

Holst, T.↗

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗

Process prediction and detection of faults using probabilistic bidirectional recurrent neural networks on real plant data

Attaining Industry 4.0 for manufacturing operations requires advanced monitoring systems and real-time data analytics of plant data, among other topics. We propose a Probabilistic Bidirectional Recurrent Network (PBRN) for industrial process monitoring for the early detection of faults. The model is based on a Gated Recurrent Unit (GRU) neural network that allows the model to retain long-term dependencies between sensor data along a time horizon, hence learning the dynamic behavior of the process. To reduce the false-positive detection rate of the model, we compel the model to learn from a highly noisy sensor reading while outputting noise-free sensor outputs. The performance of the proposed model is compared to other data-driven statistical process monitoring schemes using real plant data from an industrial Air Separations Unit (ASU) containing noisy sensor readings. We show that the model can learn from noisy data without reducing its performance. Using two different fault cases, we demonstrate the model’s ability to carry out early fault detection with average false-positive rates of 2.9% and 4.9% for both fault cases. The missed detection rates are 0.1% and 0.2%, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An Analysis of Energy Needs and Seawater Extraction Approaches for the Integration of Primary Magnesium Production with Nuclear Power

The cost of primary magnesium production is highly dominated by the energy required to produce the metal. This, along with the current geographical concentration of its production, makes the magnesium price vulnerable as energy-production paradigms begin to rapidly change. In this work, four nuclear-integration pathways for the efficient supply of heat and electricity to magnesium production from a nuclear power plant are considered. A production route using electrified metallothermic reduction is compared to a nuclear-heat-integrated electrolytic pathway. It is shown that the latter can achieve significantly lower primary energy consumption, 261 MJ/kg Mg compared with 416 MJ/kg Mg. The electrolytic pathway can also reduce emissions down to around 2.1 kgCO2/kg Mg, compared to around 34 kgCO2/kg Mg in the Pidgeon process. Beyond energy and emissions, the four pathways considered are also compared to their ability to reduce locational dependence of the production plant and reduce process complexity. This work also highlights future advances in magnesium production with a potential scheme for the efficient extraction of magnesium from seawater with coproduction of desalinated water and the potential of demand-side response from magnesium electrolysis. This is the first work to consider the full integration of primary magnesium production using nuclear cogeneration, and it is hoped the study paves the way for the analysis of economic nuclear integration within metal extraction and manufacturing.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

A Global End-Member Approach to Derive aCDOM(440) from Near-Surface Optical Measurements

This study establishes an optical inversion scheme for deriving the absorption coefficient of colored (or chromophoric, depending on the literature) dissolved organic material (CDOM) at the 440 nm wavelength, which can be applied to global water masses with near-equal efficacy. The approach uses a ratio of diffuse attenuation coefficient spectral end members, i.e., a short and long wavelength pair. The global perspective is established by sampling "extremely" clear water plus a generalized extent in turbidity and optical properties that each span three decades of dynamic range. A unique data set was collected in oceanic, coastal, and inland waters (as shallow as 0.6 m) from the North Pacific Ocean, the Arctic Ocean, Hawaii, Japan, Puerto Rico, and the east and west coasts of the United States. The data were partitioned using subjective categorizations to define a validation quality subset of conservative water masses, i.e., the inflow and outflow of properties constrain the range in the gradient of a constituent, plus 15 subcategories of water masses that were not evolving conservatively. The dependence on subcategories was confirmed with an objective methodology based on cluster analysis techniques. The latter defined five distinct classes with validation quality data present in all classes, but which also decreased in percent composition as a function of increasing class number and optical complexity. Four different algorithms based on different validation quality end members were validated with accuracies of 1.–6.2 %, wherein pairs with the largest spectral span were most accurate. Although algorithm accuracy decreased with the inclusion of more subcategories containing non-conservative water masses, changes to the algorithm fit were small when a preponderance of subcategories were included. The high accuracy for all end-member algorithms was the result of data acquisition and data processing improvements, e.g., increased vertical sampling resolution to less than 1mm and a boundary constraint to mitigate wave focusing effects, respectively. An independent evaluation with a historical database confirmed the consistency of the algorithmic approach and its application to quality assurance, e.g., to flag data outside expected ranges, identify suspect spectra, and objectively determine the in-water extrapolation interval by converging agreement for all applicable end-member algorithms. The legacy data exhibit degraded performance (as 44 % uncertainty) due to a lack of high-quality near-surface observations, especially for clear waters wherein wave-focusing effects are problematic. The novel optical approach allows the in situ estimation of an in-water constituent in keeping with the accuracy obtained in the laboratory.

Stanford B Hooker↗

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

Simple high-accuracy resolution program for convective modelling of discontinuities

For steady multidimensional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

SSME Turbopump Turbine Computations

A two-dimensional viscous code was developed to be used in the prediction of the flow in the SSME high-pressure turbopump blade passages. The rotor viscous code (RVC) employs a four-step Runge-Kutta scheme to solve the two-dimensional, thin-layer Navier-Stokes equations. The Baldwin-Lomax eddy-viscosity model is used for these turbulent flow calculations. A viable method was developed to use the relative exit conditions from an upstream blade row as the inlet conditions to the next blade row. The blade loading diagrams are compared with the meridional values obtained from an in-house quasithree-dimensional inviscid code. Periodic boundary conditions are imposed on a body-fitted C-grid computed by using the GRAPE GRids about Airfoils using Poisson's Equation (GRAPE) code. Total pressure, total temperature, and flow angle are specified at the inlet. The upstream-running Riemann invariant is extrapolated from the interior. Static pressure is specified at the exit such that mass flow is conserved from blade row to blade row, and the conservative variables are extrapolated from the interior. For viscous flows the noslip condition is imposed at the wall. The normal momentum equation gives the pressure at the wall. The density at the wall is obtained from the wall total temperature.

Jorgenson, P. G. E.↗

Lagrangian formulation of nuclear–electronic orbital Ehrenfest dynamics with real-time TDDFT for extended periodic systems

Here, we present a Lagrangian-based implementation of Ehrenfest dynamics with nuclear–electronic orbital (NEO) theory and real-time time-dependent density functional theory for extended periodic systems. In addition to a quantum dynamical treatment of electrons and selected protons, this approach allows for the classical movement of all other nuclei to be taken into account in simulations of condensed matter systems. Furthermore, we introduce a Lagrangian formulation for the traveling proton basis approach and propose new schemes to enhance its application for extended periodic systems. Validation and proof-of-principle applications are performed on electronically excited proton transfer in the o-hydroxybenzaldehyde molecule with explicit solvating water molecules. These simulations demonstrate the importance of solvation dynamics and a quantum treatment of transferring protons. This work broadens the applicability of the NEO Ehrenfest dynamics approach for studying complex heterogeneous systems in the condensed phase.

Calculus of variations↗

Smoothed particle hydrodynamics with GRAPE-1A

We describe the implementation of a smoothed particle hydrodynamics (SPH) scheme using GRAPE-1A, a special-purpose processor used for gravitational N-body simulations. The GRAPE-1A calculates the gravitational force exerted on a particle from all other particles in a system, while simultaneously making a list of the nearest neighbors of the particle. It is found that GRAPE-1A accelerates SPH calculations by direct summation by about two orders of magnitudes for a ten thousand-particle simulation. The effective speed is 80 Mflops, which is about 30 percent of the peak speed of GRAPE-1A. Also, in order to investigate the accuracy of GRAPE-SPH, some test simulations were executed. We found that the force and position errors are smaller than those due to representing a fluid by a finite number of particles. The total energy and momentum were conserved within 0.2-0.4 percent and 2-5 x 10 exp -5, respectively, in simulations with several thousand particles. We conclude that GRAPE-SPH is quite effective and sufficiently accurate for self-gravitating hydrodynamics.

Umemura, Masayuki↗

Soils' dirty little secret: Depth–based comparisons can be inadequate for quantifying changes in soil organic carbon and other mineral soil properties

Quantifying changes in soil organic carbon (SOC) stocks and other soil properties is essential for understanding how soils will respond to land management practices and global change. Although they are widely used, comparisons of SOC stocks at fixed depth (FD) intervals are subject to errors when changes in bulk density or soil organic matter occur. The equivalent soil mass (ESM) method has been recommended in lieu of FD for assessing changes in SOC stocks in mineral soils, but ESM remains underutilized for SOC stocks and has rarely been used for other soil properties. In this paper, we draw attention to the limitations of the FD method and demonstrate the advantages of the ESM approach. We provide illustrations to show that the FD approach is susceptible to errors not only for quantifying SOC stocks but also for soil mass-based properties such as SOC mass percent, C:N mass ratio, and δ 13 C. Here we describe the ESM approach and show how it mitigates the FD method limitations. Using bulk density change simulations applied to an empirical dataset from bioenergy cropping systems, we show that the ESM method provides consistently lower errors than FD when quantifying changes in SOC stocks and other soil properties. To simplify the use of ESM, we detail how the method can be integrated into sampling schemes, and we provide an example R computer script that can perform ESM calculations on large datasets. We encourage future studies, whether temporal or comparative, to utilize sampling methods that are amenable to the ESM approach. Overall, we agree with previous recommendations that ESM should be the standard method for evaluating SOC stock changes in mineral soils, but we further suggest that ESM may also be preferred for comparisons of other soil properties including mass percentages, elemental mass ratios, and stable isotope composition.

59 BASIC BIOLOGICAL SCIENCES↗

Gaussian Process Optimization of Sensitivity-Based Similarity Metrics between New Nuclear Applications and New/Existing Benchmarks [Slides]

This presentation discusses Nuclear Criticality Safety (NCS) and how designing safe, new nuclear criticality experiments requires expert judgement, which could take years of experience. Sensitivity/uncertainty (S/U) analysis can be utilized by less experienced individuals to conservatively estimate uncertainties in important parameters, such as k eff , in newly proposed nuclear experiments. The presentation poses the question of how this analysis can be performed and states that the answer lies in matching new nuclear experiments with existing benchmark experiments using similarity metrics. By increasing the criticality safety of the application in this work, higher mass limits could be used in PF-4 operations. Additionally, the presentation discusses MCNP6.2®, Whisper-1.1, the software that can be used in this analysis. Also discussed is the fact that International Criticality Safety Benchmark Evaluation Project (ICSBEP) benchmarks rarely match new nuclear applications and that there are significant differences in given set of materials and/or geometry. If there are no benchmarks that match the application, the presentation discusses the possibility of creating new benchmarks. In conclusion, this work presents a Gaussian process (GP) optimization scheme that was used to generate new benchmarks with the highest sensitivity-based similarity metrics to user-defined nuclear applications. The Gaussian process optimization successfully designed 3 new experimental benchmarks that were highly correlated to the application of interest and had k eff values near critical. Optimization over c k,i-r has shown that investigating specific isotope reactions for different applications is crucial to designing benchmark experiments. Partial contribution from Pu dominates c k similarity metric. Future work includes testing new stand-alone similarity metrics or new combinations of similarity metrics as the design criterion of this optimization – design criterion is application dependent.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Distribution System Operator with Transactive (DSO+T) Study

The Distribution System Operator with Transactive (DSO+T) study investigates the engineering and economic performance of a transactive energy retail market coordinating a high penetration of customer-side flexible energy assets. The study seeks to answer whether such an implementation is cost effective for customers, recovers sufficient revenue for DSOs, and is equally applicable and beneficial to a range of flexible asset types, renewable generation scenarios, and market assumptions. Using a highly interdisciplinary co-simulation and valuation framework, this assessment encompasses the entire electrical delivery system from bulk system generation and transmission, through the distribution system, to the modeling of individual customer buildings and flexible assets (including heating, ventilation, and air conditioning [HVAC] units, water heaters, batteries, and electric vehicles). The study exercises a transactive energy retail market coordination scheme designed to integrate with an existing day-ahead and real-time competitive wholesale electricity market. Software decision-making agents are designed for the retail market operator as well as various price-responsive flexible assets. The engineering and economic performance of the transactive energy scheme is studied for two separate flexible asset deployments: flexible loads (HVAC units and residential water heaters) and behind-the-meter batteries. The results of each transactive case are compared to a business-as-usual case. These cases are subject to two different renewable generation scenarios, a moderate renewable generation scenario, representative of current levels of renewable generation deployment, and a future high renewables scenario, including the increased deployment of rooftop solar photovoltaic and electric vehicles. The transactive coordination scheme is shown to produce effective and stable control and decrease peak loads 9–15%. The resulting annual demand flexibility provides net economic savings of $3.3–5.0B per year for a region the size of Texas. Detailed analysis shows that net benefits were seen for a range of distribution system operator, customer, and flexible asset types. Both participating customer (with transactive flexible assets) and nonparticipating customers (with nonflexible assets) see reductions in annual utility bills and net annual energy expenses in the range of 10–16%.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗