Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sigma Models”

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 523 records · Page 29

DeepUQ: Assessing the Aleatoric Uncertainties from two Deep Learning Methods

Assessing the quality of aleatoric uncertainty estimates from uncertainty quantification (UQ) deep learning methods is important in scientific contexts, where uncertainty is physically meaningful and important to characterize and interpret exactly. We systematically compare aleatoric uncertainty measured by two UQ techniques, Deep Ensembles (DE) and Deep Evidential Regression (DER). Our method focuses on both zero-dimensional (0D) and two-dimensional (2D) data, to explore how the UQ methods function for different data dimensionalities. We investigate uncertainty injected on the input and output variables and include a method to propagate uncertainty in the case of input uncertainty so that we can compare the predicted aleatoric uncertainty to the known values. We experiment with three levels of noise. The aleatoric uncertainty predicted across all models and experiments scales with the injected noise level. However, the predicted uncertainty is miscalibrated to $\rm{std}(\sigma_{\rm al})$ with the true uncertainty for half of the DE experiments and almost all of the DER experiments. The predicted uncertainty is the least accurate for both UQ methods for the 2D input uncertainty experiment and the high-noise level. While these results do not apply to more complex data, they highlight that further research on post-facto calibration for these methods would be beneficial, particularly for high-noise and high-dimensional settings.

Nevin, Rebecca↗

Stellar Mixing: I. Formalism

In this paper we use the Reynolds stress models (RSM) to derive algebraic expressions for the following variables: a) heat fluxes; b) J.l fluxes; and c) momentum fluxes. These relations, which are fully 3D, include: 1) stable and unstable stratification, represented by the Brunt-Vaislila frequency, N(exp 2) =-g/H(sub p_(del - del(sub ad))(1 - RI(sub mu)); 2) double diffusion, salt-fingers, and semi-convection, represented by the density ratio R(sub mu) = del(sub mu)/(del - del(sub ad)); 3) shear (differential rotation), represented by the mean squared shear Sigma(exp 2) or by the Richardson number, Ri =N(exp 2)Sigma(exp -2); 4) radiative losses represented by a Peclet number, Pe; 5) a complete analytical solution of the ID version of the model. In general, the model requires the solution of two differential equations for the eddy kinetic energy K and its rate of dissipation, epsilon. In the local and stationary cases, when production equals dissipation, the model equations are all algebraic.

Canuto, V .M.↗

Magnetic shaping of planetary nebulae and other stellar wind bubbles

As in the case of the solar wind, the magnetic field in the wind from a magnetized, rotating star becomes increasingly toroidal with distance from the star. The strength of the magnetic field can be characterized by sigma, the ratio of toroidal magnetic energy density to kinetic energy density in the equatorial plane of the wind. A fast wind shocks against the external medium and creates a bubble whose volume is dominated by shocked gas. The toroidal magnetic field increases in the shocked bubble and can dominate the thermal pressure. Because of the low velocities in the bubble, hydrostatic equilibrium is a good approximation and allows the calculation of the thermal and magnetic pressure in the bubble, as in the model of Begelman and Li (1992) for the Crab Nebula. The structure, which is axisymmetric and extended in the polar direction, depends on two parameters: sigma nu(sub w)/w(sub 0), where nu(sub w) is the wind velocity and w(sub 0) is the shell velocity in the polar direction, and lambda = nu(sub a)/w(sub 0), where nu(sub a) is the velocity of the slow wind. For small values of lambda, there is a cusp in the shell in the equatorial plane, i.e., there is an equatorial ring. For larger values of lambda, the maximum of the surface density moves away from the equator i.e., a double ring structure. Our models should apply to planetary nebulae, if their central stars are sufficiently magnetized; the calculated shapes do resemble the observed shapes of planetaries. In all cases, our model predicts that X-ray emission from the bubble is concentrated toward the polar axis. Finally, we briefly discuss the asymmetry of the Crab Nebula and 3C 58.

Chevalier, Roger A.↗

The Determination of Titan Gravity Field from Doppler Tracking of the Cassini Spacecraft

In its tour of the Saturnian system, the spacecraft Cassini is carrying out measurements of the gravity field of Titan, whose knowledge is crucial for constraining the internal structure of the satellite. In the five flybys devoted to gravity science, the spacecraft is tracked in X (8.4 GHz) and Ka band (32.5 GHz) from the antennas of NASA's Deep Space Network. The use of a dual frequency downlink is used to mitigate the effects of interplanetary plasma, the largest noise source affecting Doppler measurements. Variations in the wet path delay are effectively compensated by means of advanced water vapor radiometers placed close to the ground antennas. The first three flybys occurred on February 27, 2006, December 28, 2006, and June 29, 2007. Two additional flybys are planned in July 2008 and May 2010. This paper presents the estimation of the mass and quadrupole field of Titan from the first two flybys, carried out by the Cassini Radio Science Team using a short arc orbit determination. The data from the two flybys are first independently fit using a dynamical model of the spacecraft and the bodies of the Saturnian system, and then combined in a multi-arc solution. Under the assumption that the higher degree harmonics are negligible, the estimated values of the gravity parameters from the combined, multi-arc solution are GM = 8978.1337 +/- 0.0025 km(exp 3) / s(exp 2), J (sub 2) = (2.7221 +/- 0.0185) 10 (exp -5) and C (sub 22) = (1.1159 +/- 0.0040) 10 (exp -5) The excellent agreement (within 1.7 sigma) of the results from the two flybys further increases the confidence in the solution and provides an a posteriori validation of the dynamical model.

Iess, L.↗

Spiral Arms in the Asymmetrically Illuminated Disk of MWC 758 and Constraints on Giant Planets

We present the first near-IR scattered light detection of the transitional disk associated with the Herbig Ae star MWC 758 using data obtained as part of the Strategic Exploration of Exoplanets and Disks with Subaru, and 1.1 micrometer Hubble Space Telescope/NICMOS data. While submillimeter studies suggested there is a dust-depleted cavity with r = 0".35, we find scattered light as close as 0".1 (20-28 AU) from the star, with no visible cavity at H, K', or Ks . We find two small-scaled spiral structures that asymmetrically shadow the outer disk. We model one of the spirals using spiral density wave theory, and derive a disk aspect ratio of h approximately 0.18, indicating a dynamically warm disk. If the spiral pattern is excited by a perturber, we estimate its mass to be 5(exp +3)(sub -4) M(sub J), in the range where planet filtration models predict accretion continuing onto the star. Using a combination of non-redundant aperture masking data at L' and angular differential imaging with Locally Optimized Combination of Images at K' and Ks , we exclude stellar or massive brown dwarf companions within 300 mas of the Herbig Ae star, and all but planetary mass companions exterior to 0".5. We reach 5 sigma contrasts limiting companions to planetary masses, 3-4 M(sub J) at 1".0 and 2 M(sub J) at 1".55, using the COND models. Collectively, these data strengthen the case for MWC 758 already being a young planetary system.

Grady, C. A.↗

Automating Traffic Microsimulation from SYNCHRO UTDF to SUMO

Modern transportation research relies on seamlessly integrating traffic signal data with robust network representation and simulation tools. This study presents utdf2gmns, an open-source Python tool that automates conversion of the Universal Traffic Data Format, including network representation, signalized intersections, and turning volumes into the General Modeling Network Specification (GMNS) Standard. The resulting GMNS-compliant network can be converted for microsimulation in SUMO. By automatically extracting intersection control parameters and aligning them with GMNS conventions, utdf2gmns minimizes manual preprocessing and data loss. utdf2gmns also integrates with the Sigma-X engine to extract and visualize key traffic control metrics, such as phasing diagrams, turning volumes, volume-tocapacity ratios, and control delays. This streamlined workflow enables efficient scenario testing, accurate model building, and consistent data management. Validated through case studies, utdf2gmns reliably models complex urban corridors, promoting reproducibility and standardization. Documentation is available on GitHub and PyPI, supporting easy integration and community engagement.

Luo, Roy [ORNL] (ORCID:0009000312909983)↗

Superhorizon isocurvature fluctuations relax tensions

Here, we present a new class of models that have potential to alleviate tensions present in the cosmological data today. We postulate about variation in the sound horizon scale on superhorizon scales, i.e., on scales that are larger than that of the present observable low-redshift Universe (≳ 1Gpc) while at the same time smaller than the largest scales probed by the cosmic microwave background (CMB) (≲ 10Gpc). In this scenario, CMB peaks are naturally smoothed as preferred by the Planck data, while at the same time the low-redshift baryon acoustic oscillation calibration is partially decoupled from the CMB. Taking superhorizon variations in baryon fraction as an example and using approximate modeling, we find improvement in the best fit Planck power spectrum model Δχ 2 ~ 6 for 1 extra degree of freedom with the relevant extension parameter 10 3 σ b = 2.10 ± 0.60, implying about 10% variations in baryon fraction across the Universe. At the same time, S 8 drops by about 1 sigma, easing tension with weak lensing surveys. We find that the combination of Planck 2018 data, eBOSS BAO data, and Riess et al. distance ladder Hubble parameter determination produce a satisfactory fit in our model if we allow for a phantom dark energy equation of state.

79 ASTRONOMY AND ASTROPHYSICS↗

Graph neural networks for CO 2 solubility predictions in Deep Eutectic Solvents

Deep Eutectic Solvents (DESs) are a promising class of solvents for CO 2 capture. DESs are complex mixtures that can be designed to optimize CO solubility and overall capture process efficiency. However, the vast design landscape of DES mixtures makes experimental investigation prohibitive; as such, there is a need for computational models that can quickly and efficiently navigate the design space and inform data collection efforts. In this work, we propose Graph Neural Network (GNN) models for predicting CO 2 solubility for DESs; the GNN leverages a mixture graph representation that captures the molecular structure of the DES components as well as their intermolecular interactions. Here, we compare the GNN framework against alternative architectures (neural networks, graph convolution networks, and random forests) and data representations (molecular fingerprints, sigma profiles, and graphs). We show that the proposed approach offers superior predictive performance; specifically, we show that solubility can be predicted reliably directly from molecular structure (without the need of using sigma profiles as proposed in previous studies). This result is important, as obtaining sigma profiles requires expensive density functional theory computations. We also explored the ability of GNNs to predict solubility for new DES mixtures and operating conditions. We found that the model extrapolates across temperature reliably. However, we also found deficiencies in the ability of the model to predict solubility for DES mixtures, pressures, and molar ratio not included in the training sets; we show that this is due to an inherent lack of chemical diversity in datasets available in the literature. The proposed computational capabilities can thus help navigate the design space of DES and inform data collection efforts. Our models, data, and benchmarks are shared as Python code implemented in Jupyter notebooks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of tree structure on X-band microwave signature of conifers

Experimental studies are performed on some coniferous trees (Austrian pine, Nordmann spruce, and Norway spruce) to investigate the relation between the tree architecture and radar signal at X-band. For a single tree, the RCS is measured as a function of the scatterer location at 90 deg incidence. It is found that the main scatterers are the leafy branches and the difference between sigma(vv) and sigma(hh) is significant at the upper portion of the tree. At the lower portion of the tree, sigma(vv) and sigma(hh) have almost the same level. For a group of trees the angular trends of sigma(vv) and sigma(hh) are measured. It is found that the levels of sigma(vv) and sigma(hh) are of the same order, but their angular trends vary from one tree species to the other depending on the tree species structure. The interpretation of these experimental results is carried out with the help of a theoretical model which accounts for the structure of the tree. According to this theoretical study, the major scattering trend is due to the leaves, while the perturbation to the angular trend and the level difference between sigma(vv) and sigma(hh) are due to the branch orientation distributions (i.e., the tree architecture).

Mougin, Eric↗

Astrometry With the Hubble Space Telescope: Trigonometric Parallaxes of Planetary Nebula Nuclei NGC 6853, NGC 7293, ABELL 31, and DeHt 5

We present absolute parallaxes and relative proper motions for the central stars of the planetary nebulae NGC 6853 (The Dumbbell), NGC 7293 (The Helix), Abell 31, and DeHt 5. This paper details our reduction and analysis using DeHt 5 as an example. We obtain these planetary nebula nuclei (PNNi) parallaxes with astrometric data from Fine Guidance Sensors FGS 1r and FGS 3, white-light interferometers on the Hubble Space Telescope. Proper motions, spectral classifications and VJHKT2M and DDO51 photometry of the stars comprising the astrometric reference frames provide spectrophotometric estimates of reference star absolute parallaxes. Introducing these into our model as observations with error, we determine absolute parallaxes for each PNN. Weighted averaging with previous independent parallax measurements yields an average parallax precision, sigma (sub pi)/ pi = 5%. Derived distances are: d(sub NGC6853) = 405(exp +28 sub -25) pc, d(sub NGC7293) = 216(exp +14 sub -12) pc, d(sub Abell31) = 621(exp +91 sub -70) pc, and d(sub DeHt5) = 345(exp +19 sub -17) pc. These PNNi distances are all smaller than previously derived from spectroscopic analyses of the central stars. To obtain absolute magnitudes from these distances requires estimates of interstellar extinction. We average extinction measurements culled from the literature, from reddening based on PNNi intrinsic colors derived from model SEDs, and an assumption that each PNN experiences the same rate of extinction as a function of distance as do the reference stars nearest (in angular separation) to each central star. We also apply Lutz-Kelker bias corrections. The absolute magnitudes and effective temperatures permit estimates of PNNi radii through both the Stefan-Boltzmann relation and Eddington fluxes. Comparing absolute magnitudes with post-AGB models provides mass estimates. Masses cluster around 0.57 solar Mass, close to the peak of the white dwarf mass distribution. Adding a few more PNNi with well-determined distances and masses, we compare all the PNNi with cooler white dwarfs of similar mass, and confirm, as expected, that PNNi have larger radii than white dwarfs that have reached their final cooling tracks.

STARS↗

Muon $$\mathbf {g-2}$$, neutralino dark matter and stau NLSP

Abstract We explore the implications of resolving the muon $$g-2$$ g - 2 anomaly in a $$SU(4)_c \times SU(2)_L \times SU(2)_R$$ S U ( 4 ) c × S U ( 2 ) L × S U ( 2 ) R model, where the soft supersymmetry breaking scalar and gaugino masses break the left-right (LR) symmetry. A 2 $$\sigma $$ σ resolution of the anomaly requires relatively light sleptons, chargino and LSP neutralino. The stau turns out to be the NLSP of mass $$m_{\tilde{\tau }}\lesssim 400$$ m τ ~ ≲ 400 GeV, and the sleptons from the first two families can be as heavy as about 800 GeV. The chargino is also required to be lighter than about 600 GeV to accommodate the muon $$g-2$$ g - 2 solutions consistent with the dark matter relic density constraint. The dominant right-handed nature of the light slepton states suppress the sensitivity of possible signals which can be probed in Run3 experiments at the LHC. We also discuss the impact of accommodating the Higgs boson mass and the vacuum stability of the scalar potential for these solutions. Although a light stau can be compatible with the stability of the scalar potential, the Higgs boson mass constraint has a strong impact on the solutions with $$\tan \beta $$ tan β bounded from above, namely $$\tan \beta \lesssim 20$$ tan β ≲ 20 . The Higgsinos are heavier than about 4 TeV, and the LSP neutralino has the correct relic density if it is Bino-like. We identify stau–neutralino coannihilation as the dominant mechanism for realizing the desired dark matter relic density, with sneutrino–neutralino coannihiliation playing a minor role. These bino-like dark matter solutions can yield a spin-independent scattering cross-section on the order of $$10^{-13}$$ 10 - 13 pb which hopefully, can be expected to be tested in the near future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sensitivity to Dijet Resonances at Proton-Proton Colliders

A significant benchmark for discovery at a proton-proton collider is the sensitivity to a dijet resonance, X, the intermediate state of the s-channel process $pp \rightarrow X \rightarrow 2\mbox{ jets}$. To probe the highest resonance masses, hadron collider experiments have used the classic technique of searching for bumps in the mass spectrum of two individually resolved jets. In this Snowmass 2021 study, we explore the search sensitivity to multiple benchmark models of dijet resonances at current and future proton-proton colliders. We present the expected masses for $5\sigma$ discovery and 95% confidence level exclusion of diquarks, colorons, excited quarks, $W^{\prime}s$, $Z^{\prime}s$ and Randall-Sundrum gravitons, resulting from accumulation of integrated luminosities between 10 and $10^5$ fb$^{-1}$, at proton-proton colliders operating at energies $\sqrt{s}=$ 13, 14, 27, 75, 100, 150, 300 and 500 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deep HST imaging of distant weak radio and field galaxies

We present deep Hubble Space Telescope (HST) Wide-Field Camera (WFC) V- and I-band images of three distant weak radio galaxies with z = 0.311-2.390 and seven field galaxies with z = 0.131-0.58. The images were deconvolved with both the Lucy and multiresolution CLEAN methods, which yield a restoring Full Width at Half Maximum (FWHM) of less than or equal to 0.2 sec, (nearly) preserve photons and signal-to-noise ratio at low spatial frequencies, and produce consistent light profiles down to our 2 sigma surface brightness sensitivity limit of V approximately 27.2 and I approximately 25.9 mag/sq arcsec. Multi-component image modeling was used to provide deconvolution-independent estimates of structural parameters for symmetric galaxies. We present 12-band (m(sub 2750) UBVRIgriJHK) photometry for a subset of the galaxies and bootstrap the unknown FOC/48 zero point at 2750 A in three independent ways (yielding m(sub 2750) = 21.34 +/- 0.09 mag for 1.0 e(-)/s). Two radio galaxies with z = 0.311 and 0.528, as well as one field galaxy with z = 0.58, have the colors and spectra of early-type galaxies, and a(exp 1/4)-like light profiles in the HST images. The two at z greater than 0.5 have little or no color gradients in V - I and are likely giant ellipticals, while the z = 0.311 radio galaxy has a dim exponential disk and is likely an S0. Six of the seven field galaxies have light profiles that indicate (small) inner bulges following a(exp 1/4) laws and outer exponential disks, both with little or no color gradients. These are (early-type) spiral galaxies with z = 0.131-0.528. About half have faint companions or bars. One shows lumpy structure, possibly a merger. The compact narrow-line galaxy 53W002 at z = 2.390 has less than or = 30% +/- 10% of its HST V and I flux in the central kiloparsec (due to its weak Active Galactic Nucleus (AGN)). Most of its light (V approximately equal to 23.3) occurs in a symmetric envelope with a regular a(exp 1/4)-like profile of effective radius a approximately equal to 1.1 sec (approximately equal to 12 kpc for H(sub 0) = 50, q(sub 0) = 0.1. Its (HST) V - I color varies at most from approximately 0.3 mag at a approximately equal to 0.2 sec to approximately 1.2 mag at a approximately greater than 0.4 sec, and possibly to approximately greater than 2.2 mag at a approximately greater than 1.2 sec. Together with its I - K color (approximately equal to 2.5 mag for a approximately greater than 1.0 sec-2.0 sec), this is consistent with an aging stellar population approximately 0.3-0.5 Gyr old in the galaxy center (a approx. less than 2 kpc radius), and possibly approximately 0.5-1.0 Gyr old at a approximately greater than 10 kpc radius. While its outer part may thus have started to collapse at z = 2.5-4, its inner part still is aligned with its redshifted Ly(alpha) cloud and its radio axis, possibly caused by star formation associated with the radio jet, or by reflection from its AGN cone.

Windhorst, R. A.↗

The Compatibility of Friedmann Cosmological Models with Observed Properties of Gamma-Ray Bursts and a Large Hubble Constant

The distance scale to cosmic gamma-ray bursts (GRB's) is still uncertain by many orders of magnitude; however, one viable scenario places GRB's at cosmological distances, thereby permitting them to be used as tracers of the cosmological expansion over a significant range of redshifts zeta. Also, several recent measurements of the Hubble constant H(sub 0) appearing in the referred literature report values of 70-80 km/s /Mpc. Although there is significant debate regarding these measurements, we proceed here under the assumption that they are evidence of a large value for H(sub 0). This is done in order to investigate the additional constraints on cosmological models that can be obtained under this hypothesis when combined with the age of the universe and the brightness distribution of cosmological gamma-ray bursts. We show that the range of cosmological models that can be consistent with the GRB brightness distribution, a Hubble constant of 70-80 km/s/Mpc, and a minimum age of the universe of 13-15 Gyr is constrained significantly, largely independent of a wide range of assumptions regarding the evolutionary nature of the burst population. Low-density, Lambda greater than 0 cosmological models with deceleration parameter in the range -1 less than q(sub 0) less than 0 and density parameter sigma(sub 0) in the range approximately equals 0.10-0.25(Omega(sub 0) approximately equals 0.2-0.5) are strongly favored.

Horack, John M.↗

Kinetics and Thermochemistry of the Br((sup 2)P3/2) + NO2 Association Reaction

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the Br((sup 2)P3/2) + NO2 association reaction as a function of temperature (259-432 K) pressure (12.5 - 700 Torr), and buffer gas identity (He, Ar, H2, N2, CO2, CF4, SF6). The reaction is found to be in the falloff regime between third and second order over the entire range of conditions investigated. At temperatures below 350 K, the association reaction is found to be irreversible on the time scale of the experiment (approximately 30 ms). At higher temperatures reversible addition is observed, allowing equilibrium constants for BrNO2 formation and dissociation to be determined. Second- and third-law analyses of the equilibrium data are in only fair agreement and lead to the following thermochemical parameters for the association reaction: Delta-H(298) = 19.6 +/- 1.7 kcal/mol, Delta-H(0) = -18.6 +/- 2.0 kcal/mol, Delta-S(298) = 29.3 +/- 4.2 cal/mol/K, Delta-H(sub f)(sub 298)(BrNO2) = 17.0 +/-1.8 kcal/mol(uncertainties are 2 sigma estimates of absolute accuracy). The value for Delta-H(0) determined in this study has been employed to calculate k(sub 0)(sup SC), the low-pressure third-order rate coefficient in the strong collision limit, by using the method of Troe; calculated values of k(sub 0)(sup SC) are inconsistent with experimental results unless Delta-H(0) is assigned a value near the lower limit derived from analysis of the high-temperature approach to equilibrium data, i.e. delta-H(0) approximately equals -16.6 kcal/mol. A potential source of systematic error in the calculation of both k(sub 0)(sup SC) and the absolute entropy of BrNO2 results from the complete lack of knowledge of the energies and degeneracies of the electronic states of BrNO3. The procedure developed by Troe and co-workers has been employed to extrapolate experimental falloff curves to the low- and high-pressure limits. Derived values for k(sub 0)(M,298K) in units of 10(exp -31) cm(exp 6)/sq molecule/s range from 2.75 for M = He to 6.54 for M = CO2; 2 sigma uncertainties are estimated to be +/- 20%. Values for k(sub 0)(N2,T) in units of 10(exp -31) cm(exp 6)/sq molecule/s are 5.73 at 259 K, 4.61 at 298 K, and 3.21 at 346 K; the observed temperature dependence for k(sub 0)(N2,T) is consistent with the theoretical temperature dependence for Beta(sub c)k(sub 0)(sup SC). Values for k(sub infinity)(T) in units of 10(exp -11) cu cm/molecule/s are 2.86 at 259 K, 3.22 at 298 K, and 3.73 at 346 K; 2 sigma uncertainties are estimated to be a factor of 2. Approximate falloff parameters in a convenient format for atmospheric modeling are also derived.

Kreutter, K. D.↗

LIMEX '87 ice surface characteristics - Implications for C-band SAR backscatter signatures

Ice surface characterization data collected in 1987, during the Labrador Sea Ice Margin Experiment, are analyzed to estimate the changes in snow and ice properties at the onset of melt. Surface measurements were made from an ice research vessel on several days (some of which had coincident remote-sensing flights) at a number of locations in the marginal ice zone. These data are used as input parameters in a simple scattering model to simulate the effects of variations in material properties upon C-band scattering signatures. Snow moisture and large-scale surface roughness are demonstrated to have the largest effect upon HH polarization scattering cross-section sigma(HH)(O) and large differences are predicted between undeformed floe surfaces and deformed or undulating rough ice surfaces. Using a parametric approach, an approximate picture of rough and smooth ice signatures and their relative contrast are provided. The models reproduce a trend observed in synthetic aperture radar (SAR) images of increasing backscatter contrast between deformed and underformed ice over an early period of warning: such observations are consistent with the results of analysis of SAR images to date. The model also reproduces a calibrated SAR-derived signature with a reasonable degree of accuracy.

Drinkwater, Mark R.↗

Effects of heavy-element settling on solar neutrino fluxes and interior structure

We consider the effects of gravitational settling of both He and heavier elements on the predicted solar neutrino fluxes and interior sound speed and density profiles. We find that while the structural changes that result from the inclusion of both He and heavy-element settling are only slightly larger than the changes resulting from the inclusion of He settling alone, the additional increases in expected neutrino fluxes are of comparable size. Our preferred model with both He and heavy-element settling has neutrino count rates of 9.0 SNU for Cl-37 detectors and 137 SNU for Ga-71 detectors, as compared to 7.1 and 127 SNU for a comparable model without any diffusive separation, or 8.0 and 132 SNU for a model that includes He settling alone. We suggest that the correction factors by which the predicted neutrino fluxes of solar models calculated without including the effects of diffusion should be multiplied are 1.25 +/- 0.08 for Cl detectors, 1.07 +/- 0.02 for Ga detectors, and 1.28 +/- 0.09 for the B-8 flux (1 sigma errors). Comparison of internal sound speed and density profiles strongly suggests that the additional changes in calculated p-mode oscillation frequencies due to the inclusion of heavy-element settling will be small compared to the changes that result from He settling alone, especially for the higher degree modes. All models with diffusive separation give much better agreement with the observed depth of the convection zone than do nondiffusive models. The model that includes both He and heavy-element settling requires an initial He mass fraction Y = 0.280 and has a surface He abundance of Y = 0.251 at the solar age.

Proffitt, Charles R.↗

Gas dynamic simulations of galaxy formation

Results are presented from a simulation modeling the formation of a group of galaxies in a 'standard' cold, dark matter universe with delta = 1, h sub 0 = 50 km/(s(Mpc)), baryon fraction omega sub b = 0.1 and spectrum normalization sigma sub 8 = 0.6 (bias parameter b = 1.7). Initial conditions are generated within a periodic box with comoving length 16 Mpc in a manner constrained to produce a small cluster of total mass approximately 10 exp 14 solar mass. Two sets of 643 particles are used to model the dark matter and baryon fluids. Each gas particle represents 1.08 x 10 exp -8 solar mass, implying an L* galaxy is resolved by approximately 1000 particles. The system is evolved self-consistently in three dimensions using the combined N-body/hydrodynamic scheme P3MSPH up to a final redshift z = 1. Evolving to the present is prohibited by the fact that the mean density in the simulated volume is above critical and the entire volume would be going nonlinear beyond this point, We are currently analyzing another run with somewhat poorer mass resolution which was evolved to the present.

Evrard, August E.↗