Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffuse interface methods”

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 19 records

Diffuse-Interface Methods in Fluid Mechanics

The authors review the development of diffuse-interface models of hydrodynamics and their application to a wide variety of interfacial phenomena. The authors discuss the issues involved in formulating diffuse-interface models for single-component and binary fluids. Recent applications and computations using these models are discussed in each case. Further, the authors address issues including sharp-interface analyses that relate these models to the classical free-boundary problem, related computational approaches to describe interfacial phenomena, and related approaches describing fully-miscible fluids.

SURFACE DIFFUSION↗

Microgravity Crystallization of Alpha-Crustacyanin Onboard the Unmanned Carrier, EURECA

alpha-Crustacyanin, the lobster carapace astaxanthin-protein, was crystallized using the European Space Agency's (ESA) automated Protein Crystallization Facility (PCF) which flew onboard the unmanned EUropean REtrievable CArrier (EURECA). A free interface linear, liquid - liquid diffusion, method was used. Crystals grew larger and thicker in the microgravity case compared to the biggest crystals grown on earth. Video observation on EURECA revealed variations in crystal sizes through-out the reactor neatly correlated with depletion of this coloured protein from the solution. The video observations most importantly revealed no visible movement of crystals over the initial 7 weeks of the experiment, although an obvious temperature induced jump occurred at that time in a mission spanning 11 months. An important observation from this mission, over the first 7 weeks, of completely stationary crystal growth contrasts with crystal motions viewed on manned microgravity missions, even using linear liquid - liquid geometries, and much shorter flights (eg. 12 to 16 days).

Boggon, T. J.↗

Interfacial geometry and D-variation effects in two-phase systems

Numerical solutions of the governing diffusion equation for two-phase concentration dependent diffusion coefficients are examined. Solutions were also calculated for planar, cylindrical, and spherical geometries to compare the effect of interface geometries with those caused by concentration-dependent diffusion coefficients, and two methods of averaging D were considered to determine the best averaging method for different types of D-variations. The effects of interface-location criteria on mass conservation and convergence of interface location, diffusion coefficient variation in the alpha and beta-phases of a two-phase binary alloy system, effect of D(alpha) variation in a cylindrical couple on beta-phase thickness, and geometry and D-variation effects on the degree of homogenization were determined. It is concluded that typical D(alpha)-variations can have a greater influence on the kinetics of interdiffusion than the geometry.

Tenney, D. R.↗

Protein-crystal growth experiment (planned)

To evaluate the effectiveness of a microgravity environment on protein crystal growth, a system was developed using 5 cubic feet Get Away Special payload canister. In the experiment, protein (myoglobin) will be simultaneously crystallized from an aqueous solution in 16 crystallization units using three types of crystallization methods, i.e., batch, vapor diffusion, and free interface diffusion. Each unit has two compartments: one for the protein solution and the other for the ammonium sulfate solution. Compartments are separated by thick acrylic or thin stainless steel plates. Crystallization will be started by sliding out the plates, then will be periodically recorded up to 120 hours by a still camera. The temperature will be passively controlled by a phase transition thermal storage component and recorded in IC memory throughout the experiment. Microgravity environment can then be evaluated for protein crystal growth by comparing crystallization in space with that on Earth.

Fujita, S.↗

Mathematical model for the Bridgman-Stockbarger crystal growing system

In a major technical breakthrough, a computer model for Bridgman-Stockbarger crystal growth was developed. The model includes melt convection, solute effects, thermal conduction in the ampule, melt, and crystal, and the determination of the curved moving crystal-melt interface. The key to the numerical method is the use of a nonuniform computational mesh which moves with the interface, so that the interface is a mesh surface. In addition, implicit methods are used for advection and diffusion of heat, concentration, and vorticity, for interface movement, and for internal gracity waves. This allows large time-steps without loss of stability or accuracy. Numerical results are presented for the interface shape, temperature distribution, and concentration distribution, in steady-state crystl growth. Solutions are presented for two test cases using water, with two different salts in solution. The two diffusivities differ by a factor of ten, and the concentrations differ by a factor of twenty.

Roberts, G. O.↗

Intermodule ventilation studies for the Space Station

This paper examines the ability of the Space Station intermodule ventilation system to maintain centralized control of CO2 removal and O2 supply. The resulting concentration gradients that will arise are calculated by assuming steady state, ideal gas, isothermal conditions, and perfect mixing of air within and between the pressurized elements. In order to estimate the degree of mixing actually obtained for a given ventilation scheme, a program has been developed based on a potential flow solution technique. Preliminary results from this study indicate that substantial short circuiting and recirculation air flow patterns could arise if a simple duct and diffuser air exchange method at the docking port interface were employed.

Davis, Roy G.↗

Use of Plastic Capillaries for Macromolecular Crystallization

Methods of crystallization of biomolecules in plastic capillaries (Nalgene 870 PFA tubing) are presented. These crystallization methods used batch, free-interface liquid- liquid diffusion alone, or a combination with vapor diffusion. Results demonstrated growth of crystals of test proteins such as thaumatin and glucose isomerase, as well as protein studied in our laboratory such dihydrolipoamide dehydrogenase. Once the solutions were loaded in capillaries, they were stored in the tubes in frozen state at cryogenic temperatures until the desired time of activation of crystallization experiments.

Potter, Rachel R.↗

Instability Mechanisms of Thermally-Driven Interfacial Flows in Liquid-Encapsulated Crystal Growth

During the past year, a great deal of effort was focused on the enhancement and refinement of the computational tools developed as part of our previous NASA grant. In particular, the interface mollification algorithm developed earlier was extended to incorporate the effects of surface-rheological properties in order to allow the study of thermocapillary flows in the presence of surface contamination. These tools will be used in the computational component of the proposed research in the remaining years of this grant. A detailed description of the progress made in this area is provided elsewhere. Briefly, the method developed allows for the convection and diffusion of bulk-insoluble surfactants on a moving and deforming interface. The novelty of the method is its grid independence: there is no need for front tracking, surface reconstruction, body-fitted grid generation, or metric evaluations; these are all very expensive computational tasks in three dimensions. For small local radii of curvature there is a need for local grid adaption so that the smearing thickness remains a small fraction of the radius of curvature. A special Neumann boundary condition was devised and applied so that the calculated surfactant concentration has no variations normal to the interface, and it is hence truly a surface-defined quantity. The discretized governing equations are solved subsequently using a time-split integration scheme which updates the concentration and the shape successively. Results demonstrate excellent agreement between the computed and exact solutions.

Haj-Hariri, Hossein↗

Turbulent Mixing of Multiphase Flow

Thus we conduct numerical simulations of multiphase fluids stirred by two-dimensional turbulence to assess the possibility of self-similar drop size distribution in turbulence. In our turbulence simulations, we also explore the non-diffusive limit, where molecular mobility for the interface is vanishing. Special care is needed to transport the non-diffusive interface. Numerically, we use the particle level set method to evolve the interface. Instead of using the usual methods to calculate the surface tension force from the level set function, we reconstruct the interface based on phase- field modeling, and calculate the continuum surface tension forcing from the reconstructed interface.

Young, Y.-N.↗

Two-Flux Green's Function Analysis for Transient Spectral Radiation in a Composite

An analysis is developed for obtaining transient temperatures in a two-layer semitransparent composite with spectrally dependent properties. Each external boundary of the composite is subjected to radiation and convection. The two-flux radiative transfer equations are solved by deriving a Green's function. This yields the local radiative heat source needed to numerically solve the transient energy equation. An advantage of the two-flux method is that isotropic scattering is included without added complexity. The layer refractive indices are larger than one. This produces internal reflections at the boundaries and the internal interface; the reflections are assumed diffuse. Spectral results using the Green's function method are verified by comparing with numerical solutions using the exact radiative transfer equations. Transient temperature distributions are given to illustrate the effect of radiative heating on one side of a composite with external convective cooling. The protection of a material from incident radiation is illustrated by adding scattering to the layer adjacent to the radiative source.

Siegel, Robert↗

Efficient estimation of diffusion during dendritic solidification

A very efficient finite difference method has been developed to estimate the solute redistribution during solidification with diffusion in the solid. This method is validated by comparing the computed results with the results of an analytical solution derived by Kobayashi (1988) for the assumptions of a constant diffusion coefficient, a constant equilibrium partition ratio, and a parabolic rate of the advancement of the solid/liquid interface. The flexibility of the method is demonstrated by applying it to the dendritic solidification of a Pb-15 wt pct Sn alloy, for which the equilibrium partition ratio and diffusion coefficient vary substantially during solidification. The fraction eutectic at the end of solidification is also obtained by estimating the fraction solid, in greater resolution, where the concentration of solute in the interdendritic liquid reaches the eutectic composition of the alloy.

Yeum, K. S.↗

A numerical treatment of geodynamic viscous flow problems involving the advection of material interfaces

Effective numerical treatment of multicomponent viscous flow problems involving the advection of sharp interfaces between materials of differing physical properties requires correction techniques to prevent spurious diffusion and dispersion. We develop a particular algorithm, based on modern shock-capture techniques, employing a two-step nonlinear method. The first step involves the global application of a high-order upwind scheme to a hyperbolic advection equation used to model the distribution of distinct material components in a flow field. The second step is corrective and involves the application of a global filter designed to remove dispersion errors that result from the advection of discontinuities (e.g., material interfaces) by high-order, minimally dissipative schemes. The filter introduces no additional diffusion error. Nonuniform viscosity across a material interface is allowed for by the implementation of a compositionally weighted-inverse interface viscosity scheme. The combined method approaches the optimal accuracy of modern shock-capture techniques with a minimal increase in computational time and memory. A key advantage of this method is its simplicity to incorporate into preexisting codes be they finite difference, element, or volume of two or three dimensions.

Lenardic, A.↗

Preparation and Analysis of RNA Crystals

The crystallization of RiboNucleic Acids (RNA) was studied from the standpoint of mechanisms of crystal growth in three tasks: (1) preparation of high-quality crystals of oligonuclotides for X-ray diffraction, (2) finding pathways to the growth of high-quality crystals for X-ray diffraction and (3) investigation of mechanisms of action of inertial acceleration on crystal growth. In these tasks: (1) RNA crystals were prepared and studied by X-ray diffraction; (2) a pathway to high-quality crystals was discovered and characterized; a combination of kinetic and equilibrium factors could be optimized as described below; and (3) an interplay between purity and gravity was found in a combination of space and ground experiments with nucleic acids and proteins. Most significantly, the rate of concentration of precipitant and RNA can be controlled by membrane-based methods of water removal or by diffusion of multivalent cations across an interface stabilized by a membrane. Oligonucleotide solutions are electrokinetically stabilized colloids, and crystals can form by the controlled addition of multivalent cations.

Todd, Paul↗

Cellular interface morphologies in directional solidification. III - The effects of heat transfer and solid diffusivity

The shape and stability of two-dimensional finite-amplitude cellular interfaces arising during directional solidification are compared for several solidification models that account differently for latent heat released at the interface, unequal thermal conductivities of melt and solid, and solute diffusivity in the solid. Finite-element analysis and computer-implemented perturbation methods are used to analyze the families of steadily growing cellular forms that evolve from the planar state. In all models a secondary bifurcation between different families of finite-amplitude cells exists that halves the spatial wavelength of the stable interface. The quantitative location of this transition is very dependent on the details of the model. Large amounts of solute diffusion in the solid retard the growth of large-amplitude cells.

Ungar, Lyle H.↗

Shadowgraph Study of Gradient Driven Fluctuations

A fluid or fluid mixture, subjected to a vertical temperature and/or concentration gradient in a gravitational field, exhibits greatly enhanced light scattering at small angles. This effect is caused by coupling between the vertical velocity fluctuations due to thermal energy and the vertically varying refractive index. Physically, small upward or downward moving regions will be displaced into fluid having a refractive index different from that of the moving region, thus giving rise to the enhanced scattering. The scattered intensity is predicted to vary with scattering wave vector q, as q(sup -4), for sufficiently large q, but the divergence is quenched by gravity at small q. In the absence of gravity, the long wavelength fluctuations responsible for the enhanced scattering are predicted to grow until limited by the sample dimensions. It is thus of interest to measure the mean-squared amplitude of such fluctuations in the microgravity environment for comparison with existing theory and ground based measurements. The relevant wave vectors are extremely small, making traditional low-angle light scattering difficult or impossible because of stray elastically scattered light generated by optical surfaces. An alternative technique is offered by the shadowgraph method, which is normally used to visualize fluid flows, but which can also serve as a quantitative tool to measure fluctuations. A somewhat novel shadowgraph apparatus and the necessary data analysis methods will be described. The apparatus uses a spatially coherent, but temporally incoherent, light source consisting of a super-luminescent diode coupled to a single-mode optical fiber in order to achieve extremely high spatial resolution, while avoiding effects caused by interference of light reflected from the various optical surfaces that are present when using laser sources. Results obtained for a critical mixture of aniline and cyclohexane subjected to a vertical temperature gradient will be presented. The sample was confined between two horizontal parallel sapphire plates with a vertical spacing of 1 mm. The temperatures of the sapphire plates were controlled by independent circulating water loops that used Peltier devices to add or remove heat from the room air as required. For a mixture with a temperature gradient, two effects are involved in generating the vertical refractive index gradient, namely thermal expansion and the Soret effect, which generates a concentration gradient in response to the applied temperature gradient. For the aniline/cyclohexane system, the denser component (aniline) migrates toward the colder surface. Consequently, when heating from above, both effects result in the sample density decreasing with altitude and are stabilizing in the sense that no convective motion occurs regardless of the magnitude of the applied temperature gradient. The Soret effect is strong near a binary liquid critical point, and thus the dominant effect is due to the induced concentration gradient. The results clearly show the divergence at low q and the predicted gravitational quenching. Results obtained for different applied temperature gradients at varying temperature differences from the critical temperature, clearly demonstrate the predicted divergence of the thermal diffusion ratio. Thus, the more closely the critical point is approached, the smaller becomes the temperature gradient required to generate the same signal. Two different methods have been used to generate pure concentration gradients. In the first, a sample cell was filled with a single fluid, ethylene glycol, and a denser miscible fluid, water, was added from below thus establishing a sharp interface to begin the experiment. As time went on the two fluids diffused into each other, and large amplitude fluctuations were clearly observed at low q. The effects of gravitational quenching were also evident. In the second method, the aniline/cyclohexane sample was used, and after applying a vertical temperature gradient for several hours, the top and bottom temperatures were set equal and the thermal gradient died on a time scale of seconds, leaving the Soret induced concentration gradient in place. Again, large-scale fluctuations were observed and died away slowly in amplitude as diffusion destroyed the initial concentration gradient.

Cannell, David↗

An extended classical solution of the droplet growth problem

Problems of applying the classical kinetic theory to the growth of small droplets from vapor are examined. A solution for the droplet growth equation is derived which is based on the assumption of a diffusive field extending to the drop surface. The method accounts for partial thermal and mass accommodation at the interface and the kinetic limit to the mass and heat fluxes, and it avoids introducing the artifact of a discontinuity in the thermal and vapor field near the droplet. Consideration of the environmental fields in spherical geometry utilizing directional fluxes yields boundary values in terms of known parameters and a new Laplace transform integral.

Anderson, B. J.↗

Brazing method produces solid-solution bond between refractory metals

Brazing two refractory metals by diffusion bonding minimizes distortion and avoids excessive grain growth in the metals. This method requires the selection of an interface metal that forms intermediate low-melting eutectics or solid solutions with the metals to be brazed.

Source record↗

Analytical and Numerical Solution for a Solidifying Liquid Alloy Slab

Numerical and analytical solutions are presented for the temperature and concentration distributions during the solidification of a binary liquid alloy slab. The slab is taken to be of a finite depth but infinite in the horizontal direction. The solidification process is started by withdrawing a fixed amount of heat from the lower surface of the slab. The upper surface of the slab is subjected to both radiation and convective conditions. The solution gives the concentration and temperature profiles and the interface position as a function of time. Due to the smallness of the mass diffusion coefficient in the solid, the numerical solution method breaks down whenever the ratio of the diffusivities in the solid and the liquid falls below a certain value. An analytical method is developed which gives accurate solution for any value of the diffusivity ratio.

Antar, B. N.↗