Development and Validation of 2D Finite Difference Modeling of Heat Transfer, Degree of Healing, and Stress Development in FFF Printing
Explore the source record for details and available documents.
Eric Gilmer↗
SEARCH · Search NASA
Results for “Heat Transfer”
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.
Heat Transfer Enhancement in Cryogenic Quenching Process
This paper reports a heat transfer advancement in the cryogenic quenching process. An experiment was performed to evaluate the enhancement of quenching heat transfer by the use of metal tubes with low thermal conductivity coating layers. Four coating thicknesses with various coolant mass flow rates of liquid nitrogen were investigated. The results indicated that the tube inner surface coating greatly enhanced the quenching efficiency. The quenching efficiency was found to increase with increasing number of coating layers, and the efficiency also increased with decreasing mass flow rates. In general, the efficiencies cover a range between 40.6% and 80%. Comparing to the bare surface case, the percentage increase in the quenching efficiency was the minimum at 4.2% for a single coated layer at the highest flow rate and the maximum of 109.1% for four coated layers at the lowest flow rate. The coated tubes could save up to 53% in the amount of cryogen consumption.
two-phase flow↗
Final Thermal Design and Thermal Vacuum Testing of the StarBurst Instrument
Explore the source record for details and available documents.
Small Satellite↗
A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft
Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.
heat transfer↗
Cryogenic Flow Boiling in Microgravity: Effects of Reduced Gravity on Two-Phase Fluid Physics and Heat Transfer
With the growing interest in space exploration, cryogenic technologies involving two-phase flow and heat transfer are in high demand to successfully procure advanced space applications such as fuel depots and nuclear thermal propulsion (NTP) systems for deep space missions. However, the unique and extreme thermal properties of cryogenic fluids introduce distinct flow boiling fluid physics and energy transport phenomena, which differ significantly from those observed with conventional fluids. Understanding the unique two-phase physics in cryogenic flow boiling remains an ongoing challenge. Furthermore, the lack of readily available microgravity cryogenic steady-state heat transfer data hinders the assessment of gravitational effects on cryogenic flow boiling. This study aims to elucidate the gravitational effects on two-phase fluid physics and heat transfer by conducting the first-ever experimental measurement of cryogenic flow boiling performance using a steady-state heated method in a reduced gravity environment. Parabolic flight experiments were performed to acquire both heat transfer measurements and high-speed video of interfacial behaviors, under varying gravity levels (microgravity, hypergravity, Lunar gravity, and Martian gravity). The experiments involved flow boiling of liquid nitrogen (LN 2 ) with a near-saturated inlet along a circular heated tube of dimensions 8.5-mm inner diameter and 680-mm heated length. The operating parameters varied are mass velocity of 398.3 - 1342.8 kg/m2s, inlet quality of -0.08 to -0.01, and inlet pressure of 413.68 - 689.48 kPa. Captured microgravity flow patterns range from bubbly to annular, all having vapor structures that are larger than those under higher gravity levels. Under microgravity, absence of buoyancy yields symmetrical vapor structures without flow stratification, laying a physical foundation for the distinct two-phase heat transfer trends during LN 2 flow boiling in microgravity. Transient data collected during the flight parabolas exhibited decreasing heated wall temperature as the aircraft transitioned from hypergravity to microgravity phases. The temperature variation indicated an enhancement in flow boiling heat transfer with decreasing gravity levels and a reduction with increasing gravity levels. The effect of reduced gravity on cryogenic flow boiling heat transfer coefficient (HTC) is discussed based on steady state heat transfer analysis. Seminal HTC correlations are evaluated against the measured microgravity HTC data, of which one is identified for superior accuracy in predicting microgravity data. Finally, a new HTC correlation is proposed to improve accuracy of microgravity predictions, yet there still exists room for further improvement with future terrestrial flow boiling experiments at different flow orientations relative to Earth gravity.
Microgravity↗
Numerical Modeling of the Transient Chilldown Process of a Cryogenic Propellant Transfer Line
Before cryogenic fuel depots can be fully realized, efficient methods with which to chill down the spacecraft transfer line and receiver tank are required. This paper presents numerical modeling of the chilldown of a liquid hydrogen tank-to-tank propellant transfer line using the Generalized Fluid System Simulation Program (GFSSP). To compare with data from recently concluded turbulent LH2 chill down experiments, seven different cases were run across a range of inlet liquid temperatures and mass flow rates. Both trickle and pulse chill down methods were simulated. The GFSSP model qualitatively matches external skin mounted temperature readings, but large differences are shown between measured and predicted internal stream temperatures. Discrepancies are attributed to the simplified model correlation used to compute two-phase flow boiling heat transfer. Flow visualization from testing shows that the initial bottoming out of skin mounted sensors corresponds to annular flow, but that considerable time is required for the stream sensor to achieve steady state as the system moves through annular, churn, and bubbly flow. The GFSSP model does adequately well in tracking trends in the data but further work is needed to refine the two-phase flow modeling to better match observed test data.
Temperature Sensors↗
Comparison of Cryogenic Flow Boiling in Liquid Nitrogen and Liquid Hydrogen Chilldown Experiments
This paper presents a comparison between experimental results from recent liquid hydrogen (LH 2 ) transfer line chilldown experiments at high Reynolds (Re) numbers versus liquid nitrogen (LN 2 ) experiments conducted at low Re numbers. Parasitic heat leak, inner wall temperatures, inner wall heat fluxes, and heat transfer coefficients are computed to compare between the two systems. Analysis of temperature traces and flow visualization indicates that the chilldown process evolves much more rapidly at higher Re numbers due to a quick transition from vapor flow to annular liquid flow and near immediate liquid contact along the pipe walls. The lower kinematic viscosity and surface tension of LH 2 , along with reduced parasitic heat leak and higher Re numbers relative to the LN 2 experiments, causes chilldown to proceed almost immediately into the nucleate boiling regime, in comparison to low Re flows where >75% of the chilldown is spent in vapor film boiling.
Cryogenic Heat Transfer Coefficient↗
Cryogenic Spray Quenching of A Simulated Propellant Storage Tank Wall With Heat Transfer Enhancement By A Thin-Film Coating and Flow Pulsing in Microgravity
Human space exploration to the Moon, Mars, and possibly asteroids is NASA’s biggest challenge for the new millennium. One of the critical elements to this mission is the effective, sufficient, and reliable supply of cryogenic propellant fluids. Future lower-earth-orbiting (LEO) propellant fuel depots and human-carrying orbital transfer spacecraft flying to the moon and Mars will have to utilize the high thrust and high efficiency of liquid cryogenic chemical propulsion or nuclear thermal propulsion. Efficient in-space tank-to-tank propellant transfer (propellant fuel depot to orbital transfer spacecraft) of cryogenic propellants is an enabling technology for the planned Crewed Mars Surface Mission. The transfer of cryogenic propellants in space, however, has yet to be accomplished, solely due to the unavailability of cryogenic quenching heat transfer data during chilldown (quenching) and filling of the propellant receiver tank in reduced gravity and microgravity as liquid propellant cannot be stored in a required liquid state until the tank is quenched down to the liquid temperature. Therefore, highly energy efficient thermal-fluid management breakthrough concepts to conserve and minimize the cryogen consumption during propellant transfer have become the focus of research and engineering development, especially for the deep-space mission to Mars. In this paper, we introduce such concepts and demonstrate their feasibility for cryogenic storage tank chilldown in parabolic flights under a simulated space microgravity condition. In order to maximize the storage tank chilldown efficiency for the least amount of cryogen consumption, the technology adopted included cryogenic spray cooling, Teflon thin-film coating of the simulated tank surface, and spray flow pulsing. The completed flight experiments successfully demonstrated that spray cooling is the most efficient cooling method for the tank chilldown in microgravity. In microgravity, Teflon coating alone can improve the efficiency up to 72% and the efficiency can be improved up to 59% by flow pulsing alone. However, Teflon coating together with flow pulsing was found to substantially enhance the chilldown efficiency in microgravity for up to 113%.
spray↗
Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching
Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.
line chilldown↗
The NASA Turbulent Heat Flux (THX) Experiments: Summary and Lessons Learned
The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds-averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.
RANS↗
Development and Validation of a High-Vacuum Thermal Conductivity Testbed for Aerospace Interface Materials
Spacecraft thermal margins depend on the temperature penalty of installed interfaces, yet catalog conductivity omits bondline thickness, mating surfaces, preload, and environment. The Testbed for Advanced Interface Materials in Vacuum (TAIMV) addresses this engineering-data gap by converting ambient/vacuum temperature fields from a four-coupon stack into quantities usable in spacecraft thermal models. Its 17-node steady-state model resolves axial conduction, parasitic fixture paths, grease-filled contact, radiation, and ambient convection. Four monolithic checks gave axial heat-rate ratios of 0.90–1.02, and Braycote calibration self-recovery gave 2.20% casewise mean absolute percentage error (MAPE) and 2.85% maximum difference. With transferable fixture terms frozen, Krytox LVP tested cross-material transfer against a manufacturer-derived external relation. The six Krytox cases gave 2.28% MAPE, 5.10% maximum difference, and full-field residual root-mean-square values of 0.553 °C ambient and 1.306 °C vacuum. Installed-joint resistance spanned 6.39–7.15 × 10⁻⁴ m²·K/W for Braycote and 4.09–4.49 × 10⁻⁴ m²·K/W for Krytox, corresponding to 2.33–6.66 K per modeled interface over the tested heat-flux range. TAIMV therefore supplies directly usable installed-joint resistance and conductance, plus apparent installed conductivity and explicitly model-conditioned grease conductivity/contact terms.
heat transfer↗
Extending mARC II Arc-Jet Test Duration via Design and Implementation of Model System Cooling Sleeve
The mARC II is a 30 kW arc-jet facility at NASA Ames Research Center used to generate high-enthalpy flows for low-cost thermal protection system (TPS) technology development. Sustained operation of downstream instrumentation and material samples is constrained by thermal loading transmitted through the arc-jet test environment, limiting achievable run times and experimental throughput. This work presents the design, integration, and validation of a cooling sleeve implemented on the sweep arm drive motor feedthrough to mitigate thermal accumulation during testing. The addition of the cooling sleeve is a simple, robust upgrade that translates directly into enhanced facility capability by supporting longer run durations, reduced turnaround time, and higher throughput.
numerical simulations↗
Thermal Design and Thermal Vacuum Testing of the StarBurst Instrument
The StarBurst Multimessenger Pioneer is a small satellite mission serving as a wide-field gamma-ray observatory designed to capture the initial emissions of short gamma-ray bursts, electromagnetic signatures of neutron star mergers. This paper presents the final thermal design and analysis of the StarBurst Instrument, comprising the bus-to-instrument interface plate, control electronics, and twelve crystal detector units, which form the core of the mission’s science capability. The passive thermal control system design requires consideration of restrictive keep-out zones, unknown orbital parameters, and narrow temperature limits of the detectors. Also summarized is the instrument level thermal vacuum cycle test, correlated model refinements, and updated model results. Following successful completion of the instrument test campaign, the hardware was integrated with the spacecraft bus for spacecraft level testing, including additional thermal vacuum testing. The results from the spacecraft level thermal vacuum test will further inform the instrument thermal model, ensuring accurate flight temperature predictions. StarBurst launches as a secondary payload in 2027 and has a mission duration of at least one year.
StarBurst↗
A GPU-based Approach for Turbomachinery Application
Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.
Gas Turbine↗
A Machine Learning Framework for Error Compensation in Radiative Transfer Calculations
Radiative heat transfer influences the amount of heat flux transferred to the surface of the hypersonic vehicle, which is essential to evaluate the performance of thermal protection systems. The radiative heat flux is found to be computationally prohibitive while accounting for the variation in spatial, angular, and spectral domains. A new methodology has been recently developed to alleviate the cost of computation in the spectral domain by constructing flow-agnostic reduced-order models (ROMs). The developed spectral ROM databases provide grouping strategies that account for non-equilibrium absorption and emission as well as interaction between disparate species due to spectral overlap in associated radiative processes. However, the developed ROMs need to be optimized for a specific combination of interacting gas species and would need to re-calibrated in case individual species are added/omitted. In this work, we use various machine learning (ML) techniques to approximate the radiative intensities determined by a ROM optimized for a specific gas mixture. The ML model relies on the ROM databases developed for a single species which ignores any spectral overlap. Thus, radiation evaluation starts with a simple summation of radiative intensities predicted using these non-calibrated ROMs for the contributing species. The ML framework then provides a correction to account for the interplay in the frequency, i.e., emission of photons by one species and absorption by another, and yields mixture-specific radiation fields. Once trained on the individual ROM databases, the ML framework offers instantaneous corrections that serves as a time/cost effective alternative to the optimization of ROMs for a specific gas mixture. The ML framework is trained on both the high fidelity and ROM evaluated line of sight (LOS) data from Orion, Stardust, and FIRE II cases to obtain a general purpose correction model for earth re-entry scenarios when radiation contributions from both atomic nitrogen and atomic oxygen are considered. A geometric length scale parameter is used in the training process to account for errors introduced in the ROM databases as a consequence of high optical thickness. The efficacy of the ML framework is underscored through extensive analysis of train and test errors with respect to all the re-entry scenarios. The applicability of such an ML framework was further corroborated by embedding it in a state-of-the-art US3D - NERO system for determining the radiative heat flux transferred to the hypersonic vehicle surface.
Radiation↗
LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow
Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.
Secondary Flows↗
Development and Validation of a High-Vacuum Thermal Conductivity Testbed for Aerospace Interface Materials
Thermal Interface Materials (TIMs) are critical components in spacecraft thermal management systems, where thermal performance is strongly influenced by vacuum conditions, interface contact resistance, and layered metallic joint behavior. However, manufacturer-reported thermal conductivity values are often derived under idealized conditions and may not accurately represent performance within operational aerospace applications. To address this limitation, the Testbed for Advanced Interface Materials in Vacuum (TAIMV) was developed as a modular vacuum-compatible thermal conductivity characterization platform capable of evaluating aerospace-relevant TIM configurations under both ambient and high-vacuum environments. The testbed was derived from the ASTM C1044-16 guarded hot plate methodology and incorporates interchangeable layers of stainless steel coupon geometries, independently controlled main and guard heaters, embedded resistance temperature detectors (RTDs), thermocouples, multi-layer insulation (MLI), and a temperature-controlled cold plate to characterize through-thickness thermal gradients across layered interfaces. In the current configuration, interface compression is limited to the nominal contact pressure generated by the experimental stack assembly. Initial experimental campaigns were conducted at ambient pressure and below 1×10-5 torr for vacuum cases using multiple interface materials including Braycote 601EF and Krytox-based greases across a range of thermal operating conditions. In parallel, a coupled numerical Python thermal model was developed to predict temperature distribution throughout the stack while accounting for conduction, radiation, and parasitic heat transfer pathways and effective interface resistance effects. Experimental measurements and numerical predictions showed consistent thermal trends across multiple operating conditions and environmental states. Results also revealed measurable differences between ambient and vacuum thermal behavior, demonstrating the importance of interface resistance, parasitic heat transfer mechanisms, and stack geometry in determining effective thermal performance within layered thermal interfaces. The presented work establishes a foundation for future thermal model correlation efforts and expanded characterization of aerospace thermal interface materials under representative environmental conditions. Future work will focus on the integration of a load cell system to enable controlled pressure-dependent characterization of thermal interface materials under compressive loading. This capability will allow investigation of the influence of contact pressure on effective thermal conductivity, interface resistance, and thermal performance within layered aerospace thermal interfaces under representative operational conditions.
Thermal Development Testing↗
CFD-Assisted Nodal Modeling of Sloshing in a Cryogenic Propellant Tank
During autogenous pressurization, tank sloshing causes a significant increase in pressurant consumption to maintain constant ullage pressure during draining of the tank. This increase in pressurant consumption is caused by a significant increase in condensation at the liquid-vapor interface. Sloshing strongly affects the liquid side heat transfer coefficient and thereby the condensation rate. Traditionally, sloshing is modeled by CFD code using the VOF (Volume of Fluid) method to track the liquid vapor interface during sloshing. CFD calculations require a very fine grid to accurately compute the heat and mass transfer at the interface. Therefore, computations are time consuming and prohibit performing many parametric studies often needed during the design of a new system. This paper describes an alternative approach by developing a CFD-assisted nodal model to predict system parameters more economically with reasonable accuracy. In this approach, a multi-node model of tank pressurization was developed using GFSSP. A multi-node model was needed to account for stratification. The model computes heat and mass transfer at the interface to calculate the condensation rate. The liquid side heat transfer is computed using the parameters of sloshing dynamics such as frequency, wave amplitude, and interface area. The parameters of sloshing dynamics are computed by the CFD code LOCI-Stream. The model predictions were compared with test data for several cases.
Cryogenic Tank Sloshing↗