Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial forces”

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 127 records · Page 7

Ceramic microstructure and adhesion

When a ceramic is brought into contact with a ceramic, a polymer, or a metal, strong bond forces can develop between the materials. The bonding forces will depend upon the state of the surfaces, cleanliness and the fundamental properties of the two solids, both surface and bulk. Adhesion between a ceramic and another solid are discussed from a theoretical consideration of the nature of the surfaces and experimentally by relating bond forces to interface resulting from solid state contact. Surface properties of ceramics correlated with adhesion include, orientation, reconstruction and diffusion as well as the chemistry of the surface specie. Where a ceramic is in contact with a metal their interactive chemistry and bond strength is considered. Bulk properties examined include elastic and plastic behavior in the surficial regions, cohesive binding energies, crystal structures and crystallographic orientation. Materials examined with respect to interfacial adhesive interactions include silicon carbide, nickel zinc ferrite, manganese zinc ferrite, and aluminum oxide. The surfaces of the contacting solids are studied both in the atomic or molecularly clean state and in the presence of selected surface contaminants.

Buckley, D. H.↗

Gravitational effects in dendritic growth

The theories of diffusion-controlled dendritic crystallization will be reviewed briefly, along with recently published critical experiments on the kinetics and morphology of dendritic growth in pure substances. The influence of the gravitational body force on dendrite growth kinetics will be shown to be highly dependent on the growth orientation with respect to the gravity vector and on the level of the thermal supercooling. In fact, an abrupt transition occurs at a critical supercooling, above which diffusional transport dominates the growth process and below which convective transport dominates. Our most recent work on binary mixtures shows that dilute solute additions influence the crystallization process indirectly, by altering the interfacial stability, rather than by directly affecting the transport mode. Directions for future studies in this field will also be discussed.

Glicksman, M. E.↗

Mechanical Properties of Carbon Fiber Reinforced Composites Exposed to Cryogenic Conditions and Space Radiation via Simulation and Testing

As NASA missions extend beyond low Earth orbit, increasing reliance is placed on carbon fiber reinforced polymer (CFRP) composites for spacecraft structures where mass efficiency, durability, and long-term reliability are critical. In service, these materials are subjected to a combination of ultraviolet radiation, vacuum, ionizing radiation, atomic oxygen, and extreme thermal excursions under sustained mechanical loading. Flight systems such as the Boeing Starliner and SpaceX Dragon employ external composite structures that will experience these environments for extended durations. Although prior spaceflight and ground studies have reported limited changes in bulk mechanical properties, the synergistic effects of these environments on composite microstructure, particularly at the fiber matrix interphase, remain insufficiently characterized and represent a potential qualification and reliability risk. This study investigates the effects of short-term cryogenic exposure on a radiation shielding carbon epoxy composite, SC2020, as a ground-based analog for space relevant thermal extremes. The SC2020 material system has previously flown on the International Space Station under the Materials International Space Station Experiment (MISSE) program. Composite specimens were exposed to liquid nitrogen for 6 and 24 hours and evaluated using a multiscale characterization framework that combined ASTM D3039 tensile testing, Atomic Force Microscopy (AFM) based interphase analysis, and helium gas permeability measurements. Tensile testing showed no statistically significant or permanent degradation in global strength or modulus following cryogenic exposure. In contrast, AFM measurements revealed reductions in interphase modulus, weakened adhesion, and increased nanoscale heterogeneity, indicating localized degradation mechanisms not captured by conventional bulk testing. Gas permeability measurements showed a progressive increase in helium diffusion with exposure duration, consistent with micro-void formation or partial interfacial debonding. The results indicate that cryogenic exposure initiates degradation at the fiber matrix interphase while leaving global mechanical properties largely unchanged over short durations. These findings underscore the importance of multiscale diagnostics for identifying early-stage damage mechanisms that may influence long term performance and qualification margins for spaceflight composite structures. The data presented establish a cryogenic baseline for comparison with forthcoming MISSE flight exposure results and support ongoing NASA Established Program to Stimulate Competitive Research (EPSCoR) efforts aimed at improving composite qualification methodologies, risk assessment, and reliability prediction for space environments.

composite reliability↗

Experimental Investigation of Flow Orientation Effects on Cryogenic Flow Boiling

Cryogenic flow boiling physics is of paramount importance for future cryogenic space applications employing Cryogenic Fluid Management (CFM) technologies. This experimental study investigates the two-phase flow and heat transfer performance of LN 2 flow boiling across five distinct flow orientations: vertical upflow, vertical downflow, horizontal flow, 45° inclined upflow, and 45° inclined downflow. The study employed the steady- state heating method within a circular heated tube featuring an 8.5-mm inner diameter and a 680-mm heated length. High-speed video recordings were utilized to capture two-phase flow patterns and interfacial behaviors across various flow orientations. The experiments covered a wide range of operating conditions, including mass velocities ranging from 351.80 to 1572.77 kg/m 2 s and inlet pressures from 297.14 to 1032.97 kPa, primarily with near-saturated inlet subcooling. Distinct two-phase flow patterns and regime transitions were identified for each flow orientation. Symmetrical flow patterns were evident in vertical orientations, whereas non-vertical orientations exhibited asymmetric flow stratifications, primarily influenced by the buoyancy force in a terrestrial gravity environment. Bubble dynamic parameters were quantified, and bubble collision and dispersion phenomena were visualized. Analyzing heat transfer performance based on local flow boiling curves and variations in heat transfer coefficients (HTCs), it was observed that vertical upflow demonstrated the most enhanced heat transfer performance, while vertical downflow exhibited the lowest. As mass velocity exceeded 830 kg/ m 2 s, the differences in heat transfer among orientations became less distinct, emphasizing the role of flow inertia in mitigating the influence of flow orientation. A direct comparison of microgravity HTC data from a recent study by the authors against the 1-𝘨 e HTC data indicated enhanced heat transfer performance in microgravity for flow orientation configurations in terrestrial gravity. However, this enhancement progressively diminished as the heat flux increased. Distinct HTC trends were observed as mass velocity increased for different flow orientations and gravity levels.

Cryogens↗

A Review on the Coalescence of Confined Drops with a Focus on Scaling Laws for the Growth of the Liquid Bridge

The surface–tension-driven coalescence of drops has been extensively studied because of the omnipresence of the phenomenon and its significance in various natural and engineering systems. When two drops come into contact, a liquid bridge is formed between them and then grows in its lateral dimensions. As a result, the two drops merge to become a bigger drop. The growth dynamics of the bridge are governed by a balance between the driving force and the viscous and inertial resistances of involved liquids, and it is usually represented by power–law scaling relations on the temporal evolution of the bridge dimension. Such scaling laws have been well-characterized for the coalescence of unconfined or freely suspended drops. However, drops are often confined by solid or liquid surfaces and thus are a different shape from spheres, which affects their coalescence dynamics. As such, the coalescence of confined drops poses more complicated interfacial fluid dynamics challenges compared to that of unconfined drops. Although there have been several studies on the coalescence of confined drops, they have not been systematically reviewed in terms of the properties and geometry of the confining surface. Thus, we aim to review the current literature on the coalescence of confined drops in three categories: drop coalescence on a solid surface, drop coalescence on a deformable surface, and drop coalescence between two parallel surfaces with a small gap (i.e., Hele-Shaw cell), with a focus on power–law scaling relations, and to suggest challenges and outlooks for future research on the phenomena.

Sangjin Ryu↗

Thermal buoyancy and Marangoni convection in a two fluid layered system - A numerical study

Thermal buoyancy and surface tension driven convection is numerically investigated in a system with two, immiscible fluids. The geometry investigated, has an open cavity configuration with the lighter fluid situated on top of the heavier fluid, forming a stable layered system. The upper fluid meniscus and the interface are assumed to be flat and undeformable in the calculations. The governing equations and boundary/interface conditions are solved by a control volume based finite difference scheme for two pairs of immiscible fluids; the water-hexadecane system and a so called generic system. The steady state calculations show, that dramatically different flows are predicted when the interfacial tension effects are included or excluded from the system model. These differences are particularly appreciable in surface tension dominated flows, that are typical of microgravity situations. Complex flow patterns, with induced secondary flows are noticed in both the fluids. In general, the overall system heat transfer is found to increase with increases in the thermal buoyancy and surface tension effects, but the behavior of the system flow and thermal fields is not easy to characterize, when different combinations of these forces are considered.

Ramachandran, N.↗

Transient-forced convection film boiling on an isothermal flat plate

An approach for the solution of transient-forced convection film boiling on an isothermal flat plate using the boundary layer model is developed. The similarity variables are used to convert the governing partial differential equations to ordinary ones. The results of numerical solutions of these ordinary equations indicate that the transient process can be classified as one-dimensional conduction, intermediate, and the steady state regions. The time required for the one-dimensional conduction and the time necessary to attain a steady state condition are obtained. The influence of interfacial shear is seen to be negligible while the Prandtl Number and the ratio (C sub p delta T divided by h sub fg times Pr) have major influence. The use of local similarity approximations for the intermediate regime facilitates prediction of complete boundary layer growth. Using the ratio of time at any instant to the steady state time as abscissa, the curves representing the boundary layer growth can be merged into a single mean curve within 5 percent. Further, the analysis shows that the average rate of heat transfer during transient is 50 to 100 percent higher than those at steady state. The average rate of vapor convected away is 10 to 15 percent lower than at steady state while the average rate of accumulation to form the vapor layer is 1 to 14 times larger. Further, the total heat transferred during transient increases and the evaporation decreases for increasing values of C sub p delta T divided by h sub fg times Pr

Nagendra, H. R.↗

The role of material properties in adhesion

When two solid surfaces are brought into contact strong adhesive bond forces can develop between the materials. The magnitude of the forces will depend upon the state of the surfaces, cleanliness and the fundamental properties of the two solids, both surface and bulk. Adhesion between solids is addressed from a theoretical consideration of the electronic nature of the surfaces and experimentally relating bond forces to the nature of the interface resulting from solid state contact. Surface properties correlated with adhesion include, atomic or molecular orientation, reconstruction and segregation as well as the chemistry of the surface specie. Where dissimilar solids are in contact the contribution of each is considered as is the role of their interactive chemistry on bond strength. Bulk properties examined include elastic and plastic behavior in the surficial regions, cohesive binding energies, crystal structure, crystallographic orientation and state. Materials examined with respect to interfacial adhesive interactions include metals, alloys, ceramics, polymers and diamond. They are reviewed both in single and polycrystalline form. The surfaces of the contacting solids are studied both in the atomic or molecularly clean state and in the presence of selected surface contaminants.

Buckley, D. H.↗

Wave effects on ocean-ice interaction in the marginal ice zone

The effects of wave train on ice-ocean interaction in the marginal ice zone are studied through numerical modeling. A coupled two-dimensional ice-ocean model has been developed to include wave effects and wind stress for the predictions of ice edge dynamics. The sea ice model is coupled to the reduced-gravity ocean model through interfacial stresses. The main dynamic balance in the ice momentum is between water-ice stress, wind stress, and wave radiation stresses. By considering the exchange of momentum between waves and ice pack through radiation stress for decaying waves, a parametric study of the effects of wave stress and wind stress on ice edge dynamics has been performed. The numerical results show significant effects from wave action. The ice edge is sharper, and ice edge meanders form in the marginal ice zone owing to forcing by wave action and refraction of swell system after a couple of days. Upwelling at the ice edge and eddy formation can be enhanced by the nonlinear effects of wave action; wave action sharpens the ice edge and can produce ice meandering, which enhances local Ekman pumping and pycnocline anomalies. The resulting ice concentration, pycnocline changes, and flow velocity field are shown to be consistent with previous observations.

Liu, Antony K.↗

Numerical Simulation of Combustion and Extinction of a Solid Cylinder in Low-Speed Cross Flow

The combustion and extinction behavior of a diffusion flame around a solid fuel cylinder (PMMA) in low-speed forced flow in zero gravity was studied numerically using a quasi-steady gas phase model. This model includes two-dimensional continuity, full Navier Stokes' momentum, energy, and species equations with a one-step overall chemical reaction and second-order finite-rate Arrhenius kinetics. Surface radiation and Arrhenius pyrolysis kinetics are included on the solid fuel surface description and a parameter Phi, representing the percentage of gas-phase conductive heat flux going into the solid, is introduced into the interfacial energy balance boundary condition to complete the description for the quasi-steady gas-phase system. The model was solved numerically using a body-fitted coordinate transformation and the SIMPLE algorithm. The effects of varying freestream velocity and Phi were studied. These parameters have a significant effect on the flame structure and extinction limits. Two flame modes were identified: envelope flame and wake flame. Two kinds of flammability limits were found: quenching at low-flow speeds due to radiative loss and blow-off at high flow speeds due to insufficient gas residence time. A flammability map was constructed showing the existence of maximum Phi above which the solid is not flammable at any freestream velocity.

Tien, J. S.↗

Effects of G-Jitter on Interfacial Dynamics of Two Miscible Liquids: Application of MIM

We designed an experiment to examine the effects of g-jitter on mixing of two miscible liquids using the Microgravity Vibration Isolation Mount (MIM). The global bifurcation of the interface was observed with the MIM operating alternatively to either transmit the g-jitter, isolate from the g-jitter or to provide controlled vibration levels with well defined amplitude and frequency content. With the MIM in isolation mode, the interface remains stationary indicating buoyancy induced convection is negligibly small such that mixing occurs via intrinsic mass diffusion without the masking effect of vibration driven convection. Analytical and computational results are in agreement with the experimental findings. Operation of the MIM in forced mode with conditions typical of g-jitter shows that vibration induced convective flows can excite instability mechanisms such as Kelvin-Helmholtz to generate large amplitude quasi-stationary waves oriented vertically for various cases with Stokes-Reynolds number in the range of 0.003 to 0.5. The two and four mode quasi-stationary waves are also predicted with a mathematical model. Though unplanned, the effect of a primary thruster filing was captured and shown to cause a catastrophic bifurcation, enhancing local mass transport. In light of the findings, experiments planned for the International Space Station should consider the potential effects of g-jitter.

Duval, Walter M. B.↗

Separation of crack extension modes in composite delamination problems

This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.

Beuth, Jack L.↗

Effect of Residual Noncondensables on Pressurization and Pressure Control of a Zero-Boil-Off Tank in Microgravity

The Zero-Boil-Off Tank (ZBOT) Experiment is a small-scale experiment that uses a transparent ventless Dewar and a transparent simulant phase-change fluid to study sealed tank pressurization and pressure control with applications to on-surface and in-orbit storage of propellant cryogens. The experiment will be carried out under microgravity conditions aboard the International Space Station in the 2014 timeframe. This paper presents preliminary results from ZBOT's ground-based research that focuses on the effects of residual noncondensable gases in the ullage on both pressurization and pressure reduction trends in the sealed Dewar. Tank pressurization is accomplished through heating of the test cell wall in the wetted and un-wetted regions simultaneously or separately. Pressure control is established through mixing and destratification of the bulk liquid using a temperature controlled forced jet flow with different degrees of liquid jet subcooling. A Two-Dimensional axisymmetric two-phase CFD model for tank pressurization and pressure control is also presented. Numerical prediction of the model are compared to experimental 1g results to both validate the model and also indicate the effect of the noncondensable gas on evolution of pressure and temperature distributions in the ullage during pressurization and pressure control. Microgravity simulations case studies are also performed using the validated model to underscore and delineate the profound effect of the noncondensables on condensation rates and interfacial temperature distributions with serious implications for tank pressure control in reduced gravity.

Marangoni and Natural Convection↗

Oscillatory/Chaotic Thermocapillary Flow Induced by Radiant Heating

There is a continuing need to understand the fluid physics occurring under low gravity conditions in processes such as crystal growth, materials processing, and the movement of bubbles or droplets. The fluid flow in such situations is often caused by a gradient in interfacial tension. If a temperature gradient is created due to a heat source, the resulting flow is called thermocapillary flow, a special case of Marangoni Convection. In this study, an experimental investigation was conducted using silicone oil in cylindrical containers with a laser heat source at the free surface. It was desired to determine the conditions under which steady, axisymmetrical thermocapillary flow becomes unstable and oscillatory three-dimensional flow states develop. The critical Marangoni number for each observed oscillatory state was measured as a function of the container aspect ratio and the dynamic Bond number, a measure of buoyant force versus ii thermocapillary force. Various oscillatory modes were observed during three- dimensional convection, and chaotic flow was reached in one test condition. The critical Marangoni numbers are compared with those measured in previous studies, and the power spectra and phase trajectories of the instantaneous surface temperature distributions are used to characterize the routes of transitions to the chaotic flow state. Results show that only superharmonic modes appear in the routes to chaos while infinite number of subharmonic modes occur in flow transitions for pure Rayleigh convection.

DeWitt, Kenneth J.↗

Development of a 12-Thrust Chamber Kerosene /Oxygen Primary Rocket Sub-System for an Early (1964) Air-Augmented Rocket Ground-Test System

Airbreathing/Rocket combined-cycle, and specifically rocket-based combined- cycle (RBCC), propulsion systems, typically employ an internal engine flow-path installed primary rocket subsystem. To achieve acceptably short mixing lengths in effecting the "air augmentation" process, a large rocket-exhaust/air interfacial mixing surface is needed. This leads, in some engine design concepts, to a "cluster" of small rocket units, suitably arrayed in the flowpath. To support an early (1964) subscale ground-test of a specific RBCC concept, such a 12-rocket cluster was developed by NASA's Marshall Space Flight Center (MSFC). The small primary rockets used in the cluster assembly were modified versions of an existing small kerosene/oxygen water-cooled rocket engine unit routinely tested at MSFC. Following individual thrust-chamber tests and overall subsystem qualification testing, the cluster assembly was installed at the U. S. Air Force's Arnold Engineering Development Center (AEDC) for RBCC systems testing. (The results of the special air-augmented rocket testing are not covered here.) While this project was eventually successfully completed, a number of hardware integration problems were met, leading to catastrophic thrust chamber failures. The principal "lessons learned" in conducting this early primary rocket subsystem experimental effort are documented here as a basic knowledge-base contribution for the benefit of today's RBCC research and development community.

Pryor, D.↗

The Micromechanics of the Moving Contact Line

The proposed research is divided into three components concerned with molecular structure, molecular orientation, and continuum averages of discrete systems. In the experimental program, we propose exploring how changes in interfacial molecular structure generate contact line motion. Rather than rely on the electrostatic and electrokinetic fields arising from the molecules themselves, we augment their interactions by an imposed field at the solid/liquid interface. By controling the field, we can manipulate the molecular structure at the solid/liquid interface. In response to controlled changes in molecular structure, we observe the resultant contact line motion. In the analytical portion of the proposed research we seek to formulate a system of equations governing fluid motion which accounts for the orientation of fluid molecules. In preliminary work, we have focused on describing how molecular orientation affects the forces generated at the moving contact line. Ideally, as assumed above, the discrete behavior of molecules can be averaged into a continuum theory. In the numerical portion of the proposed research, we inquire whether the contact line region is, in fact, large enough to possess a well-defined average. Additionally, we ask what types of behavior distinguish discrete systems from continuum systems. Might the smallness of the contact line region, in itself, lead to behavior different from that in the bulk? Taken together, our proposed research seeks to identify and accurately account for some of the molecular dynamics of the moving contact line, and attempts to formulate a description from which one can compute the forces at the moving contact line.

Han, Minsub↗

Ice banding as a response of the coupled ice-ocean system to temporally varying winds

This study models formation of ice bands in the marginal ice zones. A one-dimensional coupled ice-ocean model is used in which the ice model is coupled to a reduced gravity ocean model through interfacial stresses. The internal ice stresses are important only at high ice concentrations (90-100 percent); otherwise, the main balance for the ice motion is between the air-ice stress and the ice-water stress, i.e., free drift. The drag coefficients were chosen so that the air-ice momentum flux is 3 times greater than the air-ocean momentum flux. Thus the Ekman transport is larger under the ice than in the open water, so that winds parallel to the ice edge, with the ice on the right, produce upwelling. The upwelling simulation was extended to include temporally varying forcing, which was chosen to vary sinusoidally with a 4-day period. This forcing resembles successive cyclone passings perpendicular to the ice edge. When the oceanic upper layer was thin, which means that the dynamics are strongly nonlinear, the ice bands were formed. The up/downwelling signals do not disappear in wind reversals because of nonlinear advection. This leads to convergences and divergences in oceanic and ice velocities that manifest themselves as ice banding. At least one wind reversal is needed to produce one ice band.

Hakkinen, S.↗

Two-layer viscous instability in a rotating couette device

A novel experiment to study the interfacial shear instability between two liquids is described. Density-matched immiscible liquids are confined between concentric cylinders such that the interface is parallel to the cylinder walls. Interfacial waves that develop because of viscosity differences between the shearing fluids are studied as a function of rotation rate and depth ratio using optical techniques. Conditions neutral stability and the most unstable wavenumber agree reasonably well with predictions from linear stability analysis of the Navier-Stokes equations. Illumination using laser sheets allows precise measurement of the interface shape. Future experiments will verify the correctness of weakly nonlinear theories that describe energy transfer and saturation of wave growth by nonlinear effects. Measurements of solitary wave shapes, that occur far above neutral stability, will be compared to similar measurements for systems that have gravity as an important force to determine how gravity effects large disturbances. These results will be used to interpret slug and annular flow data that have been obtained in other mu g studies.

Gallagher, Christopher↗