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.
An Investigation Utilizing an Electrical Analogue of Cyclic Deicing of Hollow Steel Propellers with Internal Electric Heaters
A study has been made of the heating requirements for the cyclic de?icing of hollow steel propellers fitted with two types of internal electric heaters. Solutions to the transient?teat?flow equations depicting the cyclic de?icing of propellers were obtained by use of an electrical analogy. The study showed the impracticability of using an internal tubular heater and illustrated the advantages of employing an internal shoe?type heater, which distributes the heat more evenly to the blade surface. The importance of minimizing the thermal inertia of the system was demonstrated, and the magnitude of reductions in the total energy requirement made possible through reductions in the heating period was indicated.
A Method for Calculating the Heat Required for Windshield Thermal Ice Prevention Based on Extensive Flight Tests in Natural Icing Conditions
An equation is presented for calculating the heat flow required from the surface of an internally heated windshield in order to prevent the formation of ice accretions during flight in specified icing conditions. To ascertain the validity of the equation, comparison is made between calculated values of the heat required and measured values obtained for test windshields in actual flights in icing conditions. The test windshields were internally heated and provided data applicable to two common types of windshield configurations; namely the V-type and the type installed flush with the fuselage contours. These windshields were installed on a twin-engine cargo airplane and the icing flights were conducted over a large area of the United States during the winters of 1945-46 and 1946-47. In addition to the internally heated windshield investigation, some test data were obtained for a windshield ice-prevention system in which heated air was discharged into the windshield boundary layer. The general conclusions resulting from this investigation are as follows: 1) The amount of heat required for the prevention of ice accretions on both flush- and V-type windshields during flight in specified icing conditions can be calculated with a degree of accuracy suitable for design purposes. 2) A heat flow of 2000 to 2500 Btu per hour per square foot is required for complete and continuous protection of a V-type windshield in fight at speeds up to 300 miles per hour in a moderate cumulus icing condition. For the same degree of protection and the same speed range, a value of 1000 Btu per hour per square foot suffices in a moderate stratus icing condition. 3) A heat supply of 1000 Btu per hour per square foot is adequate for a flush windshield located well aft of the fuselage stagnation region, at speeds up to 300 miles per hour, for flight in both stratus and moderate cumulus icing conditions. 4) The external air discharge system of windshield thermal ice prevention is thermally inefficient and requires a heat supply approximately 20 times that required for an internal system having the same performance.
A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean
An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.
Rheology of Diabase: Implications for Tectonics on Venus and Mars
Two important goals of our experimental investigation of the rheological behavior of diabase rocks were: (1) to determine flow laws describing their creep behavior over wide ranges of temperature, stress and strain rate and (2) to develop an understanding of the physical mechanisms by which these rocks flow under laboratory conditions. With this basis, a primary objective then was to construct constitutive equations that can be used to extrapolate from laboratory to planetary conditions. We specifically studied the rheological properties of both natural rock samples and synthetic aggregates. The former provided constraints for geologic systems, while the latter defined the relative contributions of the constituent mineral phases and avoided the influence of glass/melt found in natural samples. In addition, partially molten samples of crustal rock composition were deformed in shear to large strains (greater than 200%) important in crustal environments. The results of this research yielded essential rheological properties essential for models of crustal deformation on terrestrial planets, specifically Venus and Mars, as well as on the geodynamical evolution of these planets. Over the past three years, we also completed our investigation of the creep behavior of water ice with applications to the glaciers, ice sheets and icy satellites. Constitutive equations were determined that describe flow over a wide ranged of stress, strain rate, grain size and temperature. In the case of ice, three creep regimes were delineate. Extrapolation demonstrates that dislocation glide and grain boundary sliding processes dominate flow in ice I under planetary conditions and that diffusion creep is not an important deformation mechanism either in the laboratory or on icy satellites. These results have already been incorporated by other investigators into models describing, for example, the thickness and stability of the ice shell on Europa and to unravel long-standing discrepancies between field observations on glaciers and laboratory results.
A Computational Icing Effects Study for a Three-Dimensional Wing: Comparison with Experimental Data and Investigation of Spanwise Variation
Steady-state solutions of the Reynolds-averaged Navier-Stokes (RANS) equations were computed using the Colbalt flow solver for a constant-section, rectangular wing based on an extruded two-dimensional glaze ice shape. The one equation Spalart-Allmaras turbulence model was used. The results were compared with data obtained from a recent wind tunnel test. Computed results indicate that the steady RANS solutions do not accurately capture the recirculating region downstream of the ice accretion, even after a mesh refinement. The resulting predicted reattachment is farther downstream than indicated by the experimental data. Additionally, the solutions computed on a relatively coarse baseline mesh had detailed flow characteristics that were different from those computed on the refined mesh or the experimental data. Steady RANS solutions were also computed to investigate the effects of spanwise variation in the ice shape. The spanwise variation was obtained via a bleeding function that merged the ice shape with the clean wing using a sinusoidal spanwise variation. For these configurations, the results predicted for the extruded shape provided conservative estimates for the performance degradation of the wing. Additionally, the spanwise variation in the ice shape and the resulting differences in the flow fields did not significantly change the location of the primary reattachment.
A scheme for parameterizing ice-cloud water content in general circulation models
The optical properties of ice clouds are a primary issue for climate and climate change. Evaluating these optical properties in three-dimensional models for studying climate will require a method to calculate the ice water content of such clouds. A procedure is developed to parameterize ice water content as a function of large-scale meteorological characteristics for use in circulation models in which the ice water content is not calculated by means of a three-dimensional prognostic equation for condensed water. The technique identifies large-scale flows in which ice clouds exist and calculates their ice water content by reconstructing the trajectory associated with cloud formation. As the cloud forms, its ice content changes both by deposition of ice from water vapor and by ice removal by sedimentation. The sedimentation process is found to modify significantly the ice water content expected from deposition alone. Ice water contents predicted by the parameterization are compared with aircraft observations collected in the middle latitudes and the tropics, and show reasonable agreement over four orders-of-magnitude of ice water content. A parameterization for the sublimation of ice crystals settling into ice-subsaturated environments is also presented.
Rates of Sea‐Level Rise Are Highly Sensitive to Ice Viscosity Parameters in Model Benchmarks
Glacier flow plays a major role in current and future rates of globally averaged sea-level rise. The viscosity of glacial ice, controlling the rate of flow, decreases as stress increases and is highly sensitive to the value of the stress exponent, $n$, in the constitutive equation for viscous flow. Glaciologists and climate modelers almost exclusively assume $n=3$ when modeling ice flow and projecting sea-level rise through forward modeling. However, recent work suggests that $n\approx 4$ better fits observations, prompting the question: How sensitive are projections of sea-level rise to the value of $n$? We use an established community ice flow model and standard benchmark experiments designed as an idealized representation of Pine Island Glacier, West Antarctica. While initializing an $n=3$ model to match observations of an $n=4$ ice sheet is possible, we find that incorrectly assuming $n=3$ when in fact $n=4$ dramatically underestimates rates of sea-level rise. The scale of this error grows nonlinearly with the magnitude of the climate forcing, acting to increase projection uncertainties. Additionally, we find that models often account for this stress-dependent rheology mismatch during model initialization in a way that masks this rheological effect in the short term while leaving model outputs vulnerable to larger biases in longer-term projections. Initializations to observations of Pine Island Glacier display similar rheology-mismatch fingerprints to our idealized example.
A Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Reacting Hypersonic Flows
The necessity of enforcing conservation in computational elements or cells (element conservation) for discontinuous solutions is well understood and respected for solving conservation laws in computational fluid dynamics (CFD). In contrast, interface conservation, where the conservation across cell interfaces is enforced, is long ignored, and yet is also ruled and required by the underlying physics just like element conservation. Violation of the interface conservation across discontinuities is the root cause why an exact discontinuous solution can never be achieved in shock capturing methods. The interface conservation is examined and explored in this talk. Moving discontinuous Galerkin (MDG) finite element method with interface conservation enforcement (MDG-ICE) [1],[2] are then presented for solving compressible flow problems with discontinuities based on the observation that the interface conservation can only be satisfied, only when mesh interfaces are aligned with discontinuities. In the MDG-ICE formulation, both conservative quantities and grid geometry are considered as independent variables. A space-time DG formulation is used to solve the multi-material compressible Euler equations in the standard discontinuous solution space and the discrete grid geometry is solved using a variational formulation in a continuous space. A self-adaptive Levenberg-Marquardt method is utilized to solve the resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation. A number of numerical experiments for a variety of flow problems are conducted to assess the accuracy and performance of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to deliver the designed order of both h- and p-convergence even for discontinuous solutions, and detect all types of interfaces, via interface condition enforcement and satisfy, via grid movement, the compressible Euler equations and the associated interface condition.
An investigation of aircraft heaters XIV : an air and heat flow analysis of a ram-operated heater and duct system
Explore the source record for details and available documents.
An Investigation Utilizing an Electrical Analogue of Cyclic Deicing of a Hollow Steel Propeller with an External Blade Shoe
A study has been made of the heat requirement for the cyclic de-icing of hollow steel propellers fitted with external blade heating shoes. Solutions to the equations for the heat flow in cyclic heating of propellers were obtained, using an electrical analogy. The study showed how the energy requirement for propeller de-icing with existing blade shoes could be decreased, and illustrated the effect of blade-shoe design on the energy requirement. It was demonstrated, for example, that by increasing the heating intensity and decreasing the heating period from those currently used the energy requirement could be decreased in the order of 60 percent. ' In addition, ft was shown that heating requirements could be decreased further, by as much as 60 percent, through proper design of the shoes. The' investigation also showed the energy requirement to increase with decreasing liquid-water content and air temperature. Uncertainties as to the exact values of convective heat-transfer coefficient prevailing over the surface of the blade and ice layer resulted in uncertainties of approximately proportional magnitude in the values of required heating intensity.
Analysis of iced wings
A method for computing ice shapes along the leading edge of a wing and a method for predicting its aerodynamic performance degradation due to icing is described. Ice shapes are computed using an extension of the LEWICE code which was developed for airfoils. The aerodynamic properties of the iced wing are determined with an interactive scheme in which the solutions of the inviscid flow equations are obtained from a panel method and the solutions of the viscous flow equations are obtained from an inverse three-dimensional finite-difference boundary-layer method. A new interaction law is used to couple the inviscid and viscous flow solutions. The application of the LEWICE wing code to the calculation of ice shapes on a MS-317 swept wing shows good agreement with measurements. The interactive boundary-layer method is applied to a tapered ice wing in order to study the effect of icing on the aerodynamic properties of the wing at several angles of attack.
Analysis of iced wings
A method for computing ice shapes along the leading edge of a wing and a method for predicting its aerodynamic performance degradation due to icing is described. Ice shapes are computed using an extension of the LEWICE code which was developed for airfoils. The aerodynamic properties of the iced wing are determined with an interactive scheme in which the solutions of the inviscid flow equations are obtained from a panel method and the solutions of the viscous flow equations are obtained from an inverse three-dimensional finite-difference boundary-layer method. A new interaction law is used to couple the inviscid and viscous flow solutions. The application of the LEWICE wing code to the calculation of ice shapes on a MS-317 swept wing show good agreement with measurements. The interactive boundary layer method is applied to a tapered iced wing in order to study the effect of icing on the aerodynamic properties of the wing at several angles of attack.
Recent progress in the analysis of iced airfoils and wings
Recent work on the analysis of iced airfoils and wings is described. Ice shapes for multielement airfoils and wings are computed using an extension of the LEWICE code that was developed for single airfoils. The aerodynamic properties of the iced wing are determined with an interactive scheme in which the solutions of the inviscid flow equations are obtained from a panel method and the solutions of the viscous flow equations are obtained from an inverse three-dimensional finite-difference boundary-layer method. A new interaction law is used to couple the inviscid and viscous flow solutions. The newly developed LEWICE multielement code is amplified to a high-lift configuration to calculate the ice shapes on the slat and on the main airfoil and on a four-element airfoil. The application of the LEWICE wing code to the calculation of ice shapes on a MS-317 swept wing shows good agreement with measurements. The interactive boundary-layer method is applied to a tapered iced wing in order to study the effect of icing on the aerodynamic properties of the wing at several angles of attack.
Navier Stokes solution of the flowfield over ice accretion shapes
The numerical simulation of flow about ice accretion shapes has been accomplished by solving the Navier-Stokes equations using MacCormack's explicit finite difference scheme. The computations were performed on a CRAY-XMP computer. The influence of turbulence is taken into account by means of an algebraic eddy-viscosity model. In order to optimize the grid spacing and to achieve near orthogonality at the surface of the complex ice shapes, a hyperbolic grid generation scheme is utilized. Particular attention is given to the heat transfer process for which good agreement between the numerical and experimental results is achieved. In addition, liquid water droplet trajectories are coupled within the flowfield along with the resulting collection efficiencies using a parabolized Navier-Stokes formulation.
Modeling ice streams: Derived quantities
The model addressed is a finite-element, map-plane, time-dependent, column-averaged continuity equation solver. The key to the fitting process involves the balance between ice motion dominated by flow in internal layers, and ice motion dominated by sliding at the bed. The fitting process involves an iterative process carried out in the time domain. Beginning with the portion of the ice sheet being modeled identical to the present ice sheet with uniform flow, sliding, and fraction specified at nominal values, the model monitors each nodal point surface elevation. As the calculated surface elevation deviates from the present surface, a correction proportional to the difference is applied to selected parameter sets. This correction is in a sense that would tend to improve the fit at the particular nodal point. A calculated surface elevation that was higher than the present surface would result in an increased fraction, which would tend to lower the calculated surface (if the flow or sliding constant were being used as the fitting parameter, they would be lowered to improve the fit). This process is allowed to proceed as long as is necessary for the situation to stabilize. Typically, this takes tens of thousands of model years, but the rate is dependent on other external forcings such as the accumulation rate. The primary result is that while a typical sample of ice streams from around Antarctica can be fitted quite reasonably using only the fraction of the velocity due to sliding, a different mechanism seems to be in play along the Siple Coast, where reduced sliding constants are required to attain a reasonable fit. Flow is more strongly channelized in this region, and velocities are, in general, higher than are observed in other regions. It is unlikely that the mechanism that controls the ice movement along the Siple Coast is exactly similar to the mechanisms in the other ice streams. The concept of deformable sediments and their contribution to the fast flow along the Siple Coast may have limited applicability to other Antarctic ice streams.
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.