Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FRACTALS”

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 289 records · Page 16

Soot and Radiation Measurements in Microgravity Jet Diffusion Flames

The subject of soot formation and radiation heat transfer in microgravity jet diffusion flames is important not only for the understanding of fundamental transport processes involved but also for providing findings relevant to spacecraft fire safety and soot emissions and radiant heat loads of combustors used in air-breathing propulsion systems. Our objectives are to measure and model soot volume fraction, temperature, and radiative heat fluxes in microgravity jet diffusion flames. For this four-year project, we have successfully completed three tasks, which have resulted in new research methodologies and original results. First is the implementation of a thermophoretic soot sampling technique for measuring particle size and aggregate morphology in drop-tower and other reduced gravity experiments. In those laminar flames studied, we found that microgravity soot aggregates typically consist of more primary particles and primary particles are larger in size than those under normal gravity. Comparisons based on data obtained from limited samples show that the soot aggregate's fractal dimension varies within +/- 20% of its typical value of 1.75, with no clear trends between normal and reduced gravity conditions. Second is the development and implementation of a new imaging absorption technique. By properly expanding and spatially-filtering the laser beam to image the flame absorption on a CCD camera and applying numerical smoothing procedures, this technique is capable of measuring instantaneous full-field soot volume fractions. Results from this technique have shown the significant differences in local soot volume fraction, smoking point, and flame shape between normal and reduced gravity flames. We observed that some laminar flames become open-tipped and smoking under microgravity. The third task we completed is the development of a computer program which integrates and couples flame structure, soot formation, and flame radiation analyses together. We found good agreements between model predictions and experimental data for laminar and turbulent flames under both normal and reduced gravity. We have also tested in the laboratory the techniques of rapid-insertion fine-wire thermocouples and emission pyrometry for temperature measurements. These techniques as well as laser Doppler velocimetry and spectral radiative intensity measurement have been proposed to provide valuable data and improve the modeling analyses.

Ku, Jerry C.↗

An Analysis of Gamma-ray Burst Time Profiles from the Burst and Transient Source Experiment

This proposal requested funding to measure the durations of gamma-ray bursts (GRB) in the 4B catalog as well as to study the structure of GRB time profiles returned by the Burst And Transient Source Experiment (BATSE) on board the Compton Gamma-Ray Observatory. The duration (T90) was to be measured using the same techniques and algorithms developed by the principal investigator for the 3B data. The profile structure studies fall into the two categories of variability and fractal analyses.

Lestrade, John Patrick↗

Resistance to Rolling in the Adhesive Contact of Two Elastic Spheres

For the stability of agglomerates of micron sized particles it is of considerable importance to study the effects of tangential forces on the contact of two particles. If the particles can slide or roll easily over each other, fractal structures of these agglomerates will not be stable. We use the description of contact forces by Johnson, Kendall and Roberts, along with arguments based on the atomic structure of the surfaces in contact, in order to calculate the resistance to rolling in such a contact. It is shown that the contact reacts elastically to torque forces up to a critical bending angle. Beyond that, irreversible rolling occurs. In the elastic regime, the moment opposing the attempt to roll is proportional to the bending angle and to the pull-off force P(sub c). Young's modulus of the involved materials has hardly any influence on the results. We show that agglomerates of sub-micron sized particles will in general be quite rigid and even long chains of particles cannot be bent easily. For very small particles, the contact will rather break than allow for rolling. We further discuss dynamic properties such as the possibility of vibrations in this degree of freedom and the typical amount of rolling during a collision of two particles.

Dominik, C.↗

Resistance to Rolling in the Adhesive Contact of Two Elastic Spheres

For the stability of agglomerates of micron sized particles it is of considerable importance to study the effects of tangential forces on the contact of two particles. If the particles can slide or roll easily over each other, fractal structures of these agglomerates will not be stable. We use the description of contact forces by Johnson, Kendall and Roberts, along with arguments based on the atomic structure of the surfaces in contact, in order to calculate the resistance to rolling in such a contact. It is shown that the contact reacts elastically to torque forces up to a critical bending angle. Beyond that, irreversible rolling occurs. In the elastic regime, the moment opposing the attempt to roll is proportional to the bending angle and to the pull-off force P(sub c). Young's modulus of the involved materials has hardly any influence on the results. We show that agglomerates of sub-micron sized particles will in general be quite rigid and even long chains of particles cannot be bent easily. For very small particles, the contact will rather break than allow for rolling. We further discuss dynamic properties such as the possibility of vibrations in this degree of freedom and the typical amount of rolling during a collision of two particles.

Dominik, C.↗

Studies of the Intrinsic Complexities of Magnetotail Ion Distributions: Theory and Observations

This year we have studied the relationship between the structure seen in measured distribution functions and the detailed magnetospheric configuration. Results from our recent studies using time-dependent large-scale kinetic (LSK) calculations are used to infer the sources of the ions in the velocity distribution functions measured by a single spacecraft (Geotail). Our results strongly indicate that the different ion sources and acceleration mechanisms producing a measured distribution function can explain this structure. Moreover, individual structures within distribution functions were traced back to single sources. We also confirmed the fractal nature of ion distributions.

Ashour-Abdalla, Maha↗

Quantifying the Complexity of Flaring Active Regions

While solar physicists have a better understanding of the importance magnetic fields play in the solar heating mechanism, it is still not possible to predict whether or when an active region will flare. In recent decades, qualitative studies of the changes in active region morphology have shown that there is generally an increase in the complexity of the spatial configuration of a solar active region leading up to a flare event. In this study, we quantify the spatial structure of the region using the differential Box-Counting Method (DBC) of fractal analysis. We analyze data from NASA/Marshall Space Flight Centr's vector magnetograph from two flaring active regions: AR 6089 from June 10, 1990, which produced one M1.7 flare, and AR 6659 from June 8, 9 and 10, 1991, this data set including one C5.7 and two M(6.4 and 3.2) flare. (AR 6659 produced several other flares). Several magnetic parameters are studied, including the transverse and longitudinal magnetic field components (Bt and B1), the total field (Bmag), and the magnetic shear, which describes the non-potentiality of the field. Results are presented for the time series of magnetograms in relation to the timing of flare events.

Stark, B.↗

An optical systems analysis approach to image resampling

All types of image registration require some type of resampling, either during the registration or as a final step in the registration process. Thus the image(s) must be regridded into a spatially uniform, or angularly uniform, coordinate system with some pre-defined resolution. Frequently the ending resolution is not the resolution at which the data was observed with. The registration algorithm designer and end product user are presented with a multitude of possible resampling methods each of which modify the spatial frequency content of the data in some way. The purpose of this paper is threefold: (1) to show how an imaging system modifies the scene from an end to end optical systems analysis approach, (2) to develop a generalized resampling model, and (3) empirically apply the model to simulated radiometric scene data and tabulate the results. A Hanning windowed sinc interpolator method will be developed based upon the optical characterization of the system. It will be discussed in terms of the effects and limitations of sampling, aliasing, spectral leakage, and computational complexity. Simulated radiometric scene data will be used to demonstrate each of the algorithms. A high resolution scene will be "grown" using a fractal growth algorithm based on mid-point recursion techniques. The result scene data will be convolved with a point spread function representing the optical response. The resultant scene will be convolved with the detection systems response and subsampled to the desired resolution. The resultant data product will be subsequently resampled to the correct grid using the Hanning windowed sinc interpolator and the results and errors tabulated and discussed.

Lyon, Richard G.↗

Horizontal Radiative Fluxes in Clouds at Absorbing Wavelengths

We discuss the effect of horizontal fluxes on the accuracy of a conventional plane-parallel radiative transfer calculation for a single pixel, known as the Independent Pixel Approximation (IPA) at absorbing wavelengths. Vertically integrated horizontal fluxes can be represented as a sum of three components; each component is the IPA accuracy on a pixel-by-pixel basis for reflectance, transmittance and absorptance, respectively. We show that IPA accuracy for reflectance always improves with more absorption, while the IPA accuracy for transmittance is less sensitive to the changes in absorption: with respect to the non-absorbing case, it may first deteriorate for weak absorption and then improve again for strongly absorbing wavelengths. EPA accuracy for absorptance always deteriorates with more absorption. As a result, vertically integrated horizontal fluxes, as a sum of IPA accuracies for reflectance, transmittance and absorptance, increase with more absorption. Finally, the question of correlations between horizontal fluxes, IPA uncertainties and radiative smoothing is addressed using wavenumber spectra of radiation fields reflected from or transmitted through fractal clouds.

Marshak, A.↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars, and ongoing field programs such as DoE/ARM. These programs are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity.

Cahalan, Robert F.↗

Cloud Power Spectra-Dependence on Solar Zenith Angle and Wavelength, Implications for Cloud Optical Property Retrievals

Scale breaks (spatial scales at which power-law exponent changes occur) observed in Landsat radiances have proven to be useful indicators of radiative interactions, and have aided the development of improved techniques in the remote sensing of clouds. This work extends previous theoretical studies to absorbing wavelengths by using both Landsat Thematic Mapper (TM) observations and Monte Carlo (MC) simulations to infer the systematic dependencies of power spectral shape on cloud characteristics, illumination conditions, and wavelength. We show that MC simulations operating on a simple fractal model of horizontally inhomogeneous clouds produce power spectra that qualitatively resemble observed spectra. We also show that the decrease in the spectra power-law exponent seen at intermediate scales (referred to as "roughening") as the Sun becomes more oblique is more pronounced at absorbing wavelengths. An automated procedure designed to detect the small scale break location is unable to find systematic differences between TM Band 4 and Band 7, despite the fact that MC simulations point to systematic differences in horizontal fluxes. The effect of these qualitative characteristics of the spatial spectra on the retrieval of cloud optical properties is examined by comparing power spectra of nadir radiances with power spectra of optical properties retrieved using either traditional Independent Pixel Approximation approaches or modifications based on normalized radiance indices and the inverse Non-local Independent Pixel Approximation. Assuming that the actual cloud properties follow perfect scaling behavior at all scales, we show the improvement of the proposed retrieval modifications.

Oreopoulos, Lazaros↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars. Ongoing field programs such as Department of Energy's Atmospheric Radiation Measurement (DoE/ARM) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of adiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is beginning to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.qsfc.nasa.crov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large- scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate.

Cahalan, Robert F.↗

Cloud-Vegetation Interaction: Use of Normalized Difference Cloud Index for Estimation of Cloud Optical Thickness

A new technique to retrieve cloud optical depth for broken clouds above green vegetation using ground-based zenith radiance measurements is developed. By analogy with the Normalized Difference Vegetation Index NDVI), the Normalized Difference Cloud Index (NDCI) is defined as a ratio between the difference and the sum of two zenith radiances measured for two narrow spectral bands in the visible and near-IR regions. The very different spectral behavior of cloud liquid water drops and green vegetation is the key physics behind the NDCI. It provides extra tools to remove the radiative effects of the 3D cloud structure. Numerical calculations based on fractal clouds and real measurements of NDCI and cloud liquid water path confirm the improvements.

Marshak, A.↗

Administrative Report on Training Awards

During the tenure of this award, the recipient (David Harker) conducted areas of research which resulted in the award of a doctorate degree (August 1999) from the University of Wyoming. The primary science activity was investigation of silicate dust mineralogy in comets, particularly comet C/1995 O1 (Hale-Bopp). Determination of the dust mineralogy permits us to address an important astrophysical question of relevance to the origins and evolution of solar systems, "Do comets contain relic interstellar dust grains? Since, comets represent the frozen reservoirs of primitive proto-solar dust and ice, we can gain insight into the problem of understanding the formation of protoplanetesimals in the. early solar nebula. Mid-infrared spectrophotometry (7 - 14 micron, R approx. = 180 - 360) of Hale-Bopp was obtained with the NASA/Ames HIFOGS at four distinct epochs pre- and post- perihelion ion from 1996 October through 1997 June. These observations were conducted at the Wyoming Infrared Observatory and the NASA Infrared Telescope Facility, and were supported by funds from this training grant. The emission at mid-IR wavelengths in cometary comae arises from carbon grains, and small silicate grains which produce resonance features. Hale-Bopp had the strongest silicate feature observed from any comet to date. Theoretical calculations utilizing Mie Scattering Theory were employed to construct Synthetic cometary spectra to fit with the observed HIFOGS 10 microns spectral feature of Hale-Bopp. Our analysis suggests that the observed spectra can be modeled with the Hanner grain size distribution peaked at alpha((sub p) = 0.2 microns of fractal porous grains with porosity parameter D = 2.5. This model spectrum also fits photometry points in the 3 - 5 microns region. Comparison with the ISO SWS spectrum of Hale-Bopp obtained 1996 October reveals that the crystalline olivine grains must be at it temperature hotter than computed from Mie theory.

Woodward, Charles E.↗

"Electrostructural Phase Changes" In Charged Particulate Clouds: Planetary and Astrophysical Implications

There is empirical evidence that freely-suspended triboelectrostatically charged particulate clouds of dielectric materials undergo rapid conversion from (nominally) monodispersed "aerosols" to a system of well-defined grain aggregates after grain motion or fluid turbulence ceases within the cloud. In United States Microgravity Laboratory Space Shuttle experiments USML-1 and USML-2, it was found that ballistically-energized grain dispersions would rapidly convert into populations of filamentary aggregates after natural fluid (air) damping of grain motion. Unless continuously disrupted mechanically, it was impossible to maintain a non-aggregated state for the grain clouds of sand-size materials. Similarly, ground- based experiments with very fine dust-size material produced the same results: rapid, impulsive "collapse" of the dispersed grains into well-defined filamentary structures. In both ground-based and microgravity experiments, the chains or filaments were created by long-range dipole electrostatic forces and dipole-induced dielectric interactions, not by monopole interactions. Maintenance of the structures was assisted by short-range static boundary adhesion forces and van der Waals interactions. When the aggregate containers in the USML experiments were disturbed after aggregate formation, the quiescently disposed filaments would rearrange themselves into fractal bundles and tighter clusters as a result of enforced encounters with one another. The long-range dipole interactions that bring the grains together into aggregates are a product of randomly-distributed monopole charges on the grain surfaces. In computer simulations, it has been shown that when the force vectors of all the random charges (of both sign) on a grain are resolved mathematically by assuming Coulombic interaction between them, the net result is a dipole moment on individual grains, even though the grains are electrically neutral insofar as there is no predominance, on their surface, of one charge sign over another. The random charges of both sign derive from natural grain-to-grain interactions that produce triboelectrification via charge exchange every time grain surfaces make contact with one another. The conversion from a random distribution of grains (upon which there are randomly distributed charges) into an organization of electrostatically-ordered aggregates, can be regarded (within the framework of granular-material science) as an "electrical or Coulombic phase change" of the particulate cloud. It is not totally dissimilar from the more normal phase-change concept in which, for example, a gas with long free-path-molecules suddenly becomes a solid as a result of structural ordering of the molecules (notably, also the result of electronic forces, albeit at a different scale). In both the gas-to-solid case, and the aerosol-to-aggregate case, the same materials and charges are present before and after the phase change, but their arrangement now has a higher degree of order and a lower-energy configuration. An input of energy into the system is required to reverse the situation. The aggregates in the USML experiments were observed to undergo at least two phase changes as noted above. The point about phase changes, and by implication, the "electrostructural" reorganizations in particulate clouds, is the following: (a) they can occur very rapidly, almost spontaneously, above a critical cloud density, (b) in going from a higher energy state to a lower energy state, they convert to a denser system, (c) energy must be required to reverse the situation, implying that energy is released during the high-to-low energy phase change. In applying this information to natural particulate clouds, some inferences can be made (it is stressed that reference is still to dielectric materials attracted by dipole forces). There are several natural settings to which the USML observations apply, and to which the phase-change implications likewise apply. Dense clouds of triboelectrically-charged, kinetically-energized grains are to be found in volcanic eruptions (particularly on earth), aeolian dust storms (particularly on Mars), meteorite impact ejecta curtains (on all planets), in "immature" debris rings around planets (e.g., that from which our own Moon may have condensed), and in gravitationally collapsing protoplanetary dust/planetesimal debris disks where dielectric granules are being increasingly brought into collisional relationships with one another (increasing both charge exchange and physical proximity). It is noted that in many of these cases, the degree of electrical charging on the grains is likely to be much higher than that in the USML experiments where charging was not enhanced above the "normal", naturally encountered level for the particular materials (quartz, glass, and various silicate minerals). Application of the phase-change concept suggests that volcanic, aeolian, and impact debris clouds may, under certain circumstances, undergo rapid, impulsive, or even catastrophic collapse into a denser state that will lead to rapid precipitation or fall-out of suspended particulates. Although this idea has been suggested previously by the author , the phase change concept possibly permits some new insights into cloud-system behavior. For example, in a protoplanetary debris disk, the work of gravity may suddenly be enhanced by electrostatically driven collapse of the system when materials reach a critical intergranular spacing or grain number density. This might reduce the rate of planet formation by orders of magnitude (considering the ratio of g-forces to electrostatic ones for very small grains in close proximity), and indeed, it might drive the collapse into a system configuration that would perhaps not be created by gravity alone. Additional information is contained in the original.

Marshall, J. R.↗

Nonprotective Alumina Growth in Sulfur-Doped NiAl(Zr)

The 1200 C oxidation behavior of NiAl was examined at various levels of sulfur and zirconium dopants to test the possibility of a critical S/Zr ratio required for adhesion. Cyclic furnace testing for 200 1 -hr cycles and interrupted testing for 500 hr were used as screening tests. Pure NiAl and NiAl(Zr) with 0. 14 at.% Zr were chosen as model base compositions; they exhibited normal, slow-growing scales (3 Mg/sq cm) with excellent adhesion for the Zr-doped alloys. NiAl with about 120 ppma S exhibited a substantial weight loss (-20 Mg/sq cm) in cyclic tests and a very large weight gain (+60 Mg/sq cm) in interrupted tests. The major surface phase remained as alpha -Al2O3. Sulfur doping the NiAl(Zr) alloy caused massive weight gains of 80 - 100 Mg/sq cm, swelling, cracking, and nearly complete conversion into NiAl2O4, and alpha- Al2O3. The initial objective of determining critical S/Zr ratios for adhesion was therefore unattainable. Initiation of the catastrophic attack was examined after a 10 hr exposure, revealing a few sites of broad, raised, and cracked ridges. In cross-section, the ridges appeared as modular intrusions, with a complex, fractal, oxide-metal interface. They were primarily alumina (with occasional entrapped islands of NiAl2O4 or pure Ni metal). They possessed a unique microstructure consisting of 0.3 microns lamellae, separated by 0.1 microns open channels. This allowed for rapid growth controlled by gaseous diffusion. The microstructure is discussed in terms of SO2 evolution and a sulfur-driven de-passivation process.

Smialek, James L.↗

MOLA-derived Roughness Data Used to Predict Surface Scattering for Mars Subsurface Radar Sounding

The Mars Express orbiter, to be launched by the European Space Agency in 2003, will carry a low-frequency radar sounding instrument, MARSIS (Mars Advanced Radar for Subsurface and Ionospheric Sounding). The primary goal of MARSIS is to map the distribution of water, both solid and (if present) liquid, in the upper several km of the martian crust. Detecting discontinuities in the crust, such as an ice-water transition, presents many challenges for a Mars orbital radar sounder. One challenge that must be overcome is the presence of radar scattering (echoes) from the surface of Mars, expected to be detected by the sounder antennas at the same time as any echoes aris-ing from subsurface interfaces. As the transmitted spherical wavefront spreads within the crust of Mars,. it also interacts with surface topography at off-nadir positions, creating a "clutter" signal that can mask the subsurface echoes. The MARSIS instrument will utilize Doppler filtering to limit the off-nadir clutter in the along-track direction, and a nadir-null secondary antenna to identify strong off-nadir clutter from the cross-track direction. To evaluate the effects of off-nadir surface clutter and the capability of these schemes to reduce the clutter, it is necessary to predict the range of scattering behavior that may be expected from martian surface topography. In this paper, we utilize Mars Orbital Laser Altimeter (MOLA) data from the current Mars Global Surveyor mission to characterize the topographic roughness of a variety of martian terrain types, at scales relevant to the MARSIS clutter problem. Segments of MOLA altimetry profiles are reduced to the topographic parameters rms slope and fractal dimension, which then are used as inputs to a near-nadir radar scattering model to predict the strength of the clutter signal. Additional information is contained in the original extended abstract.

Plaut, J. J.↗

On the Implementation of a Land Cover Classification System for SAR Images Using Khoros

The Synthetic Aperture Radar (SAR) sensor is widely used to record data about the ground under all atmospheric conditions. The SAR acquired images have very good resolution which necessitates the development of a classification system that process the SAR images to extract useful information for different applications. In this work, a complete system for the land cover classification was designed and programmed using the Khoros, a data flow visual language environment, taking full advantages of the polymorphic data services that it provides. Image analysis was applied to SAR images to improve and automate the processes of recognition and classification of the different regions like mountains and lakes. Both unsupervised and supervised classification utilities were used. The unsupervised classification routines included the use of several Classification/Clustering algorithms like the K-means, ISO2, Weighted Minimum Distance, and the Localized Receptive Field (LRF) training/classifier. Different texture analysis approaches such as Invariant Moments, Fractal Dimension and Second Order statistics were implemented for supervised classification of the images. The results and conclusions for SAR image classification using the various unsupervised and supervised procedures are presented based on their accuracy and performance.

Medina Revera, Edwin J.↗

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗