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 73 records · Page 4

Low-Reynolds Number Aerodynamics of an 8.9 Percent Scale Semispan Swept Wing for Assessment of Icing Effects

Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity and Reynolds and Mach number effects on the iced-wing aerodynamics is needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing and computational flow simulations were carried out for an 8.9 percent-scale semispan wing based upon the Common Research Model airplane configuration. The wind-tunnel testing was conducted at the Wichita State University 7 by 10 ft Beech wind tunnel from Reynolds numbers of 0.8×10(exp 6) to 2.4×10(exp 6) and corresponding Mach numbers of 0.09 to 0.27. This paper presents the results of initial studies investigating the model mounting configuration, clean-wing aerodynamics and effects of artificial ice roughness. Four different model mounting configurations were considered and a circular splitter plate combined with a streamlined shroud was selected as the baseline geometry for the remainder of the experiments and computational simulations. A detailed study of the clean-wing aerodynamics and stall characteristics was made. In all cases, the flow over the outboard sections of the wing separated as the wing stalled with the inboard sections near the root maintaining attached flow. Computational flow simulations were carried out with the ONERA elsA software that solves the compressible, threedimensional RANS equations. The computations were carried out in either fully turbulent mode or with natural transition. Better agreement between the experimental and computational results was obtained when considering computations with free transition compared to turbulent solutions. These results indicate that experimental evolution of the clean wing performance coefficients were due to the effect of three-dimensional transition location and that this must be taken into account for future data analysis. This research also confirmed that artificial ice roughness created with rapid-prototype manufacturing methods can generate aerodynamic performance effects comparable to grit roughness of equivalent size when proper care is exercised in design and installation. The conclusions of this combined experimental and computational study contributed directly to the successful implementation of follow-on test campaigns with numerous artificial ice-shape configurations for this 8.9 percent scale model.

Aerodynamics↗

Low-Reynolds Number Aerodynamics of an 8.9 Percent Scale Semispan Swept Wing for Assessment of Icing Effects

Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity and Reynolds and Mach number effects on the iced-wing aerodynamics is needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing and computational flow simulations were carried out for an 8.9%-scale semispan wing based upon the Common Research Model airplane configuration. The wind-tunnel testing was conducted at the Wichita State University 7 ft x 10 ft Beech wind tunnel from Reynolds numbers of 0.8×10(exp 6) to 2.4×10(exp 6) and corresponding Mach numbers of 0.09 to 0.27. This paper presents the results of initial studies investigating the model mounting configuration, clean-wing aerodynamics and effects of artificial ice roughness. Four different model mounting configurations were considered and a circular splitter plate combined with a streamlined shroud was selected as the baseline geometry for the remainder of the experiments and computational simulations. A detailed study of the clean-wing aerodynamics and stall characteristics was made. In all cases, the flow over the outboard sections of the wing separated as the wing stalled with the inboard sections near the root maintaining attached flow. Computational flow simulations were carried out with the ONERA elsA software that solves the compressible, three-dimensional RANS equations. The computations were carried out in either fully turbulent mode or with natural transition. Better agreement between the experimental and computational results was obtained when considering computations with free transition compared to turbulent solutions. These results indicate that experimental evolution of the clean wing performance coefficients were due to the effect of three-dimensional transition location and that this must be taken into account for future data analysis. This research also confirmed that artificial ice roughness created with rapid-prototype manufacturing methods can generate aerodynamic performance effects comparable to grit roughness of equivalent size when proper care is exercised in design and installation. The conclusions of this combined experimental and computational study contributed directly to the successful implementation of follow-on test campaigns with numerous artificial ice-shape configurations for this 8.9% scale model.

Aerodynamics↗

The Comet Giacobini-Zinner magnetotail: Axial stresses and inferred near-nucleus properties

Utilizing the electron and magnetic field data from the ICE tail traversal of comet Giacobini-Zinner along with the MHD equations, a steady state, stress balance model of the cometary magnetotail was developed, and used to infer important but unmeasured ion properties within the magnetotail at ICE and upstream at the average point along each streamline where cometary ions are picked-up. The derived tailward ion flow speed at ICE is quite constant at approx. -20 to -30 km/sec across the entire tail. The flow velocity, ion temperature, density, and ion source rates upstream from the lobes (current sheet) at the average pick-up locations are approx. -75 km/sec (approx. -12), approx. 4 million K (approx. 100,000), approx. 20 cc (approx. 400), and approx. 15 cu cm/sec. Gradients in the plasma properties between the two regions are quite strong. Implications of inferred plasma properties for the near-nucleus region and for cometary magnetotail formation are examined.

Mccomas, D. J.↗

Dynamic Gas Flow Effects on the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a dynamic version of Paschen's Law that takes into account the flow of ambient gas past aerospace vehicle surfaces. However, the classic Paschen's Law does not take into account the flow of gas of an aerospace vehicle, whose surfaces may be triboelectrically charged by dust or ice crystal impingement, traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance by the electric field between the electrodes is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised Paschen equation must be a function of the mean velocity, v(sub xm), of the ambient gas and reduces to the classical version of Paschen's law when the gas mean velocity, v(sub xm) = 0. New formulations of Paschen's Law, taking into account Mach number and dynamic pressure, derived by the authors, will be discussed. These equations will be evaluated by wind tunnel experimentation later this year. Based on the results of this work, it is hoped that the safety of aerospace vehicles will be enhanced with a redefinition of electrostatic launch commit criteria. It is also possible that new products, such as new anti-static coatings, may be formulated from this data.

Paschens Law↗

Dynamic Gas Flow Effects on the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a dynamic version of Paschen's Law that takes into account the flow of ambient gas past aerospace vehicle surfaces. However, the classic Paschen's Law does not take into account the flow of gas of an aerospace vehicle, whose surfaces may be triboelectrically charged by dust or ice crystal impingement, traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance by the electric field between the electrodes is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised Paschen equation must be a function of the mean velocity, v(sub xm), of the ambient gas and reduces to the classical version of Paschen's law when the gas mean velocity, v(sub xm) = 0. New formulations of Paschen's Law, taking into account Mach number and dynamic pressure, derived by the authors, will be discussed. These equations will be evaluated by wind tunnel experimentation later this year. Based on the results of this work, it is hoped that the safety of aerospace vehicles will be enhanced with a redefinition of electrostatic launch commit criteria. It is also possible that new products, such as new anti-static coatings, may be formulated from this data.

Electrostatic Discharge↗

Modeling of Dense Water Production and Salt Transport from Alaskan Coastal Polynyas

The main significance of this paper is that a realistic, three-dimensional, high-resolution primitive equation model has been developed to study the effects of dense water formation in Arctic coastal polynyas. The model includes realistic ambient stratification, realistic bottom topography, and is forced by time-variant surface heat flux, surface salt flux, and time-dependent coastal flow. The salt and heat fluxes, and the surface ice drift, are derived from satellite observations (SSM/I and NSCAT sensors). The model is used to study the stratification, salt transport, and circulation in the vicinity of Barrow Canyon during the 1996/97 winter season. The coastal flow (Alaska coastal current), which is an extension of the Bering Sea throughflow, is formulated in the model using the wind-transport regression. The results show that for the 1996/97 winter the northeastward coastal current exports 13% to 26% of the salt produced by coastal polynyas upstream of Barrow Canyon in 20 to 30 days. The salt export occurs more rapidly during less persistent polynyas. The inclusion of ice-water stress in the model makes the coastal current slightly weaker and much wider due to the combined effects of surface drag and offshore Ekman transport.

Signorini, Sergio R.↗

Low temperature simulation of subliming boundary layer flow in Jupiter atmosphere

A low-temperature approximate simulation for the sublimation of a graphite heat shield under Jovian entry conditions is studied. A set of algebraic equations is derived to approximate the governing equation and boundary conditions, based on order-of-magnitude analysis. Characteristic quantities such as the wall temperature and the subliming velocity are predicted. Similarity parameters that are needed to simulate the most dominant phenomena of the Jovian entry flow are also given. An approximate simulation of the sublimation of the graphite heat shield is performed with an air-dry-ice model. The simulation with the air-dry-ice model may be carried out experimentally at a lower temperature of 3000 to 6000 K instead of the entry temperature of 14,000 K. The rate of graphite sublimation predicted by the present algebraic approximation agrees to the order of magnitude with extrapolated data. The limitations of the simulation method and its utility are discussed.

Chen, C. J.↗

Modeling Climate Change in the Absence of Climate Change Data. Editorial Comment

Practitioners of climate change prediction base many of their future climate scenarios on General Circulation Models (GCM's), each model with differing assumptions and parameter requirements. For representing the atmosphere, GCM's typically contain equations for calculating motion of particles, thermodynamics and radiation, and continuity of water vapor. Hydrology and heat balance are usually included for continents, and sea ice and heat balance are included for oceans. The current issue of this journal contains a paper by Van Blarcum et al. (1995) that predicts runoff from nine high-latitude rivers under a doubled CO2 atmosphere. The paper is important since river flow is an indicator variable for climate change. The authors show that precipitation will increase under the imposed perturbations and that owing to higher temperatures earlier in the year that cause the snow pack to melt sooner, runoff will also increase. They base their simulations on output from a GCM coupled with an interesting water routing scheme they have devised. Climate change models have been linked to other models to predict deforestation.

Skiles, J. W.↗

Numerical studies of unsteady two dimensional subsonic flows using the ICE method

A numerical program was developed to compute transient compressible and incompressible laminar flows in two dimensions with multicomponent mixing and chemical reaction. The algorithm used the Los Alamos Scientific Laboratory ICE (Implicit Continuous-Fluid Eulerian) method as its base. The program can compute both high and low speed compressible flows. The numerical program incorporating the stabilization techniques was quite successful in treating both old and new problems. Detailed calculations of coaxial flow very close to the entry plane were possible. The program treated complex flows such as the formation and downstream growth of a recirculation cell. An implicit solution of the species equation predicted mixing and reaction rates which compared favorably with the literature.

Wieber, P. R.↗

A Mars Pathfinder landing on a recently drained ephemeral sea: Cerberus Plains, 6 deg N, 188 deg W

Along a 500 km-wide belt extending between 202 deg and 180 deg W and lying astride the martian equator, moderately low-albedo, uncratered smooth plains exhibit low thermal inertia and potentially favorable conditions for the preservation of near-surface ice. The Cerberus Plains occupy a topographic trough as much as 2 km below the planetary datum, and the denser atmosphere at these altitudes would also favor long residence times for near-surface ice once emplaced. The plains have previously been interpreted as the result of young (late Amazonian) low viscosity lava flows or similarly youthful fluvial deposition. However, the plains are also included in maps of possibly extensive martian paleoseas or paleolakes. Ice emplaced as such seas dissipated could still be preserved under thin (a few tens of centimeters) sedimentary cover. In any case, and if a sea once existed, aqueous-born interstitial cementation, probably including hydrated iron oxides and sulfate minerals, would have been favored and is now susceptible to investigation by the Pathfinder alpha proton x-ray spectrometer and multispectral imager.

Brakenridge, G. Robert↗

Analytical determination of propeller performance degradation due to ice accretion

A computer code has been developed which is capable of computing propeller performance for clean, glaze, or rime iced propeller configurations, thereby providing a mechanism for determining the degree of performance degradation which results from a given icing encounter. The inviscid, incompressible flow field at each specified propeller radial location is first computed using the Theodorsen transformation method of conformal mapping. A droplet trajectory computation then calculates droplet impingement points and airfoil collection efficiency for each radial location, at which point several user-selectable empirical correlations are available for determining the aerodynamic penalities which arise due to the ice accretion. Propeller performance is finally computed using strip analysis for either the clean or iced propeller. In the iced mode, the differential thrust and torque coefficient equations are modified by the drag and lift coefficient increments due to ice to obtain the appropriate iced values. Comparison with available experimental propeller icing data shows good agreement in several cases. The code's capability to properly predict iced thrust coefficient, power coefficient, and propeller efficiency is shown to be dependent on the choice of empirical correlation employed as well as proper specification of radial icing extent.

Miller, T. L.↗

The Water Table Model (WTM) (v2.0.1): coupled groundwater and dynamic lake modelling

Abstract. Ice-free land comprises 26 % of the Earth's surface and holds liquid water that delineates ecosystems, affects global geochemical cycling, and modulates sea levels. However, we currently lack the capacity to simulate and predict these terrestrial water changes across the full range of relevant spatial (watershed to global) and temporal (monthly to millennial) scales. To address this knowledge gap, we present the Water Table Model (WTM), which integrates coupled components to compute dynamic lake and groundwater levels. The groundwater component solves the 2D horizontal groundwater flow equation using non-linear equation solvers from the C++ PETSc (Portable, Extensible Toolkit for Scientific Computation) library. The dynamic lake component makes use of the Fill–Spill–Merge (FSM) algorithm to move surface water into lakes, where it may evaporate or affect groundwater flow. In a proof-of-concept application, we demonstrate the continental-scale capabilities of the WTM by simulating the steady-state climate-driven water table for the present day and the Last Glacial Maximum (LGM; 21 000 calendar years before present) across the North American continent. During the LGM, North America stored an additional 14.98 cm of sea-level equivalent (SLE) in lakes and groundwater compared to the climate-driven present-day scenario. We compare the present-day result to other simulations and real-world data. Open-source code for the WTM is available on GitHub and Zenodo.

Callaghan, Kerry L. (ORCID:0000000226740838)↗

An Empirical Method Permitting Rapid Determination of the Area, Rate and Distribution of Water-Drop Impingement on an Airfoil of Arbitrary Section at Subsonic Speeds

An empirical method for the determination of the area, rate, and distribution of water-drop impingement on airfoils of arbitrary section is presented. The procedure represents an initial step toward the development of a method which is generally applicable in the design of thermal ice-prevention equipment for airplane wing and tail surfaces. Results given by the proposed empirical method are expected to be sufficiently accurate for the purpose of heated-wing design, and can be obtained from a few numerical computations once the velocity distribution over the airfoil has been determined. The empirical method presented for incompressible flow is based on results of extensive water-drop. trajectory computations for five airfoil cases which consisted of 15-percent-thick airfoils encompassing a moderate lift-coefficient range. The differential equations pertaining to the paths of the drops were solved by a differential analyzer. The method developed for incompressible flow is extended to the calculation of area and rate of impingement on straight wings in subsonic compressible flow to indicate the probable effects of compressibility for airfoils at low subsonic Mach numbers.

Bergrun, N. R.↗

A model of Comet P/Giacobini-Zinner

A computer model of Comet P/Giacobini-Zinner is presented which contains photoprocesses, gas-phase chemical kinetics, energy balance, multifluid hydrodynamics with a transition to free molecular flow, and solar wind interaction. Recently the physics for electrons in the model has been improved by including electron impact ionization and dissociation and separately accounting for electron energetics. Electron heating and cooling mechanisms include photoprocesses, recombination processes, inelastic and elastic collisions with heavy molecules, and expansion cooling. The model incorporates an internally consistent interaction of the solar wind with the coma gas using the axisymmetric ideal fluid dynamic equations. The nuclear size and composition have been chosen to make the calculations relevant to the September 11, 1985 International Cometary Explorer (ICE) encounter with Comet P/Giacobini-Zinner. Model profiles of the temperature, velocity, and number density of the electrons are in good agreement with measurements along ICE's trajectory. The results indicate the probe passed through a region of the coma at the onset of the plasma tail.

Boice, D. C.↗

Revision of Paschen's Law Relating to the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a version of Paschens law that takes into account the flow of ambient gas past electrode surfaces. Paschens law does not consider the flow of gas past an aerospace vehicle whose surfaces may be triboelectrically charged by dust or ice crystal impingement while traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance between electrode surfaces is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised theoretical model must be a function of the mean velocity vxm of the ambient gas and reduce to Paschens law when the mean velocity is zero. A new theoretical formulation of Paschens law, taking into account the Mach number and compressible dynamic pressure, derived by the authors, will be discussed. This equation has been evaluated by wind tunnel experimentation. Initial data of the baseline wind tunnel experiments show results consistent with the hypothesis. This work may enhance the safety of aerospace vehicles through a redefinition of electrostatic launch commit criteria. It is also possible for new products, such as antistatic coatings, to be formulated based on this data.

discharge↗

Revision of Paschen's Law Relating to the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a version of Paschen's law that takes into account the flow of ambient gas past electrode surfaces. Paschen's law does not consider the flow of gas past an aerospace vehicle whose surfaces may be triboelectrically charged by dust or ice crystal impingement while traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance between electrode surfaces is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised theoretical model must be a function of the mean velocity vxm of the ambient gas and reduce to Paschen's law when the mean velocity is zero. A new theoretical formulation of Paschen's law, taking into account the Mach number and compressible dynamic pressure, derived by the authors, will be discussed. This equation has been evaluated by wind tunnel experimentation. Initial data of the baseline wind tunnel experiments show results consistent with the hypothesis. This work may enhance the safety of aerospace vehicles through a redefinition of electrostatic launch commit criteria. It is also possible for new products, such as antistatic coatings, to be formulated based on this data.

electrostatic↗

Dynamic Gas Flow Effects on the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a version of Paschen's Law that takes into account the flow of ambient gas past electrode surfaces. Paschen's Law does not consider the flow of gas past an aerospace vehicle, whose surfaces may be triboelectrically charged by dust or ice crystal impingement while traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance between electrode surfaces is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised theoretical model must be a function of the mean velocity, v (sub xm), of the ambient gas and reduce to Paschen's law when the gas mean velocity, v (sub xm) equals 0. A new theoretical formulation of Paschen's Law, taking into account the Mach number and dynamic pressure, derived by the authors, will be discussed. This equation was evaluated by wind tunnel experimentation whose results were consistent with the model hypothesis.

Electrostatic↗

A study of the accuracy of neutrally buoyant bubbles used as flow tracers in air

Research has been performed to determine the accuracy of neutrally buoyant and near neutrally buoyant bubbles used as flow tracers in air. Theoretical, computational, and experimental results are presented to evaluate the dynamics of bubble trajectories and factors affecting their ability to trace flow-field streamlines. The equation of motion for a single bubble was obtained and evaluated using a computational scheme to determine the factors which affect a bubble's trajectory. A two-dimensional experiment was also conducted to experimentally determine bubble trajectories in the stagnation region of NACA 0012 airfoil at 0 deg angle of attack using a commercially available helium bubble generation system. Physical properties of the experimental bubble trajectories were estimated using the computational scheme. These properties included the density ratio and diameter of the individual bubbles. the helium bubble system was then used to visualize and document the flow field about a 30 deg swept semispan wing with simulated glaze ice. Results were compared to Navier-Stokes calculations and surface oil flow visualization. The theoretical and computational analysis have shown that neutrally buoyant bubbles will trace even the most complex flow patterns. Experimental analysis revealed that the use of bubbles to trace flow patterns should be limited to qualitative measurements unless care is taken to ensure neutral buoyancy. This is due to the difficulty in the production of neutrally buoyant bubbles.

Kerho, Michael F.↗