Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupling coefficient”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Photochemistry of Saturn's Atmosphere: Hydrocarbon Chemistry and Comparisons with ISO Observations - 1

To investigate the details of hydrocarbon photochemistry on Saturn, we have developed a one-dimensional diurnally averaged model that couples hydrocarbon and oxygen photochemistry, molecular and eddy diffusion, radiative transfer, and condensation. The model results are compared with observations from the Infrared Space Observatory (ISO) to place tighter constraints on molecular abundances, to better define Saturn's eddy diffusion coefficient profile, and to identify important chemical schemes that control the abundances of the observable hydrocarbons in Saturn's upper atmosphere. From the ISO observations, we determine that the column 12 densities of CH3, CH3C2H, and C4H2 above 10 mbar are 4 (sup +2) (sub -1.5) x 10 (exp 13) cm (sup -2), (1.1 plus or minus 0.3) x 10 (exp 15) cm (exp -2), and (1.2 plus or minus 0.3) x 10 (exp 14) cm (sup -2), respectively. The observed ISO emission features also indicate C2H2 mixing ratios of 1.2 (sup +0.9) (sub -0.6) x 10 (exp -6) at 0.3 mbar and (2.7 plus or minus 0.8) x 10 (exp -7) at 1.4 mbar, and a C2H6 mixing ratio of (9 plus or minus 2.5) x 10 (exp -6) at 0.5 mbar. Upper limits are provided for C2H4, CH2CCH2, C3H8, and C6H2 sensitivity of the model results to variations in the eddy diffusion coefficient profile, the solar flux, the CH4 photolysis branching ratios, the atomic hydrogen influx, and key reaction rates are discussed in detail. We find that C4H2 and CH3C2H are particularly good tracers of important chemical processes and physical conditions in Saturn's upper atmosphere, and C2H6 is a good tracer of the eddy diffusion coefficient in Saturn's lower stratosphere. The eddy diffusion coefficient must be smaller than approximately 3 x 10 (exp 4) sq cm s (sup -1) at pressures greater than 1 mbar in order to reproduce the C2H6 abundance inferred from ISO observations. The eddy diffusion coefficients in the upper stratosphere could be constrained by observations of CH3 radicals if the low-temperature chemistry of CH3 were better understood. We also discuss the implications of our modeling for aerosol formation in Saturn's lower stratosphere-diacetylene, butane, and water condense between approximately 1 and 300 mbar in our model and will dominate stratospheric haze formation at nonauroral latitudes. Our photochemical models will be useful for planning observational sequences and for analyzing data from the upcoming Cassini mission.

Moses, Julianne I.↗
NASA NTRS1995

Galileo Heat Shield Response Experiment Data Reduction and Analyses

The Galileo heat shield response experiment consists of 10 Analog Resistance Ablation Detectors (ARADs) embedded in the carbon phenolic heat shield. As the vehicle descends into the atmosphere of Jupiter, these ARADs will in effect measure the recession as a function of time at various locations on the heat shield. This recession data will be used to reconstruct the time-dependent shape of the vehicle. The shape data is of critical importance to the Atmospheric Reconstruction Experiment which must know mass and drag coefficient of the vehicle in order to determine the structure of the Jovian atmosphere. The data is also intrinsically useful as a database for evaluating the accuracy of our numerical codes and methodology for predicting highly radiating reentry flowfields and coupled heat shield material response. This paper will document the reduction of the data from the ARAD sensors, reconstruction of the heat shield shape history, and comparison of the actual heat shield response with previous calculations and experiments. If possible, we will also include results from new calculations which couple the shock layer flow with the transient material response.

Milos, Frank S.↗
NASA NTRS2004

Atomic Data and Spectral Line Intensities for Ar XII

Electron impact collision strengths; energy levels; oscillator strengths and spontaneous radiative decay rates are calculated for Ar XII. The configurations used are 2s(exp 2)2p(exp 3), 2s2p(exp4), 2p(exp 5), 2s22p23s, 2s(exp 2)2p(exp 2)3p, and 2s(exp 2)2p(exp 2)3d giving rise to 72 fine-structure levels in intermediate coupling. Collision strengths are calculated at five incident energies: 35,70, 105, 140 and 175 Ry. They are complemented by Coulomb-Born limits towards infinite collision energy for forbidden transitions and line strengths for optically allowed transitions. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates, statistical equilibrium equations for level populations are solved at electron densities covering the range of 10(exp 8)-101(exp 14) cm(exp -3) at an electron temperature of log T(sub e)/K = 6.4, corresponding to maximum abundance of Ar XII. Relative spectral line intensities are calculated.

Eissner, W.↗
NASA NTRS2004

Atomic Data and Spectral Line Intensities for Ar XII

Electron impact collision strengths, energy levels, oscillator strengths and spontaneous radiative decay rates are calculated for Ar XII. The configurations used are 2s(sup 2)p(sup 3), 2s(sup 2)p(sup 4), 2p(sup 5), 2s(sup 2)2p(sup 2)3s, 2s(sup 2)2p(sup 2)3p, and 2s(sup 2)2p(sup 2)3d giving rise to 72 fine-structure levels in intermediate coupling. Collision strengths are calculated at five incident energies: 35, 70, 105, 140 and 175 Ry. They are complemented by Coulomb-Born limits towards infinite collision energy for forbidden transitions and line strengths for optically allowed transitions. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates, statistical equilibrium equations for level populations are solved at electron densities covering the range of 10(exp 8)-10(exp 14)/cm(exp 3) at an electron temperature of log T(sub e)/K = 6.4, corresponding to maximum abundance of AR-XII. Relative spectral line intensities are calculated.

Eissner, W.↗
NASA NTRS2015

Small Scale Mass Flow Plug Calibration

A simple control volume model has been developed to calculate the discharge coefficient through a mass flow plug (MFP) and validated with a calibration experiment. The maximum error of the model in the operating region of the MFP is 0.54%. The model uses the MFP geometry and operating pressure and temperature to couple continuity, momentum, energy, an equation of state, and wall shear. Effects of boundary layer growth and the reduction in cross-sectional flow area are calculated using an in- integral method. A CFD calibration is shown to be of lower accuracy with a maximum error of 1.35%, and slower by a factor of 100. Effects of total pressure distortion are taken into account in the experiment. Distortion creates a loss in flow rate and can be characterized by two different distortion descriptors.

mass-flow measurements↗
NASA NTRS2020

Jig Twist Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft and High-Fidelity CFD Validation

This paper presents a jig twist optimization study of Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft using an in-house developed aero-structural analysis solver VSPAERO coupled to BEAM3D. A vortex-lattice model of the TTBW model is developed, and a transonic and viscous flow correction method is implemented in the VSPAERO model to account for transonic and viscous flow effects. A correction method for the wing-strut interference aerodynamics is developed and applied to the VSPAERO solver. Also, a structural dynamic finite-element model of the TTBW aircraft is developed. This finite-element model includes the geometric nonlinear effect due to the tension in the struts which causes a deflection-dependent nonlinear stiffness. The VSPAERO model is coupled to the corresponding finite-element model to provide a rapid aero-structural analysis. A design flight condition corresponding to Mach 0.745 at 42000 ft is selected for the TTBW aircraft jig twist optimization to reduce the drag coefficient. After the design is implemented, the drag coefficient of the twist optimized TTBW aircraft is reduced about 8 counts. At the end, a high-fidelity CFD solver FUN3D is used to validate the design.

BEAM3D↗
NASA NTRSDate not supplied

Gimbal Bearing Friction in the SLS Core Stage Thrust Vector Control System

The Space Launch System (SLS) Core Stage Thrust Vector Control (TVC) system is comprised of eight mechanical feedback Shuttle heritage Type III TVC actuators that vector the four Shuttle heritage RS-25 engines about a Shuttle heritage gimbal block/bearing. The MSFC Controls community has long regarded gimbal friction to be a negligible effect on the overall control of gimbaled RS-25 engines. This is corroborated by Space Shuttle test and flight data that does not appear to show degraded effects, nor limit cycling at the end of the shuttle flight. For this reason, friction was not expected to be a driving factor of performance and control of the reused RS-25 engines aboard the SLS. However, after test data showed a large shift in frequency behavior and a highly damped step-response in the time series, there was further investigation into what could have caused this behavior. Heritage friction models used in previous gimbal and ball bearings were evaluated such as Coulomb, Dahl and LuGre, but the single degree of freedom friction models alone were not enough to explain the behavior and shifts seen in the test data. This paper presents the additional findings and modeling efforts regarding friction on the RS-25 engines. Using the Two Actuator Operational Simulation (TAOS), the difference from modeling separate friction degrees of freedom to coupled degrees of freedom was investigated to deduce the effects of one axis’s movement on the other. Next, due to the vibration environment, a modified LuGre model has been proposed that adds an additional term to decrease the friction coefficient at low engine velocity amplitudes. Lastly, the addition of the stiffness in each half of the gimbal bearing has increased modeling fidelity by also adding the effect on the gimbal bearing bending in compliance to both the friction torque on the surface of the gimbal bearing and the actuator force that is forcing the engine in a specified direction. Through these effects, the time and frequency domain behavior seen in test can be characterized accurately.

Friction↗
NASA NTRSDate not supplied

Gimbal Bearing Friction in the SLS Core Stage Thrust Vector Control System

The Space Launch System (SLS) Core Stage Thrust Vector Control (TVC) system is comprised of eight mechanical feedback Shuttle heritage Type III TVC actuators that vector the four Shuttle heritage RS-25 engines about a Shuttle heritage gimbal block/bearing. The MSFC Controls community has long regarded gimbal friction to be a negligible effect on the overall control of gimbaled RS-25 engines. This is corroborated by Space Shuttle test and flight data that does not appear to show degraded effects, nor limit cycling at the end of the shuttle flight. For this reason, friction was not expected to be a driving factor of performance and control of the reused RS-25 engines aboard the SLS. However, after test data showed a large shift in frequency behavior and a highly damped step-response in the time series, there was further investigation into what could have caused this behavior. Heritage friction models used in previous gimbal and ball bearings were evaluated such as Coulomb, Dahl and LuGre, but the single degree of freedom friction models alone were not enough to explain the behavior and shifts seen in the test data. This paper presents the additional findings and modeling efforts regarding friction on the RS-25 engines. Using the Two Actuator Operational Simulation (TAOS), the difference from modeling separate friction degrees of freedom to coupled degrees of freedom was investigated to deduce the effects of one axis’s movement on the other. Next, due to the vibration environment, a modified LuGre model has been proposed that adds an additional term to decrease the friction coefficient at low engine velocity amplitudes. Lastly, the addition of the stiffness in each half of the gimbal bearing has increased modeling fidelity by also adding the effect on the gimbal bearing bending in compliance to both the friction torque on the surface of the gimbal bearing and the actuator force that is forcing the engine in a specified direction. Through these effects, the time and frequency domain behavior seen in test can be characterized accurately.

Friction↗
DOE OSTI2023

Revisiting the empirical particle-fluid coupling model used in DEM-CFD by high-resolution DEM-LBM-IMB simulations: A 2D perspective

The work investigates the applicability of the unresolved Computational Fluid Dynamics and Discrete Element Method (CFDDEM) technique based on empirical equations for fluid-particle coupling. We first carry out a series of representative volume element simulations using the high-resolution particle-resolved Lattice Boltzmann method and Discrete Element Method (LBMDEM) coupled by an Immersed Moving Boundary (IMB) scheme. Then, we compare the results obtained by both LBMDEM and empirical equations used in unresolved CFDDEM with analytical solutions. It is found that the existing empirical equations used in solving fluid-particle interactions in 2D CFDDEM fail to accurately calculate the hydrodynamic force applied to solid particles. The underlying reason is that the existing empirical models are obtained based on 3D experimental results and thus are not applicable to 2D problems. Based on the simulation results, a new drag coefficient model is then proposed. The estimated drag forces using the new model are compared favourably with the simulated ones, indicating the good performance of the proposed model.

42 ENGINEERING↗
DOE OSTI2022

Process-based modeling of soil nitrous oxide emissions from United States corn fields under different management and climate scenarios coupled with evaluation using regional estimates

Direct emissions of soil nitrous oxide during a growing season (N 2 O gs ) can be quantified with process-based models considering interactions between management, climate, and soil moisture when key data are available. We used an adapted “parameterized CENTURY/DAYCENT-model” ( p CENTURY) calibrated with crop growth and soil organic matter decay coefficients at the county-level for the estimation of N 2 O gs in the United States Corn Belt. Model estimated N 2 O-emissions from corn-based biofuels scenarios considering crop rotation, fertilizer inputs, tillage, and weather were compared against meta-summary of field observations from 55 studies. Both model and meta-summary ranked N 2 O gs -emissions to be corn > wheat > soybean phase while model likely underestimated cover crop N 2 O gs -emissions. The N 2 O gs -emissions and the associated emission factors (EFs) were modeled and summarized to be greater after anhydrous ammonia than urea application and from conventional tilled than non-tilled fields. Modeled and observed N 2 O gs -emissions after organic and inorganic fertilizer amendment did not differ due to high variability associated with the treatments. However, the organic fertilizer associated EFs were greater according to meta-summary data because of N input rates. Regionalized weather scenarios indicate hotspots for N 2 O gs -emissions can occur where crop N uptake is limited during dry years and in eastern states also during normal or wet seasons. The p CENTURY-derived N 2 O gs EFs (0.91 ± 0.19%) for counties investigated were only slightly lower than literature (1.07 ± 0.57%) or Tier-1 (1%) values. Our preliminary evaluation of regional soil moisture estimates showed reasonable agreement between monthly soil moisture estimates and the North American Soil Moisture Dataset during the growing season, but overestimation of soil moisture in winter-spring can influence the estimates of annual N 2 O emissions so future work is needed to calibrate soil moisture-associated model parameters. Our work provided scenario-based estimates of climate and management impacts on soil N 2 O gs -emissions together with valuable spatial insights into EFs that will be improved by more accurate information of fertilizer inputs and more temporally refined model evaluation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Combination of the W boson polarization measurements in top quark decays using ATLAS and CMS data at $\sqrt{s} =$ 8 TeV

The combination of measurements of the W boson polarization in top quark decays performed by the ATLAS and CMS collaborations is presented. The measurements are based on proton-proton collision data produced at the LHC at a centre-of-mass energy of 8 TeV, and corresponding to an integrated luminosity of about 20 fb$^{−1}$ for each experiment. The measurements used events containing one lepton and having different jet multiplicities in the final state. The results are quoted as fractions of W bosons with longitudinal (F$_{0}$), left-handed (F$_{L}$), or right-handed (F$_{R}$) polarizations. The resulting combined measurements of the polarization fractions are F$_{0}$ = 0.693 ± 0.014 and F$_{L}$ = 0.315 ± 0.011. The fraction F$_{R}$ is calculated from the unitarity constraint to be F$_{R}$ = −0.008 ± 0.007. These results are in agreement with the standard model predictions at next-to-next-to-leading order in perturbative quantum chromodynamics and represent an improvement in precision of 25 (29)% for F$_{0}$ (F$_{L}$) with respect to the most precise single measurement. A limit on anomalous right-handed vector (V$_{R}$), and left- and right-handed tensor (g$_{L}$, g$_{R}$) tWb couplings is set while fixing all others to their standard model values. The allowed regions are [−0.11, 0.16] for V$_{R}$, [−0.08, 0.05] for g$_{L}$, and [−0.04, 0.02] for g$_{R}$, at 95% confidence level. Limits on the corresponding Wilson coefficients are also derived.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

A continuous analogue of the tensor-train decomposition

We develop new approximation algorithms and data structures for representing and computing with multivariate functions using the functional tensor-train (FT), a continuous extension of the tensor-train (TT) decomposition. The FT represents functions using a tensor-train ansatz by replacing the three-dimensional TT cores with univariate matrix-valued functions. The main contribution of this paper is a framework to compute the FT that employs adaptive approximations of univariate fibers, and that is not connected to any tensorized discretization. The algorithm can be coupled with any univariate linear or nonlinear approximation procedure. We demonstrate that this approach can generate multivariate function approximations that are several orders of magnitude more accurate, for the same cost, than those based on the conventional approach of compressing the coefficient tensor of a tensor-product basis. Our approach is in the spirit of other continuous computation packages such as Chebfun, and yields an algorithm which requires the computation of “continuous” matrix factorizations such as the LU and QR decompositions of vector-valued functions. To support these developments, we describe continuous versions of an approximate maximum-volume cross approximation algorithm and of a rounding algorithm that re-approximates an FT by one of lower ranks. We demonstrate that our technique improves accuracy and robustness, compared to TT and quantics-TT approaches with fixed parameterizations, of high-dimensional integration, differentiation, and approximation of functions with local features such as discontinuities and other nonlinearities.

42 ENGINEERING↗
DOE OSTI2024

A comprehensive numerical investigation on spray models for Direct-Injection Spark-Ignition engines

Gasoline direct-injection spark-ignition (DISI) engines generate a large portion of their unburned hydrocarbon (UHC) and soot emissions during the cold-start phase. A predictive computational fluid dynamics (CFD) modeling framework can be used to understand the physical processes that characterize fuel spray evolution and fuel-film formation at cold start conditions, which can help to reduce engine-out particulate emissions. This study systematically evaluated spray submodels and developed a set of simulation best practices for physical-numerical submodels with the goal of enabling accurate simulations of liquid spray behavior in a DISI engine. Three comprehensive experimental datasets containing free-spray projected liquid volume (PLV), liquid volume fraction (LVF), and near-field X-ray radiography data were used to validate the simulation results and evaluate the spray submodels. Systematic analysis delved into injected parcel distribution, droplet collision, spray breakup, and evaporation via a detailed assessment of the relevant spray submodels. Moreover, the effects of turbulence models and the initial turbulent flow properties on the liquid spray evolution were examined. Based on extensive calibration efforts, a set of simulation best practices for the free spray was developed and validated against the PLV/LVF data. Simulation results indicated that the uniform distribution for parcel initialization, coupled with appropriate droplet collision submodels, provides an improved spray morphology compared to the cluster distribution. The findings also underscored the importance of calibrating the Kelvin-Helmholtz Rayleigh-Taylor (KH-RT) breakup model constants and droplet heat transfer coefficient scaling factor to achieve favorable agreement regarding measured liquid penetration and spray widths. In conclusion, this study marks a substantial stride towards accurately predicting fuel film evolution and soot formation within DISI engine performance.

ECN Spray G↗
DOE OSTI2022

Atmospheric Brown Carbon on the Tibetan Plateau: Regional Differences in Chemical Composition and Light Absorption Properties

In this work, fine particulate matter (PM 2.5 ) filter samples were collected from two high-altitude remote sites located in the southern (QOMS) and northern (WLG) regions of the Tibetan Plateau (TP) to explore the regional differences in brown carbon (BrC) properties. Chemical differences in BrC composition representative of these two areas were inferred from molecular-level analysis of the PM 2.5 samples using high-performance liquid chromatography coupled with photodiode array and high-resolution mass spectrometry detectors. The results show that more polar BrC chromophores were abundant in QOMS samples, while contributions from polar and less polar chromophores were comparable in WLG samples. A higher mass absorption coefficient of BrC was observed at QOMS than at WLG. Strong BrC chromophores in QOMS samples were identified as oxygenated aromatics and nitrophenol compounds, while organosulfate compounds were found in the WLG sample. The results of this study indicate regional differences in BrC chromophores and provide insights into their sources and chemical processes, which should be considered for predictive understanding and modeling of the radiative forcing of aerosol in the TP area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Incorporating the Acclimation of Photosynthesis and Leaf Respiration in the Noah‐MP Land Surface Model: Model Development and Evaluation

Abstract Realistic simulation of leaf photosynthetic and respiratory processes is needed for accurate prediction of the global carbon cycle. These two processes systematically acclimate to long‐term environmental changes by adjusting photosynthetic and respiratory traits (e.g., the maximum photosynthetic capacity at 25°C ( V cmax,25 ) and the leaf respiration rate at 25°C ( R 25 )) following increasingly well‐understood principles. While some land surface models (LSMs) now account for thermal acclimation, they do so by assigning empirical parameterizations for individual plant functional types (PFTs). Here, we have implemented an Eco‐Evolutionary Optimality (EEO)‐based scheme to represent the universal acclimation of photosynthesis and leaf respiration to multiple environmental effects, and that therefore requires no PFT‐specific parameterizations, in a standard version of the widely used LSM, Noah MP. We evaluated model performance with plant trait data from a 5‐year experiment and extensive global field measurements, and carbon flux measurements from FLUXNET2015. We show that observed R 25 and V cmax,25 vary substantially both temporally and spatially within the same PFT ( C.V. >20%). Our EEO‐based scheme captures 62% of the temporal and 70% of the spatial variations in V cmax,25 (73% and 54% of the variations in R 25 ). The standard scheme underestimates gross primary production by 10% versus 2% for the EEO‐based scheme and generates a larger spread in r (correlation coefficient) across flux sites (0.79 ± 0.16 vs. 0.84 ± 0.1, mean ± S.D.). The standard scheme greatly overestimates canopy respiration (bias: ∼200% vs. 8% for the EEO scheme), resulting in less CO 2 uptake by terrestrial ecosystems. Our approach thus simulates climate‐carbon coupling more realistically, with fewer parameters.

Ren, Yanghang [Department of Earth System Science ↗
DOE OSTI2017

Temperature dependence of the energy bandgap of multi-layer hexagonal boron nitride

The temperature dependence of the energy bandgap of hexagonal boron nitride (h-BN) has been probed via photoluminescence emission characteristics of a donor-to-acceptor pair transition in a 20-layer h-BN epilayer. The results indicate that the universal behavior of bandgap decreasing with temperature is absent in multi-layer h-BN. Below 100 K, the bandgap energy variation with temperature, E g vs. T, is dominated by the electron-phonon coupling and conforms to the common behavior of redshift with an increase in temperature. At T > 100 K, the bandgap shows an unusual blueshift with temperature, which can be attributed to the unique behavior of the in-plane thermal expansion coefficient of h-BN that becomes negative above around 60 K. Although both graphite and h-BN have negative thermal expansion coefficients in a broad temperature range, graphite has a zero energy bandgap, which makes h-BN a unique semiconductor to exhibit this unusual temperature dependence of the energy bandgap.

42 ENGINEERING↗
DOE OSTI2017

Fluctuation-enhanced electric conductivity in electrolyte solutions

Significance Using fluctuating hydrodynamics, we demonstrate that thermal fluctuations contribute to charge transport in binary electrolyte solutions. We show the existence of an enhancement, or renormalization, of the electric conductivity due to the coupling between fluctuations of charge and fluid velocity. This coupling results in nontrivial corrections to the classical Poisson–Nernst–Planck equations, which are of the order of the square root of the salt concentration and therefore significant even for dilute solutions. Our calculations predict a cation–anion cross-diffusion coefficient that is in quantitative agreement with experimental measurements. Our findings have important implications for the fields of both mesoscale hydrodynamics and electrolyte transport.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Effects of mistuning on blade torsional flutter

An analytical model for the prediction of fan blade flutter is presented and evaluated using data from NASA tests on an advanced high performance engine. For the cascade conditions appropriate to the test points studied, the aerodynamic theory cannot predict subcritical flutter. Under the assumptions of a tuned assembly, the imaginary part of the aerodynamic coefficients does indicate flutter for a limited number of interblade phase angles, but these interblade phase angles are close to those at which the acoustic resonance is predicted. Upon using the individual blade frequencies and solving the mistuned system with aerodynamic coupling only, the results show a stable system. Eigenvectors calculated for the mistuned system demonstrate the presence of several harmonics in each mistuned mode. Inclusion of both mechanical and aerodynamic coupling in the solution of the eigenproblem influences not only the frequencies but also damping in the system with a trend toward stability.

Srinivasan, A. V.↗