Search NASA⌕ Search

SEARCH · Search NASA

Results for “upwind schemes”

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 289 records · Page 16

An Initial Study of the Fundamentals of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

This presentation shows results from an initial study of the fundamental physics of ice-crystal ice accretion using the NASA Propulsion Systems Lab (PSL). Ice accretion due to the ingestion of ice-crystals is being attributed to numerous jet-engine power-loss events. The NASA PSL is an altitude jet-engine test facility which has recently added a capability to inject ice particles into the flow. NASA is evaluating whether this facility, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. The present study utilized a NACA0012 airfoil. The mixed-phase conditions were generated by partially freezing the liquid-water droplets ejected from the spray bars. This presentation shows data regarding (1) the freeze out characteristics of the cloud, (2) changes in aerothermal conditions due to the presence of the cloud, and (3) the ice accretion characteristics observed on the airfoil model. The primary variable in this test was the PSL plenum humidity which was systematically varied for two duct-exit-plane velocities (85 and 135 ms) as well as two particle size clouds (15 and 50 m MVDi). The observed clouds ranged from fully glaciated to fully liquid, where the liquid clouds were at least partially supercooled. The air total temperature decreased at the test section when the cloud was activated due to evaporation. The ice accretions observed ranged from sharp arrow-like accretions, characteristic of ice-crystal erosion, to cases with double-horn shapes, characteristic of supercooled water accretions.

upwind schemes↗

Particle Size Measurements From the First Fundamentals of Ice Crystal Icing Physics Test in the NASA Propulsion Systems Laboratory

This paper presents particle measurements by the Artium Technologies, Inc. Phase Doppler Interferometer and High Speed Imaging instruments from the first Fundamental Ice Crystal Icing Physics test conducted in the NASA Propulsion Systems Laboratory. The work focuses on humidity sweeps at a larger and a smaller median volumetric diameter. The particle size distribution, number density, and water content measured by the Phase Doppler Interferometer and High Speed Imaging instruments from the sweeps are presented and compared. The current capability for these two instruments to measure and discriminate ICI conditions is examined.

Aircraft icing↗

Evaluation of Alternative Altitude Scaling Methods for Thermal Ice Protection System in NASA Icing Research Tunnel

A test was conducted at NASA Icing Research Tunnel to evaluate altitude scaling methods for thermal ice protection system. Two new scaling methods based on Weber number were compared against a method based on Reynolds number. The results generally agreed with the previous set of tests conducted in NRCC Altitude Icing Wind Tunnel where the three methods of scaling were also tested and compared along with reference (altitude) icing conditions. In those tests, the Weber number-based scaling methods yielded results much closer to those observed at the reference icing conditions than the Reynolds number-based icing conditions. The test in the NASA IRT used a much larger, asymmetric airfoil with an ice protection system that more closely resembled designs used in commercial aircraft. Following the trends observed during the AIWT tests, the Weber number based scaling methods resulted in smaller runback ice than the Reynolds number based scaling, and the ice formed farther upstream. The results show that the new Weber number based scaling methods, particularly the Weber number with water loading scaling, continue to show promise for ice protection system development and evaluation in atmospheric icing tunnels.

upwind schemes↗

An Initial Study of the Fundamentals of Ice Crystal Icing Physics in the Nasa Propulsion Systems Laboratory

This paper presents results from an initial study of the fundamental physics of ice-crystal ice accretion using the NASA Propulsion Systems Lab (PSL). Ice accretion due to the ingestion of ice-crystals is being attributed to numerous jet-engine power-loss events. The NASA PSL is an altitude jet-engine test facility which has recently added a capability to inject ice particles into the flow. NASA is evaluating whether this facility, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. The present study utilized a NACA0012 airfoil. The mixed-phase conditions were generated by partially freezing the liquid-water droplets ejected from the spray bars. This paper presents data regarding (1) the freeze out characteristics of the cloud, (2) changes in aerothermal conditions due to the presence of the cloud, and (3) the ice accretion characteristics observed on the airfoil model. The primary variable in this test was the PSL plenum humidity which was systematically varied for two duct-exit-plane velocities (85 and 135 ms) as well as two particle size clouds (15 and 50 m MVDi). The observed clouds ranged from fully glaciated to fully liquid, where the liquid clouds were at least partially supercooled. The air total temperature decreased at the test section when the cloud was activated due to evaporation. The ice accretions observed ranged from sharp arrow-like accretions, characteristic of ice-crystal erosion, to cases with double-horn shapes, characteristic of supercooled water accretions.

Evaporative Cooling↗

Comparisons of Mixed-Phase Icing Cloud Simulations with Experiments Conducted at the NASA Propulsion Systems Laboratory

This paper builds on previous work that compares numerical simulations of mixed-phase icing clouds with experimental data. The model couples the thermal interaction between ice particles and water droplets of the icing cloud with the flowing air of an icing wind tunnel for simulation of NASA Glenn Research Centers (GRC) Propulsion Systems Laboratory (PSL). Measurements were taken during the Fundamentals of Ice Crystal Icing Physics Tests at the PSL tunnel in March 2016. The tests simulated ice-crystal and mixed-phase icing that relate to ice accretions within turbofan engines. Experimentally measured air temperature, humidity, total water content, liquid and ice water content, as well as cloud particle size, are compared with model predictions. The model showed good trend agreement with experimentally measured values, but often over-predicted aero-thermodynamic changes. This discrepancy is likely attributed to radial variations that this one-dimensional model does not address. One of the key findings of this work is that greater aero-thermodynamic changes occur when humidity conditions are low. In addition a range of mixed-phase clouds can be achieved by varying only the tunnel humidity conditions, but the range of humidities to generate a mixed-phase cloud becomes smaller when clouds are composed of smaller particles. In general, the model predicted melt fraction well, in particular with clouds composed of larger particle sizes.

Thermal model↗

Particle Size Measurements from the first Fundamentals of Ice Crystal Icing Physics Test in the NASA Propulsion Systems Laboratory

This presentation shows particle measurements by the Artium Technologies, Inc. Phase Doppler Interferometer and High Speed Imaging instruments from the first Fundamental Ice Crystal Icing Physics test conducted in the NASA Propulsion Systems Laboratory. The work focuses on humidity sweeps at a larger and a smaller median volumetric diameter. The particle size distribution, number density, and water content measured by the Phase Doppler Interferometer and High Speed Imaging instruments from the sweeps are presented and compared. The current capability for these two instruments to measure and discriminate ICI conditions is examined.

thermal scaling↗

Altitude Scaling of Thermal Ice Protection Systems in Running Wet Operation

A study into the effects of altitude on an aircraft thermal Ice Protection System (IPS) performance has been conducted by the National Research Council Canada (NRC) in collaboration with the NASA Glenn Icing Branch. The study included tests of an airfoil model, with a heated-air IPS, installed in the NRCs Altitude Icing Wind Tunnel (AIWT) at altitude and ground level conditions.

Running wet mode↗

On High-Order Upwind Methods for Advection

Scheme III (piecewise linear) and V (piecewise parabolic) of Van Leer are shown to yield identical solutions provided the initial conditions are chosen in an appropriate manner. This result is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The result also shows a key connection between the approaches of discontinuous and continuous representations.

high-order methods↗

On High-Order Upwind Methods for Advection

Scheme III (piecewise linear) and V (piecewise parabolic) of Van Leer are shown to yield identical solutions provided the initial conditions are chosen in an appropriate manner. This result is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The result also shows a key connection between the approaches of discontinuous and continuous representations.

high-order methods↗

Analysis and design of numerical schemes for gas dynamics 1: Artificial diffusion, upwind biasing, limiters and their effect on accuracy and multigrid convergence

The theory of non-oscillatory scalar schemes is developed in this paper in terms of the local extremum diminishing (LED) principle that maxima should not increase and minima should not decrease. This principle can be used for multi-dimensional problems on both structured and unstructured meshes, while it is equivalent to the total variation diminishing (TVD) principle for one-dimensional problems. A new formulation of symmetric limited positive (SLIP) schemes is presented, which can be generalized to produce schemes with arbitrary high order of accuracy in regions where the solution contains no extrema, and which can also be implemented on multi-dimensional unstructured meshes. Systems of equations lead to waves traveling with distinct speeds and possibly in opposite directions. Alternative treatments using characteristic splitting and scalar diffusive fluxes are examined, together with modification of the scalar diffusion through the addition of pressure differences to the momentum equations to produce full upwinding in supersonic flow. This convective upwind and split pressure (CUSP) scheme exhibits very rapid convergence in multigrid calculations of transonic flow, and provides excellent shock resolution at very high Mach numbers.

Jameson, Antony↗

The Linear Bicharacteristic Scheme for Computational Electromagnetics

The upwind leapfrog or Linear Bicharacteristic Scheme (LBS) has previously been implemented and demonstrated on electromagnetic wave propagation problems. This paper extends the Linear Bicharacteristic Scheme for computational electromagnetics to treat lossy dielectric and magnetic materials and perfect electrical conductors. This is accomplished by proper implementation of the LBS for homogeneous lossy dielectric and magnetic media, and treatment of perfect electrical conductors (PECs) are shown to follow directly in the limit of high conductivity. Heterogeneous media are treated through implementation of surface boundary conditions and no special extrapolations or interpolations at dielectric material boundaries are required. Results are presented for one-dimensional model problems on both uniform and nonuniform grids, and the FDTD algorithm is chosen as a convenient reference algorithm for comparison. The results demonstrate that the explicit LBS is a dissipation-free, second-order accurate algorithm which uses a smaller stencil than the FDTD algorithm, yet it has approximately one-third the phase velocity error. The LBS is also more accurate on nonuniform grids.

Beggs, John H.↗

A comparison of ENO and TVD schemes

The numerical performance of a second-order upwind-based TVD scheme is compared with that of a uniform second-order ENO scheme on shock capturing. The cases considered include flows with Mach numbers of 2.9, 5.0, and 10.0. For cases with Mach numbers of 5.0 and 10.0, the computed ENO results are inferior to the corresponding TVD results.

Chang, Shih-Hung↗

Prediction of steady and unsteady asymmetric vortical flows around cones

Steady and unsteady, supersonic asymmetric vortical flows and their passive control around circular and noncircular cones are considered in this paper. These problems are formulated by using the unsteady, compressible, single and double, thin-layer. Navier-Stokes equations. The equations are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme, either in a pseudotime stepping or in an accurate-time stepping. An implicit, approximately-factored, central-difference finite-volume scheme has also been used to validate some applications of the upwind scheme. Steady asymmetric vortical flows have been predicted by using random and controlled disturbances for circular and noncircular cones. Unsteady asymmetric vortex-shedding flows have also been predicted, for the first time, using time-accurate solutions, for circular and noncircular cones. Control of flow asymmetry have been demonstrated computationally, for the first time, by inserting a vertical fin the leeward plane of geometric symmetry.

Kandil, Osama A.↗

Application of an Upwind High Resolution Finite-Differencing Scheme and Multigrid Method in Steady-State Incompressible Flow Simulations

The analysis and design of a submarine propulsor requires the ability to predict the characteristics of both laminar and turbulent flows to a higher degree of accuracy. This report presents results of certain benchmark computations based on an upwind, high-resolution, finite-differencing Navier-Stokes solver. The purpose of the computations is to evaluate the ability, the accuracy and the performance of the solver in the simulation of detailed features of viscous flows. Features of interest include flow separation and reattachment, surface pressure and skin friction distributions. Those features are particularly relevant to the propulsor analysis. Test cases with a wide range of Reynolds numbers are selected; therefore, the effects of the convective and the diffusive terms of the solver can be evaluated separately. Test cases include flows over bluff bodies, such as circular cylinders and spheres, at various low Reynolds numbers, flows over a flat plate with and without turbulence effects, and turbulent flows over axisymmetric bodies with and without propulsor effects. Finally, to enhance the iterative solution procedure, a full approximation scheme V-cycle multigrid method is implemented. Preliminary results indicate that the method significantly reduces the computational effort.

NAVIER STOKES EQUATIONS↗

Flux-difference split scheme for turbulent transport equations

This paper describes a high accuracy upwind method to solve the mean compressible flow equations closed with two-equation turbulence models. A flux-difference splitting scheme is developed to compute the flowfield without introducing oscillations into the solution. A high-order accuracy scheme is used to predict both the mean flow variables and the turbulence variables. Results comparing two different two-equation turbulence models are presented for supersonic flow over a flat plate and a compressible, free shear-layer. The non-oscillatory behavior of the upwind scheme is demonstrated on the free shear-layer.

Morrison, J. H.↗

On the Synergy Between Numerics and Subgrid Scale Modeling in LES of Stratified Flows: Grid Convergence of a Stratocumulus-Topped Boundary Layer

The effectiveness of a linear upwinding scalar advection scheme to suppress numerical dispersion errors near sharp inversions in large-eddy simulations of a nocturnal stratocumulus-topped boundary layer is assessed. Linear upwinding is a trade-off between non-dissipative and non-linear positive definite advection schemes. It is shown that linear upwinding does not negatively impact the model's grid convergence properties and a sharp inversion free of numerical artifacts is maintained. Even though mean profiles and turbulence fluxes show good grid convergence characteristics the liquid water amount varies significantly with grid resolution. The entrainment rate is identical for all resolutions and independent of the liquid water amount. For the present stratocumulus case, the impact of cloud-top radiative cooling is negligible and turbulence is largely driven by convection emanating from the surface.

Matheou, Georgios↗

Computation and control of asymmetric vortex flow around circular cones using Navier-Stokes equations

The unsteady, compressible, thin-layer and full Navier-Stokes equations are used to numerically simulate steady and unsteady asymmetric, supersonic, locally conical flows around a 5-deg semiapex angle circular cone. The main computational scheme is the implicit, upwind, flux-difference splitting, finite-volume scheme. Comparison of asymmetric flow solutions using the thin-layer and full Navier-Stokes equations is presented and discussed. The implicit, upwind, flux-vector splitting, finite-volume scheme has also been used to solve for the unsteady asymmetric flow with vortex shedding. The unsteady-flow solution using the flux-vector splitting scheme perfectly agrees with the previously obtained solution using the flux-difference splitting scheme. Passive control of asymmetric flows has been demonstrated and studied using sharp- and round-edged, thick and thin strakes.

Kandil, Osama A.↗

Upwind and symmetric shock-capturing schemes

The development of numerical methods for hyperbolic conservation laws has been a rapidly growing area for the last ten years. Many of the fundamental concepts and state-of-the-art developments can only be found in meeting proceedings or internal reports. This review paper attempts to give an overview and a unified formulation of a class of shock-capturing methods. Special emphasis is on the construction of the basic nonlinear scalar second-order schemes and the methods of extending these nonlinear scalar schemes to nonlinear systems via the extact Riemann solver, approximate Riemann solvers, and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows is discussed. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas dynamics problems.

Yee, H. C.↗