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 109 records · Page 6

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗

Homogeneous crystal nucleation in binary metallic melts

A method for calculating the homogeneous crystal nucleation frequency in binary metallic melts is developed. The free energy of crystallization is derived from regular solution models for the liquid and solid and is used, together with model-based estimates of the interfacial tension, to calculate the nucleation frequency from the classical theory. The method can account for the composition dependence of the maximum undercooling observed in a number of experiments on small droplet dispersions. It can also be used to calculate the driving force for crystal growth and to obtain more precise estimates of the homogeneous crystal nucleation frequency in glass-forming alloys. This method, although approximate, is simple to apply, and requires only knowledge of the phase diagram and a few readily available thermodynamic quantities as input data.

Thompson, C. V.↗

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗

Li-Doped Ionic Liquid Electrolytes: From Bulk Phase to Interfacial Behavior

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries. We present an extensive computational analysis supported by experimental comparisons of the bulk and interfacial properties of a representative set of these electrolytes as a function of Li-salt doping. We begin by investigating the bulk electrolyte using quantum chemistry and ab initio molecular dynamics to elucidate the solvation structure of Li(+). MD simulations using the polarizable force field of Borodin and coworkers were then performed, from which we obtain an array of thermodynamic and transport properties. Excellent agreement is found with experiments for diffusion, ionic conductivity, and viscosity. Combining MD simulations with electronic structure computations, we computed the electrochemical window of the electrolytes across a range of Li(+)-doping levels and comment on the role of the liquid environment. Finally, we performed a suite of simulations of these Li-doped electrolytes at ideal electrified interfaces to evaluate the differential capacitance and the equilibrium Li(+) distribution in the double layer. The magnitude of differential capacitance is in good agreement with our experiments and exhibits the characteristic camel-shaped profile. In addition, the simulations reveal Li(+) to be highly localized to the second molecular layer of the double layer, which is supported by additional computations that find this layer to be a free energy minimum with respect to Li(+) translation.

molecular dynamics↗

The Transient Dendritic Solidification Experiment (TDSE)

Dendritic solidification is a common mode of solidification. It is also an important model problem in nonequilibrium physics and pattern formation physics. Current theories couple the transfer of latent heat with selection mechanisms at the interface. Measurements of succinonitrile (SCN) dendrites in microgravity show reasonable agreement between heat transfer predictions and experiment. However, data and analysis for assessing interfacial physics theories are less definitive. We are studying, and will present data on, transient effects in dendritic growth of SCN. We employ the Clapeyron pressure/melting temperature effect to make a rapid change in a sample's hydrostatic pressure, and thereby rapidly change the specimen's melting temperature, forcing the dendrite to select a new steady-state. These initial measurements show some surprising and non-intuitive effects.

Koss, M. B.↗

Dissipating Step Bunches during Crystallization under Transport Control

In studies of crystal formation by the generation and spreading of layers, equidistant step trains are considered unstable---bunches and other spatiotemporal patterns of the growth steps are viewed as ubiquitous. We provide an example to the opposite. We monitor the spatiotemporal dynamics of steps and the resulting step patterns during crystallization of the proteins ferritin and apoferritin using the atomic force microscope. The variations in step velocity and density are not correlated, indicating the lack of a long-range attraction between the steps. We show that (i) because of its coupling to bulk transport, nucleation of new layers is chaotic and occurs at the facet edges, where the interfacial supersaturation is higher; (ii) step bunches self-organize via the competition for supply from the solution; and, (iii) bunches of weakly interacting steps decay as they move along the face. Tests by numerical modeling support the conclusions about the mechanisms underlying our observations. The results from these systems suggest that during crystallization controlled by transport, with weakly or noninteracting growth steps, the stable kinetic state of the surface is an equidistant step train, and step bunches only arise during nucleation of new layers. Since nucleation only occurs at a few sites on the surface, the surface morphology may be controllably patterned or smoothened by locally controlling nucleation.

Lin, Hong↗

Data and Modeling of Dendrites Subject to A Step Change in Pressure (TDSE)

There is considerable interest in dendritic solidification because of the influence dendrites have in the determination of microstructure, and thereby in the physical properties of cast metals and alloys. Current theories and models of dendritic growth generally couple diffusion effects in the melt with the physics of the interface. Data and subsequent analysis prior of the tip growth speed and radii of thermal succinonitrile dendrites in the near-convection free, on-orbit, free-fall environment demonstrate that these theories yield predictions that are reasonably in agreement with the results of experiment. However, data and analysis for assessing the interfacial physics component of theory are not sufficiently detailed or definitive. To study fundamental aspects of dendritic interface stability, we are measuring and modeling the kinetics and morphology of dendrites as they evolve from one well-defined steady state at a pre-set supercooling, through a transient stage, to a different well-defined steady state. More specifically, we subject succinonitrile dendrites, growing under steady-state conditions, to a rapid change in pressure. This leads to a rapid change in thermal driving force from the corresponding change in both the equilibrium melting temperature due to the Clapeyron effect, and a change in the far-field temperature due to adiabatic temperature changes in the bulk liquid and solid. Subsequently, we observe transformations from a well-characterized initial state into a new steady-state. Initial data reveal that the dendrite tip velocity changes almost as fast as the pressure charges, while the tip radius changes occur more slowly, taking from 10 60 seconds depending on the size of the step change and the final supercooling. Computer modeling of this process shows both agreements and disagreements with the experimental data. In making these observations and measurements, we are gaining new understandings of interfacial dynamics and state-selection physics.

Koss, Matthew B.↗

Materials Science Experiments on the International Space Station

The Performance Goal for NASA's Microgravity Materials Science Program reads "Use microgravity to establish and improve quantitative and predictive relationships between the structure, processing and properties of materials." The advent of the International Space Station will open up a new era in Materials Science Research including the ability to perform long term and frequent experiments in microgravity. As indicated the objective is to gain a greater understanding of issues of materials science in an environment in which the force of gravity can be effectively switched off. Thus gravity related issues of convection, buoyancy and hydrostatic forces can be reduced and the science behind the structure/processing/properties relationship can more easily be understood. The specific areas of research covered within the program are (1) the study of Nucleation and Metastable States, (2) Prediction and Control of Microstructure (including pattern formation and morphological stability), (3) Phase Separation and Interfacial Stability, (4) Transport Phenomena (including process modeling and thermophysical properties measurement), and (5) Crystal Growth, and Defect Generation and Control. All classes of materials, including metals and alloys, glasses and ceramics, polymers, electronic materials (including organic and inorganic single crystals), aerogels and nanostructures, are included in these areas. The principal experimental equipment available to the materials scientist on the International Space Station (ISS) will be the Materials Science Research Facility (MSRF). Each of these systems will be accommodated in a single ISS rack, which can operate autonomously, will accommodate telescience operations, and will provide real time data to the ground. Eventual plans call for three MSRF racks, the first of which will be shared with the European Space Agency (ESA). Under international agreements, ESA and other partners will provide some of the equipment, while NASA covers launch and integration costs. The MSRF facilities will include modular components, which can be exchanged to provide inserts specifically matched to the engineering requirements of the particular Principal Investigator. To defray costs and avoid duplication of engineering effort NASA is also pursuing the possibility of using facilities provided by international partners. By this means it is anticipated that all of the types of research outlined in the previous paragraph can be done on the ISS.

Gillies, Donald C.↗

Interfacial Structure and Capacitance of Li-doped Ionic Liquid Electrolytes from Molecular Simulations

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries, supercapacitors, and hybrid energy storage devices. Though Li-salt is often present in these systems, its influence on interfacial properties is largely uncharacterized. We, thereby, present an extensive computational analysis, supported by experimental comparisons, of the properties of a representative set of these electrolytesat an ideal carbon interface as a function of Li-salt doping and voltage. We have performed polarizable molecular (MD) dynamics simulations, using the APPLEP force field, to evaluate electric double layer (EDL) capacitance and distribution of Li+ in the EDL. Differential capacitance exhibits the characteristic camel profile and is insensitive to Li-doping. Li+ localizes in the second molecular layer of the EDL, which is a result of confinement from free energy barriers associated with ion layering. Joint MDelectronic structure computations show the electrochemical window of the electrolytes to be a weak function of Li-doping. Estimates of supercapacitor specific energy are made using the computed window and capacitance. The magnitude and trends in specific energy are in good agreement with experiment.

batteries↗

In-Situ Ice Adhesion Testing using the Deformed Skin Adhesion Test

There exists a plethora of ice adhesion testing methods such as centrifuge-based methods, push/pull methods, and lap joint shear tests. However, these methods often cannot be done in-situ during an icing spray and require researcher handling/preparation prior to testing. Additionally, many of these methods test ice accretions that are not representative of the ice shapes that grow on airframes which can have unintended consequences such as edge or corner stress concentrations. In order to mitigate the issues present with other ice adhesion tests a novel, hands-free, in-situ test method imbedded into the leading edge of the airfoil has been developed and tested in the NASA Glenn Research Center Icing Research Tunnel. The mechanics of the new test method are discussed with finite element simulations being used to highlight the fundamental mechanics of this test method. The testing procedures used during the experimental test campaign are discussed along with the icing cloud test matrix. Analysis of the force-displacement results highlights the repeatability of the test method with regards to both sample behavior across several days of testing and different ice shapes. The results from first run of the day showed a significant discrepancy from subsequent runs and differential image analysis was used to determine differences in the interfacial bonding states between runs.

Ice Adhesion↗

In-Situ Ice Adhesion Testing using the Deformed Skin Adhesion Test

There exists a plethora of ice adhesion testing methods such as centrifuge-based methods, push/pull methods, and lap joint shear tests. However, these methods often cannot be done in-situ during an icing spray and require researcher handling/preparation prior to testing. Additionally, many of these methods test ice accretions that are not representative of the ice shapes that grow on airframes which can have unintended consequences such as edge or corner stress concentrations. In order to mitigate the issues present with other ice adhesion tests a novel, hands-free, in-situ test method imbedded into the leading edge of the airfoil has been developed and tested in the NASA Glenn Research Center Icing Research Tunnel. The mechanics of the new test method are discussed with finite element simulations being used to highlight the fundamental mechanics of this test method. The testing procedures used during the experimental test campaign are discussed along with the icing cloud test matrix. Analysis of the force-displacement results highlights the repeatability of the test method with regards to both sample behavior across several days of testing and different ice shapes. The results from first run of the day showed a significant discrepancy from subsequent runs and differential image analysis was used to determine differences in the interfacial bonding states between runs.

Ice Adhesion↗

A tire contact solution technique

An efficient method for calculating the contact boundary and interfacial pressure distribution was developed. This solution technique utilizes the discrete Fourier transform to establish an influence coefficient matrix for the portion of the pressurized tire surface that may be in the contact region. This matrix is used in a linear algebra algorithm to determine the contact boundary and the array of forces within the boundary that are necessary to hold the tire in equilibrium against a specified contact surface. The algorithm also determines the normal and tangential displacements of those points on the tire surface that are included in the influence coefficient matrix. Displacements within and outside the contact region are calculated. The solution technique is implemented with a finite-element tire model that is based on orthotropic, nonlinear shell of revolution elements which can respond to nonaxisymmetric loads. A sample contact solution is presented.

Tielking, J. T.↗

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗

Thermocapillary flow with evaporation and condensation at low gravity. Part 2: Deformable surface

The free surface behavior of a volatile wetting liquid at low gravity is studied using scaling and numerical techniques. An open cavity model, which was applied in part 1 to investigate fluid flow and heat transfer in non-deforming pores, is used to evaluate the influence of convection on surface morphology with length scales and subcooling/superheating limits of 1 less than or equal to D less than or equal to 10(exp 2) microns and approximately 1 K, respectively. Results show that the menisci shapes of highly wetting fluids are sensitive to thermocapillary flow and to a lesser extent the recoil force associated with evaporation and condensation. With subcooling, thermocapillarity produces a suction about the pore centerline that promotes loss of mechanical equilibrium, while condensation exerts an opposing force that under some conditions offsets this destabilizing influence. With superheating, thermocapillarity and evaporation act in the same direction and mutually foster surface stability. All of these trends are magnified by high capillary and Biot numbers, and the stronger circulation intensities associated with small contact angles. These phenomena strongly depend on the thermal and interfacial equilibrium between the liquid and vapor, and have important ramifications for systems designed to maintain a pressure differential across a porous surface.

Schmidt, G. R.↗

Effect of electronic structure of the diamond surface on the strength of the diamond-metal interface

A diamond surface undergoes a transformation in its electronic structure by a vacuum anneal at approximately 900 C. The polished surface has no electronic states in the band gap, whereas the annealed surface has both occupied and unoccupied states in the and gap and exhibits some electrical conductivity. The effect of this transformation on the strength of the diamond metal interface was investigated by measuring the static friction force of an atomically clean meta sphere on a diamond flat in ultrahigh vacuum. It was found that low friction (weak bonding) is associated with the diamond surface devoid of gap states whereas high friction (strong bonding) is associated with the diamond surface with gap states. Exposure of the annealed surface to excited hydrogen also leads to weak bonding. The interfacial bond is discussed in terms of interaction of the metal conduction band electrons with the band gap states on the diamond surface. Effects of surface electrical conductivity on the interfacial bond are also be considered.

Pepper, S. V.↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

Theory of hydrophobicity: transient cavities in molecular liquids

Observation of the size distribution of transient cavities in computer simulations of water, n-hexane, and n-dodecane under benchtop conditions shows that the sizes of cavities are more sharply defined in liquid water but the most-probable-size cavities are about the same size in each of these liquids. The calculated solvent atomic density in contact with these cavities shows that water applies more force per unit area of cavity surface than do the hydrocarbon liquids. This contact density, or "squeezing" force, reaches a maximum near cavity diameters of 2.4 angstroms. The results for liquid water are compared to the predictions of simple theories and, in addition, to results for a reference simple liquid. The numerical data for water at a range of temperatures are analyzed to extract a surface free energy contribution to the work of formation of atomic-size cavities. Comparison with the liquid-vapor interfacial tensions of the model liquids studied here indicates that the surface free energies extracted for atomic-size cavities cannot be accurately identified with the macroscopic surface tensions of the systems.

NASA Center ARC↗

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.↗