Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relaxation”

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,153 records · Page 64

Fragmentation Dynamics of Fluorene Explored Using Ultrafast XUV-Vis Pump-Probe Spectroscopy

We report on the use of extreme ultraviolet (XUV, 30.3 nm) radiation from the Free-electron LASer in Hamburg (FLASH) and visible (Vis, 405 nm) photons from an optical laser to investigate the relaxation and fragmentation dynamics of fluorene ions. The ultrashort laser pulses allow to resolve the molecular processes occurring on the femtosecond timescales. Fluorene is a prototypical small polycyclic aromatic hydrocarbon (PAH). Through their infrared emission signature, PAHs have been shown to be ubiquitous in the universe, and they are assumed to play an important role in the chemistry of the interstellar medium. Our experiments track the ionization and dissociative ionization products of fluorene through time-of-flight mass spectrometry and velocity-map imaging. Multiple processes involved in the formation of each of the fragment ions are disentangled through analysis of the ion images. The relaxation lifetimes of the excited fluorene monocation and dication obtained through the fragment formation channels are reported to be in the range of a few tens of femtoseconds to a few picoseconds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The rapid change in mental health among college students after introduction of on-campus quarantine during the 2022 Shanghai COVID-19 lockdown

Objectives Among the various impacts of disasters in terms of emotions, quarantine has been proven to result in significant increases in mental health problems. Studies of psychological resilience during outbreaks of epidemics tend to focus on long-term social quarantine. In contrast, insufficient studies have been conducted examining how rapidly negative mental health outcomes occur and how these outcomes change over time. We evaluated the time course of psychological resilience (over three different phases of quarantine) among students at Shanghai Jiao Tong University to investigate the influence of unexpected changes on college students. Methods An online survey was conducted from 5 to 7 April 2022. A structured online questionnaire was administered using a retrospective cohort trial design. Before 9 March (Period 1), individuals engaged in their usual activities without restrictions. From 9 to 23 March (Period 2), the majority of students were asked to remain in their dormitories on campus. From 24 March to early April (Period 3), restrictions were relaxed, and students were gradually allowed to participate in essential activities on campus. We quantified dynamic changes in the severity of students’ depressive symptoms over the course of these three periods. The survey consisted of five sets of self-reported questions: demographic information, lifestyle/activity restrictions, a brief mental health history, COVID-19-related background, and the Beck Depression Inventory, second edition. Results A total of 274 college students aged 18–42 years (mean = 22.34; SE = 0.24) participated in the study (58.39% undergraduate students, 41.61% graduate students; 40.51% male, 59.49% female). The proportion of students with depressive symptoms was 9.1% in Period 1, 36.1% in Period 2, and 34.67% in Period 3. Depressive symptoms increased notably with the introduction of the quarantine in Periods 2 and 3. Lower satisfaction with the food supplied and a longer duration of physical exercise per day were found to be positively associated with changes in depression severity in Periods 2 and 3. Quarantine-related psychological distress was more evident in students who were in a romantic relationship than in students who were single. Conclusion Depressive symptoms in university students rapidly increased after 2 weeks of quarantine and no perceptible reversal was observed over time. Concerning students in a relationship, ways to take physical exercise and to relax should be provided and the food supplied should be improved when young people are quarantined.

Ma, Dongni↗

MgADP Promotes Myosin Head Movement toward Actin at Low [Ca 2+ ] to Increase Force Production and Ca 2+ -Sensitivity of Contraction in Permeabilized Porcine Myocardial Strips

Myosin cross-bridges dissociate from actin following Mg 2+ -adenosine triphosphate (MgATP) binding. Myosin hydrolyses MgATP into inorganic phosphate (P i ) and Mg 2+ -adenosine diphosphate (ADP), and release of these hydrolysis products drives chemo-mechanical energy transitions within the cross-bridge cycle to power muscle contraction. Some forms of heart disease are associated with metabolic or enzymatic dysregulation of the MgATP-MgADP nucleotide pool, resulting in elevated cytosolic [MgADP] and impaired muscle relaxation. We investigated the mechanical and structural effects of increasing [MgADP] in permeabilized myocardial strips from porcine left ventricle samples. Sarcomere length was set to 2.0 µm at 28 °C, and all solutions contained 3% dextran T-500 to compress myofilament lattice spacing to near-physiological values. Under relaxing low [Ca 2+ ] conditions (pCa 8.0, where pCa = -log 10 [Ca 2+ ]), tension increased as [MgADP] increased from 0-5 mM. Complementary small-angle X-ray diffraction measurements show that the equatorial intensity ratio, I 1,1 /I 1,0 , also increased as [MgADP] increased from 0 to 5 mM, indicating myosin head movement away from the thick-filament backbone towards the thin-filament. Ca 2+ -activated force-pCa measurements show that Ca 2+ -sensitivity of contraction increased with 5 mM MgADP, compared to 0 mM MgADP. These data show that MgADP augments tension at low [Ca 2+ ] and Ca 2+ -sensitivity of contraction, suggesting that MgADP destabilizes the quasi-helically ordered myosin OFF state, thereby shifting the cross-bridge population towards the disordered myosin ON state. Together, these results indicate that MgADP enhances the probability of cross-bridge binding to actin due to enhancement of both thick and thin filament-based activation mechanisms.

60 APPLIED LIFE SCIENCES↗

The Effect of the Flux Separability Approximation on Multigroup Neutron Transport

The angular dependence of flux-weighted multigroup cross sections is commonly neglected when generating multigroup libraries. The error of this flux separability approximation is typically not isolated from other error sources due to a lack of availability of library generation and corresponding solvers that cannot relax this approximation. These errors can now be isolated and quantified with the availability of a multigroup Monte Carlo transport and multigroup library-generation capability in the OpenMC Monte Carlo transport code. This work will discuss relevant details of the OpenMC implementation, provide an example case useful for detailing the type of errors one can expect from making the flux separability approximation, and end with more realistic problems which show the impact of the approximation and highlight how it can strongly arise from an energy-dependent resonance absorption effect. Since the angle-dependence is intrinsically linked to the energy group structure, these examples also show that relaxing the flux separability approximation with angle-dependent cross sections could be used to reduce either the fine-tuning required to set a multigroup energy structure for a specific reactor type or the number of energy groups required to obtain a desired level of accuracy for a given problem. This trade-off could increase the costs of generating multigroup cross sections, and has the potential to require more memory for storing the multigroup library during the transport calculations, but it can significantly reduce the computational time required since the runtime of a discrete ordinates or method of characteristics neutron transport solver scales roughly linearly with the number of groups.

42 ENGINEERING↗

Probing Non-Equilibrium Pair-Breaking and Quasiparticle Dynamics in Nb Superconducting Resonators Under Magnetic Fields

We conducted a comprehensive study of the non-equilibrium dynamics of Cooper pair breaking, quasiparticle (QP) generation, and relaxation in niobium (Nb) cut from superconducting radio-frequency (SRF) cavities, as well as various Nb resonator films from transmon qubits. Using ultrafast pump–probe spectroscopy, we were able to isolate the superconducting coherence and pair-breaking responses. Our results reveal both similarities and notable differences in the temperature- and magnetic-field-dependent dynamics of the SRF cavity and thin-film resonator samples. Moreover, femtosecond-resolved QP generation and relaxation under an applied magnetic field reveals a clear correlation between non-equilibrium QPs and the quality factor of resonators fabricated by using different deposition methods, such as DC sputtering and high-power impulse magnetron sputtering. These findings highlight the pivotal influence of fabrication techniques on the coherence and performance of Nb-based quantum devices, which are vital for applications in superconducting qubits and high-energy superconducting radio-frequency applications.

43 PARTICLE ACCELERATORS↗

Surface Interactions and Nanoconfinement of Methane and Methane plus CO2 Revealed by High-Pressure Magic Angle Spinning NMR Spectroscopy and Molecular Dynamics

This study explores the fundamental, molecular- to microscopic-level behavior of methane gas confined into nanoporous silica proxies with different pore diameters and surface-to-volume (S/V) ratios. Surfaces and pore walls of nanoporous silica matrices are decorated with hydroxyl (-OH) groups, resembling natural heterogeneity. High-pressure MAS NMR was utilized to characterize the interactions between methane and the engineered nanoporous silica proxies under various temperature and pressure regimes. There was a change in the chemical shift position of confined methane slightly in the mixtures with nanoporous silica up to 393 K, as shown by high-pressure 13C-NMR. The 13C-NMR chemical shift of methane was changed by pressure, explained by the densification of methane inside the nanoporous silica materials. The influence of pore diameter and S/V of the nanoporous silica materials on the behaviors and dynamics of methane were studied. The presence of CO2 in mixtures of silica and methane needs analysis with caution because CO2 in a supercritical state and gaseous CO2 change the original structure of nanoporous silica and change surface area and pore volume. According to simulation, the picosecond scale dynamics of methane confined in larger pores of amorphous silica is faster. In the 4 nm pore, the diffusivity obtained from MD simulations in the pore with a higher S/V ratio is slower due to the trapping of methane molecules in adsorbed layers close to the corrugated pore surface. In contrast, relaxation measured with NMR for smaller pores (higher S/V) exhibits larger T1, indicating slower relaxation.

03 NATURAL GAS↗

Efficient SABRE-SHEATH Hyperpolarization of Potent Branched-Chain-Amino-Acid Metabolic Probe [1-13C]ketoisocaproate

Efficient 13C hyperpolarization of ketoisocaproate is demonstrated in natural isotopic abundance and [1-13C]enriched forms via SABRE-SHEATH (Signal Amplification By Reversible Exchange in SHield Enables Alignment Transfer to Heteronuclei). Parahydrogen, as the source of nuclear spin order, and ketoisocaproate undergo simultaneous chemical exchange with an Ir-IMes-based hexacoordinate complex in CD3OD. SABRE-SHEATH enables spontaneous polarization transfer from parahydrogen-derived hydrides to the 13C nucleus of transiently bound ketoisocaproate. 13C polarization values of up to 18% are achieved at the 1-13C site in 1 min in the liquid state at 30 mM substrate concentration. The efficient polarization build-up becomes possible due to favorable relaxation dynamics. Specifically, the exponential build-up time constant (14.3 ± 0.6 s) is substantially lower than the corresponding polarization decay time constant (22.8 ± 1.2 s) at the optimum polarization transfer field (0.4 microtesla) and temperature (10 °C). The experiments with natural abundance ketoisocaproate revealed polarization level on the 13C-2 site of less than 1%—i.e., one order of magnitude lower than that of the 1-13C site—which is only partially due to more-efficient relaxation dynamics in sub-microtesla fields. We rationalize the overall much lower 13C-2 polarization efficiency in part by less favorable catalyst-binding dynamics of the C-2 site. Pilot SABRE experiments at pH 4.0 (acidified sample) versus pH 6.1 (unaltered sodium [1-13C]ketoisocaproate) reveal substantial modulation of SABRE-SHEATH processes by pH, warranting future systematic pH titration studies of ketoisocaproate, as well as other structurally similar ketocarboxylate motifs including pyruvate and alpha-ketoglutarate, with the overarching goal of maximizing 13C polarization levels in these potent molecular probes. Finally, we also report on the pilot post-mortem use of HP [1-13C]ketoisocaproate in a euthanized mouse, demonstrating that SABRE-hyperpolarized 13C contrast agents hold promise for future metabolic studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of Coarse-Grained Models for Poly(4-vinylphenol) and Poly(2-vinylpyridine): Polymer Chemistries with Hydrogen Bonding

In this paper, we identify the modifications needed in a recently developed generic coarse-grained (CG) model that captured directional interactions in polymers to specifically represent two exemplary hydrogen bonding polymer chemistries—poly(4-vinylphenol) and poly(2-vinylpyridine). We use atomistically observed monomer-level structures (e.g., bond, angle and torsion distribution) and chain structures (e.g., end-to-end distance distribution and persistence length) of poly(4-vinylphenol) and poly(2-vinylpyridine) in an explicitly represented good solvent (tetrahydrofuran) to identify the appropriate modifications in the generic CG model in implicit solvent. For both chemistries, the modified CG model is developed based on atomistic simulations of a single 24-mer chain. This modified CG model is then used to simulate longer (36-mer) and shorter (18-mer and 12-mer) chain lengths and compared against the corresponding atomistic simulation results. We find that with one to two simple modifications (e.g., incorporating intra-chain attraction, torsional constraint) to the generic CG model, we are able to reproduce atomistically observed bond, angle and torsion distributions, persistence length, and end-to-end distance distribution for chain lengths ranging from 12 to 36 monomers. We also show that this modified CG model, meant to reproduce atomistic structure, does not reproduce atomistically observed chain relaxation and hydrogen bond dynamics, as expected. Simulations with the modified CG model have significantly faster chain relaxation than atomistic simulations and slower decorrelation of formed hydrogen bonds than in atomistic simulations, with no apparent dependence on chain length.

36 MATERIALS SCIENCE↗

New Determination of Fundamental Properties of Palomar 5 Using Deep DESI Imaging Data

The legacy imaging surveys for the Dark Energy Spectroscopic Instrument project provides multiple-color photometric data, which are about 2 mag deeper than those from the SDSS. In this study, we redetermine the fundamental properties for an old halo globular cluster of Palomar 5 based on these new imaging data, including structure parameters, stellar population parameters, and luminosity and mass functions. These characteristics, together with its tidal tails, are key for dynamical studies of the cluster and constraining the mass model of the Milky Way. By fitting the King model to the radial surface density profile of Palomar 5, we derive the core radius of [[Formula]], tidal radius of [[Formula]], and concentration parameter of c = 0.78 ± 0.04. We apply a Bayesian analysis method to derive the stellar population properties and get an age of 11.508 ± 0.027 Gyr, metallicity of [Fe/H] = −1.798 ± 0.014, reddening of E(B − V) = 0.0552 ± 0.0005, and distance modulus of (m−M){sub 0} = 16.835±0.006. The main-sequence luminosity and mass functions for both the cluster center and tidal tails are investigated. The luminosity and mass functions at different distances from the cluster center suggest that there is obvious spatial mass segregation. Many faint low-mass stars have been evaporated at the cluster center, and the tidal tails are enhanced by low-mass stars. Both the concentration and relaxation times suggest that Palomar 5 is a totally relaxed system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evidence for AGN-regulated Cooling in Clusters at z ∼ 1.4: A Multiwavelength View of SPT-CL J0607-4448

Abstract We present a multiwavelength analysis of the galaxy cluster SPT-CL J0607-4448 (SPT0607), which is one of the most distant clusters discovered by the South Pole Telescope at z = 1.4010 ± 0.0028. The high-redshift cluster shows clear signs of being relaxed with well-regulated feedback from the active galactic nucleus (AGN) in the brightest cluster galaxy (BCG). Using Chandra X-ray data, we construct thermodynamic profiles and determine the properties of the intracluster medium. The cool-core nature of the cluster is supported by a centrally peaked density profile and low central entropy ( K 0 = 18 − 9 + 11 keV cm 2 ), which we estimate assuming an isothermal temperature profile due to the limited spectral information given the distance to the cluster. Using the density profile and gas cooling time inferred from the X-ray data, we find a mass-cooling rate M ̇ cool = 100 − 60 + 90 M ⊙ yr −1 . From optical spectroscopy and photometry around the [O ii ] emission line, we estimate that the BCG star formation rate is SFR [ O II ] = 1.7 − 0.6 + 1.0 M ⊙ yr −1 , roughly two orders of magnitude lower than the predicted mass-cooling rate. In addition, using ATCA radio data at 2.1 GHz, we measure a radio jet power P cav = 3.2 − 1.3 + 2.1 × 10 44 erg s −1 , which is consistent with the X-ray cooling luminosity ( L cool = 1.9 − 0.5 + 0.2 × 10 44 erg s −1 within r cool = 43 kpc). These findings suggest that SPT0607 is a relaxed, cool-core cluster with AGN-regulated cooling at an epoch shortly after cluster formation, implying that the balance between cooling and feedback can be reached quickly. We discuss the implications for these findings on the evolution of AGN feedback in galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

The Santa Barbara Binary-disk Code Comparison

We have performed numerical calculations of a binary interacting with a gas disk, using 11 different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, 0.2(a 2 Ω B /v) is the kinematic viscosity. The largest discrepancies between codes resulted from the dimensionality of the setup (3D vs. 2D disks). We find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron Influence on the Parallel Proton Firehose Instability in 10-moment, Multifluid Simulations

Instabilities driven by pressure anisotropy play a critical role in modulating the energy transfer in space and astrophysical plasmas. For the first time, we simulate the evolution and saturation of the parallel proton firehose instability using a multifluid model without adding artificial viscosity. These simulations are performed using a 10-moment, multifluid model with local and gradient relaxation heat-flux closures in high-β proton–electron plasmas. When these higher-order moments are included and pressure anisotropy is permitted to develop in all species, we find that the electrons have a significant impact on the saturation of the parallel proton firehose instability, modulating the proton pressure anisotropy as the instability saturates. Even for lower β's more relevant to heliospheric plasmas, we observe a pronounced electron energization in simulations using the gradient relaxation closure. Our results indicate that resolving the electron pressure anisotropy is important to correctly describe the behavior of multispecies plasma systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Semi-Lagrangian nodal discontinuous Galerkin method for the BGK model

In this work, we propose a semi-Lagrangian (SL) nodal discontinuous Galerkin (DG) solver for the BGK equation. The BGK model was introduced by Bhatnagar, Gross, and Krook [1] as a relaxation model for the fundamental Boltzmann equation [5], which describes the kinetic dynamic of rarefied gases with a probability distribution function. The challenges of designing efficient numerical schemes for the Boltzmann equation mainly come from its high dimensionality and complicated nonlinear collision operator. The BGK model gains interests since it has much lower computational cost, due to the relatively simple structure of the relaxation operator in replacement of the collision operator, while simultaneously preserving several important physical properties, such as macroscopic quantities and dissipation of entropy.

97 MATHEMATICS AND COMPUTING↗

Protection Against Graph-Based False Data Injection Attacks on Power Systems

Graph signal processing (GSP) has emerged as a powerful tool for practical network applications, including power system monitoring. By representing power system voltages as smooth graph signals, recent research has focused on developing GSP-based methods for state estimation, attack detection, and topology identification. Included, efficient methods have been developed for detecting false data injection (FDI) attacks, which until now were perceived as non-smooth with respect to the graph Laplacian matrix. Consequently, these methods may not be effective against smooth FDI attacks. In this paper, we propose a graph FDI (GFDI) attack that minimizes the Laplacian-based graph total variation (TV) under practical constraints. In addition, we develop a low-complexity algorithm that solves the non-convex GDFI attack optimization problem using ell_1-norm relaxation, the projected gradient descent (PGD) algorithm, and the alternating direction method of multipliers (ADMM). We then propose a protection scheme that identifies the minimal set of measurements necessary to constrain the GFDI output to high graph TV, thereby enabling its detection by existing GSP-based detectors. Our numerical simulations on the IEEE-57 bus test case reveal the potential threat posed by well-designed GSP-based FDI attacks. Moreover, we demonstrate that integrating the proposed protection design with GSP-based detection can lead to significant hardware cost savings compared to previous designs of protection methods against FDI attacks.

Morgenstern, Gal↗

Toward closure between predicted and observed particle viscosity over a wide range of temperatures and relative humidity

Abstract. Atmospheric aerosols can exist in amorphous semi-solid or glassy phase states whose viscosity varies with atmospheric temperature and relative humidity. The temperature and humidity dependence of viscosity has been hypothesized to be predictable from the combination of a water–organic binary mixing rule of the glass transition temperature, a glass-transition-temperature-scaled viscosity fragility parameterization, and a water uptake parameterization. This work presents a closure study between predicted and observed viscosity for sucrose and citric acid. Viscosity and glass transition temperature as a function of water content are compiled from literature data and used to constrain the fragility parameterization. New measurements characterizing viscosity of sub-100 nm particles using the dimer relaxation method are presented. These measurements extend the available data of temperature- and humidity-dependent viscosity to −28 ∘C. Predicted relationships agree well with observations at room temperature and with measured isopleths of constant viscosity at ∼107 Pa s at temperatures warmer than −28 ∘C. Discrepancies at colder temperatures are observed for sucrose particles. Simulations with the kinetic multi-layer model of gas–particle interactions suggest that the observed deviations at colder temperature for sucrose can be attributed to kinetic limitations associated with water uptake at the timescales of the dimer relaxation experiments. Using the available information, updated equilibrium phase-state diagrams (-80∘C

54 ENVIRONMENTAL SCIENCES↗

Dataset from: Improved multispecies Dougherty collisions

The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species primitive moments are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, avoids unphysical temperatures and satisfies the H-theorem unlike an analogous Bhatnagar-Gross-Krook operator. Formulas for both a Cartesian velocity-space and a gyroaveraged operator are provided for use in Vlasov as well as long-wavelength gyrokinetic models. We present an algorithm for the discontinuous Galerkin discretization of this operator, and provide results from relaxation and Landau damping benchmarks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Higher-order phonon scattering: advancing the quantum theory of phonon linewidth, thermal conductivity and thermal radiative properties

Phonon scattering plays a central role in the quantum theory of phonon linewidth, which in turn governs important properties including infrared spectra, Raman spectra, lattice thermal conductivity, thermal radiative properties, and also signifi- cantly affects other important processes such as hot electron relaxation. Since Maradudin and Fein’s classic work in 1962, three-phonon scattering had been considered as the dominant intrinsic phonon scattering mechanism and has seen tremendous advances. However, the role of the higher-order four-phonon scattering had been persistently unclear and so was ignored. The tremendous complexity of the formalism and computational challenges stood in the way, prohibiting the direct and quantitative treatment of four-phonon scattering. In 2016, a rigorous four-phonon scattering formalism was developed, and the prediction was realized using empirical potentials. In 2017, the method was extended using first-principles calculated force constants, and the thermal conductivities of boron arsenides (BAs), Si and diamond were predicted. The predictions for BAs were later confirmed by several independent experiments. Four-phonon scattering has since been investigated in a range of materials and established as an important intrinsic scattering mechanism for thermal transport and radiative properties. Specifically, four-phonon scattering is important when the fourth-order scattering potential or phase space becomes relatively large. The former scenario includes: (i) nearly all materials when the temperature is high; (ii) strongly anharmonic (low thermal conductivity) materials, including most rocksalt compounds, halides, hydrides, chalcogenides and oxides. The latter scenario includes: (iii) materials with large acoustic–optical phonon band gaps, such as XY compounds with a large atomic mass ratio between X and Y; (iv) two- dimensional materials with reflection symmetry, such as single-layer graphene, single-layer boron nitride and carbon nanotubes; and (v) phonons with a large density of states, such as optical phonons, which are important for Raman, infrared and thermal radiative properties. Four-phonon scattering is expected to gain broad interest in various technologically important materials for thermoelectrics, thermal barrier coatings, thermal energy storage, phase change, nuclear power, ultra-high temperature ceramics, infrared spectra, Raman spectra, radiative transport, hot electron relaxation and radiative cooling. Four-phonon scattering has been, and will continue to be, established as an important intrinsic phonon scattering mechanism beyond three-phonon scattering. The prediction of four-phonon scattering will transition from a breakthrough to a new routine in the next decade.

Feng, Tianli↗

Two-Stage Gauss-Seidel Preconditioners and Smoothers for Krylov Solvers on a GPU Cluster: Preprint

Gauss-Seidel (GS) relaxation is often employed as a preconditioner for a Krylov solver or as a smoother for Algebraic Multigrid (AMG). However, the requisite sparse triangular solve is difficult to parallelize on many-core architectures such as graphics processing units (GPUs). In the present study, the performance of the sequential GS relaxation based on a triangular solve is compared with two-stage variants, replacing the direct triangular solve with a fixed number of inner Jacobi-Richardson (JR) iterations. When a small number of inner iterations is sufficient to maintain the Krylov convergence rate, the two-stage GS (GS2) often outperforms the sequential algorithm on many-core architectures. The GS2 algorithm is also compared with JR. When they perform the same number of ops for SpMV (e.g. three JR sweeps compared to two GS sweeps with one inner JR sweep), the GS2 iterations, and the Krylov solver preconditioned with GS2, may converge faster than the JR iterations. Moreover, for some problems (e.g. elasticity), it was found that JR may diverge with a damping factor of one, whereas two-stage GS may improve the convergence with more inner iterations. Finally, to study the performance of the two-stage smoother and preconditioner for a practical problem, these were applied to incompressible uid ow simulations on GPUs.

algebraic multigrid↗