Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “isotherm model parameters”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Analysis of gas storage and transport in Eagle Ford shale using pressure pulse decay measurements with He, Kr and CO 2

Detailed characterization of gas transport and storage in low-permeability and organic-rich shales is associated with an array of challenges due to the complex morphology of the pore space, representing a broad range of pore sizes, combined with the heterogeneous fabric of the shale matrix. These factors, and their interplay during gas transport and sorption, complicate a) the analysis of shale samples at the laboratory scale (~ up to a ft. In length), and b) the prediction of natural gas production and carbon sequestration potential at larger scales. Here, in this work, tandem experiments with inert (helium – He) and adsorbing (krypton - Kr and carbon dioxide - CO 2 ) gases were performed and analyzed to develop an efficient workflow for characterizing and modeling transport and sorption at the laboratory scale. In particular, pressure pulse-decay (PPD) measurements were conducted on an Eagle Ford shale core sample at room temperature using He, Kr, and CO 2 . PPD measurements with He (a non-sorbing gas) were employed to probe the overall porosity, including natural fractures, microcracks, mesopores, and micropores. A triple-porosity model (TPM) was adopted to interpret the gas transport in the shale sample: The pore space is represented by three interacting continua, including macropores (larger fractures), mesopores (including microcracks), and micropores. To facilitate the application of the TPM, a modified analytical approach is introduced to extract effective transport parameters (in terms of the characteristic time for transport at relevant porosity levels) directly from the PPD measurements with He. To validate the modeling approach, model parameters extracted from two He PPD measurements are demonstrated to provide excellent agreement with a 3rd He PPD measurement performed at a higher pressure. The effective transport parameters, extracted from the He experiments, were subsequently converted for application to the Kr and CO 2 PPD experiments, by accounting for relevant transport modes and differences in fluid properties. Excess adsorption isotherms were extracted from the equilibrium pressures of the Kr and CO 2 PPD experiments using the He pore volume as a baseline. These adsorption isotherms were then integrated into the TPM to predict combined gas transport and sorption for Kr and CO 2 , with transport coefficients translated from the He measurements. The predictions for Kr and CO 2 are demonstrated to be in excellent agreement with the experimental observations. This, in turn, demonstrates that the proposed analytical approach provides for an effective characterization of mass transfer rates in shales, that can be applied directly in a TPM representation of mass transfer and sorption.

58 GEOSCIENCES↗
DOE OSTI2020

Metamodeling Time-Temperature Creep Parameters

Abstract There exist many time-temperature parameter (TTP) models for creep rupture prediction of components including the Larson–Miller (LM), Manson–Haferd (MH), Manson–Brown (MB), Orr–Sherby–Dorn (OSD), Manson–Succop (MS), Graham–Walles (GW), Chitty–Duval (CD), Goldhoff–Sherby (GS) models. It remains a challenge to determine which model is “best”, capable of accurate interpolation and physically realistic extrapolation of creep rupture data for a given material. In this study, metamodeling is applied to create a unified TTP metamodel that combines and regresses into twelve TTP models (eight existing and four newly derived). An analysis of the mathematical problems that exist in TTP models is provided. A matlab code is written that can: (1) calibrate the material constants of any of the twelve TTP models (using the metamodel); (2) determine the most suitable stress-parameter function; (3) and report the normalized mean square error (NMSE) of rupture predictions for a given material database. Using the metamodel, and code, a design engineer can make an intelligent selection of the “best” TTP model for creep resistant design. This process is demonstrated using four isotherms of alloy P91 creep rupture data. To assess the influence of material, further validation is performed on alloys Hastelloy X, 304SS, and 316SS. It is determined that the “best” model is dependent on material type and the quality and quantity of available data.

Engineering↗
NASA NTRS2019

An Ensemble of Bayesian Neural Networks for Exoplanetary Atmospheric Retrieval

Machine learning (ML) is now used in many areas of astrophysics, from detecting exoplanets in Kepler transit signals to removing telescope systematics. Recent work demonstrated the potential of using ML algorithms for atmospheric retrieval by implementing a random forest (RF) to perform retrievals in seconds that are consistent with the traditional, computationally expensive nested-sampling retrieval method. We expand upon their approach by presenting a new ML model, plan-net, based on an ensemble of Bayesian neural networks (BNNs) that yields more accurate inferences than the RF for the same data set of synthetic transmission spectra. We demonstrate that an ensemble provides greater accuracy and more robust uncertainties than a single model. In addition to being the first to use BNNs for atmospheric retrieval, we also introduce a new loss function for BNNs that learns correlations between the model outputs. Importantly, we show that designing ML models to explicitly incorporate domain-specific knowledge both improves performance and provides additional insight by inferring the covariance of the retrieved atmospheric parameters. We apply plan-net to the Hubble Space Telescope Wide Field Camera 3 transmission spectrum for WASP-12b and retrieve an isothermal temperature and water abundance consistent with the literature. We highlight that our method is flexible and can be expanded to higher resolution spectra and a larger number of atmospheric parameters.

Adam D. Cobb↗
DOE OSTI2022

Uranium Retardation Capacity of Lithologies from the Negev Desert, Israel—Rock Characterization and Sorption Experiments

A series of batch experiments were performed to assess the uranium sorption capacity of four mineralogically distinct lithologies from the Negev Desert, Israel, to evaluate the suitability of a potential site for subsurface radioactive waste disposal. The rock specimens consisted of an organic-rich phosphorite, a bituminous marl, a chalk, and a sandstone. The sorption data for each lithology were fitted using a general composite surface complexation model (GC SCM) implemented in PHREEQC. Sorption data were also fitted by a non-mechanistic Langmuir sorption isotherm, which can be used as an alternative to the GC SCM to provide a more computationally efficient method for uranium sorption. This is because all the rocks tested have high pH/alkalinity/calcium buffering capacities that restrict groundwater chemistry variations, so that the use of a GC SCM is not advantageous. The mineralogy of the rocks points to several dominant sorption phases for uranyl (UO 2 2+ ), including apatite, organic carbon, clays, and iron-bearing phases. The surface complexation parameters based on literature values for the minerals identified overestimate the uranium sorption capacities, so that for our application, an empirical approach that makes direct use of the experimental data to estimate mineral-specific sorption parameters appears to be more practical for predicting uranium sorption.

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

Large Destabilization of (TiVNb)-Based Hydrides via (Al, Mo) Addition: Insights from Experiments and Data-Driven Models

High-entropy alloys (HEAs) represent an interesting alloying strategy that can yield exceptional performance properties needed across a variety of technology applications, including hydrogen storage. Examples include ultrahigh volumetric capacity materials (BCC alloys → FCC dihydrides) with improved thermodynamics relative to conventional high-capacity metal hydrides (like MgH 2 ), but still further destabilization is needed to reduce operating temperature and increase system-level capacity. Here, in this work, we demonstrate efficient hydride destabilization strategies by synthesizing two new Al 0.05 (TiVNb) 0.95–x Mo x (x = 0.05, 0.10) compositions. We specifically evaluate the effect of molybdenum (Mo) addition on the phase structure, microstructure, hydrogen absorption, and desorption properties. Both alloys crystallize in a bcc structure with decreasing lattice parameters as the Mo content increases. The alloys can rapidly absorb hydrogen at 25 °C with capacities of 1.78 H/M (2.79 wt %) and 1.79 H/M (2.75 wt %) with increasing Mo content. Pressure-composition isotherms suggest a two-step reaction for hydrogen absorption to a final fcc dihydride phase. The experiments demonstrate that increasing Mo content results in a significant hydride destabilization, which is consistent with predictions from a gradient boosting tree data-driven model for metal hydride thermodynamics. Furthermore, improved desorption properties with increasing Mo content and reversibility were observed by in situ synchrotron X-ray diffraction, in situ neutron diffraction, and thermal desorption spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Intrinsic kinetics of water-inhibited ultra-lean methane oxidation over PtPd-Mg/ θ -Al 2 O 3 catalyst

Here, this study develops intrinsic methane oxidation kinetics for ultra-lean methane conditions in the presence of water over a highly active and stable PtPd–Mg/θ-Al 2 O 3 catalyst. Comprehensive laboratory experiments were conducted over a wide range of methane concentrations (150–1200 ppm CH 4 ), water contents (1–5% H 2 O), and industrially relevant space velocities (80,000 ≤ GHSV ≤ 160,000 h -1 ). These systematic experiments informed a two-dimensional, axisymmetric, multiscale reactor model that was used to develop and validate methane oxidation kinetics under practically relevant conditions, including non-isothermal operation and high conversion regimes. Combined experimental and modeling results revealed significant intraparticle diffusion resistance and transport-induced reaction exotherm at elevated temperatures, which limited catalyst utilization despite high intrinsic activity. These transport effects were explicitly incorporated into the reactor model, enabling accurate estimation of intrinsic kinetic parameters without reliance on conventional effectiveness-factor corrections. The resulting kinetic model successfully captured both kinetically controlled and mass-transfer-limited regimes and reliably predicted CH 4 conversion across broad ranges of temperature, methane concentration, and water content

Heat and mass transfer limitations↗
NASA NTRS2007

Scaling Relations from Sunyaev-Zel'dovich Effect and Chandra X-ray Measurements of High-Redshift Galaxy Clusters

We present Sunyaev-Zel'dovich Effect (SZE) scaling relations for 38 massive galaxy clusters at redshifts 0.14 less than or equal to z less than or equal to 0.89, observed with both the Chandra X-ray Observatory and the centimeter-wave SZE imaging system at the BIMA and OVRO interferometric arrays. An isothermal ,Beta-model with central 100 kpc excluded from the X-ray data is used to model the intracluster medium and to measure global cluster properties. For each Cluster, we measure the X-ray spectroscopic temperature, SZE gas mass, total mass. and integrated Compton-gamma parameters within r(sub 2500). Our measurements are in agreement with the expectations based on a simple self-similar model of cluster formation and evolution. We compare the cluster properties derived from our SZE observations with and without Chandra spatial and spectral information and find them to be in good agreement: We compare our results with cosmological numerical simulations, and find that simulations that include radiative cooling, star formation and feedback match well both the slope and normalization of our SZE scaling relations.

Bonamente, Massimiliano↗
DOE OSTI2025

The Role of Water Volume Fraction on Water Adsorption in Anion Exchange Membranes

Water absorption plays a key role in the performance of polymeric anion exchange membranes. It influences important properties such as ionic conductivity and mechanical strength and alters their performance as solid electrolytes in hydrogen electrochemical devices for energy conversion. However, computational approaches that address the relationship between the polymer design and the absorption process are scarce. In this work, we introduce a simple thermodynamic model to predict the water absorption isotherms of polyelectrolyte membranes in contact with a water vapor reservoir that incorporates the specific chemical design of the polymers. The model accurately predicts the water content and macrostructural properties of polynorbornene membranes as a function of the water activity and successfully captures the effect of various polymer design parameters. The energy of pairwise attractive interactions predicted by our model provides a means to interpret the absorption process at the molecular level. The model also reveals the most significant favorable and unfavorable contributions to the free energy and indicates that their balance is solely governed by the water volume fraction, regardless of the polymer design. This universal behavior leads to important implications in the search for better ion exchange membranes.

Zaldivar, Gervasio [New York University, Brooklyn,↗
DOE OSTI2026

Cryogenic-Refined MOSFET Modeling for Oscillator, Frequency Divider, and Amplifier Designs Below 4 K

Capturing device characteristic changes at cryogenic temperatures is crucial for cryo-CMOS circuit designs. In this work, we present an isothermal cryogenic-refined modeling approach for CMOS transistors that is simple, low overhead, and easy to implement while offering the required accuracy for predicting circuit performance at the designated temperatures. Guided by die-level measurement data and circuit design principles, the model introduces corrections to only five critical parameters: threshold voltage, carrier mobility, elevated low-frequency flicker noise, dominant high-frequency shot noise, and subthreshold swing (SS). These refinements are implemented around the foundry-provided SPICE model, which is typically validated only down to about 200 K. With these adjustments, the proposed cryogenic-refined model achieves less than 5% error in both large-signal metrics (I–V characteristics) and small-signal parameters (e.g., transconductance) when compared with device measurements at deep-cryogenic temperatures. The methodology is validated in two advanced technologies: TSMC 40-nm CMOS and GlobalFoundries (GF) 45-nm RF-SOI. We further demonstrate its applicability in three representative RF circuits: a 30-GHz LC oscillator, a high-speed current-mode-logic (CML) frequency divider (FD), and a subthreshold Gb/s amplifier, all showing close agreement between simulated predictions and measurements performed at 4 and 2.5 K. Finally, we believe that the proposed approach is implementation-friendly and can significantly accelerate the development of cryo-CMOS integrated circuits.

circuit modeling↗
NASA NTRS2008

Scaling Relations from Sunyaev-Zel'dovich Effect and Chandra X-Ray Measurements of High-Redshift Galaxy Clusters

We present Sunyaev-Zel'dovich Effect (SZE) scaling relations for 38 massive galaxy clusters at redshifts 0.14 less than or equal to z less than or equal to 0.89, observed with both the Chandra X-ray Observatory and the centimeter-wave SZE imaging system at the BIMA and OVRO interferometric arrays. An isothermal Beta-model with the central 100 kpc excluded from the X-ray data is used to model the intracluster medium and to measure global cluster properties. For each cluster, we measure the X-ray spectroscopic temperature, SZE gas mass, total mass, and integrated Compton gamma-parameters within r(sub 2500). Our measurements are in agreement with the expectations based on a simple self-similar model of cluster formation and evolution. We compare the cluster properties derived from our SZE observations with and without Chandra spatial and spectral information and find them to be in good agreement. We compare our results with cosmological numerical simulations and find that simulations that include radiative cooling, star formation, and feedback match well both the slope and normalization of our SZE scaling relations.

Bonamente, Massimiliano↗
DOE OSTI2021

Investigation of high-pressure turbulent jets using direct numerical simulation

Direct numerical simulations of free round jets at a Reynolds number (Re D ) of 5000, based on jet diameter (D) and jet-exit bulk velocity (U e ), are performed to study jet turbulence characteristics at supercritical pressures. The jet consists of nitrogen (N 2 ) that is injected into N 2 at the same temperature. To understand turbulent mixing, a passive scalar is transported with the flow at unity Schmidt number. Two sets of inflow conditions that model jets issuing from either a smooth contraction nozzle (laminar inflow) or a long pipe nozzle (turbulent inflow) are considered. By changing one parameter at a time, the simulations examine the jet-flow sensitivity to the thermodynamic condition (characterized in terms of the compressibility factor ($Z$) and the normalized isothermal compressibility), inflow condition and ambient pressure (p ∞ ) spanning perfect- to real-gas conditions. The inflow affects flow statistics in the near field (containing the potential core closure and the transition region) as well as further downstream (containing fully developed flow with self-similar statistics) at both atmospheric and supercritical p ∞ . The sensitivity to inflow is larger in the transition region, where the laminar-inflow jets exhibit dominant coherent structures that produce higher mean strain rates and higher turbulent kinetic energy than in turbulent-inflow jets. Decreasing $Z$ at a fixed supercritical p ∞ enhances pressure and density fluctuations (non-dimensionalized by local mean pressure and density, respectively), but the effect on velocity fluctuations depends also on the local flow dynamics. When $Z$ is reduced, large mean strain rates in the transition region of laminar-inflow jets significantly enhance velocity fluctuations (non-dimensionalized by local mean velocity) and scalar mixing, whereas the effects are minimal in jets from turbulent inflow.

42 ENGINEERING↗
DOE OSTI2021

Analytical modeling of the evolution of the nonlinearity parameter of sensitized stainless steel

Austenitic stainless steels are subject to the precipitation of chromium carbides (M 23 C 6 ) during exposure to high temperatures, causing these alloys to be susceptible to intergranular corrosion due to chromium depletion along grain boundaries. The acoustic nonlinearity parameter, β, shows sensitivity to the formation of carbides in these alloys. The Thermo-Calc TC-PRISMA module was used to model the nucleation and growth of grain boundary M 23 C 6 carbides. The model was verified with scanning transmission electron microscopy analysis that allowed measurements of the grain boundary precipitates. The paper introduces a reduced-order model of the acoustic nonlinearity based on the formation of misfit dislocations at the interface of the grain boundary precipitate and matrix to explain the change in β during isothermal aging. Furthermore, a direct relationship between the radius of the M 23 C 6 grain boundary carbides and β was observed and verified with nonlinear ultrasound measurements on 304L and 316L stainless steels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Towards a heat- and mass-balanced kinetic model of TATB decomposition

We report TATB (1,3,5-triamino-2,4,6-trinitrobenzene) was thermally degraded by two small-scale analytical methods – simultaneous differential scanning calorimetry and thermogravimetric analysis (SDT) and a hot-stage microscope with Fourier Transform Infrared (FTIR) analysis capabilities. SDT used ramped heating, isothermal soaking, and thermal pretreatment at various conditions. The heat flow and mass loss were monitored during various treatment conditions to derive chemical decomposition kinetics and Arrhenius parameters. FTIR experiments used isothermal heating, and changes were monitored spectroscopically. Solid samples generated at specific conditions were collected from both methods and were analyzed by DMSO extraction followed by chemical speciation by optical and mass spectrometric methods. Characterization provided the following reaction insights: TATB decreases in a sigmoidal pattern in isothermally heated samples. Other soluble products gradually increase in concentration and then abruptly decline in concentration during the second exotherm, such as diamino-dinitro-benzofurazan and amino-nitro-benzodifurazan. FTIR showed gradual changes in the amino and nitro functionality, shifting positions and decreasing intensity for the first 40 min. Then the solid gradually appeared more like an amorphous C with N incorporated, similar to previous studies on thermally degraded TATB-type materials. Extracted residues (DMSO-soluble components removed) examined by FTIR showed an abrupt change in chemical composition between 40- and 45-min isothermal treatment, indicating early forming solids are different than later forming residues. A reliable mass- and energy-balanced global reaction network must include at least two autocatalytic reactions, either in parallel or series, and at least one must have an explicit initiation reaction having a low activation energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

Fracture mechanics applied to elevated temperature crack growth

Twenty-six isothermal crack growth tests were performed on Hastelloy-X tubular specimens at a variety of temperatures and strain ranges. Conditions were selected to include nominally elastic and nominally plastic conditions. A number of parameters including the stress intensity factor, strain intensity factor, J-integral, Crack Opening Displacement, and Tompkins model were examined for their ability to correlate the data. Test conditions were selected such that growth rates at a single value of the parameter were obtained at radially different crack lengths, thus exploring the geometry independence of the correlating parameter. None of the parameters were fully satisfactory. However, COD calculated from J-integral appeared to be the most successful.

Jordan, E. H.↗
DOE OSTI2023

Establishing a Design Strategy for Corrosion Resistant Structural Materials in Molten Salt Technologies

To address the gap in material design strategies for structural components in molten salt-powered technologies, the current work aimed to identify key parameters governing corrosion induced degradation of materials during exposures in molten halide salts. The applicability of the data on dissolution rates of pure metals (e.g., Cr) in a given molten salt to enable predictions of corrosion behavior of multicomponent binary and ternary Fe-based model and Ni-based commercial alloys isothermally exposed in the purified KCl-MgCl 2 in quartz capsules at 700 and 800°C was evaluated. The influence of alloy chemistry on initial dissolution rates and the consequent impact on the time-dependent corrosion induced microstructural evolution in the alloy subsurface was predicted with a coupled thermodynamic-kinetic model. Finally, the model was well able to predict the Cr depletion in all investigated alloys demonstrating the ability of the proposed approach to enable a unique design strategy for structural materials in molten salts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Unraveling the multi-step crystallization mechanism of polytetrafluoroethylene, modified polytetrafluoroethylene, and their nanocomposites with boron nitride nanobarbs: Experimental insights and theoretical analysis

The non-isothermal crystallization behavior and kinetics of polytetrafluoroethylene (PTFE) composites with boron nitride nanobarb (BNNB), a new generation nanostructure with unique surface morphology and mechanical “barbs” have been analyzed, understanding these properties is essential for their high-end applications as thermal interface materials (TIM) for microwave, 5G and microelectronic devices. The analysis of the crystallization parameters includes crystallization onset, peak and end temperatures, crystallization half-life and overall crystallinity of PTFE, modified PTFE and their BNNB composites. The results were further analyzed using theoretical models such as the combined Avrami-Ozawa model. It was found that BNNB supports crystallization in the modified PTFE but shows minimal effect on the crystallization of PTFE. Due to the limitation of the classical theoretical models used above in fully characterizing the multi-step crystallization process of PTFE, an in-depth analysis using the model-free advanced isoconversional computation was used to characterize the PTFE crystallization based on the evolution of activation energy with fractional crystallinity and for the first time with temperature. Three kinetic regions were identified in the crystallization mechanism. Here, this study investigated the molecular organization and microstructural evolution of PTFE, modified PTFE and their composites during non-isothermal crystallization using advanced X-ray scattering measurements. An insight into the changes undergone by the material's microstructural units including crystallite size and morphology, lamellar thickness and lamellar interfacial layer thickness, and crystallographic phase dynamics during non-isothermal cooling from the melt, was provided here in this work. The effect of copolymer modification of PTFE and the inclusion of pristine and functionalized BNNB (a thermally conductive and electrically insulating ceramic) are both new investigations that provide valuable knowledge for the development of materials with strong matrix-nanofiller interaction and guidance for optimizing sintering and cooling cycles, two key steps in PTFE processing that largely affect the material microstructural features. Overall, the result of the three-part investigation demonstrates that BNNB supports crystallization in the modified PTFE up to 20 wt% concentration and at low and high cooling rates typically used in the industrial processing of PTFE.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Diskoseismology: Probing accretion disks. II - G-modes, gravitational radiation reaction, and viscosity

A scalar potential is used to derive a single partial differential equation governing the oscillation of a disk. The eigenfunctions and eigenfrequencies of a variety of disk models are found to fall into two main classes which are analogous to the p-modes and g-modes in the sun. Specifically, the eigenfunctions and eigenfrequencies of isothermal disks are computed, and the way in which these results can be generalized to other disk models is indicated. The (assumed) relatively small rates of growth or damping of the modes due to various mechanisms, in particular gravitational radiation reaction and parameterized models of viscosity are also computed. It is found that for certain parameters the p-modes are unstable to gravitational radiation reaction (CFS instability), while both the p-modes and g-modes are unstable to viscosity unless highly anisotropic viscosity models are considered.

Nowak, Michael A.↗