SEARCH · Engineering Papers

Results for “flow simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Numerical simulation of transonic flows in diffusers

Numerical simulations were made of two-dimensional transonic flows in diffusers, including flow separation induced by a shock or adverse pressure gradient. The mass-averaged, time-dependent, compressible Navier-Stokes equations, simplified by the thin-layer approximation, were solved using MacCormack's hybrid method. The eddy-viscosity formulation was described by the Wilcox-Rubesin's two-equation, k-omega model. Detailed comparison of the computed results with measurements showed good agreement in all cases, including one with massive separation induced by a strong shock. The computation correctly predicted the details of a distinct lambda shock pattern, closely duplicating the configuration observed experimentally in spark-schlieren photographs.

Liou, M.-S.↗

NASA NTRS1983

Numerical simulation of turbulent trailing edge flows

Numerical simulations of the time-dependent, Reynolds-averaged, Navier-Stokes equations, employing a two-equation turbulence model, are presented and compared with measurements from a series of trailing edge experiments at transonic Mach numbers. The test flows include an asymmetric flow with no separation, an asymmetric flow with a small region of separation and a symmetric flow with a large shock-wave induced separated zone. Comparisons are made for mean surface quantities as well as for mean and fluctuating flow-field quantities. For the trailing-edge flows with little or no separation, the solutions correctly predict all the major features of the flow field. Treatment of the viscous-inviscid interaction was found to be important for predicting these test cases. Two-equation eddy-viscosity turbulence models were found to be adequate for these flows. However, for the shock-wave induced separation case, these turbulence models were inadequate to predict this flow field. Modifications of the turbulence model to correct these deficiencies are discussed.

Horstman, C. C.↗

DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗

DOE OSTI2021

Investigation of grid-based vorticity-velocity large eddy simulation off-body solvers for application to overset CFD

Accurately predicting unsteady wakes and vortex-dominated flows is essential to a wide range of engineering applications, including aircraft, rotorcraft, shipboard operations, bio-inspired unsteady flight and propulsion, wind turbines, and urban flows. While current CFD software can model the complete flow field and wake system, the computational costs incurred in high Reynolds number unsteady turbulent flow simulations often remain prohibitive for routine engineering use, particularly for applications involving moving components. Prior work has demonstrated that by adopting a vorticity-velocity formulation in a grid-based off-body flow solver (VorTran-M and VorTran-M2) one can lower these costs by several orders of magnitude when compared to conventional approaches. This paper describes the extensions made to VorTran-M2 to support turbulent flows, and associated benchmarking activity to assess its performance for problems involving strong stretching and diffusion processes, whose competing contributions to the vorticity field are core drivers of turbulent flow evolution. Predictions are presented for: (i) the Kida-Pelz problem whose inviscid form is of mathematical interest due to its apparent formation of singular flow in finite time; and (ii) the Taylor Green vortex arrangement, which has been extensively studied as a fundamental simulation challenge in the turbulent modeling community. Here, the results are used to evaluate the overall predictive ability and performance of two sub-grid scale models incorporated into VorTran-M2. Results indicate that the computational cost savings seen previously for inviscid and convection dominated problems extend to turbulent flow simulations supporting the viability of VorTran-M2 as a low cost means for accurately modeling the far-field and background flow, particularly when long duration vorticity evolution is of interest.

42 ENGINEERING↗

NASA NTRS1995

Advanced Methodology for Simulation of Complex Flows Using Structured Grid Systems

Detailed simulations of viscous flows in complicated geometries pose a significant challenge to current capabilities of Computational Fluid Dynamics (CFD). To enable routine application of CFD to this class of problems, advanced methodologies are required that employ (a) automated grid generation, (b) adaptivity, (c) accurate discretizations and efficient solvers, and (d) advanced software techniques. Each of these ingredients contributes to increased accuracy, efficiency (in terms of human effort and computer time), and/or reliability of CFD software. In the long run, methodologies employing structured grid systems will remain a viable choice for routine simulation of flows in complex geometries only if genuinely automatic grid generation techniques for structured grids can be developed and if adaptivity is employed more routinely. More research in both these areas is urgently needed.

Steinthorsson, Erlendur↗

NASA NTRS1994

A workstation based simulator for teaching compressible aerodynamics

A workstation-based interactive flow simulator has been developed to aid in the teaching of undergraduate compressible aerodynamics. By solving the equations found in NACA 1135, the simulator models three basic fluids problems encountered in supersonic flow: flow past a compression corner, flow past two wedges in series, and flow past two opposed wedges. The study can vary the geometry or flow conditions through a graphical user interface and the new conditions are calculated immediately. Various graphical formats present the results of the flow calculations to the student. The simulator includes interactive questions and answers to aid in both the use of the tool and to develop an understanding of some of the complexities of compressible aerodynamics. A series of help screens make the simulator easy to learn and use.

Benson, Thomas J.↗

DOE OSTI2023

Latest developments in the MOOSE fluid properties module

The fluid properties module in MOOSE serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, SAM~\cite{sam}, the MOOSE thermal hydraulics module, RELAP-7~\cite{relap7} and subchannel~\cite{subchannel}. It is used for coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled simulations~\cite{anl_sam_pgh}. The module offers a consistent set of interfaces to implement to create a new fluid property. There are numerous fluid properties of interest in the entirety of all fields of fluid flow simulations, and this is exacerbated by the use of different variable sets depending on the compressibility of the fluid. For single-phase fluids, the following variable sets may be used to compute fluid properties: (pressure, temperature) and (specific volume, specific internal energy). Some properties may also be computed using the (pressure, density) or the (specific volume, specific enthalpy) variable sets. In order to reduce the challenge of adding a new fluid property, properties may be implemented partially.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Rapid Simulation of Unconventional Reservoirs by Multidomain Multiresolution Modeling Based on the Diffusive Time of Flight

Summary The fast marching method (FMM)-based rapid flow simulation has been shown to accelerate simulation efficiency by orders of magnitude by transforming 3D simulation to equivalent 1D simulation using the concept of the “diffusive time of flight” (DTOF). However, the 1D transformation does not directly apply to multiwell problems. In this paper, we propose a novel DTOF-based multidomain multiresolution discretization scheme to accelerate multiwell simulation of unconventional reservoirs. Our method formulates multiwell simulation problems based on the DTOF which displays the pressure front propagation in unconventional reservoirs. The DTOF contours are used to partition the reservoir into local and shared domains. A local domain is where the flow is dominated by a single well, and the shared domain is where the fluid flow is influenced by multiple wells. The DTOF contours expand independently in local domains and interfere in the shared domain. After the partitioning, each domain is discretized using a multiresolution scheme whereby the original 3D fine mesh is preserved near the wells to account for detailed physics including gravity, and the rest of the domain is discretized into 1D mesh based on the DTOF contours to alleviate the simulation workload. The power and efficacy of our approach are demonstrated using synthetic and field-scale simulation models with different degrees of geologic and well-completion complexity. The simulation results, number of active cells, and computation time for the proposed discretization scheme are compared with the original high-fidelity 3D model for each case. The results show that the proposed method is suitable for multiwell simulation problems in unconventional reservoirs and can accelerate flow simulations by orders of magnitude with minimal loss of accuracy. The novelty of this work is the creation of DTOF-derived multiresolution discretization with local and shared domains to simplify and accelerate the calculation of subsurface flow problems, especially in unconventional reservoirs. Our workflow can be easily interfaced with commercial simulators, making it suitable for large-scale field applications.

Engineering↗

NASA NTRSDate not supplied

Turbulent Simulations of Cooling Jets in Crossflow

Turbulent flow simulations were completed of the THX-III single hole film cooling experiment performed at NASA Glenn Research Center (GRC). Two different Large-Eddy Simulation (LES) approaches were completed along with a RANS based study. The first LES approach utilizes the Flux Reconstruction (FR) based finite element code GFR, while the second approach makes use of the compact difference solver FDL3DI. Both LES codes were run on grids larger than 1.5 Billion degrees of freedom (DoF). The RANS simulations used FUN3D and included a grid convergence as a turbulence closure study. A summary of the computational approaches and key findings are presented along with the experimental measurements. Significant progress has been made towards the completion of this work, as all but one simulation have been completed. The preliminary findings have shown the LES approaches have excellent agreement with each other, and closely follow the trends of the experiment. The Reynolds stress RANS closure, while the best of the RANS results, had difficulty matching the results of the LES and experimental measurements.

High-Order Methods↗

NASA NTRSDate not supplied

Transition Prediction for the 30P30N Airfoil and the CRM-HL using FUN3D

The ability to accurately model transition over high-lift configurations can lead to improved agreement between stationary Reynolds-Averaged Navier-Stokes (RANS) computations and experimental measurements. In this work, selected transport-equation-based RANS models for transition prediction in the NASA FUN3D flow solver are used to simulate flow over the 30P30N three-element airfoil and the High-Lift Common Research Model (CRM-HL). The primary RANS-based transition model is the Langtry-Menter (LM) 𝜸-Re𝜽𝒕 model coupled with the Menter shear-stress transport model. Additionally, the existing capability for transition prediction based on linear stability correlations is used to allow the modeling of transition over multiple elements of the 30P30N configuration. This approach pairs an algebraic transition model with the Spalart-Allmaras (SA) turbulence model that relies on an imposed transition location in FUN3D and the parabolized stability equations as implemented in the LASTRAC code. Predictions based on these different transition prediction approaches and models are compared with each other, experimental data, and previous simulations by evaluating the lift and drag coefficients, the surface-pressure/skin-friction distributions, and the locations of transition onset. Results are presented for multiple angles of attack for the 30P30N airfoil and CRM-HL. For the CRM-HL, a grid resolution study is performed for the LM transition model and the SA turbulence model.

High Lift↗

DOE OSTI2020

Mixed-flow design for microfluidic printing of two-component polymer semiconductor systems

The rational creation of two-component conjugated polymer systems with high levels of phase purity in each component is challenging but crucial for realizing printed soft-matter electronics. Here, we report a mixed-flow microfluidic printing (MFMP) approach for two-component π -polymer systems that significantly elevates phase purity in bulk-heterojunction solar cells and thin-film transistors. MFMP integrates laminar and extensional flows using a specially microstructured shear blade, designed with fluid flow simulation tools to tune the flow patterns and induce shear, stretch, and pushout effects. This optimizes polymer conformation and semiconducting blend order as assessed by atomic force microscopy (AFM), transmission electron microscopy (TEM), grazing incidence wide-angle X-ray scattering (GIWAXS), resonant soft X-ray scattering (R-SoXS), photovoltaic response, and field effect mobility. For printed all-polymer (poly[(5,6-difluoro-2-octyl-2H-benzotriazole-4,7-diyl)-2,5-thiophenediyl[4,8-bis[5-(2-hexyldecyl)-2-thienyl]benzo[1,2-b:4,5-b′]dithiophene-2,6-diyl]-2,5-thiophenediyl]) [J51]:(poly{[N,N′-bis(2-octyldodecyl)naphthalene-1,4,5,8-bis(dicarboximide)-2,6-diyl]-alt-5,5′-(2,2′-bithiophene)}) [N2200]) solar cells, this approach enhances short-circuit currents and fill factors, with power conversion efficiency increasing from 5.20% for conventional blade coating to 7.80% for MFMP. Moreover, the performance of mixed polymer ambipolar [poly(3-hexylthiophene-2,5-diyl) (P3HT):N2200] and semiconducting:insulating polymer unipolar (N2200:polystyrene) transistors is similarly enhanced, underscoring versatility for two-component π -polymer systems. Mixed-flow designs offer modalities for achieving high-performance organic optoelectronics via innovative printing methodologies.

42 ENGINEERING↗

NASA NTRS1980

Modeling of vortex generated sound in solid propellant rocket motors

There is considerable evidence based on both full scale firings and cold flow simulations that hydrodynamically unstable shear flows in solid propellant rocket motors can lead to acoustic pressure fluctuations of significant amplitude. Although a comprehensive theoretical understanding of this problem does not yet exist, procedures were explored for generating useful analytical models describing the vortex shedding phenomenon and the mechanisms of coupling to the acoustic field in a rocket combustion chamber. Since combustion stability prediction procedures cannot be successful without incorporation of all acoustic gains and losses, it is clear that a vortex driving model comparable in quality to the analytical models currently employed to represent linear combustion instability must be formulated.

Flandro, G. A.↗

DOE OSTI2015

Velocity field calculation for non-orthogonal numerical grids

Computational grids containing cell faces that do not align with an orthogonal (e.g. Cartesian, cylindrical) coordinate system are routinely encountered in porous-medium numerical simulations. Such grids are referred to in this study as non-orthogonal grids because some cell faces are not orthogonal to a coordinate system plane (e.g. xy, yz or xz plane in Cartesian coordinates). Non-orthogonal grids are routinely encountered at the Savannah River Site in porous-medium flow simulations for Performance Assessments and groundwater flow modeling. Examples include grid lines that conform to the sloping roof of a waste tank or disposal unit in a 2D Performance Assessment simulation, and grid surfaces that conform to undulating stratigraphic surfaces in a 3D groundwater flow model. Particle tracking is routinely performed after a porous-medium numerical flow simulation to better understand the dynamics of the flow field and/or as an approximate indication of the trajectory and timing of advective solute transport. Particle tracks are computed by integrating the velocity field from cell to cell starting from designated seed (starting) positions. An accurate velocity field is required to attain accurate particle tracks. However, many numerical simulation codes report only the volumetric flowrate (e.g. PORFLOW) and/or flux (flowrate divided by area) crossing cell faces. For an orthogonal grid, the normal flux at a cell face is a component of the Darcy velocity vector in the coordinate system, and the pore velocity for particle tracking is attained by dividing by water content. For a non-orthogonal grid, the flux normal to a cell face that lies outside a coordinate plane is not a true component of velocity with respect to the coordinate system. Nonetheless, normal fluxes are often taken as Darcy velocity components, either naively or with accepted approximation. To enable accurate particle tracking or otherwise present an accurate depiction of the velocity field for a non-orthogonal grid, Darcy velocity components are rigorously derived in this study from normal fluxes to cell faces, which are assumed to be provided by or readily computed from porous-medium simulation code output. The normal fluxes are presumed to satisfy mass balances for every computational cell, and if so, the derived velocity fields are consistent with these mass balances. Derivations are provided for general two-dimensional quadrilateral and three-dimensional hexagonal systems, and for the commonly encountered special cases of perfectly vertical side faces in 2D and 3D and a rectangular footprint in 3D.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1981

Higher level simulations of turbulent flows

The fundamentals of large eddy simulation are considered and various approaches to this simulation are compared. The subgrid model required by large eddy simulation is discussed as well as the use of this type of simulation in the development of models for the Reynolds-averaged equations and the application of direct simulation to the testing of both subgrid scale and Reynolds-averaged models. Numerical methods used in large eddy and direct simulation are described with emphasis on special purpose methods. Topics covered include the simulation of homogeneous flow, free shear flows, the mixing layer, wakes, and wall-bounded flows including channel flow and the boundary layer. Applications of large eddy simulation in the laboratory as in meteorological and other environmental flows are examined. Directions in which the work is proceeding and what can be expected from higher levels simulated on are examined.

Ferziger, J. H.↗

DOE OSTI2016

Intercomparison of 3D pore-scale flow and solute transport simulation methods

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include methods that 1) explicitly model the three-dimensional geometry of pore spaces and 2) those that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of class 1, based on direct numerical simulation using computational fluid dynamics (CFD) codes, against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of class 1 based on the immersed-boundary method (IMB), lattice Boltzmann method (LBM), smoothed particle hydrodynamics (SPH), as well as a model of class 2 (a pore-network model or PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and nonreactive solute transport, and intercompare the model results with previously reported experimental observations. Experimental observations are limited to measured pore-scale velocities, so solute transport comparisons are made only among the various models. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations).

computational fluid dynamics↗

NASA NTRS2015

Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the isotropic turbulent flow decay, at a relatively high turbulent Mach number, show a nicely behaved spectral decay rate for medium to high wave numbers. The high-order CESE schemes offer very robust solutions even with the presence of strong shocks or widespread shocklets. The explicit formulation in conjunction with a close to unity theoretical upper Courant number bound has the potential to offer an efficient numerical framework for general compressible turbulent flow simulations with unstructured meshes.

Chang, Chau-Lyan↗

DOE OSTI2020

Numerical simulation of mold filling water blown polyurethane foams: Effects of sequential pour

Batch molding of polyurethane foams is a widely used processing technique to produce crosslinked cellular structures which form to the shape of the mold. For water-blown polyurethane foams, water is reacted to form gaseous carbon dioxide resulting in a foam which expands to fill the mold completely. Batch molding typically requires an operator to coat the surface of a mold, introducing a significant lag time during the filling stage where significant asymmetric volume change can occur. The purpose of this work is to show, through simulations, that this lag time is significant when predicting flow profiles and part quality. When the mold geometry is complex enough to force bifurcation of the flow, simulations incorporating various filling lag times predicted significantly different locations of knit or weld lines where the flow fronts meet. Current simulation techniques, which assume the filling stage occurs instantaneously, are unable to predict variations in weld line locations. The introduction of a lag time, referred to herein as sequential pour, was achieved through user-defined modules incorporated into Ansys Fluent software.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Development and validation of the cavitation-induced erosion risk assessment tool

This work presents the development of a cavitation-induced erosion risk assessment (CIERA) tool that links multiphase flow simulation predictions with the progress towards material erosion. To develop a robust erosion modeling tool, the cavitation and erosion predictions for pressurized diesel fuel flow within a channel geometry were validated over a range of Reynolds and cavitation number conditions in two different aluminum channel geometries, one featuring a rounded inlet corner and the other with a sharp inlet corner. The multiphase flow development within the channel was modeled using a compressible mixture model, where phase change was represented with the homogeneous relaxation model and the turbulent flow evolution was modeled using a dynamic structure approach for Large Eddy Simulations. To improve representation of the incubation period before material rupture over existing approaches, a physics-based metric was derived based on the cumulative energy absorbed by the solid material from repeated hydrodynamic impacts. When the average peak pressure was related to the incubation period, the incubation period and its sensitivity to changes in flow conditions were found to be overpredicted. In contrast, predictions from CIERA provided a more accurate means to qualitatively and quantitatively predict the influence of flow conditions on the incubation period before material erosion. When the predicted stored energy was related to the solid material properties to estimate the incubation period, multiphase flow simulations demonstrated accurate representation of the sensitivity of erosion severity to changes in flow conditions. The use of CIERA led to quantitative agreement of the predicted incubation period within 5% of the experimentally measured incubation period.

42 ENGINEERING↗