Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure correlations”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Computational fluid dynamics modeling of high level waste plug formation - 15228

At Hanford, an extensive network of pipelines traversing several miles is used to transfer the high level radioactive nuclear waste from tanks to the treatment facilities. During transfer operations, however, there is a potential risk for the transport lines to plug, causing significant delays, increasing operation costs and creating hazardous conditions for personnel and the environment. In this study we investigate the phenomena of plug formation in the pipelines due to solids depositing during the waste transfer process. A three-dimensional (3-D) Computational Fluid Dynamics (CFD) model has been developed using Comsol Multiphysics 4.3b by Florida International University (FIU) to evaluate the factors that can contribute to pipeline plug formation. Several virtual scenarios are developed simulating the settling mechanics as a function of particle size, solids volume fraction and solids density. The numerical results are validated with the available experimental results and the established critical velocity correlations. The 3-D numerical results were also compared with the previous 2-D numerical studies to understand the trade-off between the two studies in terms of computing speed and numerical accuracy. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗
DOE OSTI2025

Implications of Aerosol Physicochemical Properties Including Ice Nucleation at ARM Mega Sites for Improved Understanding of Microphysical Atmospheric Cloud Processes (Final Technical Report)

Continuous, long-term measurements of atmospheric ice-nucleating particles (INPs) that influence clouds and precipitation were conducted to investigate the abundance and variability of ground-level INPs across the world. Three field campaigns were organized by the DOE Atmospheric Radiation Measurement (ARM) program, including Examining INP from Southern Great Plains (ExINP-SGP, 2019), Examining INP from Eastern North Atlantic (ExINP-ENA, 2020 – 2021), and Examining INP from North Slope of Alaska (ExINP-NSA, 2021 – 2024). Additional small-scale supporting field experiments were performed in 2019 and 2021 to collect airborne particulate matter at SGP for complementary laboratory characterization of the particles’ physical and chemical properties [Aerosol-Ice Formation Closure Pilot Study (AEROICESTUDY), 2019; ExINP-SGP II, 2021]. In these studies, the PI’s team measured INP concentration with both real-time and laboratory measurements in a wide range of freezing temperatures ($T$ from 0 °C to about –30 °C). This project elucidated spatial variability and seasonality in the abundance of immersion mode active INPs across three ARM sites using a single instrument, a Portable Ice Nucleation Experiment (PINE) chamber version 03 (PINE-03 hereafter). Collocated aerosol and meteorological data were analyzed to assess the correlation between ambient INP abundance, air mass origin region, and meteorological variability. Our findings suggest very high freezing efficiency of INPs at the NSA site across the measured temperatures (ice nucleation active surface site density, $n_s(T)$, $\approx 2 \times 10^{8} - 10^{10}$ m -2 for from –16 to –31 °C), which is a factor of 10 – 1000 times greater efficiency as compared to that found in the previous mid-latitude INP measurements in autumn using the same instrument; surprisingly high INP abundance ($\ge 1 \text{ L}^{-1}$ at –25 °C) for the examined temperatures throughout the year that PINE-03 did not measure at other sites; and high INP concentration in spring, possibly related to arctic haze episodes.

58 GEOSCIENCES↗
NASA NTRS1996

Crack Development in Cross-Ply Laminates Under Uniaxial Tension

This study addresses matrix-dominated failures in carbon fiber/polymer matrix composite laminates in a cross-ply lay-up. The events of interest are interlaminar fracture in the form of transverse cracks in the 90' plies and longitudinal splitting in the 0 deg plies and interlaminar fracture in the form of 0/90 delamination. These events were observed using various nondestructive evaluation (NDE) techniques during static tensile tests. Acoustic emission (AE) x radiography, and edge view microscopy were the principal ones utilized in a real-time environment. A comparison of the NDE results with an analytical model based on the classical linear fracture mechanics concept of strain energy release rate as a criterion for crack growth was performed. The virtual crack closure theory was incorporated with a finite element model to generate strain energy release rate curves for the analytical case. Celion carbon fiber/polyimide matrix (G30-500/PMR-15) was the material tested with cross-ply lay-ups of (0(2)/90(6))s and (0(4)/90(4))s. The test specimens contained thermally induced cracks caused by the high-temperature processing. The analytical model was updated to compensate for the initial damage and to study further accumulation by taking into account the crack interactions. By correlating the experimental and analytical data, the critical energy release rates were found for the observable events of interest.

Gyekenyesi, Andrew L.↗
NASA NTRS1993

2D axisymmetric analysis of SRM ignition transient

To analyze ignition transient of Space Shuttle solid rocket motor, a transient two-dimensional numerical model based on turbulent compressible Navier-Stokes equations in a generalized coordinate system was developed. One-dimensional numerical models (Peretz et al., 1973; Han, 1992; Pardue and Han, 1992) with empirical correlations data obtained from steady turbulent boundary layer flows agrees reasonably well with test rocket data by adjusting a few parameters. However, a 1D model can not provide a physical insight into the complex multidimensional thermal fields and flowfields in the chamber and the converging-diverging rocket nozzle. As an interim step, a 2D model was developed and compared with test data. A modified version of SIMPLE algorithm was used for the numerical model, and the standard k-epsilon model with a wall function was used for turbulence closure. Transient flowfields and thermal fields in the combustion chamber and the attached nozzle were obtained for a selected rocket geometry and propellant. Transient behaviors of the flow and thermal fields were analyzed, and were found to be in good agreement with physical expectations.

Bai, S. D.↗
DOE OSTI2018

Center for Extended Magnetohydrodrynamic Modeling

This work used the NIMROD code to solve two tokamak regimes found on the DIII-D tokamak. The first regime studied is the Quiescent H-mode regime that has the advantage of enabling tokamak high performance without the deleterious effects of Edge Localized Modes (ELMs). ELMs are a possible impediment to the success of tokamaks as a burning plasma reactor and the QH mode regime provides a mechanism for avoiding it. The QH mode achieves this state due to the appearance of long wavelength instability that saturates instead of reaching an ELM state. NIMROD simulations were able to describe many of the features of a particular discharge, including the importance of shear flow in reaching the state, the role of the parallel heat flux to decorrelate the temperature and velocity fluctations, and in the correlation of density and velocity fluctations to provide the density pumpout. This work resulted in two direct publications and an invited talk. Progress towards these publications resulted in two more publications (verification and solving the Grad-Shafranov equation with scrape-off-layer currents). The second regime studied was the Giant Sawtooth regime whereby energetic particles affect the stability of the 1/1 mode. This work resulted in performance improvements of NIMROD's energetic particle algorithms, and improvements to the drift-kinetic equation infrastructure used by both continuum and energetic particle closure methods. The studies proved more computationally demanding than expected because as the 1/1 is stabilized, the near marginal conditions cause long simulations times. Publication of this work is in progress under other funding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Waterhammer Transient Simulation and Model Anchoring for the Robotic Lunar Lander Propulsion System

Waterhammer transients have the potential to adversely impact propulsion system design if not properly addressed. Waterhammer can potentially lead to system plumbing, and component damage. Multi-thruster propulsion systems also develop constructive/destructive wave interference which becomes difficult to predict without detailed models. Therefore, it is important to sufficiently characterize propulsion system waterhammer in order to develop a robust design with minimal impact to other systems. A risk reduction activity was performed at Marshall Space Flight Center to develop a tool for estimating waterhammer through the use of anchored simulation for the Robotic Lunar Lander (RLL) propulsion system design. Testing was performed to simulate waterhammer surges due to rapid valve closure and consisted of twenty-two series of waterhammer tests, resulting in more than 300 valve actuations. These tests were performed using different valve actuation schemes and three system pressures. Data from the valve characterization tests were used to anchor the models that employed MSCSoftware.EASY5 v.2010 to model transient fluid phenomena by using transient forms of mass and energy conservation. The anchoring process was performed by comparing initial model results to experimental data and then iterating the model input to match the simulation results with the experimental data. The models provide good correlation with experimental results, supporting the use of EASY5 as a tool to model fluid transients and provide a baseline for future RLL system modeling. This paper addresses tasks performed during the waterhammer risk reduction activity for the RLL propulsion system. The problem of waterhammer simulation anchoring as applied to the RLL system is discussed with results from the corresponding experimental valve tests. Important factors for waterhammer mitigation are discussed along with potential design impacts to the RLL propulsion system.

Stein, William B.↗
DOE OSTI2024

Pressure-tuned quantum criticality in the large-D antiferromagnet DTN

Strongly correlated spin systems can be driven to quantum critical points via various routes. In particular, gapped quantum antiferromagnets can undergo phase transitions into a magnetically ordered state with applied pressure or magnetic field, acting as tuning parameters. These transitions are characterized by z = 1 or z = 2 dynamical critical exponents, determined by the linear and quadratic low-energy dispersion of spin excitations, respectively. Employing high-frequency susceptibility and ultrasound techniques, we demonstrate that the tetragonal easy-plane quantum antiferromagnet NiCl 2 · 4SC(NH 2 ) 2 (aka DTN) undergoes a spin-gap closure transition at about 4.2 kbar, resulting in a pressure-induced magnetic ordering. The studies are complemented by high-pressure-electron spin-resonance measurements confirming the proposed scenario. Powder neutron diffraction measurements revealed that no lattice distortion occurs at this pressure and the high spin symmetry is preserved, establishing DTN as a perfect platform to investigate z = 1 quantum critical phenomena. The experimental observations are supported by DMRG calculations, allowing us to quantitatively describe the pressure-driven evolution of critical fields and spin-Hamiltonian parameters in DTN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Modeling the pressure-dilatation correlation

It is generally accepted that pressure dilatation, which is an additional compressibility term in turbulence transport equations, may be important for high speed flows. Recent direct simulations of homogeneous shear turbulence have given concrete evidence that the pressure dilatation is important insofar that it contributes to the reduced growth of turbulent kinetic energy due to compressibility effects. The problem of modeling pressure dilatation is addressed. A component of the pressure dilatation is isolated which exhibits temporal oscillations and, using direct numerical simulations of homogeneous shear turbulence and isotropic turbulence, show that it has a negligible contribution to the evolution of turbulent kinetic energy. Then, an analysis for the case of homogeneous turbulence is performed to obtain a model for the nonoscillatory pressure dilatation. This model algebraically relates the pressure dilatation to quantities traditionally obtained in incompressible turbulence closures. The model is validated by direct comparison with the pressure dilatation data obtained from the simulations.

Sarkar, S.↗
NASA NTRS2004

Investigation of Navier-Stokes Code Verification and Design Optimization

With rapid progress made in employing computational techniques for various complex Navier-Stokes fluid flow problems, design optimization problems traditionally based on empirical formulations and experiments are now being addressed with the aid of computational fluid dynamics (CFD). To be able to carry out an effective CFD-based optimization study, it is essential that the uncertainty and appropriate confidence limits of the CFD solutions be quantified over the chosen design space. The present dissertation investigates the issues related to code verification, surrogate model-based optimization and sensitivity evaluation. For Navier-Stokes (NS) CFD code verification a least square extrapolation (LSE) method is assessed. This method projects numerically computed NS solutions from multiple, coarser base grids onto a freer grid and improves solution accuracy by minimizing the residual of the discretized NS equations over the projected grid. In this dissertation, the finite volume (FV) formulation is focused on. The interplay between the xi concepts and the outcome of LSE, and the effects of solution gradients and singularities, nonlinear physics, and coupling of flow variables on the effectiveness of LSE are investigated. A CFD-based design optimization of a single element liquid rocket injector is conducted with surrogate models developed using response surface methodology (RSM) based on CFD solutions. The computational model consists of the NS equations, finite rate chemistry, and the k-6 turbulence closure. With the aid of these surrogate models, sensitivity and trade-off analyses are carried out for the injector design whose geometry (hydrogen flow angle, hydrogen and oxygen flow areas and oxygen post tip thickness) is optimized to attain desirable goals in performance (combustion length) and life/survivability (the maximum temperatures on the oxidizer post tip and injector face and a combustion chamber wall temperature). A preliminary multi-objective optimization study is carried out using a geometric mean approach. Following this, sensitivity analyses with the aid of variance-based non-parametric approach and partial correlation coefficients are conducted using data available from surrogate models of the objectives and the multi-objective optima to identify the contribution of the design variables to the objective variability and to analyze the variability of the design variables and the objectives. In summary the present dissertation offers insight into an improved coarse to fine grid extrapolation technique for Navier-Stokes computations and also suggests tools for a designer to conduct design optimization study and related sensitivity analyses for a given design problem.

Vaidyanathan, Rajkumar↗
DOE OSTI2016

CASL COBRA-TF improvements for PWR DNB analysis

COBRA-TF (CTF) is a thermal hydraulic (T/H) subchannel code using either three-dimensional (3D) Cartesian or subchannel coordinate formulations for two-phase fluid flow and heat transfer solutions. CTF has been improved under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program for Pressurized Water Reactor (PWR) applications, including software optimization, new closure models, pre- and post-processing and parallelization for modeling full reactor core T/H responses under normal operating and accident conditions. As a result of collaboration among CASL partners including the Westinghouse Electric Company, the Oak Ridge National Laboratory (ORNL), and the Sandia National Laboratories, additional modeling improvements were made to CTF specifically for PWR Departure from Nucleate Boiling (DNB) analysis, including a code option to evaluate fuel thermal margin in terms of DNB Ratio (DNBR) and an axial shape factor to account for effect of non-uniform axial power distribution on DNB. Multiple DNB correlations are now linked with CTF for different applications, including the Westinghouse proprietary WRB-1 correlation for fuel designs containing mixing vane grid spacers. The improved CTF code with the WRB-1 correlation (CTF/WRB-1) was validated using the DNB data from the PWR Subchannel Bundle Tests (PSBT). In addition to the comparison with the test data, the CTF/WRB-1 DNBR results and the associated local fluid conditions were compared to the results of the Westinghouse T/H design code, VIPRE-W, which is an enhanced version of the VIPRE-01 code originally developed by the Electric Power Research Institute (EPRI). The comparisons showed that CTF/WRB-1 DNBR predictions are in good agreement with the VIPRE-W results within the applicable range of the DNB correlation. A model sensitivity study was performed to confirm that the CTF void drift model had an insignificant effect on DNBR under the steam line break (SLB) low pressure condition. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Magnetic field-tuned Fermi liquid in a Kondo insulator

Abstract Kondo insulators are expected to transform into metals under a sufficiently strong magnetic field. The closure of the insulating gap stems from the coupling of a magnetic field to the electron spin, yet the required strength of the magnetic field–typically of order 100 T–means that very little is known about this insulator-metal transition. Here we show that Ce $${}_{3}$$ 3 Bi $${}_{4}$$ 4 Pd $${}_{3}$$ 3 , owing to its fortuitously small gap, provides an ideal Kondo insulator for this investigation. A metallic Fermi liquid state is established above a critical magnetic field of only $${B}_{{\rm{c}}}\approx$$ B c ≈ 11 T. A peak in the strength of electronic correlations near $${B}_{{\rm{c}}}$$ B c , which is evident in transport and susceptibility measurements, suggests that Ce $${}_{3}$$ 3 Bi $${}_{4}$$ 4 Pd $${}_{3}$$ 3 may exhibit quantum criticality analogous to that reported in Kondo insulators under pressure. Metamagnetism and the breakdown of the Kondo coupling are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Initial use of Nek5000/Cardinal to improve closure models in Pronghorn

Heat transfer coefficient closure models for pebble bed reactors are built using a data-driven approach by leveraging the capabilities of an Evolutionary Algorithm entitled Particle Swarm Optimization (PSO). In the present work, the Computational Fluid Dynamics code nekRS was used in order to collect the high-fidelity flow data for a core with 1,568 pebbles. To characterize the heat transfer, multiple concentric regions were considered to extract the physical quantities of interest, e.g./ the Reynolds number. The PSO algorithm is employed as part of an inverse problem targeting determine what are the coefficients for a Nusselt number correlation to match the collected data. Such correlation should follow any given format that is defined a priori. Finally, two correlations are proposed, one with an implicit dependence on the pebbles’ wall temperatures and another expressed as a fully explicit correlation depending on the flow conditions and the position within the core. Anyway, given the generic nature of the proposed approach, correlations following different formats could be tested. Preliminary results for the high-fidelity simulation of a fast MSR core are presented. The target Reynolds number is currently 20K, with the expectation that this will increase, pending the availability of further computational resources. These simulations will be used to inform lower fidelity models, including a coarse CFD turbulence model in Pronghorn. Additionally, they will serve as a reference for the RANS models in Nek5000/NekRS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Environmental controls on isolated convection during the Amazonian wet season

The Amazon rainforest is a vital component of the global climate system, influencing the hydrological cycle and tropical circulation. However, understanding and modeling the evolution of convection in this region remain a scientific challenge. Here, we assess the environmental conditions associated with shallow, congestus, and isolated deep convection days during the wet season (December to April), employing measurements from the Green Ocean Amazon 2014–2015 (GoAmazon2014/5) experiment and large-scale wind fields from the constrained variational analysis. Composites of deep days show moister than average conditions below 3 km early in the morning. Analyzing the water budget at the surface through observations only, we estimated the water vapor convergence term as a residual of the water balance closure. Convergence remains nearly zero during the deep days until early afternoon (13:00 LST), when it becomes a dominant factor in the water budget. At 14:00 LST, the deep days experience a robust upward large-scale vertical velocity, especially above 4 km, which supports the shallow-to-deep convective transition occurring around 16:00–17:00 LST. In contrast, shallow and congestus days exhibit drier pre-convective conditions, along with diurnal water vapor divergence and large-scale subsidence that extend from the surface to the lower free troposphere. Moreover, afternoon precipitation exhibits the strongest linear correlation (0.6) with large-scale vertical velocity, nearly double the magnitude observed for other environmental factors, even moisture, at different levels and periods of the day. Precipitation also exhibits a moderate increase with low-level wind shear, while upper-level shear has a relatively minor negative impact on convection.

Environmental Sciences & Ecology↗
NASA NTRS1993

Results of the second Round Robin on opening-load measurement conducted by ASTM Task Group E24.04.04 on crack closure measurement and analysis

A second experimental Round Robin on the measurement of the crack opening load in fatigue crack growth tests has been completed by the ASTM Task Group E24.04.04 on Crack Closure Measurement and Analysis. Fourteen laboratories participated in the testing of aluminum alloy compact tension specimens. Opening-load measurements were made at three crack lengths during constant Delta K, constant stress ratio tests by most of the participants. Four participants made opening-load measurements during threshold tests. All opening-load measurements were based on the analysis of specimens compliance behavior, where the displacement/strain was measured either at the crack mouth or the mid-height back face location. The Round Robin data were analyzed for opening load using two non-subjective analysis methods: the compliance offset and the correlation coefficient methods. The scatter in the opening load results was significantly reduced when some of the results were excluded from the analysis population based on an accept/reject criterion for raw data quality. The compliance offset and correlation coefficient opening load analysis methods produced similar results for data populations that had been screened to eliminate poor quality data.

Phillips, E. P.↗
NASA NTRS2012

Development and Applications of a Stage Stacking Procedure

The preliminary design of multistage axial compressors in gas turbine engines is typically accomplished with mean-line methods. These methods, which rely on empirical correlations, estimate compressor performance well near the design point, but may become less reliable off-design. For land-based applications of gas turbine engines, off-design performance estimates are becoming increasingly important, as turbine plant operators desire peaking or load-following capabilities and hot-day operability. The current work develops a one-dimensional stage stacking procedure, including a newly defined blockage term, which is used to estimate the off-design performance and operability range of a 13-stage axial compressor used in a power generating gas turbine engine. The new blockage term is defined to give mathematical closure on static pressure, and values of blockage are shown to collapse to curves as a function of stage inlet flow coefficient and corrected shaft speed. In addition to these blockage curves, the stage stacking procedure utilizes stage characteristics of ideal work coefficient and adiabatic efficiency. These curves are constructed using flow information extracted from computational fluid dynamics (CFD) simulations of groups of stages within the compressor. Performance estimates resulting from the stage stacking procedure are shown to match the results of CFD simulations of the entire compressor to within 1.6% in overall total pressure ratio and within 0.3 points in overall adiabatic efficiency. Utility of the stage stacking procedure is demonstrated by estimation of the minimum corrected speed which allows stable operation of the compressor. Further utility of the stage stacking procedure is demonstrated with a bleed sensitivity study, which estimates a bleed schedule to expand the compressors operating range.

Kulkarni, Sameer↗
NASA NTRS1992

The pressure-dilatation correlation in compressible flows

Simulations of simple compressible flows have been performed to enable the direct estimation of the pressure-dilatation correlation. The generally accepted belief that this correlation may be important in high-speed flows has been verified by the simulations. The pressure-dilatation correlation is theoretically investigated by considering the equation for fluctuating pressure in an arbitrary compressible flow. This leads to the isolation of a component of the pressure-dilatation that exhibits temporal oscillations on a fast time scale. Direct numerical simulations of homogeneous shear turbulence and isotropic turbulence show that this fast component has a negligible contribution to the evolution of turbulent kinetic energy. Then, an analysis for the case of homogeneous turbulence is performed to obtain a formal solution for the nonoscillatory pressure-dilatation. Simplifications lead to a model that algebraically relates the pressure-dilatation to quantities traditionally obtained in incompressible turbulence closures. The model is validated by direct comparison with the simulations.

Sarkar, S.↗
DOE OSTI2011

Electron-impact study of the S{sub 2} molecule using the R-matrix method

The present study deals with the calculation of elastic [integrated and differential cross section (DCS)], momentum-transfer, excitation, and ionization cross sections for electron impact on S{sub 2} molecules using the R-matrix method. The target states are represented by including correlations via a configuration-interaction technique. We used a double zeta plus polarization Gaussian basis set contracted as (12,8,1)/(6,4,1) for S atoms. The results of the static exchange, correlated one-state, and 20-state close-coupling approximations are presented. We have detected a stable anionic bound state {sup 2}{Pi}{sub g} of S{sub 2}{sup -} having the configuration 1{sigma}{sub g}{sup 2}{center_dot}{center_dot}{center_dot}5{sigma}{sub g}{sup 2} 1{sigma}{sub u}{sup 2}{center_dot}{center_dot}{center_dot}4{sigma}{sub u}{sup 2} 1{pi}{sub u}{sup 4} 1{pi}{sub g}{sup 4} 2{pi}{sub u}{sup 4} 2{pi}{sub g}{sup 3}. The vertical electron affinity value is 1.42 eV, which is comparable with the experimental value of 1.67{+-}0.015 eV. We detected two shape resonances, both of {sup 2}{Pi}{sub u} symmetry in the excitation cross sections of the {sup 1}{Delta}{sub g} and {sup 1}{Sigma}{sub g}{sup +} excited states. The dissociative nature of these resonances is explored by performing scattering calculations in which the S-S bond is stretched. These resonances support dissociative attachment, yielding S and S{sup -}. We have also predicted six resonances of various symmetries ({sup 2}A{sub u}, {sup 2}B{sub 1g},{sup 4}A{sub u}, {sup 4}B{sub 1g}) in the X {sup 3}{Sigma}{sub g}{sup -} {yields} B {sup 3}{Sigma}{sub u}{sup -} transition. We have calculated the DCS, in a correlated one-state model, by using the polydcs program of Sanna and Gianturco. The data from the momentum-transfer cross section, generated from DCS, are used to compute effective collision frequencies over a wide electron temperature range (200-30 000 K). The ionization cross sections are calculated in the binary-encounter Bethe model in which Hartree-Fock molecular orbitals at a self-consistent level are used to calculate kinetic and binding energies of the occupied molecular orbitals. We have included up to g-partial wave (l=4) in the scattering calculations. For this molecule we have used a Born-closure top-up procedure to account for the higher partial waves for the convergence of the cross section for the dipole-allowed excitation from the ground state. We have also evaluated the scattering length of the S{sub 2} molecule, which is equal to 2.615a{sub 0}.

74 ATOMIC AND MOLECULAR PHYSICS↗