Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase modeling”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Improbability of Post-Closure Criticality in Array of Transuranic Waste Containers after Compaction by Salt Creep at Waste Isolation Pilot Plant

Based on the rationale presented, post-closure nuclear criticality is improbable when room closure compacts containers disposing transuranic (TRU) waste emplaced at the Waste Isolation Pilot Plant (WIPP), an operating repository in bedded salt in southeastern New Mexico. As described in the original WIPP certification, a qualitative estimate of the probability of post-closure criticality in TRU waste produced from atomic energy defense activities has been low either because remote-handled TRU waste canisters are neutronically isolated by the bedded salt or because the low fissile mass in an array of contact-handled TRU waste drums cannot be compacted sufficiently by room closure from salt creep. These situations are still valid for the majority of TRU waste that is disposed at WIPP without any disposal constraints, as updated herein. This report also qualitatively evaluates the probability of criticality after disposal of TRU waste in pipe overpack containers (POC) where every POC in a shipment may have the maximum 200 fissile gram equivalent of 239Pu content. The probability of criticality for a disposal room filled with POCs is estimated during four representative phases of repository evolution: (1) a large salt block falls onto POCs in the first 20 years, (2) salt creep closes a disposal room in the first 1000 years without brine seepage and subsequent gas generation, which permits maximum compaction, (3) some brine seepage occurs into the closed room, which initiates consumption of the fiberboard (cellulose) impact absorber in the POC in the second 1000 years, and (4) full brine inundation of a room and consumption of all fiberboard thereafter. Salt-block fall in the first phase does not greatly disrupt three tiers of POCs. The compacted spacing of POCs in the later three repository conditions is calculated through high-fidelity, geomechanical modeling. Criticality evaluation of compacted 200-g 239 Pu spheres at the compacted spacing shows that neither 12-inch nor 6-inch POCs are critical after the first 1000 years, the second 1000 years, or thereafter as the sea of reflector material changes to represent the three repository conditions. Specifically, fiberboard and iron isolates 239 Pu while dry, and brine reduces the reactivity when a room is partially and fully inundated. Because POC behavior bounds behavior of other standard TRU waste containers, post-closure criticality caused by room closure compacting containers is omitted in the performance assessments for the 2019 and 2026 WIPP compliance re-certification applications to the US Environmental Protection Agency.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Properties of Hot Nuclear Matter

A fully quantitative description of the equilibrium and dynamical properties of hot nuclear matter will be needed for the interpretation of the available and forthcoming astrophysical data, providing information on the post-merger phase of a neutron star coalescence. We discuss the results of a recently developed theoretical model, based on a phenomenological nuclear Hamiltonian including two- and three-nucleon potentials, to study the temperature dependence of average and single-particle properties of nuclear matter relevant to astrophysical applications. The potential of the proposed approach for describing dissipative processes leading to the appearance of bulk viscosity in neutron star matter is also outlined.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Low data rate digital space communications

The low available transmitter power and the large frequency uncertainty constrain the data rate to be low. An all-digital communication receiver is proposed, and its feasibility is established. Although coherent systems should be used whenever practical, the noncoherent MFSK system is more suitable for very low data rates. The effect of Rician fading on the performance of MFSK receiver is studied. Fading characteristics of the Venus channel are examined based on the exponential model and available experimental data on the Venus atmosphere. Because of the requirement of high communication efficiency, three codes are evaluated and compared. The rapidly varying phase error at low data rate has great effects on the tracking loop behaviors which are examined by extensive computer study of the phase plane trajectories.

Chen, C. H.↗
NASA NTRS1984

A numerical analysis of transient planetary waves and the vertical structure in a meso-strato-troposphere model, part 1.4A

The structure of unstable planetary waves is computed by a quasi-geostrophic model extending from the surface up to 80 km by means of eigenvalue-eigenfunction techniques in spherical coordinates. Three kinds of unstable modes of distinct phase speeds and vertical structures are identified in the winter climate state: (1) the deep Green mode with its maximum amplitude in the stratosphere; (2) the deep Charney mode with its maximum amplitude in the troposphere: and (3) the shallow Charney mode which is largely confined to the troposphere. Both the Green mode and the deep Charney mode are characterized by very slow phase speeds. They are mainly supported by upward wave energy fluxes, but the local baroclinic energy conversion within the stratosphere also contributes in supporting these deep modes. The mesosphere and the troposphere are dynamically independent in the summer season decoupled by the deep stratospheric easterly. The summer mesosphere supports the easterly unstable waves 1-4. Waves 3 and 4 are identified with the observed mesospheric 2-day wave and 1.7-day wave, respectively.

Zhang, K. S.↗
NASA NTRS1989

Explicit modeling and computational load distribution for concurrent processing simulation of the space station

Two important aspects of concurrent processing under development at TRW are discussed. These are: (1) the derivation of explicit mathematical models of multibody dynamic systems, and (2) a balanced computational load distribution (BCLD) among loosely coupled computational units (processors) of a concurrent processing system. The developed methodologies are demonstrated by way of an application to the Phase 1 of the Space Station - a task being performed by TRW under NASA/JSC contract NAS9-17778. The mathematical model of the Space Station consists of three interconnected flexible bodies capable of undergoing large, rigid-body motion with respect to each other. Body 1 is the main central body and contains the pressurized modules inboard of the two Alpha gimbals. Bodies 2 and 3 are the starboard and port bodies connected to Body 1 at the Alpha gimbals and include all components on the transverse booms outboard of the Alpha gimbals (including the solar arrays). The control systems in the model maintain Body 1 in a prescribed 3-axis attitude control mode, while producing large-angle rotations of the flexible solar arrays to position them normal to the sun-line.

Gluck, R.↗
DOE OSTI2020

Integrated three-dimensional characterization of reactive phase formation and coarsening during isothermal annealing of metastable Zn–3Mg–4Al eutectic

Microstructural instabilities associated with a metastable lamellar eutectic when exposed to an elevated temperature were analyzed in a Zn-3Mg-4Al alloy as a model system. X-ray diffraction and scanning electron microscopy showed that annealing at a temperature 43 K below the equilibrium eutectic temperature converted the metastable ternary MgZn 2 /η-Zn/β-ZnAl eutectic to a stable mixture of Mg 2 Zn 11 /η-Zn/α-Al. 3D EBSD suggested that supersaturation of η-Zn with Mg and Al eased difficulties in nucleation of Mg 2 Zn 11 and α-Al at the MgZn 2 /η-Zn phase interfaces. Quasi-in-situ X-ray nanotomography revealed that a reactive diffusion process prevailed in an early stage of annealing leading to rapid consumption of MgZn 2 for growth of Mg 2 Zn 11 and concomitant thinning and pinch-off of η-Zn lamellae. The reactive diffusion process also brought about significant and quantifiable changes in the topological characteristics of the η-Zn phase at remarkably shorter time scales compared to Rayleigh instability and Ostwald ripening mechanisms. Overall, the present results provide fresh insights on the short-term annealing effects on metastable eutectics, with relevance to structural evolution in additive manufacturing.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗
NASA NTRS2020

Validation on the Europa Clipper Mission in the Formulation Phase

The Europa Clipper mission will perform a detailed investigation of Europa’s environment and habitability by placing a spacecraft into a looping orbit around Jupiter and performing repeated close Europa flybys. The Europa Clipper Project Verification and Validation (V&V) team is responsible for verifying that the Project System is compliant with the science and engineering requirements levied on it and validating that it is capable of meeting the mission objectives. Europa Clipper has developed a V&V framework to ensure that mission capabilities are realized at a reasonable confidence level prior to launch. The V&V team will coordinate the validation activities defined by this framework throughout the project lifecycle as the spacecraft evolves from design to implementation and operations. Validation of space missions is typically grouped into three categories: 1) requirements validation, 2) model validation, and 3) system validation. System validation, along with verification activities, are often the primary focus of V&V programs and these efforts peak in the implementation phase. However, the increasing reliance of exploration missions on models and simulations, and the complexity of the Europa Clipper mission, have dictated attention to requirements and model validation in the formulation phase. This paper discusses the development of the Europa Clipper V&V framework with emphasis on the validation of models and simulations, and of requirements. Approaches include the adaptation of institutional guidelines, such as those in the NASA-STD-7009 document, to create a Project Model V&V Plan and the use of tools such as DOORS NG to plan model validation and derive requirement quality metrics.

Srivastava, Priyanka↗
DOE OSTI2025

Measurement of directed flow in Au + Au collisions at $\sqrt{s{NN}}$ = 19.6 and 27 GeV with the STAR event plane detector

In heavy-ion collision experiments, the global collectivity of final-state particles can be quantified by anisotropic flow coefficients (𝑣 𝑛 ) . The first-order flow coefficient, also referred to as the directed flow (𝑣 1 ), describes the collective sideward motion of produced particles and nuclear fragments in heavy-ion collisions. It carries information on the very early stage of the collision, especially at large pseudorapidity (𝜂), where it is believed to be generated during the nuclear passage time. Directed flow therefore probes the onset of bulk collective dynamics during thermalization, providing valuable experimental guidance to models of the pre-equilibrium stage. In 2018, the Event Plane Detector (EPD) was installed in STAR and used for the Beam Energy Scan phase-II (BES-II) data taking. The combination of EPD (2.1 < |𝜂| < 5.1) and high-statistics BES-II data enables us to extend the 𝑣 1 measurement to the forward and backward 𝜂 regions. In this paper, we present the measurement of 𝑣 1 over a wide 𝜂 range in Au + Au collisions at $\sqrt{s_{NN}}$ = 19.6 and 27 GeV using the STAR EPD. The results of the analysis at $\sqrt{s_{NN}}$ = 19.6 GeV exhibit excellent consistency with the previous PHOBOS measurement, while elevating the precision of the overall measurement. The increased precision of the measurement also revealed finer structures in heavy-ion collisions, including a potential observation of the first-order event-plane decorrelation. Multiple physics models were compared to the experimental results. Only a transport model and a three-fluid hybrid model can reproduce a sizable 𝑣 1 at large 𝜂 as was observed experimentally. The model comparison also indicates 𝑣 1 at large 𝜂 might be sensitive to the QGP phase transition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Growth modeling to evaluate alternative cultivation strategies to enhance national microalgal biomass production

In this paper, we present a model-based assessment of alternative cultivation strategies on algal biomass production within the conterminous United States. Our assessment focuses on two basic cultivation strategies: 1) seasonal rotation of three representative freshwater algal strains that are well suited for warm weather, cold weather, and all-season diverse weather conditions, respectively; and 2) variation between three open pond water depths (15cm, 20cm, and 25cm). The enhanced Biomass Assessment Tool (BAT) is applied on a site-specific, hourly basis at 5,832 North American Land Data Assimilation System Phase 2 1/8° meteorological model grid cells over a 30-year period (1980-2009) to evaluate the operational strategies. Recognizing that resource management decisions may also consider scales beyond the individual site, we also conduct a regional assessment of focused on 7 representative climate zones in the conterminous United States. Results demonstrate that the spatial variability of algal productivity is largely affected by the strain-specific growth response to light and temperature that vary significantly by climate zone and latitude. Therefore, with a limited set of alternative cultivation strategies, the BAT identified the best performing combination of cultivation strategies on a site-specific basis that considerably enhances national annual biomass productivity. Specifically, the appropriate choice of seasonal strain rotation is able to significantly dampen climate-driven seasonal and spatial variability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Critical fluid dynamics in two and three dimensions

We describe a numerical method for simulating stochastic fluid dynamics near a critical point in the Ising universality class. This theory is known as model H, and is expected to govern the nonequilibrium dynamics of quantum chromodynamics (QCD) near a possible critical endpoint of the phase transition between a hadron liquid and the quark-gluon plasma. The numerical algorithm is based on a Metropolis scheme, and automatically ensures that the distribution function of the hydrodynamic variables in equilibrium is independent of the transport coefficients and only governed by the microscopic free energy. We verify dynamic scaling near the critical point of a two and three-dimensional fluid and extract the associated critical exponent z. Here, we find z≃3 in three dimensions, and z≃2 for a two-dimensional fluid. In a finite system, we observe a crossover between the mean field value z=4 and the true critical exponent z≃3 (z≃2 in d=2). This crossover is governed by the values of the correlation length and the renormalized shear viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2014

Constraints On the Emission Geometries and Spin Evolution Of Gamma-Ray Millisecond Pulsars

Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio beams and more compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio profiles using a single-altitude hollow-cone beam, with a core component when indicated by polarimetry; however, for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase, we assume that the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence intervals are determined using amaximum likelihood technique.We divide the light curves into three model classes, with gamma-ray peaks trailing (Class I), aligned (Class II), or leading (Class III) the radio peaks. Outer gap and slot gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit with pair-starved polar cap models. Distinguishing between the model classes based on typical derived parameters is difficult. We explore the evolution of the magnetic inclination angle with period and spin-down power, finding possible correlations. While the presence of significant off-peak emission can often be used as a discriminator between outer gap and slot gap models, a hybrid model may be needed.

pulsars: general↗
DOE OSTI2026

On the role of symmetry in quenching OH tunneling in 2,6-dimethylphenol

Here, we model the spectroscopy of the two methyl torsional degrees of freedom coupled to the OH torsional motion in 2,6-dimethylphenol. Recent gas-phase rotational transitions [Welsh et al ., J. Phys. Chem. Lett. 17 , 3749–3758 (2026)] have been interpreted to be the result of certain symmetry states of the methyl groups leading to the quenching of the OH torsional motion. Both a three-dimensional model and an adiabatic model, in which the OH torsion is treated as the slow mode, are developed to test these conjectures. Good agreement is found between the exact and approximate adiabatic models. The adiabatic model is developed to interpret and elucidate the key couplings leading to the quenching of the tunneling splitting.

Sibert, III, Edwin L. [Univ. of Wisconsin, Madison↗
DOE OSTI2025

Thermodynamic modeling of aqueous acetic acid, butyric acid, and lactic acid solutions

Based on the activity coefficient – fugacity coefficient approach, a rigorous thermodynamic modeling study is presented for accurate correlation of vapor-liquid equilibrium data of aqueous solutions of acetic acid (293 to 391 K), butyric acid (325 to 436 K), lactic acid (378 to 409 K), and acetic acid + butyric acid binary mixture (358 to 421 K). In addition, the pH data of the three aqueous, single carboxylic acid solutions were measured at 298 to 328 K and successfully correlated. Given that these aqueous carboxylic acid solutions exhibit various degrees of association behavior in both vapor and liquid phases, the thermodynamic models considered for this study include the Redlich-Kwong equation of state (RK-EoS) and the Hayden-O’Connell equation of state (HOC-EoS) for the vapor phase fugacity coefficients and the electrolyte non-random two-liquid model (eNRTL) and the association electrolyte non-random two-liquid model (AeNRTL) for the liquid phase activity coefficients. The combination of the HOC-EoS for the vapor phase and the AeNRTL model for the liquid phase is found to provide the best correlation results, consistent with the fact that the HOC-EoS and the AeNRTL model explicitly account for association behaviors in the vapor phase and the liquid phase, respectively.

09 BIOMASS FUELS↗
NASA NTRS2011

Gravito-Inertial Force Resolution in Perception of Synchronized Tilt and Translation

Natural movements in the sagittal plane involve pitch tilt relative to gravity combined with translation motion. The Gravito-Inertial Force (GIF) resolution hypothesis states that the resultant force on the body is perceptually resolved into tilt and translation consistently with the laws of physics. The purpose of this study was to test this hypothesis for human perception during combined tilt and translation motion. EXPERIMENTAL METHODS: Twelve subjects provided verbal reports during 0.3 Hz motion in the dark with 4 types of tilt and/or translation motion: 1) pitch tilt about an interaural axis at +/-10deg or +/-20deg, 2) fore-aft translation with acceleration equivalent to +/-10deg or +/-20deg, 3) combined "in phase" tilt and translation motion resulting in acceleration equivalent to +/-20deg, and 4) "out of phase" tilt and translation motion that maintained the resultant gravito-inertial force aligned with the longitudinal body axis. The amplitude of perceived pitch tilt and translation at the head were obtained during separate trials. MODELING METHODS: Three-dimensional mathematical modeling was performed to test the GIF-resolution hypothesis using a dynamical model. The model encoded GIF-resolution using the standard vector equation, and used an internal model of motion parameters, including gravity. Differential equations conveyed time-varying predictions. The six motion profiles were tested, resulting in predicted perceived amplitude of tilt and translation for each. RESULTS: The modeling results exhibited the same pattern as the experimental results. Most importantly, both modeling and experimental results showed greater perceived tilt during the "in phase" profile than the "out of phase" profile, and greater perceived tilt during combined "in phase" motion than during pure tilt of the same amplitude. However, the model did not predict as much perceived translation as reported by subjects during pure tilt. CONCLUSION: Human perception is consistent with the GIF-resolution hypothesis even when the gravito-inertial force vector remains aligned with the body during periodic motion. Perception is also consistent with GIF-resolution in the opposite condition, when the gravito-inertial force vector angle is enhanced by synchronized tilt and translation.

Wood, Scott J.↗
NASA NTRS1982

The CARE 3 Phase 3 Report: Test and Evaluation

CARE 3 (Computer-Aided Reliability Estimation, version three) is a computer program designed to help estimate the reliability of complex, redundant systems; although the program can model a wide variety of redundant structures, it was developed specifically for fault-tolerant avionics systems, systems distinguished by the need for extremely reliable performance since a system failure could well result in the loss of human life. CARE 3 further generalizes the class of system structures that can be modeled and greatly expands the coverage model to take into account such effects as intermittent and transient faults, latent faults, error propagation, etc. The initial test and evaluation of CARE 3 are reported.

Stiffler, J. J.↗
NASA NTRS1982

CARE 3 phase 2 report - mathematical description

CARE III (Computer-Aided Reliability Estimation, version three) a computer program designed to help estimate the reliability of complex, redundant systems is described. Although the program can model a wide variety of redundant structures, it was developed specifically for fault tolerant avionics systems. CARE III generalizes the class of system structures that can be modeled and greatly expands the coverage model to take into account such effects as intermittent and transient faults, latent faults, and error propagation.

Stiffler, J. J.↗
NASA NTRS1986

A numerical and experimental investigation of electrochemical aircraft deicing

This study was composed of three parts. The first part involved the extension of an existing transient two dimensional numerical code for an electrothermal deicer so that it would simulate the situation where a variable thickness ice layer existed at the outer surface. The Enthalpy Method was used to simulate the phase change, and Gauss-Seidel iteration was used to solve the resulting system of finite difference equations. A set of criteria were developed for determining when a variable thickness ice layer had an effect on deicer performance. The second part was the acquisition and analysis of experimental data. The test model was a section of a Bell UH-1H helicopter blade equipped with an electrothermal deicer. A total of fifty-two thermocouples were utilized to document the thermal response of the blade and decier assembly. In the deicing runs, the experimental temperature response data clearly showed when melting, shedding or refreezing occurred. The tests illustrated that the criterion for shedding in the three cases where it did occur was that the abrasion shield interface temperature was 32 to 34 F. The third part concerned the validation of a one dimensional transient thermal model of an electrothermal deicer by comparison of the predictions with the experimental data. The Enthalpy Method was found to effectively model the phase change which occurred, and the ice shedding algorithm employed in the simulation was also evaluated.

Leffel, K. L.↗