Search NASA⌕ Search

SEARCH · Search NASA

Results for “ice flow equations”

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

A hybrid Navier-Stokes/full-potential method for the prediction of iced wing aerodynamics

A hybrid method for computing compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. The resulting CPU times are less than 50 percent of the required for a full-blown Navier-Stokes analysis. Sample applications of the method to an unswept iced wing at 4 deg and 8 deg angle of attack are presented. Surface pressures are in good agreement with the measurements obtained by Bragg et al. at the University of Illinois.

Mello, O. A. F.↗

Navier-Stokes analysis of airfoils with leading edge ice accretions

A numerical analysis of the flowfield characteristics and the performance degradation of an airfoil with leading edge ice accretions was performed. The important fluid dynamic processes were identified and calculated. Among these were the leading edge separation bubble at low angles of attack, complete separation on the low pressure surface resulting in premature shell, drag rise due to the ice shape, and the effects of angle of attack on the separated flow field. Comparisons to experimental results were conducted to confirm these calculations. A computer code which solves the Navier-Stokes equations in two dimensions, ARC2D, was used to perform the calculations. A Modified Mixing Length turbulence model was developed to produce grids for several ice shape and airfoil combinations. Results indicate that the ability to predict overall performance characteristics, such as lift and drag, at low angles of attack is excellent. Transition location is important for accurately determining separation bubble shape. Details of the flowfield in and downstream of the separated regions requires some modifications. Calculations for the stalled airfoil indicate periodic shedding of vorticity that was generated aft of the ice accretion. Time averaged pressure values produce results which compare favorably with experimental information. A turbulence model which accounts for the history effects in the flow may be justified.

Potapczuk, Mark G.↗

LEWICE/E: An Euler based ice accretion code

A new version of the LEWICE ice accretion computer code was developed which calculates the ice growth on two dimensional surfaces, incorporating the effects of compressibility through the solution of the Euler equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, and plot results. This code increases the applicability of ice accretion predictions by allowing calculations at higher Mach numbers. The new elements of the code are described. Calculated results are compared to experiment for several cases, including a LEWICE example case and a thin airfoil section at a Mach number of 0.58.

Potapczuk, Mark G.↗

LEWICE/E - An Euler based ice accretion code

A new version of the LEWICE ice accretion computer code was developed which calculates the ice growth on two dimensional surfaces, incorporating the effects of compressibility through the solution of the Euler equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, and plot results. This code increases the applicability of ice accretion predictions by allowing calculations at higher Mach numbers. The new elements of the code are described. Calculated results are compared to experiment for several cases, including a LEWICE example case and a thin airfoil section at a Mach number of 0.58.

Potapczuk, Mark G.↗

Superplastic Deformation of ice: Experimental Observations

Creep experiments on fine-grained ice reveal the existence of three creep regimes: (1) a dislocation creep regime; (2) a superplastic flow regime in which grain boundary sliding is an important deformation process; and (3) a basal slip creep regime in which the strain rate is limited by basal slip. Dislocation creep in ice is likely climb-limited, is characterized by a stress exponent of 4.0, and is independent of grain size. Superplastic flow is characterized by a stress exponent of 1.8 and depends inversely on grain size to the 1.4 power. Basal slip limited creep is characterized by a stress exponent of 2.4 and is independent of grain size. A fourth creep mechanism, diffusional flow, which usually occurs at very low stresses, is inaccessible at practical laboratory strain rates even for our finest grain sizes of approximately 3 micrometers. A constitutive equation based on these experimental results that includes flow laws for these four creep mechanisms is described. This equation is in excellent agreement with published laboratory creep data for coarse-grained samples at high temperatures. Superplastic flow of ice is the rate-limiting creep mechanism over a wide range of temperatures and grain sizes at stresses less than or equal to 0.1 MPa, conditions which overlap those occurring in glaciers, ice sheets, and icy planetary interiors.

Goldsby, D. L.↗

Ice Accretion and Performance Degradation Calculations with LEWICE/NS

The LEWICE ice accretion computer code has been extended to include the solution of the two-dimensional Navier-Stokes equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, calculate aeroperformance changes, and plot results. The new elements of the code are described. Calculated results are compared to experiment for several cases, including both ice shape and drag rise.

Potapczuk, Mark G.↗

Ice accretion and performance degradation calculations with LEWICE/NS

The LEWICE ice accretion computer code has been extended to include the solution of the two-dimensional Navier-Stokes equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, calculate aeroperformance changes, and plot results. The new elements of the code are described. Calculated results are compared to experiment for several cases, including both ice shape and drag rise.

Potapczuk, Mark G.↗

RANS and DES Computations for a Wing with Ice Accretion

A computational investigation was performed to assess the effectiveness of Detached Eddy Simulation (DES) as a tool for predicting icing effects. The AVUS code, a public domain flow solver, was employed to compute solutions for an iced wing configuration using DES and steady Reynolds Averaged Navier-Stokes (RANS) equation methodologies. The wing section considered here was a business jet airfoil (GLC305) with a 22.5-minute glaze ice accretion (944-ice shape). The section was extruded to form a rectangular planform. The model was mounted between two walls so no tip effects were considered. The numerical results were validated by comparison with experimental data for the same configuration. The time-averaged DES computations showed some improvement in lift and drag results near stall when compared to steady RANS results. However, comparisons of the flow field details did not show the level of agreement suggested by the integrated quantities. More validation is needed to determine what role DES can play as part of an overall icing effects prediction strategy.

Mogili, P.↗

Xanthos-Lake Model Source Code

This repository contains the source code for Xanthos-Lake, a lake-modeling extension of the Xanthos framework that introduces a coupled lake component comprising the Xanthos-Lake Snow and Ice Model (xLSIM) and the Xanthos-Lake Water Balance Model (xLWBM). xLSIM is a basin-aware machine-learning model for lake snow, ice, and thermal conditions. It predicts monthly lake ice thickness, snow depth, snow-cover fraction, mixing-layer temperature, and lake ice fraction from meteorological forcing and lake surface-area information. It uses sequence-based deep-learning architectures, including Transformer and hybrid Long Short-Term Memory–Transformer (LSTM–Transformer) models, together with seasonal encoding, multi-lake learning, physical masking, and basin-level cryospheric and non-cryospheric classification. The training workflow uses Ray for scalable execution and includes optional Ray Tune hyperparameter optimization. Model predictions, observations, diagnostics, and feature-importance outputs are written in NetCDF. xLWBM is the water-balance component of the new lake framework. It simulates monthly lake storage, surface area, evaporation, inflow, outflow, and lake–groundwater exchange. It combines physical water-balance equations with calibrated bathymetric relationships, weir-based outlet flow, modified Penman open-water evaporation, groundwater head relaxation, Penman–Monteith snow and ice sublimation, and snow, ice, and thermal conditions supplied by xLSIM. The model calibrates lake parameters against satellite-derived surface-area data, using evaporation-based calibration where surface-area data are unavailable, and supports small, medium, and large lake classes. For large lakes, xLWBM is integrated with the managed-routing workflow so that lake storage and outflow interact directly with downstream river routing and reservoir operations. Together, xLSIM and xLWBM provide Xanthos with a coupled lake-modeling capability. xLSIM supplies the snow, ice, and thermal conditions that affect lake evaporation and snow- and ice-related water exchanges, while xLWBM translates those conditions into dynamic lake storage, surface area, evaporation, and discharge. In return, xLWBM supplies evolving lake surface area to xLSIM. This coupling enables Xanthos to represent lakes as active hydrologic components within basin-scale water-availability and routing simulations.

Machine Learning↗

Effects of icing on the aerodynamic performance of high lift airfoils

A 2D compressible Navier-Stokes solver capable of analyzing multi-element airfoils is described. The flow field is divided into multiple zones. In each zone, the governing equations are solved using an implicit finite difference scheme. The flow solver is validated through a study of the aerodynamic characteristics of a GA(W)-1 configuration, for which good quality measured surface pressure data and load data are available. The solver is next applied to a study of the effects of icing on an iced 5-element airfoil configuration, experimentally studied at NASA Lewis Research Center. It is demonstrated that the formation of ice over the leading edge slat and the main airfoil can lead to significant flow separation, and a significant loss in lift, compared to clean configurations.

Sankar, L. N.↗

Numerical investigation of the effects of icing on fixed and rotary wing aircraft

A 2-D multi-element airfoil code was modified to study the effects of icing on the aerodynamic characteristics of high lift systems. In each zone of the flow field, the solver numerically integrates the 2-D compressible Navier-Stokes equations using a time marching scheme. The surface pressure distribution is generated over a GAW 130 airfoil/flap combination for a flap setting of 25 degrees, and an angle of attack equal to 5 degrees, at a freestream Mach number equal to 0.3. A series of calculations were performed to determine the effects of small scale ice build up on the high lift characteristics of this arifoil/flap combination. The appendix summarizes this progress. Joint studies on correlation of a 3-D iced wing code with experimental data reviewed new measured laser Doppler velocimeter data in the separated region behind the leading edge ice shape. A version of the iced wing analysis using the Roe scheme was developed to evaluate the poor correlation between the computed and measured velocities in the separated region. Work on the extension of the wing-alone analysis to wing body configuration began with modifications to the 3-D iced wing analyses to accept externally generated grids and multi-block grids.

Sankar, L. N.↗

Troughs on Martian Ice Sheets: Analysis of Their Closure and Mass Balance

At the Copenhagen workshop on Martian polar processes, Ralf Greve commented that the flow regime surrounding scarps and troughs of the Martian polar ice sheets cannot be modeled using traditional "plan view" ice-sheet models. Such models are inadequate because they typically use reduced equations that embody certain simplifications applicable only to terrestrial ice sheets where the upper ice sheet surface is smooth. In response to this suggestion, we have constructed a 2-dimensional, time dependent "side view" (two spatial dimensions: one horizontal, one vertical) model of scarp closure that is designed to overcome the difficulties described by Greve. The purpose of the model is to evaluate the scales of stress variation and styles of flow closure so as to estimate errors that may be encountered by "plan view" models. We show that there may be avenues whereby the complications associated with scarp closure can be overcome in "plan view" models through appropriate parameterizations of 3-dimensional effects. Following this, we apply the flow model to simulate the evolution of a typical scarp on the North Polar Cap of Mars. Our simulations investigate: (a) the role of "radiation trapping" (see our companion abstract) in creating and maintaining "spiral-like" scarps on the ice sheet, (b) the consequences of different flowlaws and ice compositions on scarp evolution and, in particular, scarp age, and (c) the role of dust and debris in scarp evolution.

Fountain, A.↗

Numerical study of the effects of icing on viscous flow over wings

An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved using a diagonal form of an alternating-direction implicit (ADI) approximate factorization procedure. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic-based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are less than 60 percent of that required for a full-blown Navier-Stokes analysis for steady flow applications and about 60 percent of the Navier-Stokes CPU times for unsteady flows in non-vector processing machines. Applications of the method are presented for a rectangular NACA 0012 wing in low subsonic steady flow at moderate and high angles of attack, and for an F-5 wing in steady and unsteady subsonic and transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to Navier-Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier-Stokes equations can be retained with a significant savings in computational time.

Sankar, L. N.↗

Progress in development of a Navier-Stokes solver for evaluation of iced airfoil performance

A method is being developed for evaluation of the flow field behavior about an airfoil with significant ice accretion on the leading edge. The computer code, being evaluated for this purpose, solves the Navier-Stokes equations in a body-fitted curvilinear coordinate system. This requires the use of a grid generation code to transform the x-y coordinates of the physical space into xi-eta coordinates of the computational space. Evaluation of the suitability of these two codes for predicting iced airfoil performance is presently being carried out in anticipation of use in an overall icing analysis effort. Results of this evaluation to date indicate good correlation with known information on clean airfoils. Preliminary results for rime and glaze, iced airfoil shapes are also presented.

Potapczuk, M. G.↗

Structure and evolution of Uranus and Neptune

Three-layer interior models of Uranus and Neptune with central rocky cores, mantles of water, methane, and ammonia (the 'ices'), and outer envelopes primarily composed of hydrogen and helium are presented. The models incorporate a new H2O equation of state based on experimental data which is considerably 'softer' than previous H2O equations of state. Corrections for interior temperatures approximately 5000 K are included in the models, and the thermal evolution of both planets is investigated using recent heat flow measurements. It is found that the evolutionary considerations are consistent with gravitational field data in supporting models with approximately solar abundances of 'ice' and 'rock'. Evolutionary considerations indicate that initial temperatures and luminosities for Uranus and Neptune were not substantially higher than the present value. Both planets apparently have relatively small approximately 1-2 earth masses) hydrogen-helium envelopes, with Neptune's envelope smaller than Uranus'. A monotonic trend is evident among the Jovian planets: all have central rock-ice cores of approximately 15 earth masses, but with hydrogen-helium envelopes which decrease in mass from Jupiter to Saturn to Uranus to Neptune.

Hubbard, W. B.↗

Influence of the vapor flux on temperature, density, and abundance distributions in a multicomponent, porous, icy body

We calculated the vapor flux of the icy components in the surface layer of a porous, short-period, Jupiter-class comet, in order to investigate the relationship of the observed relative molecular abundances in the coma with those in the nucleus. The model assumes a body containing one major ice component (H20) and up to three minor components of higher volatility (e.g., CO, CO2, CH3OH). The body's porous structure is modeled as a bundle of tubes with a given tortuosity and initially a constant pore diameter. The mass and energy equations for the different volatiles are solved simultaneously under appropiate boundary conditions. Heat is conducted by the matrix and carried by the vapors. The one-dimensional model includes radially inward and outward flowing vapor within the body, complete depletion of less volatile ices in outer layers, and recondensation of vapor in deeper, coller layers. As a result, we obtain the temperature and abundance distribution in the nucleus and the gas flux into the interior and into the coma for each of the volatiles at various positions in the orbit. The ratio of the gas flux of minor volatiles to that of H2) in the coma varies by several orders of magnitude throughout the orbit. Thus, the relative abundances of species observed in the coma are in most cases not the same as those in the nucleus. Results also indicate that it will be impossible to determine the relative abundances of ices more volatile than water from samples taken a few meters below the surface during a comet rendezvous mission. We made calculations for a wide range of different parameters, such as porosity, pore radius, and thermal conductivity of the matrix. To introduce the model we present typical results for a dust-free comet.

Benkhoff, J.↗

HdsaLib

HdsaLib performs hyper-differential sensitivity analysis in optimization problems that are constrained by partial differential equation problems. The software aids researchers in understanding the impact of uncertainties on their optimization models. The tool has been applied for use in ice sheet modeling, fluid flow modeling, atmospheric tracer transport modeling, and subsurface flow modeling. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

van Bloemen Waanders, Bart [Sandia National Lab. (↗

Mathematical modeling of ice accretion on airfoils

The progress toward development of a computer model suitable for predicting icing behavior on airfoils over a wide range of environmental conditions and airfoils shapes is reported. The LEWICE program was formulated to solve a set of equations which describe the physical processes which occur during accretion of ice on an airfoil, including heat transfer in a time dependent mode, with the restriction that the flow must be describable by a two-dimensional flow code. Input data comprises the cloud liquid water content, mean droplet diameter, ambient air temperature, air velocity, and relative humidity. A potential flowfield around the airfoil is calculated, along with the droplet trajectories within the flowfield, followed by local values of water droplet collection efficiency at the impact points. Both glaze and rime ice conditions are reproduced, and comparisons with test results on icing of circular cylinders showed good agreement with the physical situation.

Macarthur, C. D.↗