SEARCH · Engineering Papers

Results for “implicit pressure correction”

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 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Implicit/explicit analysis of interactive phenomena in: supersonic, chemically-reacting, mixing and boundary layer problems

An efficient implicit/explicit procedure for solving the parabolized Navier-Stokes (PNS) equations to analyze supersonic, viscous, chemically-reacting flows is presented. The approach combines a fully-implicit algorithm for solving the parabolic mixing/boundary layer equations with an explicit algorithm for determining the hyperbolic pressure field. The unified model, SPLITP, contains provisions for treating bounding/embedded subsonic regions using a pressure-splitting procedure which suppresses elliptic (upstream influence) effects without recourse to sublayer approximations. Several two-equation turbulence models with compressibility and curvature-correction terms are incorporated into SPLITP. Applications to a variety of flow problems are presented which include: ducted supersonic mixing and combustion flowfields, free jets and plumes, and, curved boundary layers and wall jets.

Dash, S. M.↗

NASA NTRS1991

A Navier-Stokes Solution of Hull-Ring Wing-Thruster Interaction

Navier-Stokes simulations of high Reynolds number flow around an axisymmetric body supported in a water tunnel were made. The numerical method is based on a finite-differencing high resolution second-order accurate implicit upwind scheme. Four different configurations were investigated, these are: (1) barebody; (2) body with an operating propeller; (3) body with a ring wing; and (4) body with a ring wing and an operating propeller. Pressure and velocity components near the stern region were obtained computationally and are shown to compare favorably with the experimental data. The method correctly predicts the existence and extent of stern flow separation for the barebody and the absence of flow separation for the three other configurations with ring wing and/or propeller.

Yang, C.-I.↗

NASA NTRS2000

On the Coupling of CDISC Design Method with FPX Rotor Code

A rotor section aerodynamics design package is developed by coupling Constrained Direct Iterative Surface Curvature (CDISC) design method with the FPX rotor code. The coupling between the CDISC design and the FPX flow analysis is fully automated. The CDISC design method employs a predictor-corrector procedure iteratively to determine a surface geometry which produces a target pressure distribution, where the target pressure distributions is either pre-defined or automatically generated through flow and geometry constraints. The FPX code is an eXtended Full-Potential rotor Computational Fluid Dynamics (CFD) code, which solves the three-dimensional unsteady full-potential equation in a strong conservative form using an implicit approximate factorization finite-difference scheme with entropy and viscosity corrections. Application of the CDISC design method coupled with the FPX rotor code is made for rotor blades in hovering motions. Several design examples are presented to demonstrate the capability of the new package in rotor section design.

Hu, Hong↗

NASA NTRS1991

Numerical study of high-area-ratio H2/O2 rocket nozzles

Calculations were made for high-area-ratio H2/O2 rocket nozzles by using the RPLUS code which employs an implicit finite volume, lower-upper symmetric successive over-relaxation scheme to solve the Navier-Stokes equations and the species equations in a coupled manner. The combustion processes of hydrogen and oxygen are modeled by an 8-species and 18-step reaction mechanism, and the turbulence is simulated by the Baldwin-Lomax turbulence model with a pressure gradient correction on the van Driest's damping constant in the inner-layer model. The calculated results for two different rocket nozzles are presented as contours, profiles, and performance values. The predicted specific impulse and thrust agree well with the experimental data and those from the industry standard code. It is demonstrated that the RPLUS code can be used for analysis of rocket nozzles.

Kim, S. C.↗

NASA NTRS1982

Implicit marching solution of compressible viscous subsonic flow in planar and axisymmetric ducts

A new streamwise marching procedure was developed and coded for compressible viscous subsonic flow in planar or axisymmetric ducts with or without centerbodies. The continuity, streamwise momentum, cross-flow momentum, and energy equations are written in generalized orthogonal curvilinear coordinates. To allow the use of a marching procedure, second derivatives in the streamwise momentum equation are written as the sum of a known two dimensional imposed pressure field and an unknown one dimensional viscous correction. For turbulent flow, the Reynolds stress and heat flux terms are modeled using two-layer eddy viscosity turbulence models.

Towne, C. E.↗

NASA NTRS1982

Conservative supra-characteristics method for splitting the hyperbolic systems of gasdynamics for real and perfect gases

A conservative flux difference splitting is presented for the hyperbolic systems of gasdynamics. The stable robust method is suitable for wide application in a variety of schemes, explicit or implicit, iterative or direct, for marching in either time or space. The splitting is modeled on the local quasi one dimensional characteristics system for multi-dimensional flow similar to Chakravarthy's nonconservative split coefficient matrix method; but, as the result of maintaining global conservation, the method is able to capture sharp shocks correctly. The embedded characteristics formulation is cast in a primitive variable the volumetric internal energy (rather than the pressure) that is effective for treating real as well as perfect gases. Finally the relationship of the splitting to characteristics boundary conditions is discussed and the associated conservative matrix formulation for a computed blown wall boundary condition is developed as an example. The theoretical development employs and extends the notion of Roe of constructing stable upwind difference formulae by sending split simple one sided flux difference pieces to appropriate mesh sites. The developments are also believed to have the potential for aiding in the analysis of both existing and new conservative difference schemes.

Lombard, C. K.↗

NASA NTRS1983

Multi-grid solution of Neumann pressure problem for viscous flows using primitive variables

The multi-grid (MG) technique has been advanced for use with the Neumann boundary-value problem in clustered curvilinear orthogonal coordinates. This comprises an important step in the analysis of viscous flows using the velocity-pressure formulation of the Navier-Stokes equations. With successive over-relaxation (SOR) as the smoothing operator and with suitably formulated restriction and coarse-grid-correction operators, a 4-grid procedure enhances the efficiency of fine-grid solutions of the Neumann problem by a factor of 3 to 5, depending on the problem parameters. Thy influence of the smoothing operator is also examined by employing the alternating-direction implicit and the strongly implicit techniques instead of SOR.

Ghia, U.↗

NASA NTRS1988

Solution of the Neumann pressure problem in general orthogonal coordinates using the multigrid technique

The multigrid (MG) technique has been advanced for use with Neumann boundary-value problems in clustered curvilinear orthogonal coordinates. This comprises an important step in the analysis of incompressible flow using the velocity-pressure formulation of the Navier-Stokes equations. The finite-difference representation of the problem and the formulation of the restriction and coarse-grid correction operators are examined in detail. Maintaining consistency between these and the integral constraint associated with the Neumann problem is found to be critical for the success of the MG technique. The influence of the smoothing operator is examined by employing Gauss-Seidel, alternating-direction implicit, and strongly implicit techniques. The MG procedure enhances the efficiency of fine-grid solutions of the Neumann problem by a factor of 3 to 14, depending on the type of smoothing operator employed and the values of the problem parameters.

Ghia, U.↗

NASA NTRS1986

Viscous compressible flow direct and inverse computation and illustrations

An algorithm for laminar and turbulent viscous compressible two dimensional flows is presented. For the application of precise boundary conditions over an arbitrary body surface, a body-fitted coordinate system is used in the physical plane. A thin-layer approximation of tne Navier-Stokes equations is introduced to keep the viscous terms relatively simple. The flow field computation is performed in the transformed plane. A factorized, implicit scheme is used to facilitate the computation. Sample calculations, for Couette flow, developing pipe flow, an isolated airflow, two dimensional compressor cascade flow, and segmental compressor blade design are presented. To a certain extent, the effective use of the direct solver depends on the user's skill in setting up the gridwork, the time step size and the choice of the artificial viscosity. The design feature of the algorithm, an iterative scheme to correct geometry for a specified surface pressure distribution, works well for subsonic flows. A more elaborate correction scheme is required in treating transonic flows where local shock waves may be involved.

Yang, T. T.↗

NASA NTRS1988

A point implicit unstructured grid solver for the Euler and Navier-Stokes equations

An upwind finite element technique that uses cell centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured triangular grids. A structured grid of quadrilaterals is laid out near the body surface. For inviscid flows the method is stable at Courant numbers of over 100,000. A first order basic scheme and a higher order flux corrected transport (FCT) scheme have been implemented. This technique has been applied to the problem of predicting type III and IV shock wave interactions on a cylinder, with a view of simulating the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively generated grids. The adaptive mesh generator and the upwind Navier-Stokes solver are combined in a set of programs called LARCNESS, an acronym for Langley Adaptive Remeshing Code and Navier-Stokes Solver.

Thareja, Rajiv R.↗

NASA NTRS1989

Upwind-biased, point-implicit relaxation strategies for hypersonic flowfield simulations on supercomputers

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The relaxation strategy is well suited for computers employing either vector or parallel architectures, and the relation between computer architecture and algorithm is emphasized. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large. block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over a blunt body. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km.) over the Aeroassist Flight Experiment (AFE) configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗

NASA NTRS1989

Euler and Navier-Stokes solutions for hypersonic flows

An upwind finite-element technique that uses cell-centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured grids in two and three dimensions. A perfect gas model as well as an equilibrium air model is implemented for solving high-speed flows. A first-order basic scheme and a higher-order flux-corrected transport (FCT) scheme have been implemented. This technique has been used to predict 'Type III and IV' shock interactions on a cylinder in two dimensions and a swept cylinder in three dimensions, with a view to determine the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively-generated grids. Three-dimensional corner flow, typically encountered in engine inlets due to compression of the flow by ramps in the walls, is also modeled. This procedure is the first step in developing an integrated fluid, thermal, structural analysis capability for hypersonic flight vehicles like the National Aero-Space Plane.

Thareja, Rajiv R.↗

NASA NTRS1991

Three-dimensional blade vortex interactions

A three-dimensional time dependent Navier-Stokes analysis was applied to the rotor blade vortex interaction problem. The numerical procedure is an iterative implicit procedure using three point central differences to represent spatial derivatives. A series of calculations were made to determine the time steps, pseudo-time steps, iterations, artificial dissipation level, etc. required to maintain a nondissipative vortex. Results show the chosen method to have excellent non-dissipative properties provided the correct parameters are chosen. This study was used to set parameters for both two- and three-dimensional blade vortex interaction studies. The case considered was the interaction between a vortex and a NACA0012 airfoil. The results showed the detailed physics during the interaction including the pressure pulse propagating from the blade. The simulated flow physics was qualitatively similar to that experimentally observed. The BVI phenomena is the result of the buildup and violent collapse of the shock waves and local supersonic pockets on the blade surfaces. The resulting pressure pulse build-up appears to be centered at the blade leading edge.

Davoudzadeh, Farhad↗

DOE OSTI2024

Avoiding pitfalls in molecular simulation of vapor sorption: Example of propane and isobutane in metal–organic frameworks for adsorption cooling applications

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal–organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance.

36 MATERIALS SCIENCE↗

NASA NTRS1990

Two- and three-dimensional blade vortex interactions

A three-dimensional time dependent Navier-Stokes analysis was applied to the rotor blade vortex interaction (BVI) problem. The numerical procedure is an iterative implicit procedure using three point central differences to represent spatial derivatives. A series of calculations were made to determine the time steps, pseudo-time steps, iterations, artificial dissipation level, etc. required to maintain a nondissipative vortex. Results show the chosen method to have excellent non-dissipative properties provided the correct parameters are chosen. This study was used to set parameters for both two- and three-dimensional blade vortex interaction studies. The two-dimensional study considered the interaction between a vortex and a NACA0012 airfoil. The results showed the detailed physics during the interaction including the pressure pulse propagating from the blade. The simulated flow physics was qualitatively similar to that experimentally observed. The 2-D BVI phenomena is the result of the buildup and violent collapse of the shock waves and local supersonic pockets on the blade surfaces. The resulting pressure pulse build-up appears to be centered at the blade leading edge. The three-dimensional interaction study considered the case of a vortex at 20 deg incidence to the blade leading edge. Although the qualitative results were similar to that of the two-dimensional interaction, details clearly showed the three-dimensional nature of the interaction process.

Davoudzadeh, F.↗

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 NTRS1998

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗

NASA NTRS1990

An upwind-biased, point-implicit relaxation algorithm for viscous, compressible perfect-gas flows

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, compressible, perfect-gas flows is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Viscous terms are discretized using central differences. The relaxation strategy is well suited for computers employing either vector or parallel architectures. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. Convergence rates and grid refinement studies are conducted for Mach 5 flow through an inlet with a 10 deg compression ramp and Mach 14 flow over a 15 deg ramp. Predictions for pressure distributions, surface heating, and aerodynamics coefficients compare well with experiment data for Mach 10 flow over a blunt body.

Gnoffo, Peter A.↗