Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D Icing”

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 55 records · Page 3

Three dimensional empirical mode decomposition analysis apparatus, method and article manufacture

An apparatus and method of analysis for three-dimensional (3D) physical phenomena. The physical phenomena may include any varying 3D phenomena such as time varying polar ice flows. A repesentation of the 3D phenomena is passed through a Hilbert transform to convert the data into complex form. A spatial variable is separated from the complex representation by producing a time based covariance matrix. The temporal parts of the principal components are produced by applying Singular Value Decomposition (SVD). Based on the rapidity with which the eigenvalues decay, the first 3-10 complex principal components (CPC) are selected for Empirical Mode Decomposition into intrinsic modes. The intrinsic modes produced are filtered in order to reconstruct the spatial part of the CPC. Finally, a filtered time series may be reconstructed from the first 3-10 filtered complex principal components.

Gloersen, Per↗

Experimental Observations of Mechanisms of Pattern Formation During Freeze Casting

Empirically, many different materials compositions have been freeze cast over the years. Few have been analyzed to better understand the mechanisms by which the ice crystal growth drives structure formation, particularly the substructures on the cell walls. Collating substructural features, for example, ridges, bridges, and “jellyfish-like” substructures, found on the cell walls of freeze-cast ceramics (platelets versus particles), fibrillated biopolymers (collagen, nanocellulose), and "smooth" biopolymers (chitosan, trehalose), we capture the ice crystal growth process and analyze it both in situ in 2D and "post-mortem" in 3D to gain insights into anisotropic ice crystal growth and other self-assembly mechanisms. Combining these experimental observations with 3D Phase-field simulations of structure formation in ice-templated materials, detailed in a companion presentation, it becomes possible to gain fundamental science insights and to obtain the desired structure-property-processing correlations not only at the more traditional macroscopic, but also at the meso- and microscopic scales.

Kaiyang Yin↗

Experimental Observations of Mechanisms of Pattern Formation During Freeze Casting

Empirically, many different materials compositions have been freeze cast over the years. Few have been analyzed to better understand the mechanisms by which the ice crystal growth drives structure formation, particularly the substructures on the cell walls. Collating substructural features, for example, ridges, bridges, and “jellyfish-like” substructures, found on the cell walls of freeze-cast ceramics (platelets versus particles), fibrillated biopolymers (collagen, nanocellulose), and "smooth" biopolymers (chitosan, trehalose), we capture the ice crystal growth process and analyze it both in situ in 2D and "post-mortem" in 3D to gain insights into anisotropic ice crystal growth and other self-assembly mechanisms. Combining these experimental observations with 3D Phase-field simulations of structure formation in ice-templated materials, detailed in a companion presentation, it becomes possible to gain fundamental science insights and to obtain the desired structure-property-processing correlations not only at the more traditional macroscopic, but also at the meso- and microscopic scales.

Kaiyan Yin↗

TPSAS-NF1676L-26990-DND

The warm phase ENSO event (El Nino) of 2015/2016 produced significant perturbations to the Tropical Pacific, from sea surface to tropical tropopause and above. Although there have been 3 El Ninos during the modern "A-Train" era, the most recent El Nino started earlier in the year, was located farther west and was stronger than those of 2006/7 and 2009/10. Using A-Train observations, we show that during the 2015/16 El Nino the center of convective activity moved from the tropical western Pacific (TWP) to the central Pacific, just to the east of the dateline. This caused a record amount of cloud ice to occur near or above the environmental cold-point tropopause. Coincidently there was a large positive anomaly in tropical lower stratospheric water vapor, warming of the tropopause over the TWP, and changes to the stratospheric QBO. Using ten+ years of A-Train data and the MERRA-2 reanalysis fields, we explore the observed relationships between the 3D patterns of cloud ice, water vapor, ozone, aerosols and other trace gases in the tropical Pacific, and correlations with the ENSO and QBO cycles.

Melody A. Avery↗

13RC and and 3DRT Working Group of IRC: Status and Plans

We summarize recent activities of the International Intercomparison of 3D Radiation Codes (DRC), including development of an open source code, approximation methods for use in cloud dynamical models, and cloud reconstructions from multiple satellite platforms such as MODIS, MISR and ICESat. We also discuss the charter of the newly formed 3DRT Working Group of the International Radiation Commission, which hopes to coordinate and integrate efforts in various 3D application areas, including clouds, sea ice, and vegetation. Finally we describe plans for a special 3D session at the next International Radiation Symposium in Busan Korea, during 23-30 August 2004.

Cahalan, Robert↗

The Impact of Firn Models on Ultrawideband Brightness Temperatures in the Partially Coherent Model

The density profile of a polar ice sheet is an important parameter for the estimation of ice mass balance. Wave reflections caused by density variations are also a key uncertainty in the retrieval of ice sheet temperature profiles in Ultra-Wide band radiometry. In this paper, we examine different firn density profile models and analyze the subsurface reflections they cause using an analytical partially coherent approach. We also examine firn density profiles obtained from borehole measurements, from past UWBRAD modeling studies, from a community firn model, and from snow radar echo measurements. In previous studies, the ice sheet has been model as a 1D random medium with density variations in depth. However, horizontal density variations also exist, so that the ice sheet is a 3D random medium. Analyses using the partially coherent model show that in the presence of horizontal fluctuations, contributions from short scale variations vanish as the horizontal correlation length decreases due to the diffraction of waves.

Haokui Xu↗

Computational Icing Analysis on NASA’s SIDRM Geometry to Investigate Collection Efficiency

Computational icing analysis results were compared to experimental icing tunnel data including aerothermal (e.g. dry air) and supercooled water droplet rime-ice conditions from tests conducted in early 2022 at the NASA Icing Research Tunnel (IRT). The Simulated Intercompressor Duct Research Model (SIDRM) test article was used in this study, and its geometry represents the inter-compressor duct region of a turbofan engine. The test article’s purpose is to study the physics of supercooled water icing and ice crystal icing. This study compared three different icing codes: FENSAP-ICE (Eulerian approach), LEWICE3D (Lagrangian approach), and GlennICE (Lagrangian approach). All three icing codes were conducted on SIDRM’s complex body flow-field and compared to different experimental supercooled water rime runs. The test article instrumentation (pressure taps, thermocouples, etc.) and 3D laser scans of final ice shapes were used to compare against the different icing code simulations. The overall objectives are to understand how the icing codes handle capturing collection efficiency on the complex test article’s unheated surfaces. In the aerothermal cases, pressure tap readings matched the CFD results, but dry air CFD underpredicted thermocouple readings. Collection efficiency results from all three icing codes matched well together on the main body leading edge, main body slope, and the strut leading edges. All three icing codes underpredicted the final ice shape using a single-shot constant ice density approach, with more difficulty coming from the strut leading edge ice shape due to the swept wing like flow field. NASA’s overall goal for this effort is to develop computational icing tools to assist in the design and certification of engines for flight in icing conditions.

Supercooled icing, SIDRM, aerodynamics, computatio↗

Swept wing ice accretion modeling

An effort to develop a three-dimensional ice accretion modeling method is initiated. This first step toward creation of a complete aircraft icing simulation code builds on previously developed methods for calculating three-dimensional flowfields and particle trajectories combined with a two-dimensional ice accretion calculation along coordinate locations corresponding to streamlines. This work is intended as a demonstration of the types of calculations necessary to predict a three-dimensional ice accretion. Results of calculations using the 3D method for a MS-317 swept wing geometry are projected onto a 2D plane normal to the wing leading edge and compared to 2D results for the same geometry. These results indicate that the flowfield over the surface and the particle trajectories differed for the two calculations. This led to lower collection efficiencies, convective heat transfer coefficients, freezing fractions, and ultimately ice accumulation for the 3D calculation.

Potapczuk, M. G.↗

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m 3 . This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion, Large drops, Wind tunnel↗

Additional Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section from January 2022

In-flight icing is an important consideration that affects aircraft design, performance, certification and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of additional ice accretion testing conducted in the NASA Icing Research Tunnel in January 2022 with a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small median volumetric diameter (MVD) = 25 μm cloud and a large MVD = 110 μm cloud. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. Another series of tests explored cloud MVD variations from 50 μm to 230 μm while holding constant certain scaling parameters. The variation in ice mass and scanned ice volume across repeat conditions was approximately 50% lower than the uncertainty in the cloud MVD and liquid water content. The measured ice mass and volume calculated from the 3D scans were used to compute the mass to volume ratio that is sometimes referred to as ice density or void fraction. When the ice volume based on the ice shape maximum combined cross section was used to determine this ratio, the resulting values were in the range of 240 to 455 kg/m3. This is consistent with analogous values previously reported in the literature. The ice shape mass and volume increased with MVD from 50 to 230 μm at fixed values of the scaling parameters. The ice mass to volume ratio was approximately constant for all of the cases which showed that the accreted mass and volume increased at approximately the same rate. These results demonstrate the significance of cloud MVD on ice shape mass and volume.

Ice accretion↗

Subsurface Sounding of Mars: The Effects of Surface Roughness

The Mars Advanced Radar for Subsurface and Ionospheric Sounding (MARSIS) will conduct a global survey of Mars from the Mars Express Orbiter starting in 2004. The primary objective of the subsurface observations is to detect material interfaces in the upper several kilometers of the crust of Mars, with a particular emphasis on mapping the 3D distribution of water and ice in that portion of the crust. In order to detect subsurface interfaces, the returned echo from the subsurface must be distinguished from noise and clutter, which can arise from a variety of sources. One source of clutter is surface topography that generates backscattered energy at the same time delay as the subsurface region of interest. Surface topography can affect the detectability of subsurface features in several other ways. Surface roughness at scales comparable or somewhat smaller than the radar wavelength reduces the coherency of the wave as it passes the upper interface. Also, surface slope (tilt) at scales of the radar footprint and larger (> 5 km) affects the apparent Doppler signature of the echoes, and effectively disperses the wave transmitted into the subsurface, making processing and interpretation difficult. In this paper, we report on the roughness characteristics of Mars at these various scales as measured by the Mars Global Surveyor Laser Altimeter (MOLA), and consider the implications for achieving the subsurface sounding goals of MARSIS. Additional information is contained in the original extended abstract.

Plaut, J. J.↗

Exploration of Atmospheric Entries at Uranus & Neptune with HEEET as Heatshield TPS

One important observation from the recent Ice Giants Study sponsored by NASA was that the predicted and margined thicknesses of HEEET (new NASA TPS technology) were much greater than could be woven with the currently established loom capabilities. Since the cost of a loom upgrade would be substantial, the present work explores the entry trajectory space to determine what combinations of entry parameters would result in HEEET thicknesses that fit within the existing loom infrastructure. Toward this end, the entry trajectory space, parameterized by ballistic coefficient and entry flight path angle, was systematically explored for 45° sphere-cone geometries of 3 different radii 0.2 m, 0.3 m, and 0.4 m which covered the range from Galileo-derived probes considered in the Ice Giants Study, and a follow-on study on the possibility of using a single probe architecture (in terms of size and mass) for various destinations, including Venus, Saturn, Uranus, and Neptune. The entry velocities, latitudes, and azimuths at Uranus and Neptune used in the present work were taken from the Ice Giants Study. For each 3D OF trajectory generated by a NASA Ames in-house code, TRAJ, the material response and thickness were computed using another NASA Ames code, FIAT, along with a margins policy proposed by the HEEET project. In the present work, ballistic coefficients ranging from 200 kg/sqm to 350 kg/sqm were considered along with entry flight path angles ranging from -16° to -36° (primarily to allow deceleration loads to vary between 50 g and 200 g).

Prabhu, Dinesh K.↗

Assimilation of Sea Ice Thickness Derived from Cryosat-2 Along-Track Freeboard Measurements into the Met Office's Forecast Ocean Assimilation Model (FOAM)

The feasibility of assimilating sea ice thickness (SIT) observations derived from CryoSat-2 along-track measurements of sea ice freeboard is successfully demonstrated using a 3D-Var assimilation scheme, NEMOVAR, within the Met Office's global, coupled ocean–sea-ice model, Forecast Ocean Assimilation Model (FOAM). The CryoSat-2 Arctic freeboard measurements are produced by the Centre for Polar Observation and Modelling (CPOM) and are converted to SIT within FOAM using modelled snow depth. This is the first time along-track observations of SIT have been used in this way, with other centres assimilating gridded and temporally averaged observations. The assimilation leads to improvements in the SIT analysis and forecast fields generated by FOAM, particularly in the Canadian Arctic. Arctic-wide observation-minus-background assimilation statistics for 2015–2017 show improvements of 0.75 m mean difference and 0.41 m root-mean-square difference (RMSD) in the freeze-up period and 0.46 m mean difference and 0.33 m RMSD in the ice break-up period. Validation of the SIT analysis against independent springtime in situ SIT observations from NASA Operation IceBridge (OIB) shows improvement in the SIT analysis of 0.61 m mean difference (0.42 m RMSD) compared to a control without SIT assimilation. Similar improvements are seen in the FOAM 5 d SIT forecast. Validation of the SIT assimilation with independent Beaufort Gyre Exploration Project (BGEP) sea ice draft observations does not show an improvement, since the assimilated CryoSat-2 observations compare similarly to the model without assimilation in this region. Comparison with airborne electromagnetic induction (Air-EM) combined measurements of SIT and snow depth shows poorer results for the assimilation compared to the control, despite covering similar locations to the OIB and BGEP datasets. This may be evidence of sampling uncertainty in the matchups with the Air-EM validation dataset, owing to the limited number of observations available over the time period of interest. This may also be evidence of noise in the SIT analysis or uncertainties in the modelled snow depth, in the assimilated SIT observations, or in the data used for validation. The SIT analysis could be improved by upgrading the observation uncertainties used in the assimilation. Despite the lack of CryoSat-2 SIT observations available for assimilation over the summer due to the detrimental effect of melt ponds on retrievals, it is shown that the model is able to retain improvements to the SIT field throughout the summer months due to prior, wintertime SIT assimilation. This also results in regional improvements to the July modelled sea ice concentration (SIC) of 5 % RMSD in the European sector, due to slower melt of the thicker sea ice.

Emma K. Fiedler↗

Study of Martian Polar Caps with GISS ROCKE-3D GCM

Martian polar caps consist of both H 2 O and CO 2 ice. While H 2 O ice is mainly passive on modern Mars, it may have not been the case in recent Martian history, when its obliquity was higher, or when it was changing rapidly. The distribution of ice species in the snowpack affects its physical and thermodynamic properties. In the upper layers, it determines its albedo and thermal emissivity. Thus understanding the mutual effect between these ices and their interaction with the atmosphere is crucial for understanding the evolution of Martian polar regions. In this study, we employ a newly-developed Exotic Ices snow model coupled to the NASA Goddard Institute for Space Studies (GISS) ROCKE-3D planetary General Circulation Model (GCM) [1] to study the behavior of Martian polar caps. ROCKE-3D is a planetary GCM developed at NASA GISS as an extension of its Earth climate model, modelE [2]. It has been extensively used to simulate climate of various planets, including Mars (e.g. [3,4]). The Exotic Ices snow model was specially developed for planetary applications which involve more than one condensable in the atmosphere, in which case snow can contain multiple species of ice (CO2 and H 2 O in the Mars case). For each species of ice, the model uses their proper physical properties and phase diagram, but otherwise it treats all species of ice on an equal footing. The combined effects on albedo, thermal inertia and mutual insulation are treated accordingly. The snowpack interacts with the atmospheric dust cycle, and can accumulate a prognostic amount of dust, though the effect of dust on snow properties is not currently treated explicitly, and is prescribed. In this study, we first validate our model against the modern Martial climate, for which we use mission results from Mars Climate Sounder (atmospheric temperature and dust optical depth), SPICAM on Mars Express (atmospheric water), and Viking 2 (surface pressure). We investigate the effect of snow radiative properties on CO 2 and water cycles and the ability of our model to accurately reproduce those with minimal model tuning. We then perform simulations for several obliquities from a recent Martian past, and investigate the behavior of the Martian polar caps in such conditions.

Mars↗

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗

Effect of a simulated glaze ice shape on the aerodynamic performance of a rectangular wing

The effect of a simulated glaze-ice accretion on the flowfield of a 3D wing is studied experimentally. The model used for these tests was a semispan wing of effective aspect ratio five, mounted from the sidewall of a subsonic wind tunnel. The model has a NACA 0012 airfoil section on a rectangular untwisted planform with interchangeable leading edges to allow for testing both the baseline and the iced-wing geometry. A four-beam two-color fiberoptic laser Doppler velocimeter (LDV) was used to map the flowfield along three spanwise cuts on the model. Measurements on the centerline of the clean model compared favorably with theory and centerline measurements on the iced model compared well with measurements on a similar 2D model. The flow has the largest separation bubble at the model midspan with the smallest separation bubble occurring near the root and the wing tip.

Khodadoust, Abdi↗

An Overview of the TROPICS NASA Earth Venture Mission

The Time-Resolved Observations of Precipitation structure and storm Intensity with a Constellation of Smallsats (TROPICS) mission was selected by NASA as part of the Earth Venture-Instrument (EVI-3) program. The overarching goal for TROPICS is to provide nearly all-weather observations of 3D temperature and humidity, as well as cloud ice and precipitation horizontal structure, at high temporal resolution to conduct high-value science investigations of tropical cyclones. TROPICS will provide rapid-refresh microwave measurements (median refresh rate better than 60 min for the baseline mission) which can be used to observe the thermodynamics of the troposphere and precipitation structure for storm systems at the mesoscale and synoptic scale over the entire storm life cycle. TROPICS comprises six Cube-Sats in three low-Earth orbital planes. Each CubeSat will host a high-performance radiometer to provide temperature profiles using seven channels near the 118.75 GHz oxygen absorption line, water vapour profiles using three channels near the 183 GHz water vapour absorption line, imagery in a single channel near 90 GHz for precipitation measurements (when combined with higher-resolution water vapour channels), and a single channel near 205 GHz which is more sensitive to precipitation-sized ice particles. This observing system offers an unprecedented combination of horizontal and temporal resolution to measure environmental and inner-core conditions for tropical cyclones on a nearly global scale and is a major leap forward in the temporal resolution of several key parameters needed for assimilation into advanced data assimilation systems capable of utilizing rapid-update radiance or retrieval data.Launch readiness is currently projected for late 2019.

TROPICS↗

LEWICE3D/GlennHT Particle Analysis of the Honeywell Al502 Low Pressure Compressor

A flow and ice particle trajectory analysis was performed for the booster of the Honeywell AL502 engine. The analysis focused on two closely related conditions one of which produced a rollback and another which did not rollback during testing in the Propulsion Systems Lab at NASA Glenn Research Center. The flow analysis was generated using the NASA Glenn GlennHT flow solver and the particle analysis was generated using the NASA Glenn LEWICE3D v3.56 ice accretion software. The flow and particle analysis used a 3D steady flow, mixing plane approach to model the transport of flow and particles through the engine. The inflow conditions for the rollback case were: airspeed, 145 ms; static pressure, 33,373 Pa; static temperature, 253.3 K. The inflow conditions for the non-roll-back case were: airspeed, 153 ms; static pressure, 34,252 Pa; static temperature, 260.1 K. Both cases were subjected to an ice particle cloud with a median volume diameter of 24 microns, an ice water content of 2.0 gm3 and a relative humidity of 100 percent. The most significant difference between the rollback and non-rollback conditions was the inflow static temperature which was 6.8 K higher for the non-rollback case.

Ice Accretion Codes↗