Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermal. Parameter”

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

Towards inverse microstructure-centered materials design using generative phase-field modeling and deep variational autoencoders

The field of Integrated Computational Materials Engineering (ICME) combines a broad range of methods to study materials’ responses over a spectrum of length scales. A relatively unexplored aspect of microstructure-sensitive materials design is uncertainty propagation and quantification (UP/UQ) of materials’ microstructure, as well as establishing process-structure–property (PSP) relationships for inverse material design. In this study, an efficient UP technique built on the idea of changing probability measures and a deep generative unsupervised representative machine learning method for microstructure-based design of thermal conductivity of materials is proposed. Probability measures are used to represent microstructure space, and Wasserstein metrics are used to test the efficiency of the UP method. By using deep Variational AutoEncoder (VAE), we identify the correlations between the material/process parameters and the thermal conductivity of heterogeneous dual-phase microstructures. Through high-throughput screening, UP, and the deep-generative VAE method, PSP relationships that are too complex can be revealed by exploiting the materials’ design space with an emphasis on microstructures. As a last point, we demonstrate generative machine learning serves as a useful tool for inverse microstructure-centered materials design, and we demonstrate this by examining the inverse design of thermal conductivity in nano-structured materials. Here, the results reveal the effects of morphology, volume fraction, characteristic length scale, and the individual thermal diffusivity of phases on the thermal conductivity of dual-phase alloys. Our findings emphasize the advantages of high-throughput phase-field modeling and generative deep learning for linking PSP and inverse microstructure-centered materials design.

36 MATERIALS SCIENCE↗

Machine learning surrogates for ion energy–angle distributions in thermal and RF plasma sheaths

Ion energy–angle distributions (IEADs) at material surfaces are a critical input for plasma–material interaction (PMI) studies in fusion devices, yet they are computationally expensive to obtain using particle-in-cell (PIC) simulations. In this work, we develop a machine learning surrogate based on a deep deconvolutional neural network (DDeCNN) trained on large databases generated with the hPIC2 code. The surrogate is capable of reconstructing IEADs from sheath parameters for both thermal and radio-frequency (RF) plasmas, including cases with multiple ion species. Across thousands of test cases, the model achieves high accuracy, with over 97 % of predictions classified as good or average based on standard error metrics (MAE, MSE, L2). Even in the more challenging RF and multi-species regimes, the surrogate reliably captures the multi-peak structure of PIC results. Once trained, the surrogate produces IEADs in milliseconds on a common workstation, yielding speedups of six to seven orders of magnitude compared with running a full PIC simulation. This computational gain enables dense parameter scans and direct coupling of IEAD predictions with PMI and erosion models on whole-device scales in fusion-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural measurement of electron-phonon coupling and electronic thermal transport across a metal-semiconductor interface

Scattering of energetic charge carriers and their coupling to lattice vibrations (phonons) in dielectric materials and semiconductors are crucial processes that determine the functional limits of optoelectronics, photovoltaics, and photocatalysts. The strength of these energy exchanges is often described by the electron-phonon coupling coefficient, which is difficult to measure due to the microscopic time- and length-scales involved. In the present study, we propose an alternate means to quantify the coupling parameter along with thermal boundary resistance and electron conductivity by performing a high angular-resolution time-resolved X-ray diffraction measurement of propagating lattice deformation following laser excitation of a nanoscale, polycrystalline metal film on a semiconductor substrate. Our data present direct experimental evidence for identifying the ballistic and diffusive transport components occurring at the interface, where only the latter participates in thermal diffusion. This approach provides a robust measurement that can be applied to investigate microscopic energy transport in various solid-state materials.

36 MATERIALS SCIENCE↗

Experimental Safe Extremum Seeking for Accelerators

Here, we demonstrate the recent designs of safe extremum seeking (Safe ES) on the 1-km-long charged particle accelerator at the Los Alamos Neutron Science Center (LANSCE). Safe ES is a modification of extremum seeking (ES) which, in addition to minimizing an analytically unknown cost, also employs a safety filter based on an analytically unknown control barrier function (CBF) safety metric. Tuning is necessitated by accelerators being large complex systems, with many drifting parameters due to thermal effects and degradation. At the same time, safe operation (the maintenance of state constraints) is crucial, as damage brings astronomical costs, both financially and in operation downtime. Our measured (but analytically unknown) safety metric is the beam current. We perform multivariable Safe ES on three accelerator applications, in which we adapt 4, 6, and 3 magnet strength parameters, respectively. Two of the three applications are for validated simulation models of beamlines at LANSCE: the first for the proton radiography (pRad) beamline of 800-MeV protons for spot size tuning; the second on a high-performance code, HPSim, for tuning the low-energy beam transport (LEBT) region that contains a beam of 750-keV protons. The third is an experimental tuning of the steering magnets in the LEBT at LANSCE.

43 PARTICLE ACCELERATORS↗

Compressed Expanded Natural Graphite (CENG) Processing for PCM Composites

The use of phase change materials (PCMs) in thermal energy storage applications has received considerable attention in recent decades. Organic PCMs are popular due to their high latent heat of fusion, noncorrosive properties, and relative stability over many charge and discharge cycles. A primary limitation of these materials is their low thermal conductivity. This has led researchers to develop various methods to increase thermal conductivity by seeding PCM with or impregnating them into more conductive materials. One method is to impregnate PCMs into compressed expanded natural graphite (CENG) matrices, which can improve thermal conductivity by a factor of 100. CENG matrices have received particular interest due to their low cost, high porosity, small (nano/micro) pore size, high pore density, high thermal conductivity, and ability to be compressed into many geometries. PCM/CENG matrix composites have been extensively studied; however, the effect that CENG processing has on PCM saturation and the overall matrix thermal conductivity has not been well investigated. This processing includes four major steps including graphite intercalation, thermal shock, compression, and PCM saturation. Intercalation involves soaking graphite flakes in sulfuric and/or nitric acid to intercalate the acids between the graphene layers. The graphite flakes are then subjected to a high-temperature thermal shock, during which the intercalated acid is gasified rapidly, pushing the graphene layers apart, resulting in accordion-shaped graphite "worms". The "worms" are then compacted to a desired bulk density and then soaked with molten PCM until fully saturated. The properties of the produced CENG matrix, and its ability to allow PCM permeation, are sensitive to the processing parameters, namely, the thermal shock temperature and exposure time, as well as the matrix apparent density or porosity. Here, we study the effect of the thermal shock conditions necessary to expand intercalated graphite flakes on PCM saturation and the expanded graphite's thermal conductivity and morphology. We found that the thermal shock temperature exhibits the greatest influence. At greater shock temperatures, SEM images showed that expanded graphite worms exhibited greater density of pores, thus increasing total surface area within the matrices. Increasing thermal shock temperature yielded greater overall PCM saturation, as well as an increased rate of saturation. Improvements in PCM saturation rate and overall saturation are obtained as the shock temperature is increased. Longer exposure to thermal shock also improves initial saturation rates and is beneficial if a shortened impregnation time is needed. Thermal shock conditions did not impact thermal conductivity; however, conductivity was largely affected by matrix porosity. A local maximum in axial thermal conductivity was observed at around 83% porosity, which is similar to that observed in previous studies.

compressed expanded natural graphite↗

Dynamic Earth Energy Storage: Terawatt-year, Grid-scale Energy Storage Using Planet Earth as a Thermal Battery (GeoTES): Phase I Project (Final Report)

Grid-scale energy storage has been identified by the U.S. Department of Energy’s (DOE) Energy Storage Grand Challenge as a necessary technology to support the continued build-out of intermittent renewable energy resources required to attain a carbon-free energy future. To meet this goal, the 2018 Department of Energy Research and Innovation Act mandated the creation of a comprehensive program to accelerate the development and commercialization of next-generation energy storage technologies. One of numerous energy storage technology options is the storage of excess energy as heated geothermal brine in suitable geologic formations. This concept, known as reservoir thermal energy storage (RTES), geologic thermal energy storage (GeoTES), aquifer thermal energy storage (ATES), etc., relies on the storage of thermal energy in geologic formations for recovery and use in large-scale direct use geothermal (e.g., district heating, industrial processes, etc.) and electrical power generation applications. This thermal energy is derived from excess or waste heat from any high-temperature heat source, such as concentrated solar or from conventional thermal/nuclear generation. As such, RTES can potentially play a significant role in meeting the energy storage shortfall in the coming decades. RTES can provide energy arbitrage through both the storage and production of thermal energy stored in geologic formations for direct use applications and can serve as a source of hot fluids that can be used to generate electricity to support peak demand ramping, thus easing stress on transmission and distribution. This energy storage option has geographic benefits in that energy can be stored locally or regionally depending on the various needs/loads. RTES can also be located across an enormous geographic area, without the need for a traditional hydrothermal resource but where thermal gradients and hydrogeology allow economic exploitation of subsurface heat. The work conducted for this project includes (1) a review of lessons learned from past high-temperature RTES international projects; (2) geochemical experimental investigation and numerical simulations of potential domestic sedimentary reservoirs and (3) development of a thermo-hydrological-mechanical (THM) numerical simulation tool for optimizing formation properties and design parameters to maximize thermal energy storage performance.

15 GEOTHERMAL ENERGY↗

Effective properties of semitransparent radiative cooling materials with spectrally variable properties

Radiative cooling is a promising solid-state, non-vapor-compression technology for passive refrigeration and air conditioning. Although this phenomenon occurs naturally, achieving a significant amount of cooling to make it a technically and economically viable technology requires highly engineered, spectrally selective radiative surfaces. These characteristics make radiative cooling difficult to estimate, particularly when it is integrated with other systems such as photovoltaic panels or building envelopes. The complexity further increases when the substrate also participates in the radiative cooling (along with the radiative coating). Energy estimation is becoming increasingly critical because of the recent focus on the semitransparent radiative coatings that transmit a variety of colors to enhance the aesthetic appeal of the system. In this work, we propose an iterative method to calculate the effective radiative properties, which provides the same net radiative cooling that would be observed using the spectral properties at both the coating and substrate surfaces. Compared to traditional methods that rely on either computationally expensive full spectral analysis or methods for averaging each radiative surface parameter locally, our proposed method focuses on calculating effective properties that provide the same the net cooling effect as a full spectral analysis by accounting the emissivity, absorptivity, and transmissivity collectively, thereby providing an overall estimation error of less than 0.2%. We believe that this study will be beneficial to the engineering communities that employ complex simulation codes and require lumped solar and thermal radiation related parameters.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Hydride-based thermal energy storage

Abstract The potential and research surrounding metal hydride (MH) based thermal energy storage is discussed, focusing on next generation thermo-chemical energy storage (TCES) for concentrated solar power. The site availability model to represent the reaction mechanisms of both the forward and backward MH reaction is presented, where this model is extrapolated to a small pilot scale reactor, detailing how a TCES could function/operate in a real-world setting using a conventional shell & tube reactor approach. Further, the important parameter of effective thermal conductivity is explored using an innovative multi-scale model, to providing extensive and relevant experimental data useful for reactor and system design. Promising high temperature MH material configurations may be tuned by either destabilisation, such as using additions to Ca and Sr based hydrides, or by stabilisation, such as fluorine addition to NaH, MgH 2 , or NaMgH 3 . This versatile thermodynamic tuning is discussed, including the challenges in accurately measuring the material characteristics at elevated temperatures (500 –700 °C). Attention to scale up is explored, including generic design and prototype considerations, and an example of a novel pilot-scale pillow-plate reactor currently in development; where materials used are discussed, overall tank design scope and system integration.

25 ENERGY STORAGE↗

Role of Computational Parameters on Predicting Self-Consistent Residual Stress and Distortion during Wire Arc Additive Manufacturing

Production of three-dimensional metallic parts through integration of an articulated robot and gas metal arc welding, also known as wire arc additive manufacturing (WAAM), can produce large-scale components with moderate geometrical complexity. This technology is particularly appealing due to its high deposition rates, scalability, and cost-effective feedstock compared to other AM processes. Despite its advantages, WAAM adoption is hindered by challenges in ensuring geometric conformity without extensive distortion, defect-free structures, and consistent mechanical properties. Finite element analysis (FEA) is often used to address the challenge of geometrical conformity. As the size of parts increases, the best practices for mesh size and temporal resolution known in the literature become computationally unviable. This research examined the effects of mesh and time-step resolutions during transient FEA of a large-scale (248 layers) metallic part. The impact of computational parameters on the thermal history, displacement, and residual stress distributions were evaluated. The results showed that predicted distortion was consistent across resolutions, while time-step length significantly affected predicted thermal history, and mesh size influenced residual stress distributions. To investigate this relationship further, directionally biased meshes were considered and analyzed. The results indicated that increasing mesh resolution perpendicular to the welding path yielded stress predictions that aligned closely with higher-resolution models while offering substantial computational savings. In conclusion, the significances of this research are related to verification and validation of WAAM models for widespread industrial adoption and pragmatic guidelines for optimizing computation parameters for balancing computational efficiency and predictive accuracy of residual stress and distortion.

Solsbee, Brandon [Univ. of Tennessee, Knoxville, T↗

STAT7 v1.2 User Guide: The STAT7 Code for Statistical Propagation of Uncertainties in Steady-State Thermal Hydraulics Analysis of Plate-Fueled Reactors

The STAT7 software was developed to perform steady-state, single-phase thermal hydraulics analysis of plate-fueled reactors based on statistical propagation of uncertainties. Application of the software is for non-power research and test reactors, including conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors such as Massachusetts Institute of Technology Research Reactor. Since it can be necessary to repeat analysis during fuel reloading, STAT7 accommodates flexibility in analyzing many realistic aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model uncertainty in common fuel fabrication parameters and other key reactor operating parameters required for thermal hydraulics analyses of research and test reactors. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux, or onset of flow instability. STAT7 supports water properties based on the IAPWS-IF97 functions (The International Association for the Properties of Water and Steam Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam) in addition to the fit functions. STAT7 predicts axial profiles of fuel, cladding, and coolant temperature along a lateral stripe that runs the full length of the fuel plate from the bottom to the top. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one latera stripe of a fuel element. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each fuel element in the core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal-model-based characterization of heavy-ion-collision systems at chemical freeze-out

We investigate the chemical freeze-out in heavy-ion collisions (HICs) and the impact of the hadronic spectrum on thermal model analyses [1, 2]. Detailed knowledge of the hadronic spectrum is still an open question, which has phenomenological consequences on the study of HICs. By varying the number of resonances included in Hadron Resonance Gas (HRG) Model calculations, we can shed light on which particles may be produced. Furthermore, we study the influence of the number of states on the so-called two flavor freezeout scenario, in which strange and light particles can freeze-out separately. We consider results for the chemical freeze-out parameters obtained from thermal model fits and from calculating net-particle fluctuations. We will show the effect of using one global temperature to fit all particles and alternatively, allowing particles with and without strange quarks to freeze-out separately.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Internal short circuit and failure mechanisms of lithium-ion pouch cells under mechanical indentation abuse conditions: An experimental study

Electromechanical structural integrity and thermal stability dictate the safety performance of lithium-ion batteries. Progressive deformation and failure across microscopic and macroscopic lengths scales that are responsible for internal short circuit (ISC) in lithium-ion cells under mechanical abuse conditions remains elusive. In this study, a series of indentation tests were conducted on lithium-ion cells with different capacities up to the occurrence of ISC. The external response and internal configuration of these cells were investigated. It is discovered that cells with different capacities and state of charges exhibited different behaviors. Maximum temperature, which is often regarded as the most important parameter related to thermal runaway (TR), varied considerably due to the complicated contact configurations. X-ray computed tomography (XCT) showed that ISC was a collective result of shear band or other strain-localization modes in the electrode assembly, shear offsets in the granular coatings of electrodes, and the accompanying ductile fracture in the metal foils. We believe that the irregular strain-localization modes (kinks, cusps, and buckles), radical mismatches in mechanical properties of different layers, and geometric features of the indenter eventually lead to the tearing/puncture of cell separator at various locations. Furthermore, the results could provide useful guidance for the micromechanical modeling of lithium-ion cells.

25 ENERGY STORAGE↗

First Leptophobic Dark Matter Search from the Coherent–CAPTAIN-Mills Liquid Argon Detector

We report the first results of a search for leptophobic dark matter (DM) from the Coherent–CAPTAIN-Mills (CCM) liquid argon (LAr) detector. An engineering run with 120 photomultiplier tubes (PMTs) and 17.9 × 10 20 protons on target (POT) was performed in fall 2019 to study the characteristics of the CCM detector. The operation of this 10-ton detector was strictly light based with a threshold of 50 keV and used coherent elastic scattering off argon nuclei to detect DM. Despite only 1.5 months of accumulated luminosity, contaminated LAr, and nonoptimized shielding, CCM’s first engineering run has already achieved sensitivity to previously unexplored parameter space of light dark matter models with a baryonic vector portal. With an expected background of 115 005 events, we observe 115 005 + 16.5 events which is compatible with background expectations. For a benchmark mediator-to-DM mass ratio of m V B =m χ = 2.1, DM masses within the range 9 MeV ≲ m χ ≲ 50 MeV are excluded at 90% C. L. in the leptophobic model after applying the Feldman-Cousins test statistic. CCM’s upgraded run with 200 PMTs, filtered LAr, improved shielding, and 10 times more POT will be able to exclude the remaining thermal relic density parameter space of this model, as well as probe new parameter space of other leptophobic DM models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A combined experimental and numerical approach that eliminates the non-uniqueness associated with the Johnson-Cook parameters obtained using inverse methods

Abstract Johnson-Cook constitutive model is a commonly used material model for machining simulations. The model includes five parameters that capture the initial yield stress, strain-hardening, strain-rate hardening, and thermal softening behavior of the material. These parameters are difficult to determine using experiments since the conditions observed during machining (such as high strain-rates of the order of $$10^5$$ 10 5 /sec - $$10^6$$ 10 6 /sec) are challenging to recreate in the laboratory. To address this problem, several researchers have recently proposed inverse approaches where a combination of experiments and analytical models are used to predict the Johnson-Cook parameters. The errors between the measured cutting forces, chip thicknesses and temperatures and those predicted by analytical models are minimized and the parameters are determined. In this work, it is shown that only two of the five Johnson-Cook parameters can be determined uniquely using inverse approaches. Two different algorithms, namely, Adaptive Memory Programming for Global Optimization (AMPGO) and Particle Swarm Optimization (PSO), are used for this purpose. The extended Oxley’s model is used as the analytical tool for optimization. For determining a parameter’s value, a large range for each parameter is provided as an input to the algorithms. The algorithms converge to several different sets of values for the five Johnson-Cook parameters when all the five parameters are considered as unknown in the optimization algorithm. All of these sets, however, yield the same chip shape and cutting forces in FEM simulations. Further analyses show that only the strain-rate and thermal softening parameters can be determined uniquely and the three parameters present in the strain-hardening term of the Johnson-Cook model cannot be determined uniquely using the inverse method. A combined experimental and numerical approach is proposed to eliminate this determine all parameters uniquely.

42 ENGINEERING↗

Ab initio study of mechanical and thermal properties of GeTe-based and PbSe-based high-entropy chalcogenides

GeTe-based and PbSe-based high-entropy compounds have outstanding thermoelectric (TE) performance and crucial applications in mid and high temperatures. Recently, the optimization of TE performance of high-entropy compounds has been focused on reducing thermal conductivity by strengthening the phonon scattering process to improve TE performance. We report a first-principles investigation on nine GeTe-based high-entropy chalcogenide solid solutions constituted of eight metallic elements (Ag, Pb, Sb, Bi, Cu, Cd, Mn, and Sn) and 13 PbSe-based high-entropy chalcogenide solid solutions: Pb 0.99-y Sb 0.012 Sn y Se 1-2x Te x S x (x = 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, and y = 0) and Pb 0.99-y Sb 0.012 Sn y Se 1-2x Te x S x (y = 0.05, 0.1, 0.15, 0.2, 0.25 and x = 0.25). We have investigated the mechanical properties focusing on Debye temperature (Θ D ), thermal conductivity (κ), Grüneisen parameter (γ α ), dominant phonon wavelength (λ dom ), and melting temperature (T m ). We find that the lattice thermal conductivity is significantly reduced when GeTe is alloyed into the following compositions: Ge 0.75 Sb 0.13 Pb 0.12 Te, Ge 0.61 Ag 0.11 Sb 0.13 Pb 0.12 Bi 0.01 Te, and Ge 0.61 Ag 0.11 Sb 0.13 Pb 0.12 Mn 0.05 Bi 0.01 Te. This reduction is due to the mass increase and strain fluctuations. The results also show that Ge 0.61 Ag 0.11 Sb 0.13 Pb 0.12 Bi 0.01 Te solid solution has the lowest Young’s modulus (30.362 GPa), bulk and shear moduli (18.626 and 12.359 GPa), average sound velocity (1653.128 m/sec), Debye temperature (151.689 K), lattice thermal conductivity (0.574 W.m –1 .K –1 ), dominant phonon wavelength (0.692 Å), and melting temperature (535.91 K). Moreover, Ge 0.61 Ag 0.11 Sb 0.13 Pb 0.12 Bi 0.01 Te has the highest Grüneisen parameter with a reduced and temperature-independent lattice thermal conductivity. The positive correlation between Θ D and κ is revealed. Alloying of PbSe-based high-entropy by Sb, Sn, Te, and S atoms at the Se and Pb sites resulted in much higher shear strains resulted in the reduction of phonon velocity, a reduced Θ D , and a lower lattice thermal conductivity.

74 ATOMIC AND MOLECULAR PHYSICS↗

Observation and Modeling of the Thermal Runaway of High-Capacity Pouch Cells Due to an Internal Short Induced by an Indenter

Owing to the complexity of coupling mechanical and electrical solvers for finite element analysis (FEA) of mechanical abuse of a battery, the simulation should be performed by incorporating a local short circuit model activated based on the status of battery elements. However, the approach should also consider various other conditions of the battery simulated, such as the initial state of charge (SOC), for precise simulations of the internal short. This study proposes an approach developed with a mechanical model and single-particle battery model, along with a reaction kinetics model for thermal abuse using LS-DYNA. Several key parameters were selected, such as the activation criteria, kinetic abuse and element integration points during the internal short circuit induced by the indenter. The impact of these key parameters in the FEA on thermal runaway response for various initial SOCs from 25 to 100 % were analyzed. The FEA results obtained were compared against experimental results with various parameters such as activation criterion of internal short, kinetic abuse based on accelerating rate calorimetry (ARC) data. Finally, the root causes of the discrepancy between the FEA and experimental results are discussed regarding integration point what solver can make difference during calculation.

ENERGY STORAGE↗

Observation of Weyl exceptional rings in thermal diffusion

A non-Hermitian Weyl equation indispensably requires a three-dimensional (3D) real/synthetic space, and it is thereby perceived that a Weyl exceptional ring (WER) will not be present in thermal diffusion given its purely dissipative nature. Here, we report a recipe for establishing a 3D parameter space to imitate thermal spinor field. Two orthogonal pairs of spatiotemporally modulated advections are employed to serve as two synthetic parameter dimensions, in addition to the inherent dimension corresponding to heat exchanges. We first predict the existence of WER in our hybrid conduction–advection system and experimentally observe the WER thermal signatures verifying our theoretical prediction. When coupling two WERs of opposite topological charges, the system further exhibits surface-like and bulk topological states, manifested as stationary and continuously changing thermal processes, respectively, with good robustness. Our findings reveal the long-ignored topological nature in thermal diffusion and may empower distinct paradigms for general diffusion and dissipation controls.

42 ENGINEERING↗