Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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 595 records · Page 33

A TCN-Based Hybrid Forecasting Framework for Hours-Ahead Utility-Scale PV Forecasting

This paper presents a Temporal Convolutional Network (TCN) based hybrid PV forecasting framework for enhancing hours-ahead utility-scale PV forecasting. The hybrid framework consists of two forecasting models: a physics-based trend forecasting (TF) model and a data-driven fluctuation forecasting (FF) model. Three TCNs are integrated in the framework for: i) blending the inputs from different Numerical Weather Prediction sources for the TF model to achieve superior performance on forecasting hourly PV profiles, ii) capturing spatial-temporal correlations between detector sites and the target site in the FF model to achieve more accurate forecast of intra- hour PV power drops, and iii) reconciling TF and FF results to obtain coherent hours-ahead PV forecast with both hourly trends and intra-hour fluctuations well preserved. To automatically identify the most contributive neighboring sites for forming a detector network, a scenario-based correlation analysis method is developed, which significantly improves the capability of the FF model on capturing large power fluctuations caused by cloud movements. Here, the framework is developed, tested, and validated using actual PV data collected from 95 PV farms in North Carolina. Simulation results show that the performance of 6 hours ahead PV power forecasting is improved by 20% - 30% compared with state-of-the-art methods.

42 ENGINEERING↗

Atmospheric Profile Builder

APBuilder is a software tool designed to generate atmospheric profiles for use with AC2Dr by downloading weather model data and transforming it into 1D or 2D binary profiles. These binary profiles serve as required inputs for AC2Dr, an LLNL-developed, open source, two-dimensional numerical solver for the acoustic wave equation based on the finite difference method.

Kim, Keehoon [Lawrence Livermore National Laborato↗

Climate Influences on Capacity Expansion Planning with Application to the Western U.S

Electric power system planners utilize a variety of planning tools to inform decisions concerning generation and transmission additions to the electric grid, the need for operational changes, and to evaluate potential stressors on the system. Numerous factors contribute to the planning process including projected fuel and technology costs, policy and load profiles. There is also a growing recognition of the interdependency of the electric grid with other natural and engineered systems. Here we explore how future climate change and hydropower operability might influence decisions related to electricity capacity expansion planning and operations. To do so we assemble a multi-model framework. Specifically, water resource modeling is used to simulate climate impacts on future water supply for thermoelectric and hydropower generation. Separately, temperature impacts on electricity load are evaluated. Together, these climate factors spatially constrain a capacity expansion model that projects generation and transmission additions to the grid. The projected new capacity-builds are then evaluated on their operations, reliability, and cost under average and extreme climate conditions using production cost modeling. This coupled framework is demonstrated on the electric grid in the Western U.S., supporting capacity expansion planning by WECC, the North American Electric Reliability Corporation (NERC) regional entity responsible for reliability assurance of the Western Interconnection. This region was selected in part because the West is unique in that it has high potential for renewable penetrations and is experiencing large retirements/displacements of baseload resources, primarily coal, leading to possible operational challenges in terms of changing resource mix and the need for resource flexibility. Toward this challenge, planning scenarios encompass a range of alternative energy, climate and drought futures. In this context we explore answers to two strategic questions: 1) How does changing climate influence electricity expansion planning (generation and transmission) and future operations, including type and capacity of new builds, system reliability, cost and environmental impacts? 2) How does the representation of hydropower in the modeling framework influence the evaluation of bulk power system operations? Results indicate that climate has a measurable influence on recommendations concerning the capacity, type and location of new generation and transmission additions, with up to 17 GW additional capacity needed by 2038 to meet peak loads (~6.6% increase over capacity-builds based on historical climate). The extent of additional infrastructure needs is strongly influenced by future water availability for hydropower and the potential deployment of demand response technologies. Systems designed for future climate conditions were found to maintain high system reliability under a range of electricity and water availability scenarios (including significant drought), with minimal system curtailments. Additional capacity needs due to higher load tend to increase cumulative 20-year investment and operating costs by $\$$5-$\$$17 billion and generation costs increase by 9 to 19%. Finally, changing the representation of hydropower flexibility has a relatively small influence on capacity expansion in the Western Interconnection through 2038, but hydropower flexibility impacts generation costs to a similar extent as climate.

13 HYDRO ENERGY↗

Hydrodynamic Coupling to a Homogenized Radiation Transport Method based on Young Measures

Resolving radiation transport fields subject to opacity profiles with strong, oscillatory line structure while potentially falling under intermediate optical depth conditions presents a numerical challenge in radiation transport modeling. The Young measure-based homogenization technique formulated by Haut et al. (2017) was investigated as a candidate method for resolving radiation fields under such conditions more accurately. The method was compared against frequently-utilized mean opacity methods as the Rosseland and Planck formulations. In this work, all methods were tested through radiation slab calculations separately comprised of aluminum, copper, and krypton, each for different thermodynamic conditions. Following these offline radiation slab calculations, demonstrations shifted towards the SCEPTRE radiation transport code and, subsequently, the multiphysics ALEGRA code for approximately-coupled radiation-material simulations. Throughout all the simulations shown in this study, for a fixed computational cost, the homogenized method was observed to be more accurate than any of the solutions determined through traditional mean opacity approaches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SOWFA Large-Eddy Simulation of a Diurnal Cycle at the SWiFT Site

These are large-eddy simulation (LES) data output from the Simulator for Wind Farm Applications (SOWFA) with internal mesoscale-to-microscale coupling. A diurnal cycle was simulated over the Scaled Wind Farm Technology (SWIFT) facility near Lubbock, Texas, from 2013-11-08 12:00:00 UTC to 2013-11-09 12:00:00 UTC. Indirect profile assimilation was used to enforce nonstationary mesoscale conditions. Mesoscale forcing data were extracted from WRF numerical weather prediction outputs or derived from observations at SWiFT.

17 WIND ENERGY↗

Turbulence dissipation rate estimated from lidar observations during the LAPSE-RATE field campaign

The International Society for Atmospheric Research using Remotely-piloted Aircraft (ISARRA) hosted a flight week in July 2018 to demonstrate unmanned aircraft systems' (UASs) capabilities in sampling the atmospheric boundary layer. This week-long experiment was called the Lower Atmospheric Profiling Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE) field campaign. Numerous remotely piloted aircraft and ground-based instruments were deployed with the objective of capturing meso- and microscale phenomena in the atmospheric boundary layer. The University of Oklahoma deployed one Halo Streamline lidar, and the University of Colorado Boulder deployed two WindCube lidars. In this paper, we use data collected from these Doppler lidars to estimate turbulence dissipation rate throughout the campaign. We observe large temporal variability of turbulence dissipation close to the surface with the WindCube lidars that is not detected by the Halo Streamline. However, the Halo lidar enables estimating dissipation rate within the whole boundary layer, where a diurnal variability emerges. We also find a higher correspondence in turbulence dissipation between the WindCube lidars, which are not co-located, compared to the Halo and WindCube lidar that are co-located, suggesting a significant influence of measurement volume on the retrieved values of dissipation rate.

WIND ENERGY↗

Time dependent shear stress and temperature distribution over an insulated flat plate moving at hypersonic speed.

The laminar two-dimensional flow over a stepwise accelerated flat plate moving with hypersonic speed at zero angle of attack is analysed. The governing equations in the self-similar form are linearized and solved numerically for small times. The solutions obtained are the deviations of the velocity and the temperature profiles from those of steady state. The presented results may be used to find the first order boundary layer induced pressure on the plate.

Rodkiewicz, C. M.↗

Numerical solution of the Navier-Stokes equations for arbitrary two-dimensional airfoils

A method for numerical solution of the Navier-Stokes equations for the flow about arbitrary airfoils or other bodies is presented. This method utilizes a numerically generated curvilinear coordinate system having a coordinate line coincident with the body contour. Streamlines, velocity profiles, and pressure and force coefficients for several airfoils and an arbitrary rock are given. Potential flow solutions are also presented. The procedure capable of treating multiple-element airfoils, and potential flow results are presented.

Thames, F. C.↗

FORTRAN subroutines for out-of-core solutions of large complex linear systems

The design and usage of two main subprograms using direct methods to solve large linear complex systems, of the form Ax = b, whose coeffficient matrices are too large to be stored in core are described. The first main subprogram is for systems whose coefficient matrices are of a particular sparse structure, namely, the matrix A can be written in the form B + D, where B is a block-banded system, and D has only a few columns of nonzeros. Key elements of the algorithms used in the subprograms include: the data structure, the strategy for preserving numerical stability, the adaptability of the algorithms for dense systems as well as for block-profile systems.

Yip, E. L.↗

Turbulent solution of the Navier-Stokes equations for an inhomogenous developing shear layer

To study the nonlinear physics of inhomogeneous turbulent shear flow, the unaveraged Navier-Stokes equations are solved numerically. For initial conditions a three-dimensional cosine velocity fluctuation and a mean-velocity profile with a step are used. Although the initial conditions are nonrandom. The flow soon becomes turbulent. Concentrated turbulent energy develops near the plane where the mean velocity gradient is initially infinite. The terms in the one-point correlation equation for turbulent energy, including those for the diffusion and production of turbulence, are calculated, the diffusion terms tend to make the turbulence more homogeneous.

Deissler, R. G.↗

Kinetic theory model for the flow of a simple gas from a two-dimensional nozzle

A system of nonlinear integral equations equivalent to the Krook kinetic equation for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust two-dimensional nozzles. The method of characteristics was used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it exhausts into a vacuum. Results predict backscatter and show the effect of the inside wall boundary layer on the flowfields external to the nozzle.

Riley, B. R.↗

An observational and theoretical study of highly supercooled altocumulus

The microphysics of two altocumulus clouds sampled at a temperature of about -30 C through the NCAR King Air aircraft and coincident lidar has been characterized. The clouds are structurally similar to stratocumulus and have extensive cloudtop entrainment, a capping temperature inversion, and a dry layer above. Microphysical and radiative properties of both clouds have been modeled numerically. When entrainment effects are incorporated in the model, calculations of droplet concentration and mean diameter profiles compare favorably with the measurements. Radiative transfer calculations indicate that radiative cooling causes sufficient negative buoyancy in cloudtop parcels for producing convective instability and reproducing the observed downdraft velocities. It is shown that entrainment of warmer, drier air near a cloudtop counteracts the radiatively induced negative buoyancy in the downdrafts.

Heymsfield, Andrew J.↗

Optimization and analysis of large chemical kinetic mechanisms using the solution mapping method - Combustion of methane

A method of systematic optimization, solution mapping, as applied to a large-scale dynamic model is presented. The basis of the technique is parameterization of model responses in terms of model parameters by simple algebraic expressions. These expressions are obtained by computer experiments arranged in a factorial design. The developed parameterized responses are then used in a joint multiparameter multidata-set optimization. A brief review of the mathematical background of the technique is given. The concept of active parameters is discussed. The technique is applied to determine an optimum set of parameters for a methane combustion mechanism. Five independent responses - comprising ignition delay times, pre-ignition methyl radical concentration profiles, and laminar premixed flame velocities - were optimized with respect to thirteen reaction rate parameters. The numerical predictions of the optimized model are compared to those computed with several recent literature mechanisms. The utility of the solution mapping technique in situations where the optimum is not unique is also demonstrated.

Frenklach, Michael↗

Navier-Stokes simulation of flow through a highly contoured subsonic diffuser

The flow through a highly offset subsonic diffuser with cross-sectional profiles that varied from rectangular at the duct entrance to circular at the engine face was numerically simulated. A multizonal approach combined with a two-grid topology was used to represent both the internal and external flowfields, and an implicit, approximately-factored, partially flux-split finite-difference algorithm was used to solve the three-dimensional thin-layer Navier-Stokes equations. The computed static pressures along the inlet wall and total pressures on the engine face were compared with experimental data. In addition, the overall flowfield within the duct was examined in detail. Good agreement is shown between experiment and computations, with the limiting factor being the lack of a reliable turbulence model for internal flow problems.

Chyu, W. J.↗

Kinetic theory model for the flow of a simple gas from a three-dimensional axisymmetric nozzle

A system of nonlinear integral equations equivalent to the Krook kinetic equations for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust three-dimensional axisymmetric nozzles. The method of characteristics is used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it expands into a vacuum. Results predict backscatter and show the effect of the nozzle wall boundary layer on the external flowfields.

Riley, B. R.↗

Experimental and numerical investigation of supersonic turbulent flow in an annular duct

Experimental and numerical results are presented for developing supersonic turbulent flow in an annular duct formed by a circular centerbody and outer shroud. The experimental results are based on data taken in a new flow facility that was designed to generate a shock-free, supersonic annular flow. Numerical computations were performed using the Baldwin-Lomax turbulence model for comparison with experimentally measured profiles. The results demonstrate that computed and measured profiles are in excellent agreement, so that studies can now be conducted of shock wave/boundary layer interaction phenomena within the duct, such as those induced by changes in downstream duct geometry or by the placement of struts between the duct walls.

Williams, K. E.↗

A new methodology to determine kinetic parameters for one- and two-step chemical models

In this paper, a new methodology to determine kinetic parameters for simple chemical models and simple transport properties classically used in DNS of premixed combustion is presented. First, a one-dimensional code is utilized to performed steady unstrained laminar methane-air flame in order to verify intrinsic features of laminar flames such as burning velocity and temperature and concentration profiles. Second, the flame response to steady and unsteady strain in the opposed jet configuration is numerically investigated. It appears that for a well determined set of parameters, one- and two-step mechanisms reproduce the extinction limit of a laminar flame submitted to a steady strain. Computations with the GRI-mech mechanism (177 reactions, 39 species) and multicomponent transport properties are used to validate these simplified models. A sensitivity analysis of the preferential diffusion of heat and reactants when the Lewis number is close to unity indicates that the response of the flame to an oscillating strain is very sensitive to this number. As an application of this methodology, the interaction between a two-dimensional vortex pair and a premixed laminar flame is performed by Direct Numerical Simulation (DNS) using the one- and two-step mechanisms. Comparison with the experimental results of Samaniego et al. (1994) shows a significant improvement in the description of the interaction when the two-step model is used.

Mantel, T.↗

Percentage Contributions from Atmospheric and Surface Features to Computed Brightness Temperatures

Over the past few years, there has become an increasing interest in the use of millimeter-wave (mm-wave) and sub-millimeter-wave (submm-wave) radiometer observations to investigate the properties of ice particles in clouds. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. One methodology used since the 1960's to discern the relationship between the physical state observed and the brightness temperature (TB) is through the temperature weighting function profile. In this research, the temperature weighting function concept is exploited to analyze the sensitivity of various characteristics of the cloud profile, such as relative humidity, ice water path, liquid water path, and surface emissivity. In our numerical analysis, we compute the contribution (in Kelvin) from each of these cloud and surface characteristics, so that the sum of these various parts equals the computed TB. Furthermore, the percentage contribution from each of these characteristics is assessed. There is some intermingling/contamination of the contributions from various components due to the integrated nature of passive observations and the absorption and scattering between the vertical layers, but all in all the knowledge gained is useful. This investigation probes the sensitivity over several cloud classifications, such as cirrus, blizzards, light snow, anvil clouds, and heavy rain. The focus is on mm-wave and submm-wave frequencies, however discussions of the effects of cloud variations to frequencies as low as 10 GHz and up to 874 GHz will also be presented. The results show that nearly 60% of the TB value at 89 GHz comes from the earth's surface for even the heaviest blizzard snow rates. On the other hand, a significant percentage of the TB value comes from the snow in the cloud for 166, and 183 plus or minus 7 GHz for the heavy and medium snow rates. For submm-wave channels, there is no contribution from the surface because these channels cannot probe through clouds, nor normal water vapor amounts in clear air regions. This work is extremely valuable in physically-based retrieval algorithm development research.

Gail Skofronick-Jackson↗