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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Prediction of the ullage gas thermal stratification in a NASP vehicle propellant tank experimental simulation using FLOW-3D

As part of the National Aero-Space Plane (NASP) project, the multi-dimensional effects of gravitational force, initial tank pressure, initial ullage temperature, and heat transfer rate on the 2-D temperature profiles were studied. FLOW-3D, a commercial finite difference fluid flow model, was used for the evaluation. These effects were examined on the basis of previous liquid hydrogen experimental data with gaseous hydrogen pressurant. FLOW-3D results were compared against an existing 1-D model. In addition, the effects of mesh size and convergence criteria on the analytical results were investigated. Suggestions for future modifications and uses of FLOW-3D for modeling of a NASP tank are also presented.

Hardy, Terry L.↗
NASA NTRS2005

Upgrades of Two Computer Codes for Analysis of Turbomachinery

Major upgrades have been made in two of the programs reported in "ive Computer Codes for Analysis of Turbomachinery". The affected programs are: Swift -- a code for three-dimensional (3D) multiblock analysis; and TCGRID, which generates a 3D grid used with Swift. Originally utilizing only a central-differencing scheme for numerical solution, Swift was augmented by addition of two upwind schemes that give greater accuracy but take more computing time. Other improvements in Swift include addition of a shear-stress-transport turbulence model for better prediction of adverse pressure gradients, addition of an H-grid capability for flexibility in modeling flows in pumps and ducts, and modification to enable simultaneous modeling of hub and tip clearances. Improvements in TCGRID include modifications to enable generation of grids for more complicated flow paths and addition of an option to generate grids compatible with the ADPAC code used at NASA and in industry. For both codes, new test cases were developed and documentation was updated. Both codes were converted to Fortran 90, with dynamic memory allocation. Both codes were also modified for ease of use in both UNIX and Windows operating systems.

Chima, Rodrick V.↗
DOE OSTI2024

Solvent-Dependent Dynamics of Cellulose Nanocrystals in Process-Relevant Flow Fields

Flow-assisted alignment of anisotropic nanoparticles is a promising route for the bottom-up assembly of advanced materials with tunable properties. While aligning processes could be optimized by controlling factors such as solvent viscosity, flow deformation, and the structure of the particles themselves, it is necessary to understand the relationship between these factors and their effect on the final orientation. In this study, we investigated the flow of surface-charged cellulose nanocrystals (CNCs) with the shape of a rigid rod dispersed in water and propylene glycol (PG) in an isotropic tactoid state. In situ scanning small-angle X-ray scattering (SAXS) and rheo-optical flow-stop experiments were used to quantify the dynamics, orientation, and structure of the assigned system at the nanometer scale. The effects of both shear and extensional flow fields were revealed in a single experiment by using a flow-focusing channel geometry, which was used as a model flow for nanomaterial assembly. Due to the higher solvent viscosity, CNCs in PG showed much slower Brownian dynamics than CNCs in water and thus could be aligned at lower deformation rates. Moreover, CNCs in PG also formed a characteristic tactoid structure but with less ordering than CNCs in water owing to weaker electrostatic interactions. The results indicate that CNCs in water stay assembled in the mesoscale structure at moderate deformation rates but are broken up at higher flow rates, enhancing rotary diffusion and leading to lower overall alignment. Albeit being a study of cellulose nanoparticles, the fundamental interplay between imposed flow fields, Brownian motion, and electrostatic interactions likely apply to many other anisotropic colloidal systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantifying the Impacts of Land-Cover Change on the Hydrologic Response to Hurricane Ida in the Lower Mississippi River Basin

Abstract The Lower Mississippi River basin (LMRB) has experienced significant changes in land cover and is one of the most vulnerable regions to hurricanes in the United States. Here, we study the impacts of land-cover change on the hydrologic response to Hurricane Ida in LMRB. By using an integrated surface–subsurface hydrologic model, Energy Exascale Earth System Model (E3SM) Land Model coupled with the three-dimensional ParFlow subsurface flow model (ELM-ParFlow), we simulate the effects of land-cover change on the flood volume and peak timing induced by rainfall from Hurricane Ida. The results show that land-cover changes from 1850 to 2015, which resulted in a smoother surface and less vegetation, exacerbated both flood peak time and volume induced by Hurricane Ida. The effects of land-cover changes can be decomposed into two mechanisms: a smoother surface routes more water faster to a watershed outlet and less vegetation allows more water to contribute to surface runoff. By comparing scenarios in which the two mechanisms were isolated, we found that changes in soil moisture due to vegetation cover change have more dominant effects on floods in the southern part and changes in Manning’s coefficient have the largest effect on floods in the northern part of the LMRB. The study provides important insights into the complex relationship between land-use, land-cover, and hydrologic processes in coastal regions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1978

Magma production and migration within the moon

Partial melting is likely to have occurred throughout much of the moon due to heating during accretion and the volumetric heating of radioactive isotopes. Important problems that have received relatively little attention concern the migration of the resulting magmas to form surface or near surface volcanic rock. In the paper the basic mechanism for the migration of the magma through the lunar asthenosphere is considered. A porous flow model is proposed. The magma behaves like a liquid flowing through a porous matrix. The volume fraction of liquid present determines the saturated porosity. The differential buoyancy of the magma drives it upwards. It is shown that the per cent partial melt in the lunar interior will only slightly exceed that required to provide interconnecting porosity. Assuming that the radioactive isotopes are preferentially segregated into the magma, the time dependence of the partial melting of the lunar interior is found. It is shown that the total degree of partial melting of the deep lunar interior is likely to be between five and ten per cent.

Turcotte, D. L.↗
NASA NTRS1989

Time-dependent Navier-Stokes computations for flow-induced vibrations of vanes

Flows over two curved vane configurations were computed using a time-accurate compressible Navier-Stokes flow model. One configuration showed the presence of strong flow-induced vibrations at Strouhal numbers near 0.19 and 0.38 for bending and torsional excitation. In the other configuration, a simple modification reduced both types of response. Laminar flows were analyzed for the effects of flow-induced vibrations, and flow fields were solved for a rigid vane and a vane undergoing forced vibrations at prescribed amplitude and frequency simulating vibration response to a coupled vortex-shedding/elastic motion feedback cycle.

Liu, B. L.↗
NASA NTRS1995

Flow properties along field lines in a 3-D tilted-dipole geometry

A 3-D MHD simulation of a global, tilted-dipole solar wind flow pattern is analyzed to determine flow properties along individual magnetic field lines. In the model, flow conditions near the Sun are chosen to provide a reasonable match to the interplanetary configuration prevailing during the recent south polar passage by Ulysses, i.e., a streamer belt inclined approximately 30 deg to the solar equator and speeds ranging from 325-800 km/s. Field lines all across the stream pattern are traced from 1 to 10 AU by following the motion of marker particles embedded in the flow. It is found that those field lines threading the core of the interaction region are subject to significant latitudinal and relative longitudinal displacement over this range of heliocentric distance. Thus, observations taken at a fixed latitude in the inner solar system sample, over the course of a solar rotation, field lines which connect to a range of latitudes in the outer heliosphere. Maps of the field line displacements are presented to help visualize these connections. In addition, it is found that depending upon the location relative to the CIR structure, the radial evolution of fluid elements frozen to different field lines can deviate considerably from that of the canonical solar wind. That is, for selected subsets of field lines, large speed changes (not just at shocks) can be experienced; the density variation can be far from 1/r(exp 2), and the magnetic field intensity need not decay monotonically with distance.

Pizzo, V. J.↗
NASA NTRS1997

Investigations of Arcjet Flow Thermochemistry

An experimental and numerical investigation of the thermochemical state of arcjet flows is currently being conducted at NASA Ames Research Center. The experimental approach relies on the use of laser- and emission-spectroscopic diagnostic techniques in three regions of the flow. A fiber optic sensor is used to record spectrally resolved emission signals from the electrode package region, where the flow is most likely to be in thermochemical equilibrium. A second emission diagnostic measurement is made in the shock layer formed over a blunt-body test article placed in the stream, and a CCD camera is used to simultaneously record spectral emission from several measurement locations along the stagnation streamline. Downstream of the nozzle exit, but upstream of the test article, Laser-Induced Fluorescence (LIF) of atomic nitrogen is used to assess the nonequilibrium distribution of flow enthalpy in the free stream. Results from the measurements are compared with predictions from a two-temperature, axisymmetric flow model that solves the nozzle and shock-layer flows.

Fletcher, Douglas G.↗
NASA NTRS2017

Models of Sector Flows Under Local, Regional and Airport Weather Constraints

Recently, the ATM community has made important progress in collaborative trajectory management through the introduction of a new FAA traffic management initiative called a Collaborative Trajectory Options Program (CTOP). FAA can use CTOPs to manage air traffic under multiple constraints (manifested as flow constrained areas or FCAs) in the system, and it allows flight operators to indicate their preferences for routing and delay options. CTOPs also permits better management of the overall trajectory of flights by considering both routing and departure delay options simultaneously. However, adoption of CTOPs in airspace has been hampered by many factors that include challenges in how to identify constrained areas and how to set rates for the FCAs. Decision support tools providing assistance would be particularly helpful in effective use of CTOPs. Such DSTs tools would need models of demand and capacity in the presence of multiple constraints. This study examines different approaches to using historical data to create and validate models of maximum flows in sectors and other airspace regions in the presence of multiple constraints. A challenge in creating an empirical model of flows under multiple constraints is a lack of sufficient historical data that captures diverse situations involving combinations of multiple constraints especially those with severe weather. The approach taken here to deal with this is two-fold. First, we create a generalized sector model encompassing multiple sectors rather than individual sectors in order to increase the amount of data used for creating the model by an order of magnitude. Secondly, we decompose the problem so that the amount of data needed is reduced. This involves creating a baseline demand model plus a separate weather constrained flow reduction model and then composing these into a single integrated model. A nominal demand model is a flow model (gdem) in the presence of clear local weather. This defines the flow as a function of weather constraints in neighboring regions, airport constraints and weather in locations that can cause re-routes to the location of interest. A weather constrained flow reduction model (fwx-red) is a model of reduction in baseline counts as a function of local weather. Because the number of independent variables associated with each of the two decomposed models is smaller than that with a single model, need for amount of data is reduced. Finally, a composite model that combines these two can be represented as fwx-red (gdem(e), l) where e represents non-local constraints and l represents local weather. The approaches studied to developing these models are divided into three categories: (1) Point estimation models (2) Empirical models (3) Theoretical models. Errors in predictions of these different types of models have been estimated. In situations when there is abundant data, point estimation models tend to be very accurate. In contrast, empirical models do better than theoretical models when there is some data available. The biggest benefit of theoretical models is their general applicability in wider range situations once the degree of accuracy of these has been established.

Kulkarni, Deepak↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
NASA NTRS1988

Prediction of incidence effects on oscillating airfoil aerodynamics by a locally analytical method

A complete mathematical model is formulated to analyze the effects of mean-flow incidence angle on the unsteady aerodynamics of an oscillating airfoil in an incompressible flow field. A velocity potential formulation is utilized. The steady flow is independent of the unsteady flow field but coupled to it through the boundary conditions on the oscillating airfoil. The numerical solution technique for both the steady and unsteady flow fields is based on a locally analytical method. The flow model and solution method are then verified through the excellent correlation obtained with the Theodorsen oscillating-flat-plate and Sears transverse-gust classical solutions. The effects of mean flow incidence on the steady and oscillating airfoil aerodynamics are then investigated.

Chiang, Hsiao-Wei D.↗
NASA NTRS1984

Finite-difference fluid dynamics computer mathematical models for the design and interpretation of experiments for space flight

Numerical methods are used to design a spherical baroclinic flow model experiment of the large scale atmosphere flow for Spacelab. The dielectric simulation of radial gravity is only dominant in a low gravity environment. Computer codes are developed to study the processes at work in crystal growing systems which are also candidates for space flight. Crystalline materials rarely achieve their potential properties because of imperfections and component concentration variations. Thermosolutal convection in the liquid melt can be the cause of these imperfections. Such convection is suppressed in a low gravity environment. Two and three dimensional finite difference codes are being used for this work. Nonuniform meshes and implicit iterative methods are used. The iterative method for steady solutions is based on time stepping but has the options of different time steps for velocity and temperature and of a time step varying smoothly with position according to specified powers of the mesh spacings. This allows for more rapid convergence. The code being developed for the crystal growth studies allows for growth of the crystal as the solid-liquid interface. The moving interface is followed using finite differences; shape variations are permitted. For convenience in applying finite differences in the solid and liquid, a time dependent coordinate transformation is used to make this interface a coordinate surface.

Roberts, G. O.↗
DOE OSTI2021

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗
NASA NTRS1987

Analysis of eccentric annular incompressible seals. II - Effects of eccentricity on rotordynamic coefficients

A new analysis procedure has been presented which solves for the flow variables of an annular pressure seal in which the rotor has a large static displacement (eccentricity) from the centered position. The present paper incorporates the solutions to investigate the effect of eccentricity on the rotordynamic coefficients. The analysis begins with a set of governing equations based on a turbulent bulk-flow model and Moody's friction factor equation. Perturbations of the flow variables yields a set of zeroth- and first-order equations. After integration of the zeroth-order equations, the resulting zeroth-order flow variables are used as input in the solution of the first-order equations. Further integration of the first order pressures yields the eccentric rotordynamic coefficients. The results from this procedure compare well with available experimental and theoretical data, with accuracy just as good or slightly better than the predictions based on a finite-element model.

Nelson, C. C.↗
NASA NTRS1994

Turbulence modeling for separated flow

Two projects are described in this report. The first involves assessing turbulence models in separated flow. The second addresses the anomalous behavior of certain turbulence models in stagnation point flow. The primary motivation for developing turbulent transport models is to provide tools for computing non-equilibrium, or complex, turbulent flows. Simple flows can be analyzed using data correlations or algebraic eddy viscosities, but in more complicated flows such as a massively separated boundary layer, a more elaborate level of modeling is required. It is widely believed that at least a two-equation transport model is required in such cases. The transport equations determine the evolution of suitable velocity and time-scales of the turbulence. The present study included assessment of second-moment closures in several separated flows, including sharp edge separation; smooth wall, pressure driven separation; and unsteady vortex shedding. Flows with mean swirl are of interest for their role in enhancing mixing both by turbulent and mean motion. The swirl can have a stabilizing effect on the turbulence. An axi-symmetric extension to the INS-2D computer program was written adding the capability of computing swirling flow. High swirl can produce vortex breakdown on the centerline of the jet and it occurs in various combustors.

Durbin, Paul A.↗
NASA NTRS1986

Modeling turbulent, reacting flow

Several of the approximations or models involved in the development of a numerical combustor flow code are examined. In the first section, the importance of numerical accuracy is illustrated, and the impact that improved-accuracy schemes have on slowing convergence is demonstrated. Solution algorithms that can speed convergence are discussed and some performance features of these algorithms are illustrated. A sample calculation displaying the importance of boundary conditions on a three-dimensional numerical prediction is presented. The inaccuracy of a current turbulence model in highly turbulent (nonequilibrium) regions is described. Finally, the surprisingly good performance of a six-flux model in describing radiation heat transfer is displayed. In all the areas examined, continued research is still needed, but valuable engineering tools are available today.

Claus, Russell W.↗
NASA NTRS1980

Loss model for off-design performance analysis of radial turbines with pivoting-vane, variable-area stators

An off-design performance loss model for a radial turbine with pivoting, variable-area stators is developed through a combination of analytical modeling and experimental data analysis. A viscous loss model is used for the variation in stator loss with setting angle, and stator vane end-clearance leakage effects are predicted by a clearance flow model. The variation of rotor loss coefficient with stator setting angle is obtained by means of an analytical matching of experimental data for a rotor that was tested with six stators, having throat areas from 20 to 144% of the design area. An incidence loss model is selected to obtain best agreement with experimental data. The stator vane end-clearance leakage model predicts increasing mass flow and decreasing efficiency as a result of end-clearances, with changes becoming significantly larger with decreasing stator area.

Meitner, P. L.↗
NASA NTRS2002

Aeroelastic Analysis for Aeropropulsion Applications

Aeroelastic codes with advanced capabilities for modeling flow require substantial computational time. On the other hand, fast-running linear aeroelastic codes lack the capability to model three-dimensional, transonic, vortical, and viscous flows. The goal of this work was to develop an aeroelastic code with accurate modeling capabilities and small computational requirements.

Keith, Theo G., Jr.↗