Search NASA⌕ Search

SEARCH · Search NASA

Results for “Near Zero Scattering”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

The Evolution and Disruption of Planetary Systems

Planetary systems that encounter passing stars can experience severe orbital disruption, and the efficiency of this process is greatly enhanced when the impinging systems are binary pairs rather than single stars. Using a Monte Carlo approach, we have performed nearly half a million numerical experiments to examine the long term ramifications of planetary scattering on planetary systems. We have concluded that systems which form in dense environments such as Orion's Trapezium cluster have roughly a ten percent chance of being seriously disrupted. We have also used our programs to explore the long-term prospects for our own Solar system. Given the current interstellar environment, we have computed the odds that Earth will find its orbit seriously disrupted prior to the emergence of a runaway greenhouse effect driven by the Sun's increasing luminosity. This estimate includes both direct disruption events and scattering processes that seriously alter the orbits of the Jovian planets, which then force severe changes upon the Earth's orbit. We then explore the consequences of the Earth being thrown into deep space. The surface biosphere would rapidly shut down under conditions of zero insolation, but the Earth's radioactive heat is capable of maintaining life deep underground, and perhaps in hydrothermal vent communities, for some time to come. Although unlikely for the Earth, this scenario may be common throughout the universe, since many environments where liquid water could exist (e.g., Europa and Callisto) must derive their energy from internal (rather than external) heating.

Laughlin, Gregory↗
NASA NTRS1996

Electric Field Induced Interfacial Instabilities

The study of the interface in a charge-free, nonpolar, critical and near-critical binary fluid in the presence of an externally applied electric field is presented. At sufficiently large fields, the interface between the two phases of the binary fluid should become unstable and exhibit an undulation with a predefined wavelength on the order of the capillary length. As the critical point is approached, this wavelength is reduced, potentially approaching length-scales such as the correlation length or critical nucleation radius. At this point the critical properties of the system may be affected. In zero gravity, the interface is unstable at all long wavelengths in the presence of a field applied across it. It is conjectured that this will cause the binary fluid to break up into domains small enough to be outside the instability condition. The resulting pattern formation, and the effects on the critical properties as the domains approach the correlation length are of acute interest. With direct observation, laser light scattering, and interferometry, the phenomena can be probed to gain further understanding of interfacial instabilities and the pattern formation which results, and dimensional crossover in critical systems as the critical fluctuations in a particular direction are suppressed by external forces.

Kusner, Robert E.↗
DOE OSTI2026

ZTF-SEDm Type Ia supernova sample for Twins Embedding spectrophotometric standardization

Aims. This paper has two aims: the first aim is to build a large homogeneous spectrophotometric sample of Type Ia supernovae (SNe Ia) from the second data release of the Zwicky Transient Facility (ZTF DR2). We used the spectrum sample from the low-resolution ( R ∼ 100) SEDmachine (SEDm) Integral Field Spectrograph (IFS) that gathered 3069 spectra. This is one of the largest samples of such collections that can attempt to reproduce the Twins Embedding (TE) spectrophotometric standardization method. This is our second objective. The method was developed based on high-quality spectra from 200 SNe Ia of the Nearby Supernova factory (SNfactory) and led to an exceptionally low value of 0.073 mag for the intrinsic scatter. Methods. As the SEDm is not designed as a spectrophotometric instrument, we first improved the flux-calibration accuracy of the SN Ia spectrum sample using the ZTF photometric data, which were calibrated at the percent level. We corrected the spectra for second-order polynomials, fitted by comparing the synthetic photometry in the ZTF g , r , i filters with the light-curve (LC) data. We then applied the three steps of the TE parameterization to a subset of 783 ZTF SN spectra near maximum light while comparing results from SNfactory and ZTF. We finally analyzed the standardization methods based on the TE parameters. Results. The precision of the phase-correction model, which is the first step of the TE, is estimated at 0.01 mag in g band based on ZTF data. Despite the challenge posed by the spectrum-extraction pipeline associated with the SEDm (flux calibration, leftover host signal, low signal-to-noise ratio, and low resolution), we applied a first standardization in color based on the second step of the TE, called read between the lines (RBTL), to the ZTF sample. We reached a Hubble residual scatter of 0.153 mag, all in normalized median absolute deviation, which is to be compared to the ∼0.11 mag obtained with the SNfactory data. The SALT color and stretch standardization reaches a scatter of 0.164 mag for the same ZTF SN Ia sample, and its host steps are ∼0.1 mag and zero for RBTL. When considering the scatter due to the redshift error and flux calibration error, we estimated a RBTL scatter of ∼0.129 mag for this ZTF sample as an upper limit because we identified an additional contribution from a systematic error in color. We tested the standardization based on the nonlinear TE parameters, and, as expected from the low spectrum quality, it did not improve the overall dispersion. Conclusions. We release 1897 flux calibrated spectra of 1607 SNe Ia with an estimated photometric accuracy of 0.07 mag. We further demonstrate that some amount of spectrophotometric SN Ia standardization can be done with limited-quality spectra. The RBTL standardization is more efficient than that of SALT with one parameter less, and the resulting host steps are consistent with zero. This makes it less prone to astrophysical bias. For future spectroscopic surveys, targeting the extraction pipeline for a thorough flux calibration and good signal-to-noise ratio would enable us to compute the full TE standardization, which would further reduce the scatter in the distance estimate.

Ganot, C↗
DOE OSTI2025

XY-like incommensurate magnetic order in Ce 2 ⁢SnS 5

We report the synthesis of single crystals of Ce 2 ⁢SnS 5 through a two-stage chemical vapor transport method. The Ce 2 ⁢SnS 5 system is a member of the orthorhombic Pbam (No. 55) space group and realizes a distorted trigonal tricapped prism (TTP) crystal field around each cerium site. We characterized the sample through orientation-dependent magnetization and heat capacity measurements to probe the magnetic anisotropy in the system characteristic of XY-like anisotropic Heisenberg model behavior. Ce 2 ⁢SnS 5 furthermore enters a zero-field ordered phase under 𝑇 𝑁 =2.4 K; powder neutron diffraction measurements reveal incommensurate magnetic order near 𝑇 𝑁 . Furthermore, the system then locks into a commensurate, two-𝑞 magnetic structure below approximately 1.2 K. This commensurate structure belongs to the Shubnikov group 𝑃⁢𝑏′⁢𝑎′⁢𝑚′ ⁡(MSG 55.359) and realizes the propagation vectors $\overrightarrow{𝑞}$ = (1/3, 0, 0) and $\overrightarrow{𝑞}$ = (0, 0, 0).

Antiferromagnetism↗
NASA NTRS2016

NuSTAR Discovery of a Cyclotron Line in the Accreting X-Ray Pulsar IGR J16393-4643

The high-mass X-ray binary and accreting X-ray pulsar IGR J16393-4643 was observed by the Nuclear Spectroscope Telescope Array in the 3-79 keV energy band for a net exposure time of 50 ks. We present the results of this observation which enabled the discovery of a cyclotron resonant scattering feature with a centroid energy of -29.3(sup +1.1)(sub -1.3) keV. This allowed us to measure the magnetic field strength of the neutron star for the first time: B = (2.5 +/- 0.1) x 10(exp 12) G. The known pulsation period is now observed at 904.0+/- 0.1 s. Since 2006, the neutron star has undergone a long-term spin-up trend at a rate of P= -2 x 10(exp -8) s s(exp -1) (-0.6 s per year, or a frequency derivative of v = 3 x 10(exp -14) Hz s(exp -1)). In the power density spectrum, a break appears at the pulse frequency which separates the zero slope at low frequency from the steeper slope at high frequency. This addition of angular momentum to the neutron star could be due to the accretion of a quasi-spherical wind, or it could be caused by the transient appearance of a prograde accretion disk that is nearly in corotation with the neutron star whose magnetospheric radius is around 2 x 10(exp 8) cm.

accretion disks – gamma-ray burst: general – s↗
DOE OSTI2025

Charge-density-wave quantum critical point under pressure in 2 H -TaSe 2

The presence of a quantum-critical point (QCP) at which a nearby ordered phase is suppressed to zero temperature is often invoked to explain emergent quantum phases, e.g. superconductivity. Yet, identifying a QCP and establishing its correlation with superconductivity remains challenging. Materials featuring charge-density-wave (CDW) order and superconductivity offer a clear scenario as both states can be associated with electron-phonon coupling. Here, we uncover a CDW-QCP and demonstrate its interrelation with superconductivity in the prototypical transition-metal dichalcogenide 2 H -TaSe 2 . We determine the evolution of the CDW state up to and beyond its suppression at the critical pressure p c = 19.9(1) GPa by means of X-ray diffraction and inelastic X-ray scattering measurements providing a full crystallographic refinement of the commensurate CDW superstructure. The pressure-induced CDW-QCP in close vicinity to the maximum superconducting transition temperature. Ab-initio lattice dynamical calculations corroborate that 2 H -TaSe 2 features order-parameter fluctuation enhanced superconductivity and can serve as a paradigm to investigate superconductivity near a CDW-QCP.

critical phenomena↗
DOE OSTI2022

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Smooth splitting and zeros from on-shell recursion

We describe a new approach to understanding the origins of recently discovered “hidden zeros” and “smooth splitting” of tree-level amplitudes in Tr ϕ 3 , Non-Linear Sigma Model (NLSM), Yang-Mill-Scalar (YMS) and the special Galileon. Introducing a new type of linear shift in kinematic space we demonstrate that the mysterious splitting formulae follow from a simple contour integration argument in the style of on-shell recursion. The argument makes use of only standard notions of tree-level factorization on propagators, but assumes improved UV behavior in the form of the absence of a residue at infinity. In the case of Tr ϕ 3 and NLSM this is proven by identifying our shift as a special case of a more general construction called a g-vector shift; in the case of YMS it remains an unproven conjecture. This recursive perspective leads to numerous new results: we derive generalizations of the splitting formulae on more relaxed near-zero kinematics, including interesting new kinematic limits in which the amplitude splits into a triple-product; we also demonstrate that the uncolored special Galileon model has improved UV scaling and hence also splits. We also investigate the possible realization of hidden zeros in four dimensions. The conditions under which the dimensionality constraints are compatible with zero kinematics is investigated in detail for Tr ϕ 3 and YMS; for the latter we find they can be realized only with certain restrictions on external helicity states. The realizable 4d zeros are proven by a similar recursive argument based on BCFW and is found to generalize to a new class of intrinsically 4d “helicity zeros” present in all sectors of YM and also gravity.

effective field theories↗
DOE OSTI2025

Hidden zeros of the cosmological wavefunction

Motivated by the recent discovery of hidden zeros in particle and string amplitudes, we characterize zeros of individual graph contributions to the cosmological wavefunction of a scalar field theory. We demonstrate that these contributions factorize near these zeros for all tree graphs and provide evidence that this extends to loop graphs as well. We explicitly construct polytopal realizations of the relevant graph associahedra and show that the cosmological zeros have natural geometric and physical interpretations. As a byproduct, we establish an equivalence between the wavefunction coefficients of chain graphs and flat-space Tr(ϕ 3 ) amplitudes, enabling us to leverage the cosmological zeros to uncover the recently discovered hidden zeros of colored amplitudes.

Cosmological models↗
NASA NTRS2011

Large-Amplitude Transmitter-Associated and Lightning-Associated Whistler Waves in the Earth's Inner Plasmasphere at L less than 2

We report observations of very large amplitude whistler mode waves in the Earth fs nightside inner radiation belt enabled by the STEREO Time Domain Sampler. Amplitudes range from 30.110 mV/m (zero ]peak), 2 to 3 orders of magnitude larger than previously observed in this region. Measurements from the peak electric field detector (TDSMax) indicate that these large ]amplitude waves are prevalent throughout the plasmasphere. A detailed examination of high time resolution electric field waveforms is undertaken on a subset of these whistlers at L < 2, associated with pump waves from lightning flashes and the naval transmitter NPM in Hawaii, that become unstable after propagation through the ionosphere and grow to large amplitudes. Many of the waveforms undergo periodic polarization reversals near the lower hybrid and NPM naval transmitter frequencies. The reversals may be related to finite plasma temperature and gradients in density induced by ion cyclotron heating of the plasma at 200 Hz, the modulation frequency of the continuous ]mode NPM naval transmitter signal. Test particle simulations using the amplitudes and durations of the waves observed herein suggest that they can interact strongly with high ]energy (>100 keV) electrons on a time scale of <1 s and thus may be an important previously unaccounted for source of energization or pitch ]angle scattering in the inner radiation belt.

Breneman, A.↗
DOE OSTI2026

Frustration-driven magnetic correlations in the spin-$\frac{5}{2}$ triangular lattice antiferromagnet RbFe⁢(HPO 3 ) 2

Here, a detailed study of the structural and magnetic properties of a spin-5/2 triangular lattice antiferromagnet RbFe⁢(HPO 3 ) 2 is presented using x-ray diffraction, magnetization, heat capacity, and 31 P nuclear magnetic resonance (NMR) experiments on a polycrystalline sample. The crystal structure features an equilateral triangular lattice of Fe 3+ ions. The thermodynamic measurements reveal the onset of a magnetic long-range order at 𝑇 N⁢1 ≃7.8K in zero field, followed by another low-temperature field-induced ordering at 𝑇 N⁢2 in higher fields. The transition at 𝑇 N⁢1 is further confirmed from the NMR spin lattice relaxation measurements. The value of the frustration ratio (𝑓≃7) implies moderate spin frustration in the compound. The 31 P NMR spectra exhibit two distinct spectral lines corresponding to two inequivalent phosphorus sites (P1 and P2), consistent with the crystal structure. The P1 site is strongly coupled with an isotropic hyperfine coupling of 𝐴$^{iso}_{hf}$ = 0.55⁢(2)⁢ T/𝜇 B while the P2 site is weakly coupled with 𝐴$^{iso}_{hf}$ = 0.25⁢(3)⁢ T/𝜇 B with the Fe 3+ ions. The magnetic susceptibility and NMR shift data are described well assuming a spin-5/2 isotropic triangular lattice antiferromagnetic model with an average exchange coupling of 𝐽/𝑘 B = 2.8⁢(2)⁢K. Below 𝑇 N⁢1 , the spectra evolve into a nearly rectangular powder pattern, indicating a commensurate antiferromagnetic type order. The 31 P spin-lattice relaxation rate well below 𝑇N⁢1 follows a 𝑇 3 temperature dependence, implying a two-magnon Raman scattering mechanism in the ordered state. Three well-defined phase regimes are clearly ascertained in the 𝐻−𝑇 phase diagram, reflecting a weak magnetic anisotropy in the compound.

Nagpal, V. [Indian Institute of Science Education ↗
NASA NTRS1995

Multiwavelength monitoring of the BL Lacertae object PKS 2155-304. 4: Multiwavelength analysis

Simultaneous X-ray, ultraviolet, optical, infrared, and radio monitoring data were used to test and constrain models of continuum emission from the BL Lacertae object PKS 2155-304. Intensively sampled ultraviolet and soft X-ray light curves showed a clear temporal correlation with the X-rays leading the ultraviolet by 2-3 hr. This lag was found to be significantly different from zero after an exhaustive comparison of four different techniques for measuring temporal correlations. Variations in the ultraviolet trough optical wave bands were also strongly correlated, with no measurable lag down to limiting timescales of approximately less than 1-2 hr. This strong correlation extends to the near-infrared, but the less intensive sampling precludes measurement of any lag beyomnd an upper limit of approximately less than 1 day. These lags and limits of the order of hours are much shorter than most rapid observed single-band variations. Because of the very sparse radio sampling, it was not possible to measure quantitatively the correlation and lag with shorter wavelengths, but the data do suggest that the radio may lag the optical/ultraviolet by approximately 1 week, with longer delays and weaker variations to longer radio wavelengths. The epoch-folding Q(exp 2) statistic was used to test for periodicity, and no evidence for strict or quasi-periodicity was found in any of the light curves. Because they lead the lower frequencies, the soft X-rays (approximately less than 1 keV) cannot arise from synchrotron self-Compton scattering. These results also rule out the accretion disk model, which predicts a measurable lag between ultraviolet/optical wavelength bands and a correlation between hardness and brightness, neither of which were seen. They are consistent with the entire radio through X-ray continuum arising from direct synchrotron emission from a relativistic jet. However, the tapered jet model, in which the X-ray emission is produced closer in, has problems explaining the magnitude of the ultraviolet/X-ray lag, because the X-ray-emitting electrons have very short lifetimes (t(sub 1/2) much less than 1 s). The result that the lag is much smaller than the variability timescale suggests instead that the radiation may be produced in a flattened region such as a shock front.

Edelson, R.↗
NASA NTRS1995

GHRS observations of cool, low-gravity star. 2: Flow and turbulent velocities in the outer atmosphere of gamma CRU CIS (M3.4 III)

The Goddard High Resoulution Spectrograph (GHRS) on the Hubble Space Telescope (HST) has been used to obtain medium (R = 20,000) and high (R = 85,000) resoultion UV spectra of chromosphere emission features for the M3.4 III star gamma Cru. Small Science Aperture (SSA) G270M and Echelle-B spectra of selected regions in the 2300-2850 A range were obtained to determine the kinematics of the chromosphere using lines of C2), Fe2, Co2, Si1/2), Ni2, Mn2, and Mg2. Profiles of C2) (UV 0.01) lines and fluorescently excited lines of low optical depth indicate average turbulent velocities (Doppler FWHM) of 30.2 +/- 1.3 and 28.8 +/- 1.3 km/s, respectively. The fluorescent emission lines (mean RV = 21.3 +/- 0.9 km/s) and the wings of the emission components of Fe2 lines (mean RV = 22.8 +/- 0.4 km/s) are approximately at rest relative to the radial velocity of the star (21 km/s), while the C2) lines show a modest inflow (mean RV = 23.1 +/- 0.9 km/s). The more opaque lines of Fe2 and Mg2 exhibit complex profiles resulting from line formation in an optically thick, extended expanding atmosphere. The emission wings of these lines are broadened by multiple scattering, and they are centered near the photospheric radial velocity. Closer to line center, these strong lines show a strong blueshifted self-absorption feature (already seen in IUE data), indicative of formation in an expanding chromosphere, and a previously unseen dip in the profiles on the red side of line center. The absorption components, when extracted using simple Gaussian fits, show strong correlations with the relative optical depths of the lines. The derived absorption flow velocities converge to the photospheric velocity as one examines spectra features formed deeper in the atmosphere. The blueward abosrption velocity increases in magnitude from about 7 to 14 km/s with increasing line optical depth - the strong absorptions directly map the acceleration of the outflowing stellar wind, while the interpretation of the weaker redshifted abosrptions is more ambiguous, indicating either an inflow of material or formation in an extended, spherically expanding outflow. The Mg2 and Fe2 profiles, taken together, imply that the wind speed decreases between the atmospheric layers where the Mg2 and Fe2 self-absorption components are formed. Interstellar absorptions are seen in the resonance lines of Mg2 (UV 1) and Fe2 (UV1) with zero-volt lower levels, at about -3 km/s, consistent with models of the interstellar medium in the direction of gamma Cru. Finally, we have detected the Mg2 'satellite lines' seen in solar spectra obtained above the limb. In gamma Cru these lines are probably fluorescently excited by H Ly beta.

Carpenter, Kenneth G.↗
DOE OSTI2024

Noncollinear polar magnet Fe 2 (SeO 3 ) 3 (H 2 O) 3 with inequivalent Fe 3+ sites

The emergence of novel magnetic states becomes more likely when the inversion symmetry of the crystal field, relative to the center between two spins, is broken. We propose that placing magnetic spins in inequivalent sites in a polar lattice can promote a realization of nontrivial magnetic states and associated magnetic properties. To test our hypothesis, we study Fe 2 (SeO 3 )(H 2 O) 3 as a model system that displays two distinct Fe(1) and Fe(2) magnetic sites in a polar structure (R3c space group). At low fields μ 0 H≤ 0.06 T, the material undergoes an antiferromagnetic ordering with T N1 = 77 K and a second transition at T N2 ≈ 4 K. At μ 0 H≥ 0.06 T and 74 K ≤T≤ 76 K, a positive entropy change of ∼0.12 mJ mol −1 K −1 can be associated with a metamagnetic transition to possibly nontrivial spin states. At zero field, Fe(1) is nearly fully ordered at T≈ 25 K, while Fe(2) features magnetic frustration down to T = 4 K. The magnetic ground state, a result corroborated by single-crystal neutron diffraction and 57 Fe Mössbauer spectroscopy, is a noncollinear antiparallel arrangement of ferrimagnetic Fe(1)–Fe(2) dimers along the c-axis. The results demonstrate that placing distinct magnetic sites in a polar crystal lattice can enable a new pathway to modifying spin, orbital, and lattice degrees of freedom for unconventional magnetism.

Crystal lattices↗
DOE OSTI2021

Randomized Algorithms for Scientific Computing (RASC)

Randomized algorithms have propelled advances in artificial intelligence (AI) and represent a foundational research area in advancing AI for Science. Future advancements in DOE Office of Science priority areas such as climate science, astrophysics, fusion, advanced materials, combustion, and quantum computing all require randomized algorithms for surmounting challenges of complexity, robustness, and scalability. Advances in data collection and numerical simulation have changed the dynamics of scientific research and motivate the need for randomized algorithms. For instance, advances in imaging technologies such as X-ray ptychography, electron microscopy, electron energy loss spectroscopy, or adaptive optics lattice light-sheet microscopy collect hyperspectral imaging and scattering data in terabytes, at breakneck speed enabled by state-of-the-art detectors. The data collection is exceptionally fast compared with its analysis. Likewise, advances in high-performance architectures have made exascale computing a reality and changed the economies of scientific computing in the process. Floating-point operations that create data are essentially free in comparison with data movement. Thus far, most approaches have focused on creating faster hardware. Ironically, this faster hardware has exacerbated the problem by making data still easier to create. Under such an onslaught, scientists often resort to heuristic deterministic sampling schemes (e.g., low-precision arithmetic, sampling every nth element) and sacrifice potentially valuable accuracy. Dramatically better results can be achieved via randomized algorithms, reducing the data size as much as or more than naive deterministic subsampling can achieve, while retaining the high accuracy of computing on the full data set. By randomized algorithms we mean those algorithms that employ some form of randomness in internal algorithmic decisions to accelerate time to solution, increase scalability, or improve reliability. Examples include matrix sketching for solving large-scale least-squares problems (see Figure 1) and stochastic gradient descent for training machine learning models. We are not recommending heuristic methods but rather randomized algorithms that have certificates of correctness and probabilistic guarantees of optimality and near-optimality. Such approaches can be useful beyond acceleration, for example, in understanding how to avoid measure zero worst-case scenarios that plague methods such as QR matrix factorization.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Higher Order, Hybrid BEM/FEM Methods Applied to Antenna Modeling

In this presentation, the authors address topics relevant to higher order modeling using hybrid BEM/FEM formulations. The first of these is the limitation on convergence rates imposed by geometric modeling errors in the analysis of scattering by a dielectric sphere. The second topic is the application of an Incomplete LU Threshold (ILUT) preconditioner to solve the linear system resulting from the BEM/FEM formulation. The final tOpic is the application of the higher order BEM/FEM formulation to antenna modeling problems. The authors have previously presented work on the benefits of higher order modeling. To achieve these benefits, special attention is required in the integration of singular and near-singular terms arising in the surface integral equation. Several methods for handling these terms have been presented. It is also well known that achieving ~he high rates of convergence afforded by higher order bases may als'o require the employment of higher order geometry models. A number of publications have described the use of quadratic elements to model curved surfaces. The authors have shown in an EFIE formulation, applied to scattering by a PEC .sphere, that quadratic order elements may be insufficient to prevent the domination of modeling errors. In fact, on a PEC sphere with radius r = 0.58 Lambda(sub 0), a quartic order geometry representation was required to obtain a convergence benefi.t from quadratic bases when compared to the convergence rate achieved with linear bases. Initial trials indicate that, for a dielectric sphere of the same radius, - requirements on the geometry model are not as severe as for the PEC sphere. The authors will present convergence results for higher order bases as a function of the geometry model order in the hybrid BEM/FEM formulation applied to dielectric spheres. It is well known that the system matrix resulting from the hybrid BEM/FEM formulation is ill -conditioned. For many real applications, a good preconditioner is required to obtain usable convergence from an iterative solver. The authors have examined the use of an Incomplete LU Threshold (ILUT) preconditioner . to solver linear systems stemming from higher order BEM/FEM formulations in 2D scattering problems. Although the resulting preconditioner provided aD excellent approximation to the system inverse, its size in terms of non-zero entries represented only a modest improvement when compared with the fill-in associated with a sparse direct solver. Furthermore, the fill-in of the preconditioner could not be substantially reduced without the occurrence of instabilities. In addition to the results for these 2D problems, the authors will present iterative solution data from the application of the ILUT preconditioner to 3D problems.

Fink, P. W.↗
NASA NTRS2003

Concentration and Velocity Gradients in Fluidized Beds

In this work we focus on the height dependence of particle concentration, average velocity components, fluctuations in these velocities and, with the flow turned off, the sedimentation velocity. The latter quantities are measured using Particle Imaging Velocimetry (PIV). The PIV technique uses a 1-megapixel camera to capture two time-displaced images of particles in the bed. The depth of field of the imaging system is approximately 0.5 cm. The camera images a region with characteristic length of 2.6 cm for the small particles and 4.7 cm. for the large particles. The local direction of particle flow is determined by calculating the correlation function for sub-regions of 32 x 32 pixels. The velocity vector map is created from this correlation function using the time between images (we use 15 to 30 ms). The software is sensitive variations of 1/64th of a pixel. We produce velocity maps at various heights, each consisting of 3844 velocities. We break this map into three vertical zones for increased height information. The concentration profile is measured using an expanded (1 cm diameter) linearly polarized HeNe Laser incident on the fluidized bed. A COHU camera (gamma=1, AGC off) with a lens and a polarizer images the transmitted linearly polarized light to minimize the effects of multiply scattered light. The intensity profile (640 X 480 pixels) is well described by a Gaussian fit and the height of the Gaussian is used to characterize the concentration. This value is compared to the heights found for known concentrations. The sedimentation velocity is estimated using by imaging a region near the bottom of the bed and using PIV to measure the velocity as a function of time. With a nearly uniform concentration profile, the time can be converted to height information. The stable fluidized beds are made from large pseudo-monodisperse particles (silica spheres with radii (250-300) microns and (425-500) microns) dispersed in a glycerin/water mix. The Peclet number is sufficiently large that Brownian motion of the particles can be ignored and the Reynolds number sufficiently small that particle inertia is negligible. A packed particle bed is used to randomize and disperse the flowing fluid introduced by a peristaltic pump. The bed itself is a rectangular glass cell 8 cm wide (x), 0.8 cm deep and a height of 30.5 cm (z). The depth of field of the camera is approximately 0.5 cm so depth information is averaged. Over flow fluid is returned to the reservoir making a closed loop system. In these experiments the particles form a sediment approximately 5.7 cm high with the pump off and expand to 22 cm with the pump on. For the smaller particles the pump velocity is .5 millimeters per second and 1.1 millimeters per second for the large particles. At this concentration the bed has a very well defined top where particle concentration rapidly drops to zero.

McClymer, James P.↗