Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter compression”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

TE/TM Simulations of Interferometric Measurements

Interferometric synthetic aperture radar (IFSAR) measurements at X-, C-, L-, and P-band are used to derive ground topography at meter level resolution. Interpretation of the derived topography requires attention due to the complex interaction of the radar signal with ground cover. The presence of penetrable surfaces such as vegetation, and tree canopies poses a challenge since the depth of penetration depends on a number of parameters such as the operating radar frequency, polarization, incident angle, as well as terrain structure. The dependence of the reconstructed topography on polarization may lead to the characterization of the ground cover. Simulation of interferometric measurements is useful for interpretation of the derived topography (B. Houshmand, Proceedings of URSI, 314, 1997). In this talk , time domain simulations for interferometric measurement for TE- and TM- polarization are presented. Time domain simulation includes the effects of the surface material property as well geometry comparable the radar signal wavelength (B. Houshmand, Proceedings of the URSI, 25, 1998). The IFSAR simulation is carried out in two steps. First, the forward scattering data is generated based on full wave analysis. Next, the electromagnetic information is inverted to generate surface topography. This inversion is based on the well known IFSAR processing technique which is composed of signal compression, and formation of an interferogram. The full wave forward scattering data is generated by the scattered-field formulation of the FDTD algorithm. The simulation is carried out by exciting the computational domain by a radar signal. The scattered field is then computed and translated to the receiving interferometric antennas using the time-domain Huygen's principle. The inversion process starts by compressing the time-domain data. The range compressed data from both receivers are then coregistered to form an interferogram. The resulting interferogram is then related to the ground topography using the radar imaging geometry. In this talk, the simulation results are compared with the C-band TM IFSAR derived topography, and the TE/TM SAR images at L-Band.

Houshmand, Bijan↗
DOE OSTI2019

Hot deformation behavior and flow stress modeling of a Ni-based superalloy

The deformation behavior of a novel Ni-based superalloy was investigated using isothermal compression on a Gleeble system at temperatures between 1050 and 1210 °C with strain rates between 0.001 and 0.1 s –1 . Flow-stress curves and electron backscatter diffraction maps were employed to experimentally identify the various flow mechanisms operative during deformation. Deformation at temperatures below 1130 °C presented strong work hardening with limited restoration during dynamic softening leading to partially recrystallized microstructures. Increasing the deformation temperature to and above 1130°C enhanced the driving force for dislocation and grain boundary mobility thereby enabling dynamic recovery (DRV) and dynamic recrystallization (DRX) mechanisms to better operate. Here, the influence of the strain rate was more evident during deformation at these temperatures. Increasing the strain rate from 0.001 s –1 to 0.1 s –1 resulted in a transition in dominant softening mechanism from DRV to DRX. Flow stress modeling using the Zener-Hollomon parameter was performed to obtain the activation energy and the constitutive equation for hot deformation of the alloy. Strong changes in flow behavior affected the accuracy of the flow stress model, and thus, the model was used alternatively to identify deformation parameters associated with various flow regimes. In doing so, the activation energy and the other equation constants were obtained for each deformation mechanism observed experimentally.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Toward first principles-based simulations of dense hydrogen

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of vapor-liquid interactions in condensing ejectors

Ejectors are compact mechanical devices that utilize the expansion of a high-pressure primary fluid to entrain and compress a low-pressure secondary fluid by means of momentum transfer between the two streams of fluid. Condensing ejectors feature both momentum and heat transfer through the interaction between the vapor and liquid streams. The goal of this paper is to develop and validate a one-dimensional slug flow model to simulate the vapor–liquid interactions in two types of condensing ejectors – one with primary liquid and secondary vapor flows (Type I), and the other with primary vapor and secondary liquid flows (Type II). Control volume analysis of the mass, momentum, and energy balance in each phase and across the interface was conducted for the ejectors. The slug flow models for both types of ejectors are validated with published results in the literature. The friction coefficient on the inner wall of the ejectors and the interfacial heat transfer coefficient are identified as controlling parameters for the ejector performance in terms of pressure and temperature distributions along the axis. The detailed parametric study shows that the liquid inlet velocity and mixing tube diameter have a significant impact on the performance of the Type I ejectors, and the performance of the Type II ejectors is mainly controlled by the secondary liquid flow inlet temperature.

42 ENGINEERING↗
DOE OSTI2020

Different bonding type along each crystallographic axis: Computational study of poly($\textit{p}$-phenylene terephthalamide)

Poly($\textit{p}$-phenylene terephthalamide) (PPTA) exhibits van der Waals (vdW) bonding along its $\textit{a}$ axis, hydrogen bonding along its $\textit{b}$ axis, and covalent bonding along its $\textit{c}$ axis. In this work, we explore the structural and mechanical properties of PPTA using density functional theory with various functionals including LDA, PBE, PBE+rVV10L, SCAN, and SCAN+rVV10, compared with available experiments. The hierarchy of nonempirical semilocal functionals LDA, PBE, and SCAN (not fitted to any multicenter bonded system) includes differing amounts of intermediate-range vdW interaction. rVV10 is the long-range vdW correction. (rVV10L differs from rVV10 only in the value of a range parameter.) Among the tested functionals, SCAN shows the best performance for the lattice parameters of PPTA along the two crystal directions involving vdW or hydrogen-bond interaction. The equilibrium lattice constants obtained by SCAN and PBE+rVV10L are closest to experimental data, while SCAN+rVV10 slightly overbinds the system. We study the mechanical response of PPTA by applying strain along three lattice directions. Due to the inclusion of vdW interaction, SCAN, PBE+rVV10L, and SCAN+rVV10 all exhibit correct bonding strain-energy curves. On the contrary, PBE strongly underestimates the vdW interaction needed to resist uniaxial stretching along the $\textit{a}$ axis. The Young's modulus and yield strength of PPTA are computed and compared with previous results. The experimental values are much smaller than the computed ones, mainly due to the fact that the PPTA fiber samples used for measurements are mechanically weaker than the perfect molecular crystals considered in the simulations. Interestingly, when a compressive uniaxial stress of 25 Kbar is applied along the $\textit{b}$ axis, a structural phase transition, in which the hydrogen bonds reform along one diagonal of the $\textit{ab}$ rectangle, is predicted by SCAN.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗
DOE OSTI2021

Pressure Exchanger for Energy Recovery in a Trans-Critical CO2 Refrigeration System

Trans-critical CO2 vapor compression (VC) refrigeration cycles require a high compression ratio, which is associated with high expansion losses. To recover these expansion losses, a pressure exchange process between the low- and high-pressure sides of the VC cycle is proposed and examined in this study. The proposed pressure exchange system is an open type constant volume process where the high- and low-pressure flows mix inside the system. This prototype is inspired by the pressure exchangers used in reverse-osmosis (RO) desalination systems. In this system, a 2D model was generated and modeled using the computational fluid dynamics (CFD) technique. The numerical model ignored any losses due to leakage or hydraulic friction and the process is considered adiabatic. For the modeling, it was assumed that the inlet conditions for the two pressure exchanger flows are similar to the flow conditions at the evaporator and gas cooler outlets in a VC cycle. Two parameters are examined to test the validity of the system and understand their effect on the performance, including the inlet flow rate represented by the inlet velocity and the process time represented by the speed of rotation. A total of nine cases were simulated and analyzed in this study.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

A Hydro-based MCMC Analysis of SNR 0509–67.5: Revealing the Explosion Properties from Fluid Discontinuities Alone

Using Hubble Space Telescope/ACS Hα images of SNR 0509–67.5 taken ~10 yr apart, we measure the forward shock (FS) proper motions (PMs) at 231 rim locations. The average shock radius and velocity are 3.66 ± 0.036 pc and 6315 ± 310 km s –1 . Hydrodynamic simulations, recast as similarity solutions, provide models for the supernova remnant’s expansion into a uniform ambient medium. These are coupled to a Markov chain Monte Carlo (MCMC) analysis to determine explosion parameters, constrained by the FS measurements. For our baseline model, the MCMC posteriors yield an age of 315.5 ± 1.8 yr, a dynamical explosion center at 5 h 09 m 31. s 16, $-67^\circ 31^{\prime} 17\buildrel{\prime\prime}\over{.} 1$ and ambient medium densities at each azimuth ranging over 3.7–8.0 × 10 –25 g cm –3 . The age uncertainty can be an order of magnitude larger when considering other models, or subsets of the data. We detect stellar PMs corresponding to speeds in the Large Magellanic Cloud ≥ 770 km s –1 . Five stars in the remnant show measurable PMs but none are moving radially from the dynamical center. There are four stars 1."4 from the center, including three faint, previously unidentified ones. Using coronal [Fe XIV ] λ5303 emission as a proxy for the reverse shock location, we constrain the explosion energy (for a compression factor of 4) to a value of E = (1.30 ± 0.41) × 10 51 erg for the first time from shock kinematics alone. Higher compression factors (7 or more) are strongly disfavored based on multiple criteria, arguing for inefficient particle acceleration in the Balmer shocks of SNR 0509–67.5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The Impact of Alfvénic Shear Flow on Magnetic Reconnection and Turbulence

Magnetic reconnection is a fundamental and omnipresent energy conversion process in plasma physics. Novel observations of fields and particles from Parker Solar Probe (PSP) have shown the absence of reconnection in a large number of current sheets in the near-Sun solar wind. Using near-Sun observations from PSP encounters 4–11 (2020 January–2022 March), we investigate whether reconnection onset might be suppressed by velocity shear. We compare estimates of the tearing mode growth rate in the presence of shear flow for time periods identified as containing reconnecting current sheets versus nonreconnecting times, finding systematically larger growth rates for reconnection periods. Upon examination of the parameters associated with reconnection onset, we find that 85% of the reconnection events are embedded in slow, non-Alfvénic wind streams. We compare with fast, slow non-Alfvénic, and slow Alfvénic streams, finding that the growth rate is suppressed in highly Alfvénic fast and slow wind, and reconnection is not seen in these wind types, as would be expected from our theoretical expressions. These wind streams have strong Alfvénic flow shear, consistent with the idea of reconnection suppression by such flows. This could help explain the frequent absence of reconnection events in the highly Alfvénic, near-Sun solar wind observed by PSP. Finally, we find a steepening of both the trace and magnitude magnetic field spectra within reconnection periods in comparison to ambient wind. We tie this to the dynamics of relatively balanced turbulence within these reconnection periods and the potential generation of compressible fluctuations.

slow solar wind↗
NASA NTRS1994

Variability of measured sonic boom signatures. Volume 1: Technical report

Sonic boom signatures from two databases, the BOOMFILE and the XB-70, were analyzed in terms of C-weighted sound exposure level (CSEL), A-weighted sound exposure level (ASEL), and Stevens Mark VII perceived level (PLdB), as well as the more traditional peak positive overpressure and rise time. The variability of these parameters due to propagation through atmosphere was analyzed for different aircraft Mach number and altitude groups. The low Mach number/low altitude group had significantly greater variation in rise time, overpressure, and loudness level than the high Mach number/high altitude group. The loudness of measured booms were found to have a variation of up to 25 dB relative to the loudness of boom predicted for a non-turbulent atmosphere. This is due primarily to the steeper ray paths of the high Mach number/high altitude group and the corresponding shorter distances traveled by these rays through the lower atmosphere resulting in reduced refraction effects. The general trend of decreased overpressure and loudness level with increasing lateral distance was also seen. Sonic boom signatures from early morning flights had less variation in rise time and overpressure than afternoon flights because of reduced turbulence. Measures of asymmetry (difference between compression and expansion portion of the signature) showed that the variability in Delta loudness level was greater than the variability in Delta overpressure due to the large influence of turbulence on rise time. Lastly, analysis of data within 50 percent of lateral cutoff showed that the mean value for overpressure and loudness level was independent of time of day but that the frequency with which it occurred was greater in the morning. This is a clear indicator of increased turbulence in the afternoon.

Elmer, K. R.↗
NASA NTRS1994

Variability of measured sonic boom signatures. Volume 2: Data report

Sonic boom signatures from two databases, the BOOMFILE and the XB-70, were analyzed in terms of C-weighted sound exposure level (CSEL), A-weighted sound exposure level (ASEL), and Stevens Mark VII perceived level (PLdB), as well as the more traditional peak positive overpressure and rise time. The variability of these parameters due to propagation through atmosphere was analyzed for different aircraft Mach number and altitude groups. The low Mach number/low altitude group had significantly greater variation in rise time, overpressure, and loudness level than the high Mach number/high altitude group. The loudness of measured booms were found to have a variation of up to 25 dB relative to the loudness of boom predicted for a non-turbulent atmosphere. This is due primarily to the steeper ray paths of the high Mach number/high altitude group and the corresponding shorter distances traveled by these rays through the lower atmosphere resulting in reduced refraction effects. The general trend of decreased overpressure and loudness level with increasing lateral distance was also seen. Sonic boom signatures from early morning flights had less variation in rise time and overpressure than afternoon flights because of reduced turbulence. Measures of asymmetry (difference between compression and expansion portion of the signature) showed that the variability in Delta loudness level was greater than the variability in Delta overpressure due to the large influence of turbulence on rise time. Lastly, analysis of data within 50 percent of lateral cutoff showed that the mean value for overpressure and loudness level was independent of time of day but that the frequency with which it occurred was greater in the morning. This is a clear indicator of increased turbulence in the afternoon.

Elmer, K. R.↗
NASA NTRS1994

The Uranian aurora and its relationship to the magnetosphere

About 32 h of Voyager Ultraviolet Spectrometer (UVS) observations of Uranus H2 band airglow emission (875 less than or equal to lambda less than or equal to 1115 A) have been analyzed using the singular value decomposition (SVD) approach to inversion, producing an intensity map showing aurora at both magnetic poles. An H Lyman alpha aurora may also be present but is difficult to separate from scattered solar and local interstellar medium components. SVD analysis of variance shows that the intensity estimate is significantly larger than the error estimate over both Uranographic poles and part of the equatorial region, fortuitously including both magnetic polar regions. The Goddard Space Flight Center Q(sub 3) magnetic field model correctly predicts that the aurora should be larger in area and emit more power at the weaker N magnetic pole than at the stronger S magnetic pole. However, the auroral emissions are quite localized in magnetic longitude and so do not form complete auroral ovals. The brightest auroral emission at each magnetic pole is confined to a range of approximately 90 deg of magnetic longitude centered on the magnetotail direction, at moderate magnetic L parameter (5 less than or equal to L less than or equal to 10), but some emission at each pole is distributed over a range of more than 180 deg of longitude. The magnetic longitudes of the aurora are completely inconsistent with the 'windshield wiper' effect for either ions or electrons, indicating that some other effect, such as rapid depletion of the population of precipitating particles of highly localized strong pitch-angle diffusion, may be acting to localize emission. The low apparent L of the precipitating particles indicates that their energies may be less than or equal to 10 keV. Hence magnetospheric convection is likely to be important, and thus particles exciting the aurora may not remain on constant L shells. The precipitating particles may be a relatively low-energy population at high L that is heated to aurora-exciting energy by adiabatic compression during convection to low L. We estimate that the total auroral power output at H Lyman alpha and shorter wavelengths is about 3 x 10(exp 9) to 7 x 10(exp 9) W, requiring about 10 times that much power for excitation.

Herbert, Floyd↗
DOE OSTI2018

Microstructural packing- and rheology-based binder selection and characterization for Ultra-high Performance Concrete (UHPC)

This paper reports strategies to design the binder phase of ultra-high performance concretes (UHPC) from commonly available cement replacement (fly ash, slag, microsilica, metakaolin) and fine filler (limestone) materials. A packing algorithm is used to extract the number density, mean centroidal distance, and coordination number of the microstructure. Similarly, rheological studies on the pastes provide yield stress, plastic viscosity, and min-slump spread. The selection criteria involves using the three microstructural and three rheological parameters individually or in combination to define packing and flow coefficients. The selection criteria is flexible enough to allow users modify the constraints depending on the application. In this study, only ternary and quaternary blends, with a total cement replacement of 30% by mass are selected for further characterization. The highly efficient microstructural packing in these mixtures and better workability that facilitated dispersion of particles to enhance the reactivity results in beneficial pore structure and mechanical properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Open-Loop Control of Adjustable Tuned Mass Dampers for Floating Wind Turbine Platforms

In this study, we demonstrate the ability of adjustable tuned mass dampers (TMDs) to reduce the platform motion of floating offshore wind turbines (FOWTs). The TMDs are located in the hull and provide control authority by varying the amount of water ballast and compressed air in a reservoir over time. The optimal TMD settings depend on the wind and wave conditions of the FOWT. We present an open-loop control scheme that changes the TMD natural frequency based on sea state and a simple method for estimating the peak wave period and significant wave height. The performance of this open-loop control is compared to ideal control of the TMDs (set knowing the exact sea state), a TMD with a constant natural frequency targeting the design load case with the greatest platform motion, and a baseline set of simulations of the platform without TMDs. The constant natural frequency case performs nearly as effectively as the actively controlled case, suggesting that the TMDs can be simply designed for extreme load conditions and provide consistent control over a range of environments; however, the methods used to find the optimal TMD parameters and estimate sea state statistics could be used in other aspects of FOWT control and design.

artificial intelligence↗
DOE OSTI2018

Advanced 3D characterisation of iodine induced stress corrosion cracks in zirconium alloys

Highlights: • Macroscopic crack direction is affected by metallurgical state. • High resolution XCT reveals complex crack morphology. • 3D-EBSD used to fully characterise attacked grain boundaries. • Evidence of iodine segregation ahead of a crack tip. - Abstract: In this study, a range of advanced techniques have been used to characterise iodine-induced stress corrosion cracks, and investigate their interaction with microstructure in zirconium alloys used as cladding material in nuclear power reactors. C-ring samples, machined from cladding tubes in a recrystallised and cold worked condition, were subjected to compressive stress while submerged in an iodine-ethanol solution to simulate stress corrosion cracks arising from pellet cladding interaction. The complex morphology of these stress corrosion cracks was imaged in 3D using X-ray computed tomography (XCT), and subsequent analysis using both 2D and 3D electron backscatter diffraction (EBSD) allowed for the interaction of these cracks with microstructure to be investigated. Intergranular cracking modes were observed to dominate in the recrystallised microstructure, resulting in cracking in the radial direction with variations on the order of the grain size. Crack propagation in the cold worked microstructure was observed to be correlated with the direction of the elongated grains, resulting in a more tangential crack path with less local variation and a microstructure potentially more susceptible to iodine-induced stress corrosion cracking. 3D-EBSD using serial sectioning allows for the determination of the full five-parameter grain boundary description, and opens up the possibility of a statistical analysis of susceptible boundaries to SCC. Nanoscale secondary ion mass spectrometry (NanoSIMS) was also used to map the chemical distribution in the crack tip region and provides possible evidence of iodine segregation ahead of the crack tip. The results provide new insight into crack propagation mechanisms in recrystallised and cold worked microstructures, and pave the way for a more mechanistic understanding of crack propagation through advanced characterisation of iodine-induced stress corrosion cracks in zirconium alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Scale Separation Effects on Simulations of Plasma Turbulence

Abstract Understanding plasma turbulence requires a synthesis of experiments, observations, theory, and simulations. In the case of kinetic plasmas such as the solar wind, the lack of collisions renders the fluid closures such as viscosity meaningless and one needs to resort to higher-order fluid models or kinetic models. Typically, the computational expense in such models is managed by simulating artificial values of certain parameters such as the ratio of the Alfvén speed to the speed of light ( v A / c ) or the relative mass ratio of ions and electrons ( m i / m e ). Although, typically care is taken to use values as close as possible to realistic values within the computational constraints, these artificial values could potentially introduce unphysical effects. These unphysical effects could be significant at sub-ion scales, where kinetic effects are the most important. In this paper, we use the 10-moment fluid model in the Gkeyll framework to perform controlled numerical experiments, systematically varying the ion–electron mass ratio from a small value down to the realistic proton–electron mass ratio. We show that the unphysical mass ratio has a significant effect on the kinetic range dynamics as well as the heating of both plasma species. The dissipative process for both ions and electrons becomes more compressive in nature, although the ions remain nearly incompressible in all cases. The electrons move from being dominated by incompressive viscous-like heating/dissipation to very compressive heating/dissipation dominated by compressions/rarefactions. While the heating change is significant for the electrons, a mass ratio of m i / m e ∼ 250 captures the asymptotic behavior of electron heating.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

the role of shock waves in modulation of galactic cosmic rays

The understanding of modulation of the galactic cosmic rays has considerably progressed by the exploration by space probes of major heliospheric structures, such as the Corotating Interaction Regions, the neutral sheet, and the compression regions of intense heliospheric magnetic fields. Also relevant in this context were the detections in the outer heliosphere of long lasting Forbush type decreases of cosmic ray intensity. The results of recent theoretical studies on the changes in intensity and energy, at different location from the Sun, induced by the passage of shocks across the heliosphere are presented. In this version of the research, the simplest cases of modulation of uGV and 2GV particles by single or several shocks during periods of positive and negative solar field polarity are reviewed. The results of the theoretical aspects of the search is reported. The comparison of the theoretical predictions with space probe data allows conclusions to be drawn on the role of shocks on the modulation on both the 11 and 22 year galactic cosmic ray cycles in the outer heliosphere and on the plausibility of the models and parameters used.

Gall, R.↗