Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW MODELS”

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 649 records · Page 36

Flap Side-Edge Noise: Acoustic Analysis of Sen's Model

The two-dimensional flap side-edge flow model developed by Sen is analyzed to reveal the noise production potential of the proposed mechanism. The model assumes that a vortex will form at the equilibrium position off the side edge of the flap. The vortex is then perturbed away from the equilibrium position by incoming turbulence causing it to oscillate and thus radiate sound. The noise field is calculated three-dimensionally by taking the flap to have a finite chord. Spectra and directivity of the farfield sound are presented. In addition, the effect of retarded time differences is evaluated. The parameters in the model are related to typical aircraft parameters and noise reduction possibilities are proposed.

Hardin, Jay C.↗

A Review of Computational Models for the Flow of Milled Biomass Part I: Discrete-Particle Models

Biomass is a renewable and sustainable energy resource. Current design of biomass handling and feeding equipment leverage both experiments and numerical modeling. This paper reviews the state-of-the-art discrete element methods (DEM) for the flow of milled biomass (Part I), accompanied by a comprehensive review on continuum-based computational models (Part II). The present review on DEM is primarily focused on the features and suitability of various particle shape models for different types of milled biomass because particle shape is the predominant attribute controlling the flow behavior of complex-shaped granular material. The general strengths and weaknesses in the applicability of those models for the milled biomass modeling are summarized. In particular, comments are provided to balance the numerical model capabilities and the computational cost for the development of DEM models. To our best knowledge, this is the first-of-its-kind review on DEM specifically for biomass. Our study indicates that the current DEM models require further development, calibration, and validation based on a deep understanding of biomass particle contact mechanics and experimental data support before they can be reliably used for predictive simulations in handling and feeding systems.

42 ENGINEERING↗

Centrifugal and Axial Pump Design and Off-Design Performance Prediction

A meanline pump-flow modeling method has been developed to provide a fast capability for modeling pumps of cryogenic rocket engines. Based on this method, a meanline pump-flow code PUMPA was written that can predict the performance of pumps at off-design operating conditions, given the loss of the diffusion system at the design point. The design-point rotor efficiency and slip factors are obtained from empirical correlations to rotor-specific speed and geometry. The pump code can model axial, inducer, mixed-flow, and centrifugal pumps and can model multistage pumps in series. The rapid input setup and computer run time for this meanline pump flow code make it an effective analysis and conceptual design tool. The map-generation capabilities of the code provide the information needed for interfacing with a rocket engine system modeling code. The off-design and multistage modeling capabilities of PUMPA permit the user to do parametric design space exploration of candidate pump configurations and to provide head-flow maps for engine system evaluation.

Veres, Joseph P.↗

Artificial dissipation models for hypersonic external flow

Four artificial dissipation models which augment central difference schemes were examined for hypersonic external flows. The models were a first and third order dissipation model, a directionally scaled first and third order dissipation model, a flux limited dissipation model, and a flux difference split dissipation model. Each model was implemented in the lower-upper symmetric-Gauss-Seidel (LU-SGS) algorithm to solve the full Navier-Stokes equations. The latter two models can be regarded as total variation diminishing (TVD) schemes. Test results for model problems showed that the flux limited dissipation model was robust enough to predict a high speed blunt body flow with strong shock and expansion waves. The flux difference split dissipation model was capable of shock capturing with higher resolution, but was less robust. First and third order dissipation models turned out to be neither accurate nor robust enough for high Mach number flow computations.

Yoon, Seokkwan↗

On the Thermal Model of Transverse Flow of Unidirectional Materials

The thermal model for transverse heat flow of having single filament in a unit cell is extended. In this model, we proposed that two circular filaments in a unit cell of square packing array and obtained the transverse thermal conductivity of an unidirectional material.

Tai, Hsiang↗

Longitudinal Stability and Control Characteristics of a Semispan Wind-Tunnel Model of the XF7U-1 Airplane and a Comparison with Complete-Model Wind-Tunnel Tests and Semispan-Model Wing-Flow Tests

An investigation was conducted on an 0.08-scale semispan model of the Chance Vought XF7U-1 airplane in the Langley high-speed 7- by 10-foot tunnel in the Mach number range from 0.40 to 0.97. The results are compared with those obtained with an 0.08-scale sting-mounted complete model tested in the same tunnel and with an 0.026-scale semispan model tested by the wing-flow method. The lift-curve slopes obtained for the 0.08-scale semispan model and the 0.026-scale wing-flow model were in good agreement but both were generally lower than the values obtained for the sting model. The results of an unpublished investigation have shown that tunnel-wall boundary-layer and strut-leakage effects can came the difference noted between the lift-curve slopes of the sting and the semispan data. Fair agreement was obtained among the data of the three models as regard the variation of pitching-moment coefficients with lift coefficient. The agreement between the complete and the semispan models was more favorable with the vertical fine on, because the wall-boundary-layer and strut leakage effects were less severe. In the Mach number range between 0.94 and 0.97, ailavator-control reversal was indicated in the wing-flow data near zero lift; Whereas, these same trends were indicated in the larger scale semispan data at somewhat higher lift coefficients.

Goodson, Kenneth W.↗

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗

Solar wind eddies and the heliospheric current sheet

Ulysses has collected data between 1 and 5 AU during, and just following solar maximum, when the heliospheric current sheet (HCS) can be thought of as reaching its maximum tilt and being subject to the maximum amount of turbulence in the solar wind. The Ulysses solar wind plasma instrument measures the vector velocity and can be used to estimate the flow speed and direction in turbulent 'eddies' in the solar wind that are a fraction of an astronomical unit in size and last (have either a turnover or dynamical interaction time of) several hours to more than a day. Here, in a simple exercise, these solar wind eddies at the HCS are characterized using Ulysses data. This character is then used to define a model flow field with eddies that is imposed on an ideal HCS to estimate how the HCS will be deformed by the flow. This model inherently results in the complexity of the HCS increasing with heliocentric distance, but the result is a measure of the degree to which the observed change in complexity is a measure of the importance of solar wind flows in deforming the HCS. By comparison with randomly selected intervals not located on the HCS, it appears that eddies on the HCS are similar to those elsewhere at this time during the solar cycle, as is the resultant deformation of the interplanetary magnetic field (IMF). The IMF deformation is analogous to what is often termed the 'random walk' of interplanetary magnetic field lines.

Suess, S. T.↗

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗

Trace Chemistry

The goals of the trace chemistry group were to identify the processes relevant to aerosol and aerosol precursor formation occurring within aircraft gas turbine engines; that is, within the combustor, turbine, and nozzle. The topics of discussion focused on whether the chemistry of aerosol formation is homogeneous or heterogeneous; what species are important for aerosol and aerosol precursor formation; what modeling/theoretical activities to pursue; what experiments to carry out that both support modeling activities and elucidate fundamental processes; and the role of particulates in aerosol and aerosol precursor formation. The consensus of the group was that attention should be focused on SO2, SO3, and aerosols. Of immediate concern is the measurement of the concentration of the species SO3, SO2, H2SO4 OH, HO2, H2O2, O, NO, NO2, HONO, HNO3, CO, and CO2 and particulates in various engines, both those currently in use and those in development. The recommendation was that concentration measurements should be made at both the combustor exit and the engine exit. At each location the above species were classified into one of four categories of decreasing importance, Priority I through IV, as follows: Combustor exit: Priority I species - SO3:SO2 ratio, SO3, SO2, and particulates; Priority II species: OH and O; Priority III species - NO and NO2; and Priority IV species - CO and CO2. For the Engine exit: Priority I species - SO3:SO2 ratio, SO3, SO2,H2SO4, and particulates; Priority II species: OH,HO2, H2O2, and O; Priority III species - NO, NO2, HONO, and HNO3; and Priority IV species - CO and CO2. Table I summarizes the anticipated concentration range of each of these species. For particulate matter, the quantities of interest are the number density, size distribution, and composition. In order to provide data for validating multidimensional reacting flow models, it would be desirable to make 2-D, time-resolved measurements of the concentrations of the above species and, in addition, of the pressure, temperature, and velocity. A near term goal of the experimental program should be to confirm the nonlinear effects of sulfur speciation, and if present, to provide an explanation for them. It is also desirable to examine if the particulate matter retains any sulfur. The recommendation is to examine the effects on SOx production of variations in fuel-bound sulfur and aromatic content (which may affect the amount of particulates formed). These experiments should help us to understand if there is a coupling between particulate formation and SO, concentration. Similarly, any coupling with NOx can be examined either by introducing NOx into the combustion air or by using fuel-bound nitrogen. Also of immediate urgency is the need to establish and validate a detailed mechanism for sulfur oxidation/aerosol formation, whose chemistry is concluded to be homogeneous, because there is not enough surface area for heterogeneous effects. It is envisaged that this work will involve both experimental and theoretical programs. The experimental work will require, in addition to the measurements described above, fundamental studies in devices such as flow reactors and shock tubes. Complementing this effort should be modeling and theoretical activities. One impediment to the successful modeling of sulfur oxidation is the lack of reliable data for thermodynamic and transport properties for several species, such as aqueous nitric acid, sulfur oxides, and sulfuric acid. Quantum mechanical calculations are recommended as a convenient means of deriving values for these properties. Such calculations would also help establish rate constants for several important reactions for which experimental measurements are inherently fraught with uncertainty. Efforts to implement sufficiently detailed chemistry into computational fluid dynamic codes should be continued. Zero- and one-dimensional flow models are also useful vehicles for elucidating the minimal set of species and reactions that must be included in two- and three-dimensional modeling studies.

Radhakrishnan, Krishnan↗

Calibrating uncertain parameters in melt pool simulations of additive manufacturing

Melt pool scale numerical modeling of additive manufacturing (AM) processes can provide predictive capabilities and theoretical insight into the process-property-structure-performance relationships for AM parts. Despite capabilities of numerical models to solve complex multi-physics problems, it is often important to consider a tradeoff between detailed physics and computational cost. Therefore, sources of uncertainty in both experimental conditions and the parameters needed for modeling require models to be validated against empirical evidence. Here, a method is proposed to calibrate uncertain parameters used in continuum-scale melt pool models for powder bed fusion (PBF) AM. Both a simplified heat transfer model and a heat transfer and fluid flow model were investigated. A surrogate model and Markov chain-based optimization algorithm calibrated melt pool geometry for models within experimental variation of the target melt pool width and depth from the NIST AM-Bench 2018-02 dataset. The melt pool temperature distributions, solidification parameters, and simulated multi-layer solidification microstructures were compared between the two models. Similar results from both models indicate that calibrated, lower fidelity numerical models may be used in place of higher fidelity models to generate melt pool solidification data. Finally, these calibrated models therefore enable lower computational cost melt pool simulations without a noticeable decrease in simulation accuracy for grain-scale microstructure simulations.

36 MATERIALS SCIENCE↗

Deposit geometry and oxygen concentration spatial variations due to composition change in printed functionally graded components

Functionally graded materials (FGMs) with site specific chemical composition are commonly manufactured by directed energy deposition (DED). Although previous work fabricated an FGM with a compositional variation between a ferritic and austenitic alloy, difficulties arose due to variations in deposit shape with composition change. This problem also occurs for FGMs in literature; however, unlike other cases, the thermophysical properties of these two alloys were similar throughout the build. Here, we investigate the role of chemical composition and surface-active elements on deposit geometry during the manufacture of FGMs by laser DED. Single-track experiments for the relevant FGM compositions are analyzed with results from a well-tested, three-dimensional, transient numerical heat transfer and fluid flow model and thermodynamic calculations. Experiments showed deposit shape varied as a function of composition for constant laser power and scan speed. Thermodynamic analysis indicated that the oxygen solubility in the fusion zone varied significantly for each composition used for the FGM. Numerical modeling revealed that the change in fluid flow caused by Marangoni convection due to dissolved oxygen in the fusion zone were mainly responsible for the changes in deposit shape observed in experiments. Finally, because oxygen can be introduced into the fusion zone through the feedstock as well as the surrounding atmosphere, these findings elucidate a previously unconsidered aspect of process control during DED fabrication of FGMs.

36 MATERIALS SCIENCE↗

Coupling of CTF and RELAP5-3D Within an Enhanced Fidelity Nuclear Power Plant Simulator

A robust and accurate multiphysics engineering simulator is being developed to model the core behavior and system response of pressurized water reactors. This simulator relies on the NESTLE and CTF computer codes to model the neutronics and thermal hydraulics (TH), respectively, inside the core on a nodal scale and on the Reactor Excursion and Leak Analysis Program—Three Dimensional (RELAP5-3D) to model the entire nuclear steam supply system. The RELAP5-3D model includes highly detailed nodalization and multidimensional flow modeling throughout the vessel. Previously, pin-resolved data generated via the Virtual Environment for Reactor Analysis core simulator were used to improve the accuracy of the NESTLE core predictions. The engineering simulator being developed as part of this work uses the 3KEYMASTER platform to couple the enhanced NESTLE model to a nodal-fidelity CTF model to balance run time with accuracy; NESTLE provides node-dependent powers to CTF, and CTF provides node-dependent coolant densities and fuel temperatures to NESTLE.An overlapping domain approach is used for the core TH in which RELAP5-3D provides core boundary conditions based on the system response and CTF provides a node-dependent coolant heating rate to the RELAP5-3D core solution. In the preliminary TH demonstration discussed in this paper, CTF and RELAP5-3D provided similar steady-state core predictions, indicating the hydraulic compatibility between the codes, as well as reasonable and expected behavior under hypothetical transient conditions. This provides an initial step in ongoing efforts toward a robust, multiscale TH/neutronics engineering simulator capability.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Data-Driven Multi-agent Deep Reinforcement Learning for Distribution System Decentralized Voltage Control with High Penetration of PVs

This paper proposes a novel model-free/data-driven centralized training and decentralized execution multi-agent deep reinforcement learning (MADRL) framework for distribution system voltage control with high penetration of PVs. The proposed MADRL can coordinate both the real and reactive power control of PVs with existing static var compensators and battery storage systems. Unlike the existing DRL-based voltage control methods, our proposed method does not rely on a system model during both the training and execution stages. This is achieved by developing a new interaction scheme between the surrogate modeling of the original system and the multi-agent soft actor critic (MASAC) MADRL algorithm. In particular, the sparse pseudo-Gaussian process with a few-shots of measurements is utilized to construct the surrogate model of the original environment, i.e., power flow model. This is a data-driven process and no model parameters are needed. Furthermore, the MASAC enabled MADRL allows to achieve better scalability by dividing the original system into different voltage control regions with the aid of real and reactive power sensitivities to voltage, where each region is treated as an agent. This also serves as the foundation for the centralized training and decentralized execution, thus significantly reducing the communication requirements as only local measurements are required for control. Comparative results with other alternatives on the IEEE 123-nodes and 342-nodes systems demonstrate the superiority of the proposed method.

14 SOLAR ENERGY↗

ParFlow

ParFlow is an open-source, modular, parallel watershed flow model. It includes fully-integrated overland flow, the ability to simulate complex topography, geology and heterogeneity and coupled land-surface processes including the land-energy budget, biogeochemistry and snow (via CLM). It is multi-platform and runs with a common I/O structure from laptop to supercomputer. ParFlow is the result of a long, multi-institutional development history and is now a collaborative effort between CSM, LLNL, UniBonn and UCB. ParFlow has been coupled to the mesoscale, meteorological code ARPS and the NCAR code WRF.

Smith, Steven↗

ParFlow

ParFlow is an open-source, modular, parallel watershed flow model. It includes fully-integrated overland flow, the ability to simulate complex topography, geology and heterogeneity and coupled land-surface processes including the land-energy budget, biogeochemistry and snow (via CLM). It is multi-platform and runs with a common I/O structure from laptop to supercomputer. ParFlow is the result of a long, multi-institutional development history and is now a collaborative effort between CSM, LLNL, UniBonn and UCB. ParFlow has been coupled to the mesoscale, meteorological code ARPS and the NCAR code WRF.

Smith, Steven↗