Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak estimation”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗
DOE OSTI2023

Late-time post-merger modeling of a compact binary: effects of relativity, r-process heating, and treatment of transport

Detectable electromagnetic counterparts to gravitational waves from compact binary mergers can be produced by outflows from the black hole-accretion disk remnant during the first 10 s after the merger. Two-dimensional axisymmetric simulations with effective viscosity remain an efficient and informative way to model this late-time post-merger evolution. In addition to the inherent approximations of axisymmetry and modeling turbulent angular momentum transport by a viscosity, previous simulations often make other simplifications related to the treatment of the equation of state and turbulent transport effects. In this paper, we test the effect of these modeling choices. By evolving with the same viscosity the exact post-merger initial configuration previously evolved in Newtonian viscous hydrodynamics, we find that the Newtonian treatment provides a good estimate of the disk ejecta mass but underestimates the outflow velocity. We find that the inclusion of heavy nuclei causes a notable increase in ejecta mass. An approximate inclusion of r-process effects has a comparatively smaller effect, except for its designed effect on the composition. Diffusion of composition and entropy, modeling turbulent transport effects, has the overall effect of reducing ejecta mass and giving it a speed with lower average and more tightly-peaked distribution. Also, we find significant acceleration of outflow even at distances beyond 10 000 km, so that thermal wind velocities only asymptote beyond this radius and at higher values than often reported.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Relevance of E × B drifts for particle and heat transport in divertors

Radial electric fields up to ~4 kV m –1 are observed in the boundary between the private flux region (PFR) and the scrape-off layer (SOL) driving E × B drifts between the inner and outer targets at speeds up to 2.8 km s –1 in the Tokamak à configuration variable divertor. The resulting E × B fluxes, located in a narrow region ($\Delta {\rho _\Psi } < 0.012$ in normalized radius or $\Delta $ R – R sep <4 mm mapped to the outer midplane) are equivalent to around 20% of the total heat and particle flux to the divertor targets (inner + outer). At the peak E r , the E × B poloidal transport is equivalent to parallel flows with M ∥ ~ 3. In the snowflake divertor with a second X-point in the outer SOL, the drifts in the PFR-SOL boundary were equivalent to around 30% of the total heat and particle flux to the divertor targets and cover a region ~50% wider than in the single null ($\Delta {\rho _\Psi }$ ~ 0.018, $\Delta $R – R sep ~ 6 mm). Furthermore, the location of the PFR-SOL boundary drift shifts radially in the E ∥ × B direction when reversing the toroidal field direction. Peaks in density and electron pressure have been identified near the primary X-point along with large gradients in density, temperature, and potential, the latter resulting in a local electric field ~2.7 kV m –1 which drives a drift (1.9 km s –1 ) upwards towards the closed flux surfaces. Floating potential (V f ) magnitudes up to 75 V (~2 kT e ) were measured, indicating that V f and parallel currents should not be neglected when estimating plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Milky Way Satellite Census -- II. Galaxy-Halo Connection Constraints Including the Impact of the Large Magellanic Cloud

The population of Milky Way (MW) satellites contains the faintest known galaxies, and thus provides essential insight into galaxy formation and dark matter microphysics. Here, we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high-Galactic-latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6.5\pm 1.5$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from $\textit{Gaia}$ proper motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $2.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Corresponding Standard Reference Material Data used in Partial Least Squares Regression Models for Sugar Composition Estimates in Biomass in: Economic Impact of Yield and Composition Variation in Bioenergy Crops: Populus trichocarpa

Corresponding Standard Reference Material Data used in Partial Least Squares Regression Models for Sugar Composition Estimates in Biomass in: Economic Impact of Yield and Composition Variation in Bioenergy Crops: Populus trichocarpa (for corresponding manuscript: DOI: 10.1002/bbb.2148) PDF Files: Images of 1H NMR spectra for neutralized 2-stage acid hydrolysates of 4 NIST Standard Reference Material biomass samples (Monterey Pine 8493, Sugarcane Bagasse 8491, Wheat Straw 8494, and Eastern Cottonwood/Poplar 8492) and 2 Center for Bioenergy Innovation reference biomass samples (Poplar - Populus trichocarpa and Switchgrass - Panicum Virgatum). Suppression of the water peak was achieved using a NOESY-1D with presaturation, a recycle delay of 5 s, and a total of 64 scans. Spectra were acquired at 298 K and processed with automatic phase correction, baseline correction, and chemical shift referencing to TSP-d4. Images show all 1H data from 10 to 1ppm with inset spectra of region of interest (4.0 to 3.1 ppm). Text Files: Spectra for neutralized 2-stage acid hydrolysates of 4 NIST Standard Reference Material biomass samples (Monterey Pine 8493, Sugarcane Bagasse 8491, Wheat Straw 8494, and Eastern Cottonwood/Poplar 8492) and 2 Center for Bioenergy Innovation reference biomass samples (Poplar - Populus trichocarpa and Switchgrass - Panicum Virgatum) were converted into text files for plotting. Files contain 8192 points of raw spectral data from 12.23 to -2.78 ppm. The text file contains 4 columns of data and includes: Point number, Intensity, Hz, and ppm. Xcel Spreadsheet: HPLC measured monomeric sugar concentrations and bucketed 1H NMR data used to build monomeric sugar composition prediction models. Sugar composition in biomass determined from HPLC analyses are given in mg sugar/mg of biomass. Spectral bucketing was performed using Bruker’s AMIX software. Spectra were divided into 0.005 ppm buckets in the region of 3.10– 4.15 ppm for a total of 210 buckets. Headers for the bucketed data are the chemical shift in ppm of the center of the bucket. Bucketed data was used to build partial least squares models for subsequent predictions in The Unscrambler v. 10.5(CAMO A/S, Trondheim, Norway). The formation of methanol during hydrolysis interferes with the quantitative NMR analysis of sugars, so the methanol peak centered at 3.37 ppm and spanning four buckets (3.2925 – 3.2775 ppm) was set to zero for all spectra.

09 BIOMASS FUELS↗
NASA NTRS1976

Multistream and Monte Carlo calculations of the sun's aureole

The scattered UV intensities in the sun's aureole are calculated using both a multi-stream scattering method and a Monte Carlo approach. Angular distributions for both Rayleigh and aerosol scatterings are obtained with a realistic atmospheric model. Moderate- and sharply-forward-peaked phase functions of aerosol scattering, corresponding to realistic analytic size distributions, are incorporated. Results from the two independent calculations are in reasonable agreement for a realistic atmospheric model. The results indicate that the scattered UV intensity in the sun's aureole is about four orders smaller than the transmitted intensity while the scattered intensity for pure Rayleigh scattering is only 1-10 millionths times that of the transmitted sun. Based on these intensity ratios, we estimate that the integrated scattered contributions from the aureole to a well-collimated sun photometer of acceptance aperture 2-3 deg are below about 1% of the direct contribution.

Furman, D. R.↗
NASA NTRS1999

Satellite Altimetry, Ocean Circulation, and Data Assimilation

Ocean circulation is a critical factor in determining the Earth's climate. Satellite altimetry has been proven a powerful technique for measuring the height of the sea surface for the study of global ocean circulation dynamics. A major objective of my research is to investigate the utility of altimeter data for ocean circulation studies. The 6 years' data record of TOPEX/POSEIDON have been analyzed to study the spatial and temporal characteristics of large-scale ocean variability. A major result obtained in 1998 is the discovery of large-scale oscillations in sea level with a period of 25 days in the Argentine Basin of the South Atlantic Ocean (see diagram). They exhibit a dipole pattern with counterclockwise rotational propagation around the Zapiola Rise (centered at 45S and 317E), a small seamount in the abyssal plain of the basin. The peak-to-trough amplitude is about 10 cm over a distance of 500-1000 km. The amplitude of these oscillations has large seasonal-to-interannual variations. The period and rotational characteristics of these oscillations are remarkably similar to the observations made by two current meters deployed near the ocean bottom in the region. What TOPEX/POSEIDON has detected apparently are manifestations of the movement of the entire water column (barotropic motion). The resultant transport variation is estimated to be about 50 x 10(exp 6) cubic M/S, which is about 50% of the total water transport in the region. Preliminary calculations suggest that these oscillations are topographically trapped waves. A numerical model of the South Atlantic is used to investigate the nature of and causes for these waves. A very important property of sea surface height is that it is directly related to the surface geostrophic velocity, which is related to deep ocean circulation through the density field. Therefore altimetry observations are not only useful for determining the surface circulation but also for revealing information about the deep ocean. Another objective of my research is to investigate the extent of information one can infer about the deep ocean from altimetric observations. The approach is data assimilation by ocean models. The objective is to develop a modeling/data assimilation system that will produce estimates of the 3-dimensional state of the ocean for the entire duration of the TOPEX/POSEIDON Mission.

Fu, Lee-Lueng↗
NASA NTRS2007

Hinode Observations of the Onset Stage of a Solar Filament Eruption

We used Hinode X-Ray Telescope (XRT) and Solar Optical Telescope (SOT) filtergraph (FG) Stokes-V magnetogram observations, to study the early onset of a solar eruption that includes an erupting filament that we observe in TRACE EUV images. The filament undergoes a slow rise for at least 20min prior to its fast eruption and strong soft X-ray (SXR) flaring; such slow rises have been previously reported, and the new Hinode data elucidate the physical processes occurring during this period. XRT images show that during the slow-rise phase, an SXR sigmoid forms from apparent reconnection low in the sheared core field traced by the filament, and there is a low-level intensity peak in both EUV and SXRs during the slow rise. MDI and SOT FG Stokes-V magnetograms show that the pre-emption filament is along a neutral line between opposing-polarity enhanced network cells, and the SOT magnetograms show that these opposing fields are flowing together and canceling for at least six hours prior to eruption. From the MDI data we measured the canceling network fields to be approx. 40 G, and we estimated that approx. 10(exp 19)Mx of flux canceled during the five hours prior to eruption; this is only approx.5% of the total flux spanned by the eruption and flare, but apparently its tether-cutting cancellation was enough to destabilize the sigmoid field holding the filament and resulted in that field's eruption.

Sterling, Alphonse C.↗
NASA NTRS2016

MAVEN Observations of Magnetic Flux Ropes with a Strong Field Amplitude in the Martian Magnetosheath During the ICME Passage on 8 March 2015

We present initial results of strong field amplitude flux ropes observed by Mars Atmosphere and Volatile EvolutioN (MAVEN) mission around Mars during the interplanetary coronal mass ejection (ICME) passage on 8 March 2015. The observed durations were shorter than 5 s and the magnetic field magnitudes peaked above 80 nT, which is a few times stronger than those usually seen in the magnetosheath barrier. These are the first unique observations that MAVEN detected such flux ropes with a strong field at high altitudes (greater than 5000 km). Across these structures, MAVEN coincidentally measured planetary heavy ions with energies higher than a few keV. The spatial properties inferred from the Grad-Shafranov equation suggest that the speed of the structure can be estimated at least an order of magnitude faster than those previously reported quiet-time counterparts. Hence, the space weather event like the ICME passage can be responsible for generating the observed strong field, fast-traveling flux ropes.

8 March 2015↗
NASA NTRS2020

Tidal Disruptions of Main-sequence Stars. IV. Relativistic Effects and Dependence on Black Hole Mass

Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass( MM105 105BH7)and stellar mass (MM0.33) as larger MBH leads to stronger relativistic effects. For fixed Må, as MBH increases, the ratio of the maximum pericenter distance yielding full disruptions (t) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for= ́MM510BH7, while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of 2 correction at= ́MM510BH7. We provide approximate formulae that express the relativistic corrections of both t and the energy width relative to their Newtonian approximate estimates. For partial disruptions, we find that the fractional remnant mass for a given ratio of the pericenter to t is higher for larger MBH. These results have several implications. As MBH increases above~M107, the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross-section ratio for partial to complete disruptions depends only weakly on MBH. The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. ForMM10BH7, the MBH-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics.

Taeho Ryu↗
NASA NTRSDate not supplied

Nanoflare Heating Frequency of an X-ray Bright Point Observed by MaGIXS

Nanoflares have been considered to be one of the most likely candidates for heating the solar corona to multi-million kelvin temperatures. Individual nanoflares are difficult to detect with today's instruments, but their presence may be established by comparing simulated nanoflare-heated plasma emissions to observed emissions. We present a simulation of emission from an X-ray Bright Point (XBP) detected by the MaGIXS, as well as simultaneous observations from SDO/AIA and Hinode/XRT. To simulate the XBP loops, we utilize the HYDRAD code. The length and magnetic field strength of these loops are determined using potential field extrapolation of SDO/HMI's observed photospheric magnetogram. Each loop is considered to be heated by random nanoflares, the amplitude and frequency of which are governed by the length of the loop and the strength of the magnetic field. The simulated outputs are used to estimate the intensity of spectrally pure maps of Fe-18, Fe-17, Ne-9, O-8, O-9, Ne-9, and so on, which is then compared to the intensity determined from MaGIXS observations. In addition, we derived the intensity maps obtained by AIA and XRT and compared them to the observed data. The composite distribution of the delay time of the nanoflares for which the simulated loops morphology and intensities match with observation shows a peak at 200s-500s, indicating that most of the nanoflares have a high/intermediate frequency.

coronal heating↗
DOE OSTI2026

Local Air Quality Impacts of a Peak-Shaving Diesel Generator Unit in Waynesville, North Carolina

The Waynesville Electric Department in Waynesville, North Carolina, operates a 2,000-kW diesel generator with a 4% capacity factor, emitting less than 1 ton annually of PM2.5, VOCs, carbon monoxide (CO), and sulfur dioxide (SO2), and approximately 1.8 tons of NO. The National Laboratory of the Rockies (NLR) conducted high-resolution, near-source air-quality modeling to evaluate the generator's impacts. Results indicate minimal effects on ambient pollutant concentrations: annual PM2.5 increases are below 0.0004 microgrm/m3 for local census tracts, while CO and SO2 levels increase by less than 0.00004 ppm and 0.0002 ppb, respectively. Estimated annual premature mortality attributable to PM2.5 exposure is less than 0.001 cases across Haywood County, with an associated economic impact of approximately $12,000. This report complements a separate technoeconomic analysis to inform Waynesville's investment strategies for demand reduction.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Assessment of thermal and radiation induced creep in the dual cooled lead lithium blanket

The creep in recently designed blanket for the fusion nuclear science facility is described by decomposing it into thermal, irradiation and cavity swelling creep. Initial estimates of thermal creep using elastic high temperature rules showed that the recent heat transfer improvements at the first wall (FW) is feasible, and the FW can withstand heat flux of 0.35 MW/m 2 at 3000 hrs without rupture. A viscoplastic model that is based on Norton's power law and relies on Multiphysics coupling of solid mechanics and heat transfer modules is used to capture the inelastic deformation in the blanket. Results showed that the design peak heat flux of 0.25 MW/m 2 produced maximum thermal creep strain of 0.45 % and the relaxation of the thermal stress at the FW. Irradiation creep is prescribed to be proportional to the displacement damage dose and applied stress. The displacements from irradiation creep radially decrease from the FW to the back wall, but the maximum deformation is found at the back wall that is connected to the Helium manifold due to the high stresses at the region. In conclusion, cavity swelling creep is dominant at steady state and combined with radiation creep to produce displacement of about 8 mm at the FW.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Evaluation of the absolute single-photon detection efficiency of HRPPD

Pixelated High Rate Picosecond Photon Detectors (HRPPDs) by Incom Inc. are promising photosensors for use in Ring Imaging CHerenkov (RICH) detectors, where a high gain, sub-mm position resolution and sub-100ps timing resolution are required in a single photon mode. Quantum Efficiency (QE) has been measured for the first batch of EIC HRPPDs both by Incom and EIC research groups at Jefferson Lab and Brookhaven Lab, with peak values at ∼365 nm typically exceeding 30%. Furthermore, in this study, we present a first direct measurement of Photon Detection Efficiency being equal to (17.1 ± 0.1 [stat] ± 0.3 [sys])% at 398.6 nm, for a pixel near the center of HRPPD, in a photoelectron pulse counting mode using a picosecond diode laser. HRPPD QE at the same spot and at the same wavelength was evaluated to be (24.4 ± 0.1 [stat] ± 0.3 [sys])%, leading to a Collection Efficiency estimate of (70.3 ± 1.6)%, which is consistent with the expectations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Valued peaks: Sustainable water allocation for small hydropower plants in an era of explicit ecological needs

Optimizing hydropower operations to balance economic profitability and support functioning ecosystem services is integral to river management policy. In this article, we propose a dynamic, constrained optimization framework for small hydropower plants (SHPs) to evaluate trade-offs between economic profitability and socio-ecological requirements. Specifically, we examine the balance between short-term losses in hydropower generation and the potential for compensatory benefits in the form of revenue from recreational ecosystem services, irrespective of the direct beneficiary. Our framework integrates a fish habitat model, a hydropower optimization model, and a recreational ecosystem service estimate to evaluate different environmental flow scenarios. The optimization process gives three outflow release scenarios, informed by previous streamflow realisations (dam inflow), and designed environmental flow constraints. The framework is applied and tested for the river Kuusinkijoki in North-eastern Finland, which is a habitat for migratory brown trout and grayling populations. We show that the revenue loss due to the environmental flow constraints arises through a reduction in revenue per generated energy unit and through a reduction in turbine efficiency. Additionally, the simulation results reveal that all the designed environmental flow constraints cannot be met simultaneously. Under the environmental flow scenario with both minimum flow and flow ramping rate constraints, the annual hydropower revenue decreases by 16.5 %. An annual increase of 8 % in recreational fishing visits offsets the revenue loss. In conclusion, the developed framework provides knowledge of the costs and benefits of hydropower environmental flow constraints and guides the prioritizing process of environmental measures.

13 HYDRO ENERGY↗
DOE OSTI2023

Ly α emission in low-redshift most metal-deficient compact star-forming galaxies

ABSTRACT We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of nine most metal-deficient compact star-forming galaxies with oxygen abundances 12 + log(O/H) = 6.97–7.23, redshifts z = 0.02811–0.13320, and stellar masses M⋆ ≤ 107 M⊙. We aim to study the properties of Ly α emission in these extremely metal-deficient objects. We find that all nine galaxies are Ly α emitters (LAEs). We examine various relations between the Ly α escape fraction fesc(Ly α) and other characteristics – such as absolute UV magnitude, oxygen abundance, O32 ratio, stellar mass, Lyman-alpha luminosity, and equivalent width EW(Ly α), ionizing photon production efficiency ξion and velocity separation Vsep between the two peaks of the Ly α profile – of a large sample of LAEs, including our lowest metallicity galaxies and other objects from the literature. We find a relatively tight correlation between fesc(Ly α) and two characteristics, EW(Ly α) and Vsep, whereas no correlation is found between fesc(Ly α) and the oxygen abundance. We also find a relatively tight relation between the Ly α and LyC escape fractions. We propose to use the latter relation to estimate indirectly the escaping ionizing radiation in LAEs, when direct measurements of LyC emission are not possible. We show that the global properties of low-z LAEs are very similar to those of z > 6 galaxies. They are thus ideal local proxies for studying physical processes during the epoch of reionization of the Universe.

Izotov, Y. I. (ORCID:0000000214166082)↗
NASA NTRS1992

Estimation of the nucleation rate by differential scanning calorimetry

A realistic computer model is presented for calculating the time-dependent volume fraction transformed during the devitrification of glasses, assuming the classical theory of nucleation and continuous growth. Time- and cluster-dependent nucleation rates are calculated by modeling directly the evolving cluster distribution. Statistical overlap in the volume fraction transformed is taken into account using the standard Johnson-Mehl-Avrami formalism. Devitrification behavior under isothermal and nonisothermal conditions is described. The model is used to demonstrate that the recent suggestion by Ray and Day (1990) that nonisothermal DSC studies can be used to determine the temperature for the peak nucleation rate, is qualitatively correct for lithium disilicate, the glass investigated.

Kelton, Kenneth F.↗
NASA NTRS1994

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗