Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Volume Scheme”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗
NASA NTRSDate not supplied

WMLES for the Fifth High-Lift PredictionWorkshop Cases Using FUN3D

This paper presents solution assessments and grid convergence studies for the test cases outlined in the Fifth High-Lift Prediction Workshop (HLPW-5), focusing on the high-lift Common Research Models (CRM-HL). The study utilizes a wall-modeled large-eddy simulation (WMLES) methodology developed in the unstructured-grid, node-centered flow solver FUN3D. The second-order accurate simulations conducted in this study utilize a finite-volume spatial discretization and an implicit temporal scheme. Large-scale turbulent features are resolved away from the wall, with small-scale effects captured by the Vreman subgrid-scale model. An equilibrium wall function uses the first grid point off the wall serving as the critical interface between the wall model and the large-eddy simulation region, thus requiring careful placement in grid design. WMLES solutions are assessed for HLPW-5 cases, including a clean wing-body configuration and geometry-buildup configurations corresponding to the 5.1\% ONERA CRM-HL model. Grid-convergence studies are systematically conducted using uniformly refined grids. Moreover, simulation results and grid sensitivity are presented for the NASA 5.2\% CRM-HL configuration at both moderate and flight-scale Reynolds numbers. Overall, WMLES results are satisfactory and agree well with available experimental data, especially on sufficiently fine grids.

high-lift aerodynamics↗
DOE OSTI2024

Finite domain solution of a KGD hydraulic fracture in the viscosity-dominated regime

This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future.

Domain-based method↗
NASA NTRS1995

A comparative study of upwind and MacCormack schemes for CAA benchmark problems

In this study, upwind schemes and MacCormack schemes are evaluated as to their suitability for aeroacoustic applications. The governing equations are cast in a curvilinear coordinate system and discretized using finite volume concepts. A flux splitting procedure is used for the upwind schemes, where the signals crossing the cell faces are grouped into two categories: signals that bring information from outside into the cell, and signals that leave the cell. These signals may be computed in several ways, with the desired spatial and temporal accuracy achieved by choosing appropriate interpolating polynomials. The classical MacCormack schemes employed here are fourth order accurate in time and space. Results for categories 1, 4, and 6 of the workshop's benchmark problems are presented. Comparisons are also made with the exact solutions, where available. The main conclusions of this study are finally presented.

Viswanathan, K.↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
DOE OSTI2020

Moments of nucleon isovector structure functions in 2+1+1-flavor QCD

We present results on the isovector momentum fraction, $\langle x \rangle_{u–d}$ , helicity moment, $\langle x \rangle_{Δ u–Δd}$ , and the transversity moment, $\langle x \rangle_{δu–δd}$, of the nucleon obtained using nine ensembles of gauge configurations generated by the MILC Collaboration using 2 + 1 + 1 -flavors of dynamical highly improved staggered quarks. The correlation functions are calculated using the Wilson-Clover action, and the renormalization of the three operators is carried out nonperturbatively on the lattice in the RI'–MOM scheme. The data have been collected at lattice spacings a ≈ 0.15 , 0.12, 0.09, and 0.06 fm and $M_π$ ≈ 310 , 220, and 135 MeV, which are used to obtain the physical values using a simultaneous chiral-continuum-finite-volume fit. The final results, in the $\overline{\text{MS}}$ scheme at 2 GeV, are $\langle$ x $\rangle$ u – d = 0.173 ( 14 ) ( 07 ) , $\langle x \rangle_{Δu–Δd}$ = 0.213 ( 15 ) ( 22 ) , and $\langle x \rangle_{δu–δd}$ = 0.208 ( 19 ) ( 24 ) , where the first error is the overall analysis uncertainty and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with other recent lattice calculations and phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1981

Numerical method for gas dynamics combining characteristic and conservation concepts

An efficient implicit numerical method that solves the compressible Navier-Stokes equations in arbitrary curvilinear coordinates by the finite-volume technique is presented. An intrinsically dissipative difference scheme and a fully implicit treatment of boundary conditions, based on characteristic and conservation concepts, are used to improve stability and accuracy. Efficiency is achieved by using a diagonal form of the implicit algorithm and spatially varying time-steps. Comparisons of various schemes and methods are presented for one- and two-dimensional flows, including transonic separated flow past a thick circular-arc airfoil in a channel. The new method is equal to or better than a version of MacCormack's hybrid method in accuracy and it converges to a steady state up to an order of magnitude faster.

Coakley, T. J.↗
NASA NTRS1986

An Euler aerodynamic method for leading-edge vortex flow simulation

The current capabilities and the future plans for a three dimensional Euler Aerodynamic Method are described. The basic solution algorithm is based on the finite volume, Runge-Kutta pseudo-time-stepping scheme of FLO-57. Several modifications to improve accuracy and computational efficiency were incorporated and others are being investigated. The computer code is used to analyze a cropped delta wing at 0.6 Mach number and an arrow wing at 0.85 Mach number. Computed aerodynamic parameters are compared with experimental data. In all cases, the configuration is impulsively started and no Kutta condition is applied at sharp edges. The results indicate that with additional development and validation, the present method will be a useful tool for engineering analysis of high speed aircraft.

Raj, P.↗
NASA NTRS1988

Investigation of supersonic chemically reacting and radiating channel flow

The 2-D time-dependent Navier-Stokes equations are used to investigate supersonic flows undergoing finite rate chemical reaction and radiation interaction for a hydrogen-air system. The explicit multistage finite volume technique of Jameson is used to advance the governing equations in time until convergence is achieved. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The multidimensional radiative transfer equations for a nongray model are provided for a general configuration and then reduced for a planar geometry. Both pseudo-gray and nongray models are used to represent the absorption-emission characteristics of the participating species. The supersonic inviscid and viscous, nonreacting flows are solved by employing the finite volume technique of Jameson and the unsplit finite difference scheme of MacCormack. The specified problem considered is of the flow in a channel with a 10 deg compression-expansion ramp. The calculated results are compared with those of an upwind scheme. The problem of chemically reacting and radiating flows are solved for the flow of premixed hydrogen-air through a channel with parallel boundaries, and a channel with a compression corner. Results obtained for specific conditions indicate that the radiative interaction can have a significant influence on the entire flow field.

Mani, Mortaza↗
NASA NTRS1987

Calculation of wall and free turbulent-shear flows at supersonic speeds

Supersonic turbulent flows are simulated numerically by solving the Reynolds-averaged full Navier-Stokes equations by an implicit finite-volume method. This flux-vector-split upwind scheme uses approximate factorization followed by line-Gauss-Seidel relaxations. The effects of turbulence are modeled by two eddy viscosity models. High-Reynolds-number form of the k-epsilon model is coupled with a wall-function to avoid excessive refinement of the grid in the low-Reynolds number regions. The k-epsilon equations are solved by the explicit-implicit MacCormack method. The algebraic Baldwin-Lomax model is also being used as an inexpensive alternative for the cases which do not experience massive separation. Several examples of two-dimensional solutions are given to illustrate both wall and free turbulent shear flows which include fluid dynamic phenomena, such as shocks, boundary layers, shear layers, wakes, separations and recirculations. The results compared with experimental data show good agreeent.

Baysal, O.↗
NASA NTRS1989

Comparison of two- and three-dimensional Navier-Stokes solutions with NASA experimental data for CAST-10 airfoil

The two-dimensional (2-D) and three-dimensional Navier-Stokes equations are solved for flow over a NAE CAST-10 airfoil model. Recently developed finite-volume codes that apply a multistage time stepping scheme in conjunction with steady state acceleration techniques are used to solve the equations. Two-dimensional results are shown for flow conditions uncorrected and corrected for wind tunnel wall interference effects. Predicted surface pressures from 3-D simulations are compared with those from 2-D calculations. The focus of the 3-D computations is the influence of the sidewall boundary layers. Topological features of the 3-D flow fields are indicated. Lift and drag results are compared with experimental measurements.

Swanson, R. Charles↗
NASA NTRS1991

Time accurate Navier-Stokes computations of solid rocket motor internal field joint cavity flows

The low-speed flow of combustion gases in the Space Shuttle Solid Rocket Motor (SRM) following ignition are investigated with a time-accurate CFD technique. The flow of the combustion gases is examined in the forward joint section of the SRM at 80 s after ignition with attention given to the amplification of existing 'organ-pipe' modes of 15, 30, and 45 Hz resulting from frequency lock-in with vortex-shedding modes. The CFD code employed is a multizone compressible Navier-Stokes flow solver using a finite-volume, upwind-biased, total-variation-diminishing scheme. The vortex-generated oscillations generated by the inhibitor and the cavity at the forward field joint are described in terms of amplitude and frequency and related to the flow dynamics in the SRM. Oscillations are shown to occur spontaneously when the trailing edge of the cavity at the joint is sharp and can be amplified if the unburned portion of the inhibitor protrudes into the flow.

Farr, Rebecca A.↗
NASA NTRS1993

Conical Euler analysis and active roll suppression for unsteady vortical flows about rolling delta wings

A conical Euler code was developed to study unsteady vortex-dominated flows about rolling, highly swept delta wings undergoing either forced motions or free-to-roll motions that include active roll suppression. The flow solver of the code involves a multistage, Runge-Kutta time-stepping scheme that uses a cell-centered, finite-volume, spatial discretization of the Euler equations on an unstructured grid of triangles. The code allows for the additional analysis of the free to-roll case by simultaneously integrating in time the rigid-body equation of motion with the governing flow equations. Results are presented for a delta wing with a 75 deg swept, sharp leading edge at a free-stream Mach number of 1.2 and at 10 deg, 20 deg, and 30 deg angle of attack alpha. At the lower angles of attack (10 and 20 deg), forced-harmonic analyses indicate that the rolling-moment coefficients provide a positive damping, which is verified by free-to-roll calculations. In contrast, at the higher angle of attack (30 deg), a forced-harmonic analysis indicates that the rolling-moment coefficient provides negative damping at the small roll amplitudes. A free-to-roll calculation for this case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation, which is characteristic of highly swept delta wings. This limit cycle oscillation may be actively suppressed through the use of a rate-feedback control law and antisymmetrically deflected leading-edge flaps. Descriptions of the conical Euler flow solver and the free-to roll analysis are included in this report. Results are presented that demonstrate how the systematic analysis of the forced response of the delta wing can be used to predict the stable, neutrally stable, and unstable free response of the delta wing. These results also give insight into the flow physics associated with unsteady vortical flows about delta wings undergoing forced motions and free-to-roll motions, including the active suppression of the wing-rock type phenomenon. The conical Euler methodology developed is directly extend able to three-dimensional calculations.

Lee-Rausch, Elizabeth M.↗
NASA NTRS1993

Calculations of turbulent separated flows

A numerical study of incompressible turbulent separated flows is carried out by using two-equation turbulence models of the K-epsilon type. On the basis of realizability analysis, a new formulation of the eddy-viscosity is proposed which ensures the positiveness of turbulent normal stresses - a realizability condition that most existing two-equation turbulence models are unable to satisfy. The present model is applied to calculate two backward-facing step flows. Calculations with the standard K-epsilon model and a recently developed RNG-based K-epsilon model are also made for comparison. The calculations are performed with a finite-volume method. A second-order accurate differencing scheme and sufficiently fine grids are used to ensure the numerical accuracy of solutions. The calculated results are compared with the experimental data for both mean and turbulent quantities. The comparison shows that the present model performs quite well for separated flows.

Zhu, J.↗
NASA NTRS1995

Moving and adaptive grid methods for compressible flows

This paper describes adaptive grid methods developed specifically for compressible flow computations. The basic flow solver is a finite-volume implementation of Roe's flux difference splitting scheme or arbitrarily moving unstructured triangular meshes. The grid adaptation is performed according to geometric and flow requirements. Some results are included to illustrate the potential of the methodology.

Trepanier, Jean-Yves↗
NASA NTRS1996

Assessment of an Unstructured-Grid Method for Predicting 3-D Turbulent Viscous Flows

A method Is presented for solving turbulent flow problems on three-dimensional unstructured grids. Spatial discretization Is accomplished by a cell-centered finite-volume formulation using an accurate lin- ear reconstruction scheme and upwind flux differencing. Time is advanced by an implicit backward- Euler time-stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one-equation model, which is coupled with a wall function to reduce the number of cells in the sublayer region of the boundary layer. A systematic assessment of the method is presented to devise guidelines for more strategic application of the technology to complex problems. The assessment includes the accuracy In predictions of skin-friction coefficient, law-of-the-wall behavior, and surface pressure for a flat-plate turbulent boundary layer, and for the ONERA M6 wing under a high Reynolds number, transonic, separated flow condition.

Frink, Neal T.↗
NASA NTRS1986

Comparison of Euler and full potential marching techniques for flows over complex configurations

Two recently developed aerodynamic prediction techniques based on the steady full potential equation and the unsteady Euler equations have been applied to a variety of three-dimensional supersonic flow problems exhibiting embedded subsonic regions. Both techniques utilize planar Gauss-Seidel relaxation in the marching direction and approximate factorization in the cross-flow plane. A conservative switching scheme and flux bias technique are employed in the full potential method to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. A new unified approach with finite volume, high accuracy (up to third order) Total Variation Diminishing formulation (based on Roe's scheme) is used in the Euler solver. In the supersonic regions of the flow an 'infinitely large' time step is employed, and a finite time step is applied in the subsonic regions of the flow to reach the steady-state as a time-asymptote. Numerical solutions are obtained for a number of complex configurations, including: (1) an elliptic waverider, (2) a realistic fighter configuration, (3) the Space Shuttle, and (4) a Shuttle-like configuration. Both the Full Potential and Euler numerical results are in good agreement with available experimental data.

Szema, K. Y.↗