Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generative diffusion”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Effects of vacuum-ultraviolet irradiation on copper penetration into low-k dielectrics under bias-temperature stress

The effects of vacuum-ultraviolet (VUV) irradiation on copper penetration into non-porous low-k dielectrics under bias-temperature stress (BTS) were investigated. By employing x-ray photoelectron spectroscopy depth-profile measurements on both as-deposited and VUV-irradiated SiCOH/Cu stacks, it was found that under the same BTS conditions, the diffusion depth of Cu into the VUV-irradiated SiCOH is higher than that of as-deposited SiCOH. On the other hand, under the same temperature-annealing stress (TS) without electric bias, the Cu distribution profiles in the VUV-irradiated SiCOH were same with that for the as-deposited SiCOH. The experiments suggest that in as-deposited SiCOH, the diffused Cu exists primarily in the atomic state, while in VUV-irradiated SiCOH, the diffused Cu is oxidized by the hydroxyl ions (OH{sup −}) generated from VUV irradiation and exists in the ionic state. The mechanisms for metal diffusion and ion injection in VUV irradiated low-k dielectrics are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Recovering doping profiles in semiconductor devices with the Boltzmann-Poisson model

We investigate numerically an inverse problem related to the Boltzmann-Poisson system of equations for transport of electrons in semiconductor devices. The objective of the (ill-posed) inverse problem is to recover the doping profile of a device, presented as a source function in the mathematical model, from its current-voltage characteristics. To reduce the degree of ill-posedness of the inverse problem, we proposed to parameterize the unknown doping profile function to limit the number of unknowns in the inverse problem. We showed by numerical examples that the reconstruction of a few low moments of the doping profile is possible when relatively accurate time-dependent or time-independent measurements are available, even though the later reconstruction is less accurate than the former. We also compare reconstructions from the Boltzmann-Poisson (BP) model to those from the classical drift-diffusion-Poisson (DDP) model, assuming that measurements are generated with the BP model. We show that the two type of reconstructions can be significantly different in regimes where drift-diffusion-Poisson equation fails to model the physics accurately. However, when noise presented in measured data is high, no difference in the reconstructions can be observed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Isochrones for old (>5 Gyr) stars and stellar populations. I. Models for –2.4 ≤ [Fe/H] ≤+0.6, 0.25 ≤ Y ≤ 0.33, and –0.4 ≤ [α/Fe] ≤+0.4

Canonical grids of stellar evolutionary sequences have been computed for the helium mass-fraction abundances Y = 0.25, 0.29, and 0.33, and for iron abundances that vary from –2.4 to +0.4 (in 0.2 dex increments) when [α/Fe] =+0.4, or for the ranges –2.0 ≤ [Fe/H] ≤+0.6, –1.8 ≤ [Fe/H] ≤+0.6 when [α/Fe] =0.0 and –0.4, respectively. The grids, which consist of tracks for masses from 0.12 M{sub ⊙} to 1.1-1.5 M{sub ⊙} (depending on the metallicity) are based on up-to-date physics, including the gravitational settling of helium (but not metals diffusion). Interpolation software is provided to generate isochrones for arbitrary ages between ≈5 and 15 Gyr and any values of Y, [α/Fe], and [Fe/H] within the aformentioned ranges. Comparisons of isochrones with published color-magnitude diagrams (CMDs) for the open clusters M67 ([Fe/H] ≈0.0) and NGC 6791 ([Fe/H] ≈0.3) and for four of the metal-poor globular clusters (47 Tuc, M3, M5, and M92) indicate that the models for the observed metallicities do a reasonably good job of reproducing the locations and slopes of the cluster main sequences and giant branches. The same conclusion is reached from a consideration of plots of nearby subdwarfs that have accurate Hipparcos parallaxes and metallicities in the range –2.0 ≲ [Fe/H] ≲ –1.0 on various CMDs and on the (log T {sub eff}, M{sub V} ) diagram. A relatively hot temperature scale similar to that derived in recent calibrations of the infrared flux method is favored by both the isochrones and the adopted color transformations, which are based on the latest MARCS model atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Molecular Dynamics-based Simulations of Bulk/Interfacial Structures and Diffusion Behaviors in Nuclear Waste Glasses

This NEUP Project aimed to generate accurate atomic structural models of nuclear waste glasses by using large-scale molecular dynamics-based computer simulations and to use these models to investigate self-diffusion behaviors, interfacial structures, and hydrated gel structures formed during dissolution of these glasses. The goal was to obtain realistic and accurate short and medium range structures of these complex oxide glasses, to provide a mechanistic understanding of the dissolution behaviors, and to generate reliable information with predictive power in designing nuclear waste glasses for long-term geological storage. Looking back of the research accomplishments of this project, most of the scientific goals initially proposed have been achieved through intensive research in the three and a half year period of the project. This project has also generated a wealth of scientific data and vibrant discussions with various groups through collaborations within and outside of this project. Throughout the project one book chapter and 14 peer reviewed journal publications have been generated (including one under review) and 16 presentations (including 8 invited talks) have been made to disseminate the results of this project in national and international conference. Furthermore, this project has trained several outstanding graduate students and young researchers for future workforce in nuclear related field, especially on nuclear waste immobilization. One postdoc and four PhD students have been fully or partially supported through the project with intensive training in the field material science and engineering with expertise on glass science and nuclear waste disposal

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1954

Investigation of Conical Subsonic Diffusers for Ram-jet Engines

The efficiency of a 30 degree conical diffuser was improved as much as 20 percent and separation was eliminated by the use of vortex generators. The use of splitter cones gave only small efficiency gains, but relatively uniform diffuser-outlet velocity profiles were obtained with the better designs. A configuration which incorporated both vortex generators and a splitter cone gave efficiencies higher than those obtained with any other splitter-cone configuration and also higher than those obtained when the same vortex generators were used without the splitter cone. Moderate diffuser-efficiency increases were obtained with guide vanes; however, in all cases combustion occurred in the vane wakes.

Farley, John M↗
DOE OSTI2025

Physically constrained 3D diffusion for inverse design of fiber-reinforced polymer composite materials

Designing fiber-reinforced polymer composites (FRPCs) with a tailored nonlinear stress-strain response is crucial for applications such as energy absorption in crash structures, flexible robotics, and impact-resistant protective gear. However, the inherent complexities of composite materials and the multitude of parameters involved, render traditional design and optimization methods inadequate for achieving effective inverse design of composites. In this paper, we present an AI-based inverse design framework that effectively and efficiently generates FRPCs with targeted nonlinear stress-strain responses. We introduce a physically constrained diffusion model (PC3D_Diffusion) capable of managing the complexities of composite materials and producing detailed, high-quality designs. We propose a loss-guided, learning-free approach to generate physically feasible microstructure designs by explicitly enforcing physical constraints during the generation process. For training purposes, 1.35 million FRPC samples were created, and their corresponding stress-strain curves were computed using established physics-based computational models. The results show that PC3D_Diffusion consistently generates high-quality designs with tailored mechanical behaviors, while guaranteeing compliance with the physical constraints. PC3D_Diffusion advances FRPC inverse design and may facilitate the inverse design of other 3D materials, offering potential applications in industries reliant on materials with custom mechanical properties.

Xu, Pei [Clemson Univ., SC (United States)]↗
DOE OSTI2023

High‐Throughput Evaluation of Hardening Coefficients of Eight Alloying Elements in Magnesium

Liquid–solid diffusion couples (LSDCs) are employed to generate a composition gradient in the single‐phase hexagonal closed‐packed (hcp) solid solution with compositions up to the solubility limit of various solutes in Mg. Nanoindentation scanning across the composition gradient in LSDCs allows effective evaluation of composition‐dependent hardness of eight alloying elements (Al, Ca, Ce, Gd, Li, Sn, Y, and Zn) in the hcp Mg phase. The hardening coefficients, an indicator of the potency of solid‐solution hardening, are evaluated from the measured composition‐hardness data and correlated with various materials properties such as atomic radius, shear modulus, and elastic modulus of the solutes. The rank of hardening potency of Al, Gd, Sn, Y, and Zn measured by nanoindentation is in good agreement with that measured by microindentation reported in the literature. The hardening coefficient (potency) from the strongest to the weakest is Ce > Ca > Y ≈ Gd > Zn > Al ≈ Sn > Li in Mg‐based hcp binary solid solutions. The hardening coefficient is found to be closely correlated with the strengthening potency.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mesoscale modeling and semi-analytical approach for the microstructure-aware effective thermal conductivity of porous polygranular materials

Here we established a comprehensive modeling approach for investigating the microstructure-aware effective thermal conductivity ($κ_{eff}$) for porous microstructures containing solid particles and gaseous pores. Our approach combines the mesoscale computational modeling framework and the semi-analytical method, allowing for efficient prediction of $κ_{eff}$ for realistic porous microstructures, while considering complicated microstructural thermal conduction pathways effectively in the prediction. We used the diffuse-interface mesoscale computational model to generate extensive simulated $κ_{eff}$ data for realistic digital representations of microstructures with wide ranges of porosity ($f_p$), thermal conductivity of the gas phase ($κ_g$), and thermal conductivity of the solid phase ($κ_s$). From the simulated data, we identified two property variation regimes for $κ_{eff}$: (1) a slow $κ_{eff}$ increase for $κ_s ~ κ_g$; and (2) a faster $κ_{eff}$ increase for $κ_s \gg κ_g$. To capture the key features of the relationship between the microstructure and $κ_{eff}$, we derived a semi-analytical model by introducing structure and intensification factors. The two new factors incorporate the calibrated effective contribution of the solid volume with $κ_s$ and additional interfacial effects into the prediction of $κ_{eff}$, respectively, allowing for consideration of parallel, serial, and interfacial conduction mechanisms effectively. Using the selected simulation data, we quantified key model parameters within the semi-analytical model and verified that the parameterized model exhibits excellent agreement with simulated $κ_{eff}$ for the entire range of the parameter space.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Spatially resolved probing of electrochemical reactions via energy discovery platforms

The electrochemical reactivity of solid surfaces underpins functionality of a broad spectrum of materials and devices ranging from energy storage and conversion, to sensors and catalytic devices. The surface electrochemistry is, however, a complex process, controlled by the interplay of charge generation, field-controlled and diffusion-controlled transport. Here we explore the fundamental mechanisms of electrochemical reactivity on nanocrystalline ceria, using the synergy of nanofabricated devices and time-resolved Kelvin probe force microscopy (tr-KPFM), an approach we refer to as energy discovery platform. Through tr-KPFM, the surface potential mapping in both the space and time domains and current variation over time are obtained, enabling analysis of local ionic and electronic transport and their dynamic behavior on the 10 ms to 10 s scale. Based on their different responses in the time domain, conduction mechanisms can be separated and identified in a variety of environmental conditions, such as humidity and temperature. The theoretical modeling of ion transport through finite element method allows for creation of a minimal model consistent with observed phenomena, and establishing of the dynamic characteristics of the process, including mobility and diffusivity of charged species. Furthermore, the future potential of the energy discovery platforms is also discussed.

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

Hydrate-filled Fracture Formation at Keathley Canyon 151, Gulf of Mexico, and Implications for Non-vent Sites

Near-vertical hydrate-filled fractures are found in subseafloor marine muds, at both advective methane vent sites and at sites without obvious methane and fluid advection (non-vent sites). At non-vent sites, the mechanisms that transport methane to the fractures and control how hydrate-filled fractures form are not well understood. However, these mechanisms are important to establish, as most of Earth's natural gas hydrate is likely bound in marine mud systems. In this work, we focus on understanding the origin of hydrate and how fracture form at non-vent sites by examining previously hydrate-bearing fractures in conventional cores taken from Keathley Canyon 151, U.S. northern Gulf of Mexico, drilled by the Gas Hydrate Joint Industry Project in 2005. We combine information from well logs, sediment cores, and science party results and add new X-ray computed tomography of archival sections and scanning electron microcopy of core samples to develop a conceptual model. We propose that locally generated microbial methane is transported via diffusion from small pores in marine mud into biomineralized burrows with larger pore size in a process called short-range migration. Hydrate forms in burrows once the methane diffuses into them and the dissolved methane concentration exceeds the solubility threshold. When hydrate fills a burrow, heave from additional hydrate growth places stress on the burrow edges, expands the fracture, and creates additional void space in which methane can diffuse and continue forming hydrate. Fractures slowly propagate in the direction of maximum principal stress.

58 GEOSCIENCES↗
DOE OSTI2021

Preliminary Analysis of Source Physics Experiment Explosion-Triggered Microseismicity Using the Back-Projection Method

A series of four chemical explosions were detonated in a deep borehole within the Yucca Flat Dry Alluvium Geology (DAG) at the Nevada National Security Site between 2018 and 2019. The two larger chemical explosions of 50 tons (DAG-2) and 10 tons (DAG-4) TNT equivalent yield triggered energetic aftershock sequences numbering 1392 and 347 microearthquakes, respectively, within the first 10 days. No significant aftershock activity was observed for the two smaller 1-ton explosions (DAG-1 and DAG-3). Here, we used a back-projection method based on travel-time migration and stacking of signal-to-noise ratio traces to detect, associate and locate aftershocks from a subset of 22-geophones within a larger 2 urn:x-wiley:21699313:media:jgrb54910:jgrb54910-math-0001 2 km seismic array surrounding the borehole. The aftershocks located within 300 m of the borehole and the depths were above the working points of 300 and 50 m depths of DAG-2 and DAG-4, respectively, ruling out triggering slip on geologic faults or disturbances beneath neighboring collapse craters. DAG-2 and DAG-4 aftershocks decayed at similar rates, with power-law exponents of p = 1.48 and p = 1.49, respectively. These decay rates are comparable to aftershocks sequences triggered by earthquakes and historical nuclear explosions at Yucca Flat. A smooth power-law aftershock decay within the first 10 days suggests a triggering mechanism from explosion generated stress relaxation due to the diffusion of high gas pressures in the cavity and radial fractures. A more random and episodic aftershock rate would be expected due to cavity collapse or falling rubble in chimney formation.

58 GEOSCIENCES↗
DOE OSTI2018

In situ study of nucleation and growth dynamics of Au nanoparticles on MoS 2 nanoflakes

Two-dimensional (2D) substrates decorated with metal nanoparticles offer new opportunities to achieve high-performance catalytic behavior. However, little is known on how the substrates control the nucleation and growth processes of the nanoparticles. This paper presents the visualization of dynamic nucleation and growth processes of gold nanoparticles on ultrathin MoS 2 nanoflakes by in situ liquid-cell transmission electron microscopy (TEM). The galvanic displacement resulting in Au nuclei formation on MoS 2 was observed in real time inside the liquid cell. We found that the growth mechanism of Au particles on pristine MoS 2 is in between diffusion-limited and reaction-limited, possibly due to the presence of electrochemical Ostwald ripening. A larger size distribution and more orientation variation is observed for the Au particles along the MoS 2 edge than on the interior. Differing from pristine MoS 2 , sulfur vacancies on MoS 2 induce Au particle diffusion and coalescence during the growth process. Density functional theory (DFT) calculations show that the size difference is because the exposed molybdenum atoms at the edge with dangling bonds can strongly interact with Au atoms, whereas sulfur atoms on the MoS 2 interior have no dangling bonds and weakly interact with gold atoms. In addition, S vacancies on MoS 2 generate strong nucleation centers that can promote diffusion and coalescence of Au nanoparticles. Furthermore, the present work provides key insights into the role of 2D materials in controlling the size and orientation of noble metal nanoparticles vital to the design of next generation catalysts.

2D materials↗
DOE OSTI2014

Bulk ion acceleration and particle heating during magnetic reconnection in a laboratory plasma

Bulk ion acceleration and particle heating during magnetic reconnection are studied in the collisionless plasma of the Magnetic Reconnection Experiment (MRX). The plasma is in the two-fluid regime, where the motion of the ions is decoupled from that of the electrons within the ion diffusion region. The reconnection process studied here is quasi-symmetric since plasma parameters such as the magnitude of the reconnecting magnetic field, the plasma density, and temperature are compatible on each side of the current sheet. Our experimental data show that the in-plane (Hall) electric field plays a key role in ion heating and acceleration. The electrostatic potential that produces the in-plane electric field is established by electrons that are accelerated near the electron diffusion region. The in-plane profile of this electrostatic potential shows a “well” structure along the direction normal to the reconnection current sheet. This well becomes deeper and wider downstream as its boundary expands along the separatrices where the in-plane electric field is strongest. Since the in-plane electric field is 3–4 times larger than the out-of-plane reconnection electric field, it is the primary source of energy for the unmagnetized ions. With regard to ion acceleration, the Hall electric field causes ions near separatrices to be ballistically accelerated toward the outflow direction. Ion heating occurs as the accelerated ions travel into the high pressure downstream region. This downstream ion heating cannot be explained by classical, unmagnetized transport theory; instead, we conclude that ions are heated by re-magnetization of ions in the reconnection exhaust and collisions. Two-dimensional (2-D) simulations with the global geometry similar to MRX demonstrate downstream ion thermalization by the above mechanisms. Electrons are also significantly heated during reconnection. The electron temperature sharply increases across the separatrices and peaks just outside of the electron diffusion region. Unlike ions, electrons acquire energy mostly from the reconnection electric field, and the energy gain is localized near the X-point. However, the increase in the electron bulk flow energy remains negligible. These observations support the assertion that efficient electron heating mechanisms exist around the electron diffusion region and that the heat generated there is quickly transported along the magnetic field due to the high parallel thermal conductivity of electrons. Classical Ohmic dissipation based on the perpendicular Spitzer resistivity is too small to balance the measured heat flux, indicating the presence of anomalous electron heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Driven diffusion against electrostatic or effective energy barrier across α-hemolysin

We analyze the translocation of a charged particle across an α-Hemolysin (αHL) pore in the framework of a driven diffusion over an extended energy barrier generated by the electrical charges of the αHL. A one-dimensional electrostatic potential is extracted from the full 3D solution of the Poisson’s equation. We characterize the particle transport under the action of a constant forcing by studying the statistics of the translocation time. We derive an analytical expression of translocation time average that compares well with the results from Brownian dynamic simulations of driven particles over the electrostatic potential. Moreover, we show that the translocation time distributions can be perfectly described by a simple theory which replaces the true barrier by an equivalent structureless square barrier. Remarkably, our approach maintains its accuracy also for low-applied voltage regimes where the usual inverse-Gaussian approximation fails. Finally, we discuss how the comparison between the simulated time distributions and their theoretical prediction results to be greatly simplified when using the notion of the empirical Laplace transform technique.

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

Point defect states in Sb-doped germanium

Defect states in n-type Sb-doped germanium were investigated by deep-level transient spectroscopy. Cobalt-60 gamma rays were used to generate isolated vacancies and interstitials which diffuse and react with impurities in the material to form four defect states (E{sub 37}, E{sub 30}, E{sub 22}, and E{sub 21}) in the upper half of the bandgap. Irradiations at 77 K and 300 K as well as isothermal anneals were performed to characterize the relationships between the four observable defects. E{sub 37} is assigned to the Sb donor-vacancy associate (E-center) and is the only vacancy containing defect giving an estimate of 2 × 10{sup 11 }cm{sup −3} Mrad{sup −1} for the uncorrelated vacancy-interstitial pair introduction rate. The remaining three defect states are interstitial associates and transform among one another. Conversion ratios between E{sub 22}, E{sub 21}, and E{sub 30} indicate that E{sub 22} likely contains two interstitials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2012

DYNAMICS OF THE SOLAR MAGNETIC BRIGHT POINTS DERIVED FROM THEIR HORIZONTAL MOTIONS

The subarcsecond bright points (BPs) associated with the small-scale magnetic fields in the lower solar atmosphere are advected by the evolution of the photospheric granules. We measure various quantities related to the horizontal motions of the BPs observed in two wavelengths, including the velocity autocorrelation function. A 1 hr time sequence of wideband H{alpha} observations conducted at the Swedish 1 m Solar Telescope (SST) and a 4 hr Hinode G-band time sequence observed with the Solar Optical Telescope are used in this work. We follow 97 SST and 212 Hinode BPs with 3800 and 1950 individual velocity measurements, respectively. For its high cadence of 5 s as compared to 30 s for Hinode data, we emphasize more the results from SST data. The BP positional uncertainty achieved by SST is as low as 3 km. The position errors contribute 0.75 km{sup 2} s{sup -2} to the variance of the observed velocities. The raw and corrected velocity measurements in both directions, i.e., (v{sub x} , v{sub y} ), have Gaussian distributions with standard deviations of (1.32, 1.22) and (1.00, 0.86) km s{sup -1}, respectively. The BP motions have correlation times of about 22-30 s. We construct the power spectrum of the horizontal motions as a function of frequency, a quantity that is useful and relevant to the studies of generation of Alfven waves. Photospheric turbulent diffusion at timescales less than 200 s is found to satisfy a power law with an index of 1.59.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

CMB spectral distortions as solutions to the Boltzmann equations

We propose to re-interpret the cosmic microwave background spectral distortions as solutions to the Boltzmann equation. This approach makes it possible to solve the second order Boltzmann equation explicitly, with the spectral y distortion and the momentum independent second order temperature perturbation, while generation of μ distortion cannot be explained even at second order in this framework. We also extend our method to higher order Boltzmann equations systematically and find new type spectral distortions, assuming that the collision term is linear in the photon distribution functions, namely, in the Thomson scattering limit. As an example, we concretely construct solutions to the cubic order Boltzmann equation and show that the equations are closed with additional three parameters composed of a cubic order temperature perturbation and two cubic order spectral distortions. The linear Sunyaev-Zel'dovich effect whose momentum dependence is different from the usual y distortion is also discussed in the presence of the next leading order Kompaneets terms, and we show that higher order spectral distortions are also generated as a result of the diffusion process in a framework of higher order Boltzmann equations. The method may be applicable to a wider class of problems and has potential to give a general prescription to non-equilibrium physics.

79 ASTRONOMY AND ASTROPHYSICS↗