Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equation recovery”

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 91 records · Page 5

An assessment of gravity model improvements using TOPEX/Poseidon TDRSS observations

The contribution of TOPEX/Poseidon (T/P) TDRSS data to geopotential model recovery is assessed. Simulated TDRSS one-way and Bilateration Ranging Transponder System (BRTS) observations have been generated and orbitally reduced to form normal equations for geopotential parameters. These normals have been combined with those of the latest prelaunch T/P gravity model solution using data from over 30 satellites. A study of the resulting solution error covariance shows that TDRSS can make important contributions to geopotential recovery, especially for improving T/P specific effects like those arising from orbital resonance. It is argued that future effort is desirable both to establish TDRSS orbit determination limits in a reference frame compatible with that used for the precise laser/DORIS orbits, and the reduction of these TDRSS data for geopotential recovery.

Putney, B. H.↗

An improved reconstruction method for cosmological density fields

This paper proposes some improvements to existing reconstruction methods for recovering the initial linear density and velocity fields of the universe from the present large-scale density distribution. We derive the Eulerian continuity equation in the Zel'dovich approximation and show that, by applying this equation, we can trace the evolution of the gravitational potential of the universe more exactly than is possible with previous approaches based on the Zel'dovich-Bernoulli equation. The improved reconstruction method is tested using N-body simulations. When the Zel'dovich-Bernoulli equation describes the formation of filaments, then the Zel'dovich continuity equation also follows the clustering of clumps inside the filaments. Our reconstruction method recovers the true initial gravitational potential with an rms error about 3 times smaller than previous methods. We examine the recovery of the initial distribution of Fourier components and find the scale at which the recovered phases are scrambled with respect their true initial values. Integrating the Zel'dovich continuity equation back in time, we can improve the spatial resolution of the reconstruction by a factor of about 2.

Gramann, Mirt↗

Conceptual and Mathematical Foundation for the FE/NETL CO 2 Prophet Model for Simulating CO 2 Enhanced Oil Recovery, Version 2

The United States (U.S.) Department of Energy (DOE) Office of Fossil Energy (FE) at the National Energy Technology Laboratory (NETL) has developed the FE/NETL CO 2 Prophet Model, which is Version 2 of CO 2 Prophet. This document describes the mathematical foundation for the FE/NETL CO 2 Prophet Model. CO 2 Prophet was originally developed in the 1990s by Texaco Exploration and Production and Technology Department for DOE. The FE/NETL CO 2 Prophet Model is an oil reservoir simulator that is suitable for simulating water floods and supercritical carbon dioxide (CO 2 ) enhanced oil recovery (EOR). The FE/NETL CO 2 Prophet Model uses a number of assumptions to simplify the equations describing the flow of oil, water (or brine), and CO 2 in the oil reservoir.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

User’s Manual for StrmtbFlow, the Stream Tube Multiphase Flow Part of the FE/NETL CO 2 Prophet Model, Version 2

The United States (U.S.) Department of Energy (DOE) Office of Fossil Energy (FE) at the National Energy Technology Laboratory (NETL) has developed the FE/NETL CO 2 Prophet Model, which is an updated version of the CO 2 Prophet. This document is a user’s manual for StrmtbFlow, which is one part of the FE/NETL CO 2 Prophet Model. CO 2 Prophet was originally developed in the 1990s by Texaco Exploration and Production Technology Department for DOE. The FE/NETL CO 2 Prophet Model is an oil reservoir simulator that is suitable for simulating water floods and supercritical carbon dioxide (CO 2 ) enhanced oil recovery (EOR). The FE/NETL CO 2 Prophet Model uses a number of assumptions to simplify the equations describing the flow of oil, water (or brine), and CO 2 in the oil reservoir.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

User’s Manual for StrmtbGen, the Stream Tube Generating Part of the FE/NETL CO 2 Prophet Model, Version 2

The United States (U.S.) Department of Energy (DOE) Office of Fossil Energy (FE) at the National Energy Technology Laboratory (NETL) has developed the FE/NETL CO 2 Prophet Model, which is Version 2 of CO 2 Prophet. This document is a User’s Manual for StrmtbGen, which is one part of the FE/NETL CO 2 Prophet Model. CO 2 Prophet was originally developed in the 1990s by Texaco Exploration and Production Technology Department for DOE. The FE/NETL CO 2 Prophet Model is an oil reservoir simulator that is suitable for simulating water floods and supercritical carbon dioxide (CO 2 ) enhanced oil recovery (EOR). The FE/NETL CO 2 Prophet Model uses a number of assumptions to simplify the equations describing the flow of oil, water (or brine), and CO 2 in the oil reservoir.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Anisotropic constitutive modeling for nickel base single crystal superalloy Rene N4 at 982 C

A back stress/drag stress constitutive model based on a crystallographic approach to model single crystal anisotropy is presented. Experimental results demonstrated the need for the back stress variable in the inelastic flow equations. Experimental findings suggested that back stress is orientation dependent and controls both strain hardening and recovery characteristics. Due to the observed stable fatigue loops at 1800 F, drag stress is considered constant for this temperature. The constitutive model operated with constraints determined only from tensile data was extensively tested from simple tensile and fatigue to complicated strain hold tests. The model predicted very well under those conditions.

Stouffer, D. C.↗

Experimental Determination of the Recovery Factor and Analytical Solution of the Conical Flow Field for a 20 deg Included Angle Cone at Mach Numbers of 4.6 and 6.0 and Stagnation Temperatures to 2600 degree R

The local recovery factor was determined experimentally along the surface of a thin-walled 20 deg included angle cone for Mach numbers near 6.0 at stagnation temperatures between 1200 deg R and 2600 deg R. In addition, a similar cone configuration was tested at Mach numbers near 4.5 at stagnation temperatures of approximately 612 deg R. The local Reynolds number based on flow properties at the edge of the boundary layer ranged between 0.1 x 10(exp 4) and 3.5 x 10(exp 4) for tests at temperatures above 1200 deg R and between 6 x 10(exp 4) and 25 x 10(exp 4) for tests at temperatures near 612 deg R. The results indicated, generally, that the recovery factor can be predicted satisfactorily using the square root of the Prandtl number. No conclusion could be made as to the necessity of evaluating the Prandtl number at a reference temperature given by an empirical equation, as opposed to evaluating the Prandtl number at the wall temperature or static temperature of the gas at the cone surface. For the tests at temperatures above 1200 deg R (indicated herein as the tests conducted in the slip-flow region), two definite trends in the recovery data were observed - one of increasing recovery factor with decreasing stagnation pressure, which was associated with slip-flow effects and one of decreasing recovery factor with increasing temperature. The true cause of the latter trend could not be ascertained, but it was shown that this trend was not appreciably altered by the sources of error of the magnitude considered herein. The real-gas equations of state were used to determine accurately the local stream properties at the outer edge of the boundary layer of the cone. Included in the report, therefore, is a general solution for the conical flow of a real gas using the Beattie-Bridgeman equation of state. The largest effect of temperature was seen to be in the terms which were dependent upon the internal energy of the gas. The pressure and hence the pressure drag terms were unaffected.

Pfyl, Frank A.↗

Probabilistic neural networks for fluid flow surrogate modeling and data recovery

We consider the use of probabilistic neural networks for fluid flow surrogate modeling and data recovery. This framework is constructed by assuming that the target variables are sampled from a Gaussian distribution conditioned on the inputs. Consequently, the overall formulation sets up a procedure to predict the hyperparameters of this distribution which are then used to compute an objective function given training data. We demonstrate that this framework has the ability to provide for prediction confidence intervals based on the assumption of a probabilistic posterior, given an appropriate model architecture and adequate training data. The applicability of the present framework to cases with noisy measurements and limited observations is also assessed. To demonstrate the capabilities of this framework, we consider canonical regression problems of fluid dynamics from the viewpoint of reduced-order modeling and spatial data recovery for four canonical data sets. The examples considered in this study arise from (i) the shallow-water equations, (ii) a two-dimensional cylinder flow, (iii) the wake of a NACA0012 airfoil with a Gurney flap, and (iv) the NOAA sea surface temperature data set. Furthermore, the present results indicate that the probabilistic neural network not only produces a machine-learning-based fluid flow surrogate model but also systematically quantifies the uncertainty therein to assist with model interpretability.

42 ENGINEERING↗

Optimization Framework to Assess the Demand Response Capacity of a Water Distribution System

As large electricity consumers, water distribution system (WDS) pumping stations have the potential to become meaningful participants in demand response (DR) programs. The authors propose an optimization framework for assessing the DR capacity of a WDS and identifying the optimal bidding strategy for maximizing WDS revenue in the DR spot market. The proposed mixed integer linear programming (MILP) model overcomes computational constraints of previous DR optimization models by adopting a preprocessing procedure to minimize the number of binary variables and implementing a convex relaxation technique to linearize the hydraulic equations. The proposed MILP model also explicitly accounts for varying levels of risk tolerance of WDS operators by varying the recovery period over which pumping returns to business-as-usual operation. The optimization framework is implemented on a skeletonized 48-node WDS model that includes 7 pumps, 6 tanks, and 39 pipes. Using a simulated DR event and water consumption profile, the authors derive the optimal DR supply curves (i.e., compensation price versus load curtailment quantity) and revenue potential of the WDS under six scenarios for DR participation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

Some implications for cyclic plastic and viscoplastic equations based on nonproportional loading experiments

From cyclic, strain controlled, nonproportional tests on type 304 stainless steel and Hastelloy-X, the following statements may be made: (1) A dynamic recovery term is essential to properly model the backstress evolution. (2) From analysis of Hastelloy-X data obtained at 649 C, the inelastic strain rate appears to be a satisfactory directional index for direct hardening, but the backstress appears to be an inappropriate directional index of dynamic recovery. (3) Sinusoidal, 90 deg out-of-phase axial torsional tests can be very useful in aiding determination of backstress evolution functions, including both directional indices and scalar hardening functions, by virtue of the associated approximately constant magnitudes of overstress, inelastic strain rate, and effective stress. Such tests have previously been associated with the study of nonproportional hardening effects but have more far ranging applications.

Mcdowell, D. L.↗

Inverse point source location with the Helmholtz equation on a bounded domain

The problem of recovering acoustic sources, more specifically monopoles, from point-wise measurements of the corresponding acoustic pressure at a limited number of frequencies is addressed. To this purpose, a family of sparse optimization problems in measure space in combination with the Helmholtz equation on a bounded domain is considered in this study. A weighted norm with unbounded weight near the observation points is incorporated into the formulation. Optimality conditions and conditions for recovery in the small noise case are discussed, which motivates concrete choices of the weight. The numerical realization is based on an accelerated conditional gradient method in measure space and a finite element discretization.

97 MATHEMATICS AND COMPUTING↗

Performance study of a novel dew point evaporative cooler in the climate of central Europe using building simulation tools

This paper presents the performance research of a hybrid air conditioning system (System II) equipped with novel dew point evaporative cooler (DPIEC) for a retail building located in temperate climatic conditions. The new cooler has a dedicated, adaptable structure that maximizes its performance in the climate conditions of central Europe. The unit can operate as a typical regenerative air cooler when outdoor air is drier than indoor air, it can operate as a counter-flow heat recovery unit when indoor conditions are more dry than outdoor and it can operate as a combination of both arrangements. The year-round, hourly-stepped building energy simulations were carried-out using a prognostic tool of energy consumption. The novel dew point evaporative cooler application potential was established using the black-box model based on regression equations was established based on experimental tests on a prototype of novel DPIEC. The hybrid system operation was compared to the typical air handling unit equipped with a standard energy wheel for heat recovery (System I). Presented research shows that DPIEC can cover about 95% of total cooling loads while the energy wheel covers only 9% of total cooling load. System II allows us to save 65% of seasonal electricity consumption comparing to System I. It comes from the fact that during most of the cooling season DPIEC is in operation and it covers most of the cooling loads without the need to switch the chilled water system.

Air conditioning↗

Mapping Incidence and Prevalence Peak Data for SIR Modeling Applications

Infectious disease modeling and forecasting have played a key role in helping assess and respond to epidemics and pandemics. Recent work has leveraged data on disease peak infection and peak hospital incidence to fit compartmental models for the purpose of forecasting and describing the dynamics of a disease outbreak. Incorporating these data can greatly stabilize a compartmental model fit on early observations, where slight perturbations in the data may lead to model fits that forecast wildly unrealistic peak infection. We introduce a new method for incorporating historic data on the value and time of peak incidence of hospitalization into the fit for a Susceptible-Infectious-Recovered (SIR) model by formulating the relationship between an SIR model’s starting parameters and peak incidence as a system of two equations that can be solved computationally. We demonstrate how to calculate SIR parameter estimates – which describe disease dynamics such as transmission and recovery rates – using this method, and determine that there is a noticeable loss in accuracy whenever prevalence data is misspecified as incidence data. To exhibit the modeling potential, we update the Dirichlet-Beta State Space modeling framework to use hospital incidence data, as this framework was previously formulated to incorporate only data on total infections. This approach is assessed for practicality in terms of accuracy and speed of computation via simulation.

97 MATHEMATICS AND COMPUTING↗

NASTRAN finite element idealization study

The investigation of the effects of variations of mesh refinement and mesh pattern were conducted using a basic rectangular mesh pattern. When employing the constant strain TRMEM element, the basic rectangular pattern was subdivided into triangles. This subdivision employs two different triangular patterns and allows results to be obtained which demonstrate the effect of modelling bias. Errors in tip deflection, direct stress, and shearing stress as a function of mesh size and element aspect ratio were obtained as well as mid-span stress distributions. All problems were solved on an IBM 360/95 computer using MacNeal-Schwendler Version MSC-38 Rigid Format-1 of the NASTRAN computer program. While NASTRAN uses double precision arithmetic for the solution of the global equations for displacements, subsequent computations to obtain element stresses are carried out in single precision. This suggests that some improvement in stress recovery might be expected when using the higher precision DCD machines.

Case, W. R.↗

Viscoplasticity based on overstress with a differential growth law for the equilibrium stress

Two coupled, nonlinear differential equations are proposed for the modeling of the elastic and rate (time) dependent inelastic behavior of structural metals in the absence of recovery and aging. The structure of the model is close to the unified theories but contains essential differences. It is shown that the model reproduces almost elastic regions upon initial loading and in the unloading regions of the hysteresis loop. Under loading, unloading and reloading in strain control the model simulated the experimentally observed sharp transition from nearly elastic to inelastic behavior. When a formulation akin to existing unified theories is adopted the almost elastic regions reduce the points and the transition upon reloading is very gradual. For different formulations the behavior under sudden in(de)creases of the strain rate by two orders of magnitude is simulated by numerical experiments and differences are noted. The model represents cyclically neutral behavior and contains three constants and two positive, decreasing functions. The determination of constants and functions from monotonic loading with strain rate changes and relaxation periods is described.

Krempl, E.↗

Viscoplasticity based on overstress with a differential growth law for the equilibrium stress

Two coupled, nonlinear differential equations are proposed for the modeling of the elastic and rate (time) dependent inelastic behavior of structural metals in the absence of recovery and aging. The structure of the model is close to the unified theories but contains essential differences. It is shown that the model reproduces almost elastic regions upon initial loading and in the unloading regions of the hysteresis loop. Under loading, unloading and reloading in strain control the model simulated the experimentally observed sharp transition from nearly elastic to inelastic behavior. When a formulation akin to existing unified theories is adopted the almost elastic regions reduce the points and the transition upon reloading is very gradual. For different formulations the behavior under sudden in(de)creases of the strain rate by two orders of magnitude is simulated by numerical experiments and differences are noted. The model represents cyclically neutral behavior and contains three constants and two positive, decreasing functions. The determination of constants and functions from monotonic loading with strain rate changes and relaxation periods is described.

Krempl, E.↗

Vortex Generators in a Streamline-Traced, External-Compression Supersonic Inlet

Vortex generators within a streamline-traced, external-compression supersonic inlet for Mach 1.66 were investigated to determine their ability to increase total pressure recovery and reduce total pressure distortion. The vortex generators studied were rectangular vanes arranged in counter-rotating and co-rotating arrays. The vane geometric factors of interest included height, length, spacing, angle-of-incidence, and positions upstream and downstream of the inlet terminal shock. The flow through the inlet was simulated numerically through the solution of the steady-state, Reynolds-averaged Navier-Stokes equations on multi-block, structured grids using the Wind-US flow solver. The vanes were simulated using a vortex generator model. The inlet performance was characterized by the inlet total pressure recovery and the radial and circumferential total pressure distortion indices at the engine face. Design of experiments and statistical analysis methods were applied to quantify the effect of the geometric factors of the vanes and search for optimal vane arrays. Co-rotating vane arrays with negative angles-of-incidence positioned on the supersonic diffuser were effective in sweeping low-momentum flow from the top toward the sides of the subsonic diffuser. This distributed the low-momentum flow more evenly about the circumference of the subsonic diffuser and reduced distortion. Co-rotating vane arrays with negative angles-of-incidence or counter-rotating vane arrays positioned downstream of the terminal shock were effective in mixing higher-momentum flow with lower-momentum flow to increase recovery and decrease distortion. A strategy of combining a co-rotating vane array on the supersonic diffuser with a counter-rotating vane array on the subsonic diffuser was effective in increasing recovery and reducing distortion.

computational fluid dynamics↗