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 1,027 records · Page 57

Direct simulation of typical AOTV entry flows

The direct simulation Monte Carlo method is applied to two flows that are relevant to the projected aeroassisted orbital transfer vehicle. The first is the axially symmetric flow about a 70 degree blunted cone shield with a cylindrical afterbody, and the second is the one-dimensional flow along the stagnation streamline of the first flow. The shock standoff distance for the stagnation streamline study is obtained from the full axially symmetric calculation. Advantage is then taken of the speed of the one-dimensional calculation to study the sensitivity of the stagnation point heat transfer to a number of physical variables. These include the chemical reaction rate coefficients, the surface catalycity, the electric field associated with ambipolar diffusion, and the nonequilibrium thermal radiation from the shock front.

Bird, G. A.↗

Wind Plant Flow Physics and Power Performance in Complex Environments: Cooperative Research and Development (Final Report)

Cornell University will partner with NLR on the topic of wind farm wake effects to improve understanding of interactions between complex atmospheric flows, terrain, and wind turbine wakes and plant efficiency. Wind plant flow simulation tools will also be validated. The work performed will help improve wind farm modeling by analyzing data, applying models, designing and performing experiments to acquire additional wind farm data, and develop better models.

17 WIND ENERGY↗

Multigrid Strategies for Viscous Flow Solvers on Anisotropic Unstructured Meshes

Unstructured multigrid techniques for relieving the stiffness associated with high-Reynolds number viscous flow simulations on extremely stretched grids are investigated. One approach consists of employing a semi-coarsening or directional-coarsening technique, based on the directions of strong coupling within the mesh, in order to construct more optimal coarse grid levels. An alternate approach is developed which employs directional implicit smoothing with regular fully coarsened multigrid levels. The directional implicit smoothing is obtained by constructing implicit lines in the unstructured mesh based on the directions of strong coupling. Both approaches yield large increases in convergence rates over the traditional explicit full-coarsening multigrid algorithm. However, maximum benefits are achieved by combining the two approaches in a coupled manner into a single algorithm. An order of magnitude increase in convergence rate over the traditional explicit full-coarsening algorithm is demonstrated, and convergence rates for high-Reynolds number viscous flows which are independent of the grid aspect ratio are obtained. Further acceleration is provided by incorporating low-Mach-number preconditioning techniques, and a Newton-GMRES strategy which employs the multigrid scheme as a preconditioner. The compounding effects of these various techniques on speed of convergence is documented through several example test cases.

Movriplis, Dimitri J.↗

Estimating Spatial Distribution Impacts of Rooftops Solar PV on Dynamic Hosting Capacity Evaluation for a Real Distribution Feeder

This paper models the sensitivity of dynamic hosting capacity analysis to the spatial distribution of distributed PV scenarios, random, close and far from substation deployments, at various penetration levels on a real distribution feeder, in a quasi-static time series (QSTS) power flow simulation using high resolution time-aware metrics. These spatial distribution scenarios were chosen to capture a wide range of potential operation impacts of these deployments using a set of thermaland voltage-based metrics such as instantaneous violations and moving averages for both thermal and voltage constraints. This study uses actual load and PV data up-scaled from 1-hour time step to 1-minute resolution to fully characterize the interaction between the daily changes in load and PV output, and their impacts on distribution system operations.

distribution system↗

Characterizing the Impact of Fractured Caprock Heterogeneity on Supercritical CO 2 Injection

We present a set of multiphase flow simulations where supercritical CO 2 (scCO 2 ) displaces water at hydrostatic conditions within three-dimensional discrete fracture networks that represent paths for potential leakage through caprock above CO 2 storage reservoirs. The simulations are performed to characterize and compare the relative impact of hydraulic and structural heterogeneity in fractured media on the initial movement of scCO 2 through these caprock formations. In one scenario, intrinsic fracture permeabilities are varied stochastically within a fixed network structure. In another scenario, we generate multiple independent, identically distributed network realizations with varying fracture network densities to explore a wide range of geometric and topological configurations. Analysis of the simulations indicates that network structure, specifically connectivity and the presence of hanging fractures, plays a larger role in controlling the displacement of water by scCO 2 than variations in local hydraulic properties. Finally, we identify active surface area of the network as a single-phase feature that could provide a lower bound on the percentage of the network surface area reached by scCO 2 .

54 ENVIRONMENTAL SCIENCES↗

Wind-US Users Guide Version 3.0

Wind-US is a computational platform which may be used to numerically solve various sets of equations governing physical phenomena. Currently, the code supports the solution of the Euler and Navier-Stokes equations of fluid mechanics, along with supporting equation sets governing turbulent and chemically reacting flows. Wind-US is a product of the NPARC Alliance, a partnership between the NASA Glenn Research Center (GRC) and the Arnold Engineering Development Complex (AEDC) dedicated to the establishment of a national, applications-oriented flow simulation capability. The Boeing Company has also been closely associated with the Alliance since its inception, and represents the interests of the NPARC User's Association. The "Wind-US User's Guide" describes the operation and use of Wind-US, including: a basic tutorial; the physical and numerical models that are used; the boundary conditions; monitoring convergence; the files that are read and/or written; parallel execution; and a complete list of input keywords and test options. For current information about Wind-US and the NPARC Alliance, please see the Wind-US home page at http://www.grc.nasa.gov/WWW/winddocs/ and the NPARC Alliance home page at http://www.grc.nasa.gov/WWW/wind/. This manual describes the operation and use of Wind-US, a computational platform which may be used to numerically solve various sets of equations governing physical phenomena. Wind-US represents a merger of the capabilities of four CFD codes - NASTD (a structured grid flow solver developed at McDonnell Douglas, now part of Boeing), NPARC (the original NPARC Alliance structured grid flow solver), NXAIR (an AEDC structured grid code used primarily for store separation analysis), and ICAT (an unstructured grid flow solver developed at the Rockwell Science Center and Boeing).

Navier-Stokes equation↗

Optimal gas cooling with minimal aero-opticaldistortion for next generation high-average-power lasers (Progress report, Year 2)

Next-generation lasers that combine high-peak-power ultrashort pulses with high repetition rates have the potential to efficiently deliver the high-average-power outputs required for inertial fusion energy generation, laser-based particle acceleration, and to generate compact secondary x- and 𝛾-ray radiation sources, among others. The combination of long gain length with high surface-to-volume ratio to remove residual heat makes gas-cooled multi-slab amplifier heads an inherently aperture- and average power scalable technology that enables this new class of solid-state lasers, and their thermal management becomes the most pressing technical challenge. As the index-of refraction depends on the gas density, heat removal from the slab surfaces unavoidably entails scattering of propagated light from regions of high temperature fluctuations. The requirement to maximize convective heat extraction from the gain medium under the constraint of minimum optical aberration of the beam renders amplifier slab cooling a technological challenge at the intersection of aero-optics, fluid-mechanics and heat transfer. The goal of this research is to leverage computational fluid dynamics and aero optics simulations of different levels of fidelity to devise a numerical modeling and optimization framework for gas-cooled multi-slab amplifiers. The project will be conducted in close collaboration with the Advanced Photon Technologies Group at Lawrence Livermore National Laboratory, where a gas-cooled test stand for laser materials and thermal management structures will be developed and installed. The final goal of the collaborative effort is to manufacture and benchmark a physical prototype of an numerically optimized amplifier head design. Measurement data on solid-state and thermo-optical distortions will be used for the validation of the numerical models. The key technical challenge in this multi-objective design problem is that of balancing optimal heat removal with minimal optical distortion. To facilitate computationally efficient design optimization and, at the same time, deepen our basic understanding of the interaction of light with turbulence, a complementary study that combines high and low-fidelity flow simulations and aero-optical models is proposed. Along the lines of these two objectives and solution strategies, the project is organized into two thrusts.

42 ENGINEERING↗

Heat‐Generating Nuclear Waste in Salt: Field Testing and Simulation

Core Ideas A field‐scale experiment and numerical simulations confirm salt backfill behavior. Simulations closely match temperature around and under the piled salt backfill. Results indicate limited dissolution–precipitation reactions around the heat source. Alteration of backfill is unlikely if the drift is allowed to dry before emplacement. Investigations relating to in‐drift disposal of heat‐generating nuclear waste in salt have raised questions about heat–brine interactions in the unsaturated run‐of‐mine (RoM) salt pile used as backfill. These interactions have the potential to change the structure of the RoM salt surrounding the canister, possibly altering long‐term containment of the source. An experiment is in progress at the Waste Isolation Pilot Plant (WIPP), New Mexico, in which a heated canister was placed on the floor of an open drift, covered in a pile of RoM salt, and energized with 1000 W. Temperature in the RoM salt pile had stabilized after about 15 d, allowing evaluation of the heat‐up period of the ongoing experiment. Using a multiphase porous flow simulator that has been modified to handle salt‐specific coupled processes, we examined coupled thermal–hydrological–chemical behavior in the RoM salt pile. Our simulations suggest that for the relatively dry cases examined, porosity changes within RoM salt in a generic salt repository are likely to be minor in the period between waste emplacement and plastic closure of the drift. The primary sensitivity for porosity change is to the early moisture content of the RoM salt used to cover the canister. Secondary influences include moisture availability from the disturbed rock zone (DRZ) surrounding the drift and the capillary pressure ratio between the DRZ and the RoM salt. Early changes in porosity and permeability may be affected by moisture content, but this was not observed in the test. Such changes would be most likely to occur when using damp RoM salt or if waste is emplaced in a drift immediately following opening of the drift before evaporative dewatering of the drift walls occurs.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Unstructured grid large eddy simulation of wall bounded turbulent flows

Historically, large eddy simulations (LES) have been restricted to simple geometries where spectral or finite difference methods have dominated due to their efficient use of structured grids. Structured grids, however, not only have difficulty representing complex domains and adapting to complicated flow features, but also are rather inefficient for simulating flows at high Reynolds numbers. The lack of efficiency stems from the need to resolve the viscous sub layer which requires very fine resolution in all three directions near the wall. Structured grids make use of a stretching to reduce the normal grid spacing but must carry the fine resolution in the streamwise and spanwise directions throughout the domain. The unnecessarily fine grid for much of the domain leads to disturbingly high grid estimates. Chapman (1979), and later Moin & Jimenez (1993), pointed out that, in order to advance the technology to airfoils at flight Reynolds numbers, structured grids must be abandoned in lieu of what are known as nested or unstructured grids. We illustrate the ability of an unstructured mesh to refine only the near wall region. Note the large number of points near the wall (where the fine vortical features need better resolution) and the coarseness in all directions away from the wall (where the scales are much larger). The important difference between this approach and the usual structured grid stretching is that the number of elements used to discretize the spanwise and streamwise features of the flow is reduced in each successive layer coming off the wall. This is due to the fact that the elements not only grow in the normal direction, but in the other directions as well. This greatly reduces the total number of points or elements required for a given Reynolds number flow.

Jansen, Kenneth↗

The Effect of Anode Position on Operation of a 25-A Class Hollow Cathode

The effect of the distance between a hollow cathode and a cylindrical anode on cathode operation is investigated for two anode geometries. Neutral flow simulations demonstrate that the anode diameter and distance from the cathode exit can elevate the downstream pressure as much as two orders of magnitude above what the cathode experiences while operating within a Hall thruster. Based on the results of this modeling, two axially-segmented cylindrical molybdenum anodes were constructed: a 64-mm diameter one that replicated the anode geometry used in recent NASA hollow cathode development testing and a larger 254-mm diameter one designed to reduce the neutral pressure in front of the cathode to thruster-like values. For each anode design, cathode performance was characterized for varying anode/cathode distance using metrics such as discharge voltage and oscillation magnitudes, and the ion voltage spectra were characterized using a radially-positioned retarding potential analyzer. It was found that as local neutral pressure decreased, discharge voltage and high-voltage ion content in the plume increased. For the 254-mm diameter anode, an ion voltage tail in excess of 200 V was found for nominal cathode flow rates. The implications of these results for standalone hollow cathode development tests are discussed.

Hall, Scott J.↗

Multiphase flow studies for microscale hydrodynamics in the structured packed column

Post combustion carbon capture by solvent absorption in a structured packed column is a promising technology for mitigating greenhouse gas emissions. Computational fluid dynamic (CFD) modelling of such a column is a challenging multiscale problem due to the range in length and time scales. The microscale hydrodynamics play a key role in the overall column efficiency with the interfacial area significantly influencing the mass transfer between the gas and liquid phases. In this context, multiphase flow simulations using the volume of fluid (VOF) method in a representative elementary unit (REU) of the packed column can provide fundamental insights into the microscale hydrodynamics, such as, interfacial area and liquid holdup. The present paper systematically examines the impact of various factors (e.g., physical properties and contact angle) on the interfacial area. The results are compared with existing correlations and a scaling analysis is also performed. The solvent physical properties are characterized by the Kapitza number (Ka), a dimensionless number that depends only on fluid properties. So that the Ka number decreases with increasing viscosity. At a fixed liquid load, the interfacial area and liquid holdup are observed to increase with decreasing Ka number. The impact of contact angle (i.e., solid surface characteristics) is effectively investigated by modifying the wall boundary conditions. The interfacial area and liquid holdup are found to decrease with increasing contact angle. Subsequently, a phenomenological correlation for interfacial area is proposed that includes the impact of these parameters. Finally, this correlation may be used to predict the interfacial area for gas-liquid flow in a structured packing for rivulet to fully wetted flow regimes.

42 ENGINEERING↗

Hawaii Play Fairway Analysis: Oahu Groundwater Flow Model

Groundwater flow model for the island of Oahu. Data is from the following sources: Rotzoll, K., A.I. El-Kadi. 2007. Numerical Ground-Water Flow Simulation for Red Hill Fuel Storage Facilities, NAVFAC Pacific, Oahu, Hawaii - Prepared TEC, Inc. Water Resources Research Center, University of Hawaii, Honolulu.; Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report - Volume VII - Island of Oahu Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008.; and Whittier, R. and A.I. El-Kadi. 2009. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems - Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. December 2009.

15 GEOTHERMAL ENERGY↗

DRAGON Grid: A Three-Dimensional Hybrid Grid Generation Code Developed

Because grid generation can consume 70 percent of the total analysis time for a typical three-dimensional viscous flow simulation for a practical engineering device, payoffs from research and development could reduce costs and increase throughputs considerably. In this study, researchers at the NASA Glenn Research Center at Lewis Field developed a new hybrid grid approach with the advantages of flexibility, high-quality grids suitable for an accurate resolution of viscous regions, and a low memory requirement. These advantages will, in turn, reduce analysis time and increase accuracy. They result from an innovative combination of structured and unstructured grids to represent the geometry and the computation domain. The present approach makes use of the respective strengths of both the structured and unstructured grid methods, while minimizing their weaknesses. First, the Chimera grid generates high-quality, mostly orthogonal meshes around individual components. This process is flexible and can be done easily. Normally, these individual grids are required overlap each other so that the solution on one grid can communicate with another. However, when this communication is carried out via a nonconservative interpolation procedure, a spurious solution can result. Current research is aimed at entirely eliminating this undesired interpolation by directly replacing arbitrary grid overlapping with a nonstructured grid called a DRAGON grid, which uses the same set of conservation laws over the entire region, thus ensuring conservation everywhere. The DRAGON grid is shown for a typical film-cooled turbine vane with 33 holes and 3 plenum compartments. There are structured grids around each geometrical entity and unstructured grids connecting them. In fiscal year 1999, Glenn researchers developed and tested the three-dimensional DRAGON grid-generation tools. A flow solver suitable for the DRAGON grid has been developed, and a series of validation tests are underway.

Liou, Meng-Sing↗

Numerical simulation of the flow about the F-18 HARV at high angle of attack

As part of NASA's High Alpha Technology Program, research has been aimed at developing and extending numerical methods to accurately predict the high Reynolds number flow about the NASA F-18 High Alpha Research Vehicle (HARV) at large angles of attack. The HARV aircraft is equipped with a bidirectional thrust vectoring unit which enables stable, controlled flight through 70 deg angle of attack. Currently, high-fidelity numerical solutions for the flow about the HARV have been obtained at alpha = 30 deg, and validated against flight-test data. It is planned to simulate the flow about the HARV through alpha = 60 deg, and obtain solutions of the same quality as those at the lower angles of attack. This report presents the status of work aimed at extending the HARV computations to the extreme angle of attack range.

Murman, Scott M.↗

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗

Large-scale large eddy simulation of nuclear reactor flows: Issues and perspectives

Numerical simulation has been an intrinsic part of nuclear engineering research since its inception. In recent years a transition is occurring toward predictive, first-principle-based tools such as computational fluid dynamics. Even with the advent of petascale computing, however, such tools still have significant limitations. In the present work some of these issues, and in particular the presence of massive multiscale separation, are discussed, as well as some of the research conducted to mitigate them. Petascale simulations at high fidelity (large eddy simulation/direct numerical simulation) were conducted with the massively parallel spectral element code Nek5000 on a series of representative problems. These simulations shed light on the requirements of several types of simulation: (1) axial flow around fuel rods, with particular attention to wall effects; (2) natural convection in the primary vessel; and (3) flow in a rod bundle in the presence of spacing devices. Finally, the focus of the work presented here is on the lessons learned and the requirements to perform these simulations at exascale. Additional physical insight gained from these simulations is also emphasized.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Discrete-vortex simulation of pulsating flow on a turbulent leading-edge separation bubble

Studies are made of the turbulent separation bubble in a two-dimensional semi-infinite blunt plate aligned to a uniform free stream with a pulsating component. The discrete-vortex method is applied to simulate this flow situation because this approach is effective for representing the unsteady motions of the turbulent shear layer and the effect of viscosity near the solid surface. The numerical simulation provides reasonable predictions when compared with the experimental results. A particular frequency with a minimum reattachment is related to the drag reduction. The most effective frequency is dependent on the amplified shedding frequency. The turbulent flow structure is scrutinized. This includes the time-mean and fluctuations of the velocity and the surface pressure, together with correlations between the fluctuating components. A comparison between the pulsating flow and the non-pulsating flow at the particular frequency of the minimum reattachment length of the separation bubble suggests that the large-scale vortical structure is associated with the shedding frequency and the flow instabilities.

Sung, Hyung Jin↗