SEARCH · Engineering Papers

Results for “grid integration study”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Conical Euler analysis and active roll suppression for unsteady vortical flows about rolling delta wings

A conical Euler code was developed to study unsteady vortex-dominated flows about rolling, highly swept delta wings undergoing either forced motions or free-to-roll motions that include active roll suppression. The flow solver of the code involves a multistage, Runge-Kutta time-stepping scheme that uses a cell-centered, finite-volume, spatial discretization of the Euler equations on an unstructured grid of triangles. The code allows for the additional analysis of the free to-roll case by simultaneously integrating in time the rigid-body equation of motion with the governing flow equations. Results are presented for a delta wing with a 75 deg swept, sharp leading edge at a free-stream Mach number of 1.2 and at 10 deg, 20 deg, and 30 deg angle of attack alpha. At the lower angles of attack (10 and 20 deg), forced-harmonic analyses indicate that the rolling-moment coefficients provide a positive damping, which is verified by free-to-roll calculations. In contrast, at the higher angle of attack (30 deg), a forced-harmonic analysis indicates that the rolling-moment coefficient provides negative damping at the small roll amplitudes. A free-to-roll calculation for this case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation, which is characteristic of highly swept delta wings. This limit cycle oscillation may be actively suppressed through the use of a rate-feedback control law and antisymmetrically deflected leading-edge flaps. Descriptions of the conical Euler flow solver and the free-to roll analysis are included in this report. Results are presented that demonstrate how the systematic analysis of the forced response of the delta wing can be used to predict the stable, neutrally stable, and unstable free response of the delta wing. These results also give insight into the flow physics associated with unsteady vortical flows about delta wings undergoing forced motions and free-to-roll motions, including the active suppression of the wing-rock type phenomenon. The conical Euler methodology developed is directly extend able to three-dimensional calculations.

Lee-Rausch, Elizabeth M.↗

NASA NTRS2005

Brayton-Cycle Power-Conversion Unit Tested With Ion Thruster

Nuclear electric propulsion has been identified as an enabling technology for future NASA space science missions, such as the Jupiter Icy Moons Orbiter (JIMO) now under study. An important element of the nuclear electric propulsion spacecraft is the power conversion system, which converts the reactor heat to electrical power for use by the ion propulsion system and other spacecraft loads. The electrical integration of the power converter and ion thruster represents a key technical challenge in making nuclear electric propulsion technology possible. This technical hurdle was addressed extensively on December 1, 2003, when a closed- Brayton-cycle power-conversion unit was tested with a gridded ion thruster at the NASA Glenn Research Center. The test demonstrated end-to-end power throughput and marked the first-ever coupling of a Brayton turbo alternator and a gridded ion thruster, both of which are candidates for use on JIMO-type missions. The testing was conducted at Glenn's Vacuum Facility 6, where the Brayton unit was installed in the 3-m-diameter vacuum test port and the ion thruster was installed in the 7.6-m-diameter main chamber.

Hervol, David S.↗

DOE OSTI2024

Optimizing Grid Regulation With Gravity Storage Systems: A Comparative Analysis With Different Motor Inertias: Preprint

The integration of renewable energy sources into power grids necessitates solutions for grid support and stability during fluctuations in electricity generation and demand. Gravity energy storage systems (GESS) are emerging as a promising technology for managing the balance between the energy supply and demand. However, their capacity to optimize energy flow and offer voltage and frequency regulation amid generation-demand imbalances is less reported. This paper investigates the control of GESS for optimizing energy flow during voltage and frequency regulation. The study evaluates the regulation capabilities of GESS with different motor inertias during a Texas grid event: one with a high-speed, low-inertia motor and another with a low-speed, high-inertia motor. Results indicate that both GESS systems provide fast frequency response by converting potential energy into kinetic energy and vice versa, with a response time of 1.5 seconds from the frequency variation. This aligns with grid requirements for primary frequency response from traditional synchronous generators and motors with large inertia. Furthermore, a GESS with a high-inertia motor may be able operate over a broader range of frequency variations whereas a low-inertia system may be limited by thermal constraints of the motor.

frequency regulation↗

DOE OSTI2024

Optimizing Grid Regulation with Gravity Energy Storage Systems: A Comparative Analysis with Different Motor Inertias

The integration of renewable energy sources into power grids necessitates solutions for grid support and stability during fluctuations in electricity generation and demand. Gravity energy storage systems (GESS) are emerging as a promising technology for managing the balance between energy supply and demand. However, their capacity to optimize energy flow and offer voltage and frequency regulation amid imbalances in generation and demand is less reported. This paper investigates the control of GESS for optimizing energy flow during voltage and frequency regulation. The study evaluates the regulation capabilities of GESS with different motor inertias during a Texas grid event: one with a high-speed, low-inertia motor and another with a low-speed, high-inertia motor. Results indicate that both GESS scenarios provide fast frequency response by converting potential energy into kinetic energy and vice versa, with a response time of 1.5 s from the frequency variation. This aligns with grid requirements for primary frequency response from traditional synchronous generators and motors with large inertia. Furthermore, a GESS based on a highinertia motor may be able to operate over a broader range of frequency variations, whereas a low-inertia system may be limited by thermal constraints of the motor.

ENERGY STORAGE,POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Asymptotically compatible reproducing kernel collocation and meshfree integration for the peridynamic Navier equation

Here, we study reproducing kernel (RK) collocation method for peridynamic Navier equation. In the first part, we apply a linear RK approximation to both displacement and dilatation, and then back-substitute dilatation and solve the peridynamic Navier equation in a pure displacement form. The RK collocation scheme converges to the nonlocal limit for a fixed nonlocal interaction length and also to the local limit as nonlocal interactions vanish. The stability is shown by comparing the collocation scheme with the standard Galerkin scheme using Fourier analysis. In the second part, we apply the RK collocation to the quasi-discrete peridynamic Navier equation and show its convergence to the correct local limit when the ratio between the nonlocal length scale and the discretization parameter is fixed. The analysis is carried out on a special family of rectilinear Cartesian grids for the RK collocation method with a designated kernel with finite support. We assume the Lamé parameters satisfy λ ≥ μ to avoid extra assumptions on the nonlocal kernel. Finally, numerical experiments are conducted to validate the theoretical results.

42 ENGINEERING↗

DOE OSTI2023

An evaluation of the demand response potential of integrated dynamic window and HVAC systems

Demand response (DR) increases the flexibility and reliability of the electricity grid as use of intermittent renewable energy sources increases. HVAC and envelope DR measures present the largest aggregate energy and peak demand savings potential of all commercial building end uses because their net demand savings occur during critical peak demand periods. Controllable envelope measures include switchable electrochromic windows, operable window attachments such as outdoor louvers, roller shades, and awnings, as well as other innovative facade technologies that can modulate both solar heat gain and daylight admission over a broad solar-optical range. This study evaluated the technical potential of DR-enabled dynamic windows to reduce critical peak demand for a prototypical medium office building situated in all 16 U.S. climates. Model predictive control (MPC) algorithms were designed to minimize electricity cost in daylit perimeter office zones through control of an electrochromic window with and without HVAC thermostat setpoint control. Conventional and time-of-use rates were used to shape the degree of DR. Median annual peak demand savings with window and thermostat control across all climate zones were 24.3 kW (4.4 W/m 2 ) per building or 15.9 W/m 2 for non-north perimeter zones. Resource adequacy at the whole building level was estimated to be 13.1 to 43.4 $/kW per year over the 30-year life of the installation. Co-benefits were increased energy efficiency, and reduced electricity cost and emissions. Visual and thermal comfort requirements were met at all times. Dynamic facades controlled by MPC have substantial technical potential for DR across all U.S. climates and warrant serious consideration for inclusion in DR portfolios.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS2013

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: 1) a k-epsilon turbulence model computation on a 6.8 million point grid using wall functions, 2) a k-epsilon turbulence model computation on a 14 million point grid integrating to the wall, and 3) a k-omega turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-epsilon and k-omega computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-epsilon computations which extended from 70% to 100% span at the exit rating plane, whereas the k-omega computation had reversed flow from 95% to 100% span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-epsilon case resulted in an under-prediction in adiabatic efficiency of 8.3 points, whereas the k-omega case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗

NASA NTRS2013

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: (1) a k-ε turbulence model computation on a 6.8 million point grid using wall functions, (2) a k-ε turbulence model computation on a 14 million point grid integrating to the wall, and (3) a k-ω turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-ε and k-ω computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-ε computations which extended from 70 to 100 percent span at the exit rating plane, whereas the k-ω computation had reversed flow from 95 to 100 percent span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-ε case resulted in an underprediction in adiabatic efficiency of 8.3 points, whereas the k-ω case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗

DOE OSTI2022

What you get is not always what you see—pitfalls in solar array assessment using overhead imagery

Effective integration planning for small, distributed solar photovoltaic (PV) arrays into electric power grids requires access to high quality data: the location and power capacity of individual solar PV arrays. Unfortunately, national databases of small-scale solar PV do not exist; those that do are limited in their spatial resolution, typically aggregated up to state or national levels. While several promising approaches for solar PV detection have been published, strategies for evaluating the performance of these models are often highly heterogeneous from study to study. The resulting comparison of these methods for practical applications for energy assessments becomes challenging and may imply that the reported performance evaluations overly optimistic. The heterogeneity comes in many forms, each of which we explore in this work: the degree of diversity of the locations and sensors (e.g. different satellites, aerial photography) from which the training and validation data originate, the validation of ground truth (manual annotation of imagery vs known solar PV locations), the level of spatial aggregation (e.g. array-level vs regional estimates), and inconsistencies in the training and validation datasets (e.g. different datasets are used for each study and those data are not always made accessible). For each, we discuss emerging practices from the literature to address them or suggest directions of future research. As part of our investigation, we evaluate solar PV identification performance in two large regions: the entire state of Connecticut and the city of San Diego, CA. In Connecticut, we also use 33,114 known parcel-level solar PV installations from Berkeley Lab’s Tracking the Sun dataset to evaluate parcel-level performance and evaluate capacity estimates using 169 municipalities. We also make our code (which we call SolarMapper), pre-trained models, training data, and predictions publicly available and provide a web portal for interactively inspecting each prediction that was made. Here our findings suggest that traditional performance evaluation of the automated identification of solar PV from satellite imagery may be optimistic due to common limitations in the validation process. The takeaways from this work are intended to inform and catalyze the large-scale practical application of automated solar PV assessment techniques by energy researchers and professionals.

14 SOLAR ENERGY↗

DOE OSTI2022

Methods for estimating X-ray machine output through measurement and simulation

We report ball grid arrays are increasingly being applied in the electronics industry and may require X-ray inspection to ensure the integrity and correct placement of solder pins. However, as the architecture of integrated circuits continues to narrow while simultaneously growing more complex, the risk of electronic failure due to radiation damage increases. While medical X-ray devices have been held to high standards and are repeatedly shown to be well characterized, devices used for electronic inspection are often lacking detailed characterization. This study presents unique methods to solve for important properties in X-ray inspection devices such as source to object distance and energy spectrum. This information can then be applied to Monte Carlo models to achieve better overall dose estimates to electronics, which will lead to superior manufactured products. Since X-ray devices can vary greatly in source characteristics, this work investigates spectral measurement and Monte Carlo representation of three X-ray devices. For a Philips SRO 33 100 medical diagnostic device, the spectral output followed expected trends given by the prediction software SpekCalc and Spektr. For the Dage XD7500NT, direct measurement showed a spectral artifact that through the use of Gafchromic films, was shown to be a contributing effect in the dose output. For the Rad Source RS1800, a high powered irradiation device, direct spectral measurement was not achieved. However, a Monte Carlo model using an assumed spectra was found to match ion chamber measurements to a high degree.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

From Bricks to Clicks: Mapping the White Space in Building Innovation

It is a critical national imperative to transform the buildings sector, yet innovation is impeded by deployment failures that leave promising technologies stranded. Conventional market reports and techno-economic analysis provide an insufficient understanding of markets and resource allocation for emerging building technologies. They omit crucial commercialization factors such as ecosystem maturity and adoption friction, where the coordinated participation of a network of suppliers, contractors, financiers, regulators, and integrators is required to scale solutions. This study addresses these gaps by introducing an evaluation framework grounded in front-line data from six years of the DOE's IMPEL incubator, comprising experience from 300 building-sector innovators and the adjacent, complex ecosystem. Our methodology synthesizes top-down market analysis with bottom-up, practitioner-level data across five megatrends: (M1) Affordable materials and industrialized construction; (M2) Healthy and efficient mechanical systems; (M3) Intelligent building operations; (M4) Buildings as grid assets; and (M5) High-density power and cooling for data centers and therein identify twelve "white space" technology opportunities. Next, we develop a multi-criteria scoring rubric to rank these opportunities based on parameters, i.e., Affordability, Quality of Life, Reliability, and Security, yielding composite ‘Demand’ and ‘Maturity’ indices. Our results indicate that the most significant white spaces may not be incremental products but a new class of ‘Ecosystem Enablers’, such as logistics platforms, orchestration layers, and automated compliance software that solve structural deployment gaps. This paper summarizes this transparent, evidence-based, practitioner-informed evaluation framework for policymakers and investors to re-evaluate policy and resource allocation and unlock scalable market transformation.

Singh, Reshma↗

DOE OSTI2024

Improved Transient Stability Analysis of Multi-loop Droop-Controlled Grid Forming Inverters with Current Reference Limiter

Grid-forming control is regarded as a promising technology for achieving high levels of renewable energy integration in future power grids. To fulfill this role, grid-forming devices is expected to possess fault ride through capability. This paper focuses on the transient stability of a droop control-based grid-forming inverter that utilizes current reference saturation as a current limiting strategy. The objective of this work is to investigate the fault recovery process of the grid-forming inverter and analyze the criterion for transitioning the inverter out of the current limiting mode. The theoretical analysis in this study is validated through time-domain simulations, which compare the results obtained using the circuit current limiting method and the d-axis priority-based current limiting method. Various operational conditions are considered to comprehensively evaluate the performance of the grid-forming inverter during a three-phase short circuit scenario. The simulation results confirm the effectiveness of the proposed theoretical approach.

Lyu, Xue↗

DOE OSTI2023

Expanded verification and validation studies of hypersonic aerodynamics with multiple physics-fidelity models

Hypersonic aerothermodynamics is an important domain of modern multiphysics simulation. The Multi-Fidelity Toolkit is a simulation tool being developed at Sandia National Laboratories to predict aerodynamic properties for compressible flows from a range of physics fidelities and computational speeds. These models include the Reynolds-averaged Navier–Stokes (RANS) equations, the Euler equations with momentum-energy integral technique (MEIT), and modified Newtonian aerodynamics with flat-plate boundary layer (MNA+FPBL) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities, as well as model validation with uncertainty considerations. Code verification activities on the MNA+FPBL model build on previous work by focusing on the viscous portion of the model. Viscous quantities of interest are compared against those from an analytical solution for flat-plate, inclined-plate, and cone geometries. The code verification methodology for the MEIT model is also presented. Test setup and results of code verification tests on the laminar and turbulent models within MEIT are shown. Solution-verification activities include grid-refinement studies on simulations that model the HIFiRE-1 wind tunnel experiments. These experiments are used for validation of all model fidelities. A thorough validation comparison with prediction error and uncertainty is also presented. Three additional HIFiRE-1 experimental runs are simulated in this study, and the solution verification and validation work examines the effects of the associated parameter changes on model performance. Finally, a study is presented that compares the computational costs and fidelities from each of the different models.

42 ENGINEERING↗

NASA NTRS2017

Lessons Learned from Inlet Integration Analysis of NASA's Low Boom Flight Demonstrator

In 2016, NASA's Aeronautics Research Mission Directorate announced the New Aviation Horizons Initiative with a goal of designing/building several X-Planes, including a Low Boom Flight Demonstrator (LBFD). That same year, NASA awarded a contract to Lockheed Martin (LM) to advance the LBFD concept through preliminary design. Several configurations of the LBFD aircraft were analyzed by both LM engineers and NASA researchers. This presentation focuses on some of the CFD simulations that were run by NASA Glenn researchers. NASA's FUN3D V13.1 code was used for all adjoint-based grid refinement studies and Spalart-Allmaras turbulence model was used during adaptation. It was found that adjoint-based grid adaptation did not accurately capture inlet performance for high speed top-aft-mounted propulsion.

Adjoint-based Grid Adaptatio↗

DOE OSTI2020

The Cost of Floating Offshore Wind Energy in California Between 2019 and 2032

California’s energy planning is centered around meeting the emissions reduction and renewable energy requirements of Senate Bill 350 by 2030. However, state power system planning is expected to eventually address California’s requirement to achieve 100% of total retail electricity sales from renewable energy and zero-carbon resources by 2045, as mandated by Senate Bill 100. To comply with these directives, California needs to investigate the further development of energy efficiency, storage, and a diverse range of renewable energy, zero-carbon emission, and transmission resources, including offshore wind. Wind resources off the coast of California have the potential to generate a significant portion of the state’s electric energy as it moves toward a zero-carbon economy and can help diversify its energy mix. Floating offshore wind technology, which is suitable for the deep waters along the California coast, is currently in a precommercial phase, with approximately 84 megawatts (MW) installed worldwide at the end of 2019. Globally there are over 7,000 MW in planning and permitting phases of development, with the first commercial-scale projects expected to be operational in 2024. This study provides site-specific cost and performance data for floating offshore wind to inform California’s long-term energy planning. The identification of new resources to meet California’s policy goals at least cost is part of the Integrated Resource Planning (IRP) process, which is coordinated by the California Public Utilities Commission (CPUC). In 2019–2020 IRP modeling, offshore wind was included for the first time as a candidate resource in some sensitivity cases (CPUC 2019). The data and information presented in this report can be used to update offshore wind inputs in future IRP cycles. The authors conducted a geospatial cost analysis over portions of the offshore wind resource area of California. The analyzed spatial domain includes sites with a mean wind speed of at least 7 meters per second and water depths between 40 meters (m) and 1,300 m. Costs and energy production vary across this analysis domain. We calculated these parameters on a grid layout with over 750 sites, with each site representing a 1,000-MW commercial offshore wind power plant. Levelized cost of energy (LCOE) was calculated at each site over the analysis domain. The resulting variation in LCOE across the analysis domain is illustrated through heat maps in this report. Five study areas were selected within the analysis domain where more detailed cost analysis was conducted and cost parameters, such as annual energy production, capital cost expenditures (CapEx), operational cost expenditures (OpEx), and net capacity factors are reported. These five study areas include Morro Bay, Diablo Canyon, Humboldt, Cape Mendocino and Del Norte.

17 WIND ENERGY↗

DOE OSTI2023

Distribution Grid Impacts of Community Solar [Slides]

Community solar (CS) projects often face uncertain interconnection costs and fees associated with distribution grid infrastructure upgrades required to connect the project. These costs can determine the economic viability of a CS project, but they are difficult to assess. Cost uncertainty can discourage new projects and prevent communities from accessing the benefits of community solar projects. At the same time, CS deployment strategies hold potential to defer or avoid some distribution costs due to new loads. To mitigate CS interconnection costs, it is important to find least-cost combinations of distribution system infrastructure solutions (transformer upgrades, reconductoring, voltage regulators, storage), and to understand how location of CS projects within a feeder impact distribution grid upgrade costs. This study aims to quantify CS impacts on the distribution grid and provide policy and regulatory insights and CS deployment strategies to address them. It is the first analysis that has systematically studied the technical impacts of community solar projects on a wide range of distribution feeders using state-of-the-art optimization and power flow tools. The analysis employs Berkeley Lab’s novel Least-cost Optimal Distribution Grid Expansion (LODGE) model, a deterministic version of the REPAIR model, that optimally upgrades hundreds or even thousands of distribution circuits or feeders. This is the first application of the LODGE model. LODGE finds the least-cost portfolio of traditional distribution system upgrades to integrate CS in combination with alternative solutions, such as utility-owned storage and downsizing CS capacity. Working with a set of least-cost solutions per feeder allows us to benchmark, compare and find techno-economic trends in CS interconnection.

14 SOLAR ENERGY↗

DOE OSTI2026

WINDPROF: Merged Best-Estimate Wind Profile Data – Nantucket (WFIP3 Campaign)

WINDPROF provides 10-minute wind and turbulence profiles, integrating Doppler lidars, wind profiling radars, and sonic anemometers across Northeast U.S. coastal/offshore sites during the WFIP3 campaign. Key data include wind speed, direction, vertical velocity, and turbulence parameters, with standardized quality control (e.g., instrument-specific thresholds and inter-instrument validation). Profiles are interpolated to a height grid (20 m spacing below 100 m; 30 m above) and include comprehensive uncertainty estimates. The Nantucket dataset covers February 2024–September 2025, offering reproducible methods for atmospheric research, model validation, and wind energy studies.

17 WIND ENERGY↗

NASA NTRS2002

A Mechanistic Study of Nucleate Boiling Under Microgravity Conditions

The overall objective of this work is to study nucleate boiling heat transfer under microgravity conditions in such a way that while providing basic knowledge of the phenomena, it also leads to development of simulation models and correlations that can be used as design tools for a wide range of gravity levels. In the study a building block type of approach is used and both pool and low velocity flow boiling are investigated. Starting with experiments using a single bubble, the complexity of the experiments is increased to two or three inline bubbles, to five bubbles placed on a two-dimensional grid. Finally, experiments are conducted where a large number of prescribed cavities nucleate on the heater and when a commercial surface is used. So far experiments have been conducted at earth normal gravity and in the reduced gravity environment of the KC-135 aircraft whereas experiments on the space station are planned. Modeling/complete numerical simulation of the boiling process is an integral part of the total effort. Experiments conducted with single bubbles formed on a nucleation site microfabricated on a polished silicon wafer show that for gravity levels (g) varying from 1.5g(sub e) to 0.01g(sub e), the bubble diameter at departure varies approximately as (g(sub e)/g)(exp 1/2) and the growth period as (g(sub e)/g). When bubbles merge either inline or in a plane, the bubble diameter at departure is found to be smaller than that obtained for a single bubble and shows a weaker dependence on the level of gravity. The possible reason is that as the bubbles merge they create fluid circulation around the bubbles, which in turn induces a lift force that is responsible for the earlier departure of the bubbles. The verification of this proposition is being sought through numerical simulations. There is a merger of two inline, three inline, and several bubbles in a plane in the low gravity environment of the KC-135 aircraft. After merger and before departure, a mushroom type of bubble with several stems attached to the heater surface is clearly evident. Local heat fluxes during growth and departure of a single bubble were also measured. It was found that during most of the growth period of the bubble, generally the wall heat flux decreased with time because of the increased dry area under the bubble. However, the heat flux increased rapidly just prior to departure of the bubble because of the transient conduction into the cold liquid rushing to fill the space vacated by the bubble as the bubble base shrinks. The measured heat fluxes at various radial locations are found to be in qualitative agreement with the numerical predictions. Single bubble studies at earth normal gravity have also been performed on surfaces oriented at different angles to the gravitational acceleration with flow parallel to the surface. It is found that in all cases the bubbles slide along the surface before lift-off from the surface. The lift force generated as a result of the relative motion between the sliding bubbles and the imposed flow is found to play an important role when the normal force due to buoyancy is reduced. An experimental apparatus for the study of the bubble behavior with imposed flow under reduced gravity conditions has been developed and will soon be employed for experiments in the KC-135 aircraft.

Dhir, V. K.↗