Engineering Papers⌕ Search

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.

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

At least 127 records · Page 7

Methods and Tools to allow molecular flow simulations to be coupled to higher level continuum descriptions of flows in porous/fractured media and aerosol/dust dynamics

The purpose of this project was to develop methods and tools that will aid in safety evaluation of nuclear fuels and licensing of nuclear reactors relating to accidents.The objectives were to develop more detailed and faster computations of fission product transport and aerosol evolution as they generally relate to nuclear fuel and/or nuclear reactor accidents. The two tasks in the project related to molecular transport in nuclear fuel and aerosol transport in reactor vessel and containment. For both the tasks, explorations of coupling of Direct Simulation Monte Carlo with Navier-Stokes solvers or the Sectional method were not successful. However, Mesh free methods for the Direct Simulation Monte Carlo method were successfully explored.These explorations permit applications to porous and fractured media, and arbitrary geometries.The computations were carried out in Mathematica and are fully parallelized. The project has resulted in new computational tools (algorithms and programs) that will improve the fidelity of computations to actual physics, chemistry and transport of fission products in the nuclear fuel and aerosol in reactor primary and secondary containments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assembly flow simulation of a radar

A discrete event simulation model has been developed to predict the assembly flow time of a new radar product. The simulation was the key tool employed to identify flow constraints. The radar, production facility, and equipment complement were designed, arranged, and selected to provide the most manufacturable assembly possible. A goal was to reduce the assembly and testing cycle time from twenty-six weeks. A computer software simulation package (SLAM 2) was utilized as the foundation for simulating the assembly flow time. FORTRAN subroutines were incorporated into the software to deal with unique flow circumstances that were not accommodated by the software. Detailed information relating to the assembly operations was provided by a team selected from the engineering, manufacturing management, inspection, and production assembly staff. The simulation verified that it would be possible to achieve the cycle time goal of six weeks. Equipment and manpower constraints were identified during the simulation process and adjusted as required to achieve the flow with a given monthly production requirement. The simulation is being maintained as a planning tool to be used to identify constraints in the event that monthly output is increased. 'What-if' studies have been conducted to identify the cost of reducing constraints caused by increases in output requirement.

Rutherford, W. C.↗

Validating compositional fluid flow simulations using 4D seismic interpretation and vice versa in the SECARB Early Test—A critical review

We report on strengths and weaknesses of 4D seismic interpretation as a technique for monitoring carbon dioxide that was injected as part of a large scale test associated with commercial enhanced oil recovery (EOR) at Cranfield Field, Mississippi, USA. The goals of the monitoring effort are 1) to make measurements to verify that the CO 2 is contained in the reservoir according to operational designs and model predictions, and 2) that if there are deviations, to provide data which can be used to update the earth models and determine if any mitigation is needed. Our work uses a compositional numerical simulation to model CO 2 flow in the reservoir and compares the results with estimates of changes in seismic properties between a pre-injection and a survey after more than 2 million metric tons injected. The complicated physics of the problem in Cranfield field present challenges to seismic interpretations. These findings reflect partial agreement between the results of the numerical simulation and the seismic interpretations. Possible causes of discrepancies among the fluid flow model and multiple interpretations of the same seismic data sets performed with different work flows illuminates the types of uncertainties that should be considered to achieve the goals of monitoring.

4D seismic↗

Model equation for simulating flows in multistage turbomachinery

A steady, three-dimensional average-passage equation system is derived for use in simulating multistage turbomachinery flows. These equations describe a steady, viscous flow that is periodic from blade passage to blade passage. From this system of equations, various reduced forms can be derived for use in simulating the three-dimensional flow field within multistage machinery. It is suggested that a properly scaled form of the average-passage equation system would provide an improved mathematical model for simulating the flow in multistage machines at design and, in particular, at off-design conditions.

Adamczyk, J. J.↗

Model equation for simulating flows in multistage turbomachinery

A steady, three-dimensional average-passage equation system is derived for use in simulating multistage turbomachinery flows. These equations describe a steady, viscous flow that is periodic from blade passage to blade passage. From this system of equations, various reduced forms can be derived for use in simulating the three-dimensional flow field within multistage machinery. It is suggested that a properly scaled form of the averaged-passage equation system would provide an improved mathematical model for simulating the flow in multistage machines at design and, in particular, at off-design conditions.

Adamczyk, J. J.↗

Lip leakage flow simulation for the Gravity Probe B gas spinup

The lip leakage flow for the Gravity Probe B (GP-B) gas spinup system is investigated using a particle simulation code on the Connection Machine (PSiCM). Particle simulation is employed because the flow conditions are in the transition regime between continuum and free molecule where particle methods are of greatest use. The dominant flow is Couette in nature and the simulation is first validated through comparison to theoretical results for Couette flow in the transition regime. A GP-B type geometry is then simulated and results are presented for two conditions, those corresponding to near the inlet and near the outlet of the spinup channel. Comparison to experiment is not made because experimental data is not yet available.

Dagum, Leonardo↗

Model equations for simulating flows in multistage turbomachinery

A steady, three dimensional average-passage equation system was derived. The purpose was to simulate multistage turbomachinery flows. These equations describe a steady, viscous flow that is periodic from blade passage to blade passage. Moreover, these equations have a closure problem that is similar to that of the Reynolds-average Navier-Stokes equations. A scaled form of the average-passage equation system could provide an improved mathematical model for simulating the flow in the design and in the off-design conditions of a multistage machine.

Adamczyk, John J.↗

Flow Simulation of N2B Hybrid Wing Body Configuration

The N2B hybrid wing body aircraft was conceptually designed to meet environmental and performance goals for the N+2 generation transport set by the subsonic fixed wing project. In this study, flow fields around the N2B configuration is simulated using a Reynolds-averaged Navier-Stokes flow solver using unstructured meshes. Boundary conditions at engine fan face and nozzle exhaust planes are provided by response surfaces of the NPSS thermodynamic engine cycle model. The present flow simulations reveal challenging design issues arising from boundary layer ingestion offset inlet and nacelle-airframe interference. The N2B configuration can be a good test bed for application of multidisciplinary design optimization technology.

Kim, Hyoungjin↗

Three-dimensional rarefied flow simulations for the Aeroassist Flight Experiment vehicle

Hypersonic rarefied flow about the Aeroassist Flight Experiment vehicle has been investigated using a three-dimensional direct simulation Monte Carlo method. Calculations are performed for the transitional flows encountered during the vehicle's atmospheric entry for altitudes of 110 and 100 km with an entry velocity of 9.9 km/s. The simulations are performed using a five-species reacting gas model that account for rotational and vibrational internal energies. The solutions indicate that dissociation is important at altitudes of 110 km and below. Results are presented for surface pressures, convective heating, flowfield structure, and aerodynamic coefficient variations with altitude.

Celenligil, M. Cevdet↗

Incompressible and anelastic flow simulations on numerically generated grids

In the numerical simulation of incompressible and anelastic flows, it is necessary to solve an elliptic equation at each time step. When the boundaries of such flows are nonrectangular, it may be advantageous to solve the equations on a new, numerically generated coordinate grid, in which the property of orthogonality has been preserved. Flow equations in general curvilinear coordinates maintaining the conservative form are given for both anelastic models using the momentum equations, and for incompressible models, using the vorticity equation. The general problem of grid-generation in two dimensions is presented, and a quasi-conformal transformation technique is discussed in detail. Some examples of grids generated by this technique are exhibited. Three examples of the flow of a stratified fluid over obstacles are presented, in which the grid-generation permits some new results to be obtained.

Sharman, R. D.↗

Turbulent Vortex-Flow Simulation Over a 65 deg Sharp and Blunt Leading-Edge Delta Wing at Subsonic Speeds

Turbulent thin-layer, Reynolds-Averaged Navier-Stokes solutions, based on a multi-block structured grid, are presented for a 65 deg delta wing having either a sharp leading edge (SLE) or blunt leading edge (BLE) geometry. The primary objective of the study is to assess the prediction capability of the method for simulating the leading-edge flow separation and the ensuing vortex flow characteristics. Computational results are obtained for two angles of attack of approximately 13 and 20 deg, at free-stream Mach number of 0.40 and Reynolds number of 6 million based on the wing mean aerodynamic chord. The effects of two turbulence models of Baldwin-Lomax with Degani-Schiff (BL/DS) and the Spalart-Allmaras (SA) on the numerical results are also discussed. The computations also explore the effects of two numerical flux-splitting schemes, i.e., flux difference splitting (fds) and flux vector splitting (fvs), on the solution development and convergence characteristics. The resulting trends in solution sensitivity to grid resolution for the selected leading-edge geometries, angles of attack, turbulence models and flux splitting schemes are also presented. The validity of the numerical results is evaluated against a unique set of experimental wind-tunnel data that was obtained in the National Transonic Facility at the NASA Langley Research Center.

Ghaffari, Farhad↗

Vortex methods for two- and three-dimensional flow simulations

The point vortex and vortex blob methods for two dimensional flows are presented. Several results are discussed concerning the numerical analysis of the latter scheme, e.g., the preservation of globally conserved quantities and the analysis of the spatial discretization error resulting from the convection of fixed blobs of vorticity. An application to the two dimensional mixing layer is briefly described. The contour dynamics method is also discussed. The simulation of three dimensional flows with vortex methods is discussed. A natural way to represent the vorticity is in the form of closed tubes of filaments of vorticity, although other schemes are examined. Applications to aircraft trailing vortices and to a turbulent spot in a laminar boundary layer are presented. Hybrid schemes that use an Eulerian mesh to solve the Poisson equation for the velocity field are discussed. The goal of these schemes is to avoid the high cost of the Biot-Savart integration if many vortex elements are used while enjoying most of the advantages of pure Lagrangian schemes.

Leonard, A.↗

Insights into Waterflooding in Hydrocarbon-Bearing Nanochannels of Varying Cross Sections from Mesoscopic Multiphase Flow Simulations

Waterflooding is one of the geotechnique used to recover fuel sources from nanoporous geological formations. The scientific understanding of the process that involves the multiphase flow of nanoconfined fluids, however, has lagged, mainly due to the complex nanopore geometries and chemical compositions. To enable benchmarked flow of nanoconfined fluids, the architected geomaterials, such as synthesized mesoporous silica with tunable pore shapes and surface chemical properties, are used for designing and conducting experiments and simulations. This work uses a modified many-body dissipative particle dynamics (mDPD) model with accurately calibrated parameters to perform parametric flow simulations for studying the influences of waterflooding driven power, pore shape, surface roughness, and surface wettability on the multiphase flow in heptane-saturated silica nanochannels. Remarkably, up to 80\% reduction in the effective permeability is found for water-driven heptane flow in a baseline 4.5 nm-wide slit channel, when compared with the Hagen–Poiseuille equation. In the 4.5 nm-wide channels with architected surface roughness, the flow rate is found either higher or lower than the baseline case, depending on the shape and size of cross-sections. High wettability of the solid surface to water is essential for achieving high recovery of heptane, regardless of surface roughness. When the solid surface is less wetting or non-wetting to water, the existence of an optimal waterflooding driven power is found to allow for the highest possible recovery. A detailed analysis on the evolution of the transient water-heptane interface in those nanochannels is presented to elucidate the underlying mechanisms that impact or dictate the multiphase flow behaviors.

58 GEOSCIENCES↗

Determining the impacts of venoarterial extracorporeal membrane oxygenation on cerebral oxygenation using a one-dimensional blood flow simulator

Venoarterial extracorporeal membrane oxygenation (VA-ECMO) is a mechanical system that provides rapid and short-term support for patients with cardiac failure. In many patients, pulmonary function is also impaired, resulting in poorly-oxygenated cardiac outflow competing against well-oxygenated VA-ECMO outflow, a condition known as North-South syndrome. North-South syndrome is a primary concern because of its potential to cause cerebral hypoxia, which has a critical influence on neurological complications often seen in this patient population. In order to reduce ischemic neurological complications, it is important to understand how clinical decisions regarding VA-ECMO parameters influence blood oxygenation. In our work, we studied the impacts of flow rate and cannulation site on oxygenation using a one-dimensional (1D) model to simulate blood flow. Our model was initially tested by comparing blood flow results to those observed from experimental work in VA-ECMO patients. The 1D model was combined with a two-phase flow model to simulate oxygenation. Additionally, the influence of various other clinician-tunable parameters on oxygenation in the common carotid arteries (CCAs) were tested, including, blood viscosity, cannula position within the insertion artery, heart rate, and systemic vascular resistance (SVR), as well as geometrical changes such as arterial radius and length. Our results indicated that blood oxygenation to the brain strongly depended on the cannula insertion site and the VA-ECMO flow rate with a weaker but potentially significant dependence on arterial radius. During femoral cannulation, VA-ECMO flow rates greater than ~4.9 L/min were needed to perfuse the CCAs. However, axillary and central cannulation began to perfuse the CCAs at significantly lower flow (~1 L/min). These results may help explain the incidence of cerebral hypoxia in this patient population and the common need to change cannulation strategies during treatment to address this clinical problem. While this work describes patient-averaged results, determining these relationships between VA-ECMO parameters and cerebral hypoxia is an important step towards future work to develop patient-specific models that clinicians can use to improve outcomes.

1D blood flow↗

A three-dimensional incompressible flow simulation method and its application to the Space Shuttle main engine. II Turbulent flow

An implicit finite difference code cast in general curvilinear coordinates is further developed for three-dimensional incompressible turbulent flows. The code is based on the method of pseudocompressibility and utilizes the Beam and Warming implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is further extended to include composite-grids to overcome the excessive computer memory required for solving turbulent flows in complex three-dimensional geometries. A simple turbulence model is proposed for internal flows. The code is being used for the Space Shuttle Main Engine (SSME) internal flow analyses.

Chang, J. L. C.↗

Transforming Energy Through Computational Excellence: Bringing Low Mach Number Reactive Flow Simulations at the Exascale

PeleLMeX's unique capabilities are allowing for reactive flow modeling at unprecedented scales and a reasonable time and cost. The code is currently being extended to tackle more practical, design-oriented simulations by implementing Large Eddy Simulation and data-driven chemical models, providing a fast but accurate tool for engineers considering the emergence of GPU-accelerated platforms. These extensions are critical for enabling the physical insight required to design the next generation of combustion devices as a key component of a renewable energy future.

MATHEMATICS AND COMPUTING↗