Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scattering geometry”

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 577 records · Page 32

A new spectrally-invariant approach to the remote sensing of inhomogeneous clouds

Current operational satellite retrievals of cloud optical and microphysical properties go back to the Nakajima-King technique developed at the end of the 1980s. This technique is based on library calculations for plane-parallel homogeneous clouds. It often works well for overcast skies but leads to substantial errors for inhomogeneous and broken cloud fields where the plane-parallel-geometry assumption is no longer valid. The basic concept of a new technique was first introduced in nuclear physics to quantify the critical conditions of reactors. Based on the eigenvalues of the radiative transfer equation, their approach provides a powerful means to parameterize the structure of 3D media. This parameterization was later successfully applied to relate surface reflectance spectra to 3D canopy structure and is known as the spectrally-invariant approximation. The proposed approach adapts this technique to the remote sensing of cloud properties such as droplet single scattering albedo and the average number of scatterings (which are the fundamental parameters in radiative transfer theory), with an emphasis on quantifying the associated errors and uncertainties. This retrieval is free from the plane-parallel homogeneous cloud assumption.

clouds↗

Electronic structure of BaO/W cathode surfaces

The local electronic structure of the emissive layer of barium dispenser thermionic cathodes is investigated theoretically using the relativistic scattered-wave approach. The interaction of Ba and O with W, Os, and W-Os alloy surfaces is studied with atomic clusters modeling different absorption environments representative of B- and M-type cathodes. Ba is found to be strongly oxidized, while O and the metal substrate are in a reduced chemical state. The presence of O enhances the surface dipole and Ba binding energy relative to Ba on W. Model results for W-Os alloy substrates show only relatively small changes in Ba and O for identical geometries, but very large charge redistributions inside the substrate, which are attributed to the electronegativity difference between Os and W. If Os is present in the surface layer, the charge transfer from Ba to the substrate and the Ba binding energy increase relative to W. Explanations are offered for the improved electron emission from alloy surfaces and the different emission enhancement for different alloy substrates.

Muller, Wolfgang↗

Wave properties near the subsolar magnetopause - Pc 1 waves in the sheath transition layer

We study the waves in the frequency range of Pc 1 in the sheath transition layer of the magnetopause from the ISEE 1 and 2 observations. The waves are enhanced in the sheath transition layer, although they are scattered into the magnetosheath when the outer edge of the sheath transition layer is not sharp. The wave frequency is proportional to, and equal to, about 44 percent of the ion gyrofrequency. The waves are left-handed polarized for southward IMF, but linearly polarized for northward IMF. The direction of maximum variation is perpendicular to both the background field and the gradients of the field and density for northward IMF; for southward IMF, the waves are more turbulent. Wave generation mechanisms may depend on IMF orientations rather than the shock geometry. To investigate the free energy generating the waves for northward IMF, a method is developed combining the measurements from the fast plasma experiment and Lepedea to obtain a high time resolution estimate of the temperature anisotropy for strongly northward IMF. The estimated ion temperature anisotropy is enhanced, up to a factor of 2, within the sheath transition layer for northward IMF.

Song, P.↗

Paths and patterns: the biology and physics of swimming bacterial populations

The velocity distribution of swimming micro-organisms depends on directional cues supplied by the environment. Directional swimming within a bounded space results in the accumulation of organisms near one or more surfaces. Gravity, gradients of chemical concentration and illumination affect the motile behaviour of individual swimmers. Concentrated populations of organisms scatter and absorb light or consume molecules, such as oxygen. When supply is one-sided, consumption creates gradients; the presence of the population alters the intensity and the symmetry of the environmental cues. Patterns of cues interact dynamically with patterns of the consumer population. In suspensions, spatial variations in the concentration of organisms are equivalent to variations of mean mass density of the fluid. When organisms accumulate in one region whilst moving away from another region, the force of gravity causes convection that translocates both organisms and dissolved substances. The geometry of the resulting concentration-convection patterns has features that are remarkably reproducible. Of interest for biology are (1) the long-range organisation achieved by organisms that do not communicate, and (2) that the entire system, consisting of fluid, cells, directional supply of consumables, boundaries and gravity, generates a dynamic that improves the organisms' habitat by enhancing transport and mixing. Velocity distributions of the bacterium Bacillus subtilis have been measured within the milieu of the spatially and temporally varying oxygen concentration which they themselves create. These distributions of swimming speed and direction are the fundamental ingredients required for a quantitative mathematical treatment of the patterns. The quantitative measurement of swimming behaviour also contributes to our understanding of aerotaxis of individual cells.

Non-NASA Center↗

Reanalysis of Clementine Bistatic Radar Data from the Lunar South Pole

On 9 April 1994 the Clementine spacecraft high-gain antenna was aimed toward the Moon's surface and the resulting 13-cm wavelength radio echoes were received on Earth. Using these data, we have found that the lunar surface generally follows a Lambertian bistatic scattering function sigma(sub 0) = K(sub D)cos(theta(sub i) with K(sub D) approx. 0.003 for the opposite (expected) sense of circular polarization and K(sub D) approx. 0.001 for the same (unexpected) sense. But there are important deviations-of up to 50% in some parts of the echo spectrum-from this simple form. Based on an earlier analysis of these same data, Nozette et al. claimed detection of an enhancement in echoes with right circular polarization from regions near the South Pole in a near-backscatter geometry. Such behavior would be consistent with presence of perhaps large quantities of water ice near the Pole. We have been unable to reproduce that result. Although we find weak suggestions of enhanced echoes at the time of South Pole backscatter, similar features are present at earlier and later times, adjacent frequencies, and in left circular polarization. If enhanced backscatter is present, it is not unique to the South Pole; if not unique to the Pole, then ice is less likely as an explanation for the enhancement.

Simpson, Richard A.↗

Development of 3D electromagnetic modeling tools for airborne vehicles

The main goal of this project is to develop methodologies for scattering by airborne composite vehicles. Although our primary focus continues to be the development of a general purpose code for analyzing the entire structure as a single unit, a number of other tasks are also pursued in parallel with this effort. These tasks are important in testing the overall approach and in developing suitable models for materials coatings, junctions and, more generally, in assessing the effectiveness of the various parts comprising the final code. Here, we briefly discuss our progress on the five different tasks which were pursued during this period. Our progress on each of these tasks is described in the detailed reports (listed at the end of this report) and the memoranda included. The first task described below is, of course, the core of this project and deals with the development of the overall code. Undoubtedly, it is the outcome of the research which was funded by NASA-Ames and the Navy over the past three years. During this year we developed the first finite element code for scattering by structures of arbitrary shape and composition. The code employs a new absorbing boundary condition which allows termination of the finite element mesh only 0.3 lambda from the outer surface of the target. This leads to a remarkable reduction of the mesh size and is a unique feature of the code. Other unique features of this code include capabilities to model resistive sheets, impedance sheets and anisotropic materials. This last capability is the latest feature of the code and is still under development. The code has been extensively validated for a number of composite geometries and some examples are given. The validation of the code is still in progress for anisotropic and larger non-metallic geometries and cavities. The developed finite element code is based on a Galerkin's formulation and employs edge-based tetrahedral elements for discretizing the dielectric sections and the region between the target and the outer mesh termination boundary (ATB). This boundary is placed in conformity with the target's outer surface, thus resulting in additional reduction of the unknown count.

Volakis, John L.↗

Radiative Heat Transfer in Finite Cylindrical Enclosures with Nonhomogeneous Participating Media

Results of a numerical solution for radiative heat transfer in homogeneous and nonhomogeneous participating media are presented. The geometry of interest is a finite axisymmetric cylindrical enclosure. The integral formulation for radiative transport is solved by the YIX method. A three-dimensional solution scheme is applied to two-dimensional axisymmetric geometry to simplify kernel calculations and to avoid difficulties associated with treating boundary conditions. As part of the effort to improve modeling capabilities for turbulent jet diffusion flames, predicted distributions for flame temperature and soot volume fraction are used to calculate radiative heat transfer from soot particles in such flames. It is shown that the nonhomogeneity of radiative property has very significant effects. The peak value of the divergence of radiative heat flux could be underestimated by 2 factor of 7 if a mean homogeneous radiative property is used. Since recent studies have shown that scattering by soot agglomerates is significant in flames, the effect of magnitude of scattering is also investigated and found to be nonnegligible.

Hsu, Pei-Feng↗

Laser Brazing of High Temperature Braze Alloy

The Space Shuttle Main Engine (SSME) consists of 1080 conical tubes, which are furnace brazed themselves, manifolds, and surrounding structural jacket making almost four miles of braze joints. Subsequent furnace braze cycles are performed due to localized braze voids between the coolant tubes. SSME nozzle experiences extremely high heat flux (180 mW/sq m) during hot fire. Braze voids between coolant tubes may result in hot combustion gas escape causing jacket bulges. The nozzle can be disqualified for flight or result in mission failure if the braze voids exceed the limits. Localized braze processes were considered to eliminate braze voids, however, damage to the parent materials often prohibited use of such process. Being the only manned flight reusable rocket engine, it has stringent requirement on the braze process. Poor braze quality or damage to the parent materials limits the nozzle service life. The objective of this study was to develop a laser brazing process to provide quality, localized braze joints without adverse affect on the parent materials. Gold (Au-Cu-Ni-Pd-Mn) based high temperature braze alloys were used in both powder and wire form. Thin section iron base superalloy A286 tube was used as substrate materials. Different Laser Systems including CO2 (10.6 micrometers, 1kW), ND:YAG (1.06 micrometers, 4kW). and direct diode laser (808nm. 150W) were investigated for brazing process. The laser process variables including wavelength. laser power, travel speed and angle of inclination were optimized according to bead geometry and braze alloy wetting at minimum heat input level, The properties of laser brazing were compared to that of furnace brazing. Microhardness profiles were used for braze joint property comparison between laser and furnace brazing. The cooling rate of laser brazing was compared to furnace brazing based on secondary dendritic arm spacing, Both optical and Scanning Electron Microscope (SEM) were used to evaluate the microstructures of the braze materials and tube substrate. Metallography of the laser braze joint was compared to the furnace braze. SEM Energy Disperse X-Ray Spectra (EDX) and back scattered imaging were used to analyze braze alloy segregation. Although all of the laser systems, CO2, ND:YAG, and direct diode laser produced good braze joint, the direct diode laser was selected for its system simplicity, compactness and portability. Excellent laser and braze alloy coupling is observed with powder alloy compared to braze alloy wire. Good wetting is found with different gold based braze alloys. The laser brazing process can be optimized so that the adverse affect on the parent materials can be eliminated. Metallography of the laser braze joint has shown that quality braze joint was produced with laser brazing process. Penetration of the laser braze to the substrate is at neglectable level. Zero penetration is observed. Microstructure examinations shown that no observable changes of the microstructure (grain structure and precipitation) in the HAZ area between laser braze and furnace braze. Wide gaps can be laser brazed with single pass for up to 0.024 inches. Finer dendritic structure is observed in laser brazing compared with equiaxial and coarser grain of the furnace brazing microstructure. Greater segregation is also found in the furnace braze. Higher hardness of the laser braze joint comparing to furnace braze is observed due to the fast cooling rate and Finer microstructure in the laser brazing. Laser braze joint properties meet or exceed the furnace joint properties. Direct diode laser for thin section tube brazing with high temperature braze alloys have been successfully demonstrated. The laser's high energy density and precise control has shown significant advantages in reducing process heat input to the substrates and provide high quality braze joints comparing to other localized braze process such as torch, TIG, and MPTA processes. Significant cost savings can be realized particularly with localized braze comparing to a full furnace braze cycle.

Gao, Y. P.↗

Canadian Led X-ray Polarimeter Mission CXP

We propose a Canadian-led X-ray Polarimetry Mission (CXP), to include a scattering X-ray Polarimeter and sensitive All-Sky X-ray Monitor (ASXM). Polarimetry would provide a new observational window on black holes, neutron stars, accretion disks and jets, and the ASXM would offer sensitive monitoring of the volatile X-ray sky. The envisioned polarimeter consists of a hollow scattering beryllium cone surrounded by an annular proportional counter, in a simple and elegant design that is reliable and low-risk. It would be sensitive in the 6-30 keV band to approx. 3% polarization in approx. 30 Galactic sources and 2 AGN in a baseline 1-yr mission, and have sensitivity greater than 10 times that of the previous X-ray polarimeter flown (NASA's OSO-8, 1975-78) for most sources. This X-ray polarimeter would tackle questions like, Do black holes spin?, How do pulsars pulse?, What is the geometry of the magnetic field in accreting neutron stars? Where and how are jets produced in microquasars and AGN?, What are the geometries of many of the most famous accretion-disk systems in the sky? This will be done using a novel and until-now unexploited technique that will greatly broaden the available observational phase space of compact objects by adding to timing and spectroscopy observations of polarization fraction and position angle as a function of energy. The All-Sky X-ray Monitor would scan for transients, both as potential targets for the polarimeter but also as a service to the worldwide astronomical community. The entire CXP mission could be flown for $40- 60M CDN, according to estimates by ComDev International, and could be built entirely in Canada. It would fall well within the CSA's SmallSat envelope and would empower the growing and dynamic Canadian High-Energy Astrophysics community with world leadership in a potentially high impact niche area.

Kaspi, V.↗

Propeller noise prediction

Analytic propeller noise prediction involves a sequence of computations culminating in the application of acoustic equations. The prediction sequence currently used by NASA in its ANOPP (aircraft noise prediction) program is described. The elements of the sequence are called program modules. The first group of modules analyzes the propeller geometry, the aerodynamics, including both potential and boundary layer flow, the propeller performance, and the surface loading distribution. This group of modules is based entirely on aerodynamic strip theory. The next group of modules deals with the actual noise prediction, based on data from the first group. Deterministic predictions of periodic thickness and loading noise are made using Farassat's time-domain methods. Broadband noise is predicted by the semi-empirical Schlinker-Amiet method. Near-field predictions of fuselage surface pressures include the effects of boundary layer refraction and (for a cylinder) scattering. Far-field predictions include atmospheric and ground effects. Experimental data from subsonic and transonic propellers are compared and NASA's future direction is propeller noise technology development are indicated.

Zorumski, W. E.↗

Model images of radio halos around supernova remnants

I present model calculations of profiles and two-dimensional images of the radio synchrotron emission of young supernova remnants, concentrating on observable effects of relativistic eletrons diffusing upstream of the shock wave. If the preshock electron scattering mean free path is sufficiently long, observable synchrotron halos outside the bulk of the radio emission can potentially result; their absence can constrain the mean free path from above. If scattering is primarily due, as expected, to Alfven waves with amplitude detla(B), the halo is expected to extend a distance of order r(sub g)c(delta(B)/B)(exp 2)/v(sub s) beyond the shock, where r(sub g) is the gyroradius of the electrons emitting at the observed frequency, B is the upstream magnetic field strength, v(sub s) is the shock velocity, and the amplitude delta(B) refers to wave with wavelength comparable to r(sub g), of order 10(exp 13) cm for typical supernova-remnant parameters. However, the detailed geometry of the halo varies with the assumptions about particle acceleration in the shock wave. I present an atlas of model profiles and images as a function of preshock diffusion length, of aspect angle between the magnetic field and the line of sight, and of other relevant parameters.

Reynolds, Stephen P.↗

A Probabilistic Design Method Applied to Smart Composite Structures

A probabilistic design method is described and demonstrated using a smart composite wing. Probabilistic structural design incorporates naturally occurring uncertainties including those in constituent (fiber/matrix) material properties, fabrication variables, structure geometry and control-related parameters. Probabilistic sensitivity factors are computed to identify those parameters that have a great influence on a specific structural reliability. Two performance criteria are used to demonstrate this design methodology. The first criterion requires that the actuated angle at the wing tip be bounded by upper and lower limits at a specified reliability. The second criterion requires that the probability of ply damage due to random impact load be smaller than an assigned value. When the relationship between reliability improvement and the sensitivity factors is assessed, the results show that a reduction in the scatter of the random variable with the largest sensitivity factor (absolute value) provides the lowest failure probability. An increase in the mean of the random variable with a negative sensitivity factor will reduce the failure probability. Therefore, the design can be improved by controlling or selecting distribution parameters associated with random variables. This can be implemented during the manufacturing process to obtain maximum benefit with minimum alterations.

Shiao, Michael C.↗

Highly Significant Detection of X-Ray Polarization from the Brightest Accreting Neutron Star Sco X-1

The Imaging X-ray Polarimetry Explorer measured with high significance the X-ray polarization of the brightest Z-source, Sco X-1, resulting in the nominal 2–8 keV energy band in a polarization degree of 1.0% ± 0.2% and a polarization angle of 8° ± 6° at a 90% confidence level. This observation was strictly simultaneous with observations performed by NICER, NuSTAR, and Insight-HXMT, which allowed for a precise characterization of its broadband spectrum from soft to hard X-rays. The source has been observed mainly in its soft state, with short periods of flaring. We also observed low-frequency quasiperiodic oscillations. From a spectropolarimetric analysis, we associate a polarization to the accretion disk at <3.2% at 90% confidence level, compatible with expectations for an electron scattering dominated optically thick atmosphere at the Sco X-1 inclination of ∼44°; for the higher-energy Comptonized component, we obtain a polarization of 1.3% ± 0.4%, in agreement with expectations for a slab of Thomson optical depth of ∼7 and an electron temperature of ∼3 keV. A polarization rotation with respect to previous observations by OSO-8 and PolarLight, and also with respect to the radio-jet position angle, is observed. This result may indicate a variation of the polarization with the source state that can be related to relativistic precession or a change in the corona geometry with the accretion flow.

Fabio La Monaca↗

A low-frequency synthesis array in earth orbit

The scientific objectives and design concept of a space-based VLBI array for high-resolution astronomical observations at 1-30 MHz are discussed. The types of investigations calling for such an array include radio spectroscopy of individual objects, measurement of the effects of scattering and refraction by the interplanetary medium (IPM) and the ISM, mapping the distribution of low-energy cosmic-ray electrons, and determining the extent of the Galactic halo. Consideration is given to the limitations imposed on an LF VLBI array by the ionosphere, the IPM, and the ISM; the calibration advantages offered by circular polar orbits of slightly differing ascending-node longitude for the array satellites; and collection of the IF data streams from the array satellites by one 'master' satellite prior to transmission to the ground. It is shown that determination of the three-dimensional array geometry by means of intersatellite radio links is feasible if there are at least seven spacecraft in the array.

Jones, Dayton L.↗

An Examination of the Spatial Distribution of Carbon Dioxide and Systematic Errors

The industrial period and modern age is characterized by combustion of coal, oil, and natural gas for primary energy and transportation leading to rising levels of atmospheric of CO2. This increase, which is being carefully measured, has ramifications throughout the biological world. Through remote sensing, it is possible to measure how many molecules of CO2 lie in a defined column of air. However, other gases and particles are present in the atmosphere, such as aerosols and water, which make such measurements more complicated1. Understanding the detailed geometry and path length of the observation is vital to computing the concentration of CO2. Comparing these satellite readings with ground-truth data (TCCON) the systematic errors arising from these sources can be assessed. Once the error is understood, it can be scaled for in the retrieval algorithms to create a set of data, which is closer to the TCCON measurements1. Using this process, the algorithms are being developed to reduce bias, within.1% worldwide of the true value. At this stage, the accuracy is within 1%, but through correcting small errors contained in the algorithms, such as accounting for the scattering of sunlight, the desired accuracy can be achieved.

Carbon Dioxide (CO2)↗

Freely Suspended Smectic Filaments and the Structure of the B7 Phase of MHOBOW

Our recent discovery of the spontaneous formation of chiral domains in fluid smectic phases of achiral bow-shaped molecules opens up a wide variety of possibilities for new liquid crystal phases and phenomena. The basic, spontaneously chiral layer structure of the highest temperature fluid smectic phases, the B2 and B7, are shown. One of the most intriguing aspects of this structure is the plethora of possible phases coming from different stacking sequences of the polar ordering and tilt directions. The four possibilities of next-nearest neighbor alternation are shown. In the original material studied, NOBOW, the ground states found are antiferroelectric, either the racemic SmC(sub S)P(sub A) or the chiral SmC(sub A)P(sub A). We are currently studying MHOBOW, synthesized by D. Walba which, by virtue of its methyl hexyloxy tail has a tendency to form anticlinic layer interfaces, in the hope of finding a phase with a ferroelectric ground state, either SmC(sub A)P(sub S) or SmC(sub S)P(sub A), which can be obtained in NOBOW only by applying a field. Preliminary observations of MHO-BOW have made its study, from the point of view of understanding novel LC structures, extremely high priority. The following truly remarkable characteristics have been revealed: (i) The smectic phase grows out of the isotropic in the form of helical ribbons. The resulting planar aligned textures of focal conics with layers normal to glass plates exhibit bizarre modulations, including stripes and checker-boards. These have also been seen in other materials suggesting that this is a new phase (tentatively called B7), which is a fluid smectic with some kind of in-layer structure. (ii) It is virtually impossible to make freely suspended films of MHOBOW. Rather it makes the freely suspended filaments which preliminary x-ray scattering experiments reveal to have the nested cylinder layer structure indicated; (iii) The powder x-ray diffraction exhibits four resolution-limited smectic layering peaks, very close in layer spacing, which vary continuously with T. This is further evidence for a more complex three dimensional structure than NOBOW, which has a typical single layering reflection. (iv) The x-ray structure factor of the layering peak of the filaments is extraordinarily complex and rich. Varying in qL (the scattering vector component along the filament axis) from a double slit-like pattern to modulated layer-like patterns, as qH (the scattering vector component normal to the filament axis) is varied over the range where the four powder peaks are located. These results suggest some kind of mosaic structure, perhaps with different layer spacings corresponding to the different stacking sequences. Recent x-ray diffraction experiments show that the peaks are modulated in intensity upon translation along a filament, in domains of several hundred microns dimension. These preliminary experiments suggest that the B7 is a fluid smectic with extremely unusual and fascinating structures. Of all of the many hundreds of fluid smectic materials we have attempted to study in the freely suspended film geometry over the years, only a few have failed to form films, and none showed any great tendency to form filaments, although this clearly should be a possible freely suspended smectic LC morphology. On several occasions in the past we have intentionally tried to make filaments from a variety of smectics without success. Thus the smectic filament formation property makes the B7 phase unique. It seems quite likely that the stability of filaments is related to the in-plane structure. The filaments exhibit other interesting structural and optical features. They are birefringent with a local optic axis which is oblique and which can vary continuously along filament and which can be manipulated with an electric field applied normal to the fiber, as if the field were causing a rotation of the optic axis about the fiber axis. Rapid displacement of the ends of the fiber toward one another causes a macroscopic helixing at low T and causes thick regions to transiently appear at high T, a 1D analog of island formation on a rapidly compressed film.

Clark, N.↗

Convective flow effects on protein crystal growth

A high-resolution microscopic interferometric setup for the monitoring of protein morphologies has been developed. Growth or dissolution of a crystal can be resolved with a long-term depth resolution of 200 A and a lateral resolution of 2 microns. This capability of simultaneously monitoring the interfacial displacement with high local depth resolution has yielded several novel results. We have found with lysozyme that (1) the normal growth rate is oscillatory, and (2) depending on the impurity content of the solution, the growth step density is either greater or lower at the periphery of a facet than in its center. The repartitioning of Na plus and Cl minus ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed, to interpret the large body of data in unified way. The results strongly suggest that (1) the ion to lysozyne ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter electron microscopy results appear to confirm this finding, which could have far-reaching consequences for x-ray diffraction studies. A computational model for diffusive-convective transport in protein crystallization has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies and our kinetics data for the growth of lysozyme. The results show that even in the small cell employed, protein concentration nonuniformities and gravity-driven solutal convection can be significant. The calculated convection velocities are of the same order to magnitude as those found in earlier experiments. As expected, convective transport, i.e., at Og, lysozyme crystal growth remains kinetically limited. The salt distribution in the crystal is predicted to be non-uniform at both 1g and 0g, as a consequence of protein depletion in the solution. Static and dynamic light scattering studies in undersaturated and supersaturated solutions have been performed. Diffusivities in undersaturated solutions, were found to vary with lysozyme concentrations. Depending on the salt concentration, the diffusivities either increase or decrease. Interestingly, the corresponding static scattering intensities behave oppositely, Our current analysis indicates that these changes are inconsistent with aggregation in undersaturated solutions. However, the data are compatible with concentration-dependent changes of the interactions between protein and salt.

Rosenberger, Franz↗

A 3D High-Order Unstructured Finite-Volume Algorithm for Solving Maxwell's Equations

A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Liu, Yen↗