Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Contact”

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,117 records · Page 62

An interface-aware sub-scale dynamics multi-material cell model for solids with void closure and opening at all speeds

Here, we present a multi-material cell model (closure model) for demanding arbitrary Lagrangian-Eulerian (ALE) simulations of fluids and solids. It is based on the interface-aware sub-scale dynamics (IASSD) approach which utilizes the exact material interface geometry within the computational cell to calculate internal material interactions. Our formulation of the closure model also aims to improve the accuracy in low-speed impact events. Voids are used to represent ambient vacuum and internal free boundaries of the distinct materials. Void regions can close and open at contact surfaces, allowing a transition from contact physics to free motion in vacuum. The coupling of void closure and opening with a new formulation of the IASSD model for solids is tested on several one- and two-dimensional numerical examples, ranging from gas expansion in vacuum to planar and round object impacts at various speeds.

42 ENGINEERING↗

Predictive dynamic wetting, fluid–structure interaction simulations for braze run-out

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. Here, the quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

36 MATERIALS SCIENCE↗

Influence of Initial fabric and water-wetting on particle fracture and force chain evolution in natural silica sand

Particle fracture has a significant influence on the engineering behavior of granular materials. However, the combined influence of the initial fabric and wetting on particle fracture and force chain evolution in granular materials, particularly in sands, remains insufficiently explored in the current literature. Here, in this paper, one-dimensional (1D) confined compression experiments were conducted on dry and wet specimens, and particle-level fracture was captured using three-dimensional (3D) in-situ Synchrotron Micro-Computed Tomography (SMT). A quantitative assessment of particle fracture revealed that wet specimens, which have disturbed fabric, exhibited a higher fracture percentage in comparison to dry specimens. Complementary 3D finite element (FE) simulations were performed using dry and wet (only solid fabric was considered and the influence of pore water is not considered) sand assemblies to assess interparticle contact forces and particle-scale stress distributions within the specimens. Force chain analysis was subsequently conducted, encompassing quantification of interparticle forces, characterization of force network structures, and monitoring the dynamic evolution of force chains under different strain levels. The results show that the specimens with disturbed fabric led to a more dynamic and less persistent force network, more fabric instability, and thus more reorganization of force chain structures. In addition, the frequent rearrangement of the force network in the presence of water (with reduced inter-particle contact friction) likely exacerbates localized stress concentrations, promoting failure in previously unengaged or weakly connected particles. The results reported in this paper offer a new insight into how the initial fabric and wetting cause different fracture behavior. The findings can also pave the way for more in-depth future investigations into the mechanics governing particle fracture in granular assemblies.

Finite element analysis↗

Revisiting the π → π* transition of the nitrite ion at the air/water interface: A combined experimental and theoretical study

Broadband deep ultraviolet electronic sum frequency generation spectroscopy allows measurement of electronic |χ (2) | 2 -spectra at aqueous interfaces, providing considerable improvement over deep ultraviolet electronic second harmonic generation spectroscopy that yields error-prone, pointwise spectra. Reexamination of the π → π* transition of nitrite at the air/water interface reveals that the interfacial |χ (2) | 2 -spectrum is strikingly similar to the bulk absorption spectrum. Molecular dynamics simulations and SOS-CIS(D 0 ) electronic structure calculations provide no evidence for the previously reported contact ion pair-induced red-shift of the π → π* transition at the air/water interface. Finally, our results thus revise the previous description for nitrite interfacial adsorption as a contact ion pair.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Superhydrophobic composite asymmetric electrospun membrane for sustainable vacuum assisted air gap membrane distillation

Membrane distillation (MD) relies on superhydrophobic membranes that must achieve the competing goals of high liquid entry pressure (LEP) and high permeability. Here, in this study, we developed the first superhydrophobic composite electrospun asymmetric membrane using fluorinated graphite (FG). We also provided among the first comparisons of the relative benefits of membrane asymmetry and composite coatings via 3D electrospinning and spraying. FG particles increased the roughness and reduced the zeta potential of the membranes. Superhydrophobic composite membranes had a contact angle of 157° (sliding angle of 7°) and an LEP of 327 kPa. Long-term MD testing revealed that FG particles could increase MD's sustainability while decreasing membrane wetting. Vacuum-assisted air gap membrane distillation (V-AGMD) showed that the composite asymmetric membrane outperformed the non-asymmetric version (composite nanofiber) membrane by 21 % in average water flux and 25 % in average thermal efficiency because of the asymmetric structure of the membrane. The numerical model developed for the V-AGMD system also successfully predicted the water flux generated experimentally from the composite nanofiber and composite asymmetric membranes. Our results confirmed that electrospun asymmetric FG membranes with high contact angle, low sliding angle, high permeability, and high gained output ratio (GOR) would be promising candidates for sustainable MD.

36 MATERIALS SCIENCE↗

Interface stability of LiCl-rich argyrodite Li 6 PS 5 Cl with propylene carbonate boosts high-performance lithium batteries

In solid-state lithium (Li)-metal batteries (SSLMBs), sulfide-based Li-ion conductors represent one of the most popular solid electrolytes (SEs). However, the development of sulfide-based SSLMBs is significantly hampered by interfacial issues including large solid-solid contact resistance and serious side reactions at the SE/electrode interface. To address these issues, here, we demonstrated a simple and efficient strategy by using LiCl-rich argyrodite Li 6 PS 5 Cl (Li 6 PS 5 Cl-LiCl) SE and trace amount of propylene carbonate (PC) at the SE/electrode interface to facilitate the formation of stable and robust solid electrolyte interphase (SEI). Additionally, the formed SEI not only serves as a buffer layer to passivate the interfacial reactions and suppress Li dendrite growth, but also acts as a bridge for Li-ion conduction to reduce the contact resistance. As a result, the Li||Li symmetric cell exhibited long-term electrochemical cycling stability over 1000 h at a current density of 0.2 mA cm -2 . Furthermore, the assembled Li||Li 4 Ti 5 O 12 (LTO) batteries delivered a specific capacity of 175 mAh g–1 at 0.2 C, and remained an excellent specific capacity of 116 mAh g –1 after 200 cycles at a high current rate of 1 C. These features indicate a feasible strategy to enhance the interfacial property of high-performance SSLMBs.

25 ENERGY STORAGE↗

Morphology engineering of iridium electrodes via modifying titanium substrates with controllable pillar structures for highly efficient oxygen evolution reaction

Nowadays, Ti is the well-chosen anode substrate material for proton exchange membrane electrolyzer cells (PEMECs) and modifications of the substrate surfaces are essential for the fabrication of highly efficient electrodes. Herein, we introduce the morphology engineering of Ir/Ti electrodes with different acid treatments of hydrochloric acid (HCl) and oxalic acid (OA), and the comparative benefits of these two acid treatment methods are studied from the aspects of their impacts on the morphology, interfacial contact resistance (ICR), and oxygen evolution reaction (OER) performances of resultant electrodes. Notably, compared to the flat surface from OA treatment, Ti substrates with the pillar structure could be successfully achieved via HCl etching. Further, the HCl and the oxalic acid (OA) treatments would reduce the interfacial contact resistance (ICR) to 15.2% and 5.5% of the pristine Ti substrate at 1.38 MPa, respectively. By iridium electrodeposition on Ti substrates, Ir/Ti electrodes with different catalyst loadings are fabricated. With similar low-loadings of about 0.05 mgIr cm –2 , an Ir/Ti electrode with HCl treated substrate exhibits a lower overpotential of ~283 mV at 10 mA cm –2 current density than 305 mV from OA treatment, due to the boosted reaction areas from HCl treatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic interactions and stability issues at the copper-graphene interface in air and in alkaline solution under electrochemical control

A micro-electrochemical cell is sealed with a polymer-free single-layer graphene (SLG) membrane to monitor the stability of Cu nanoparticles (NPs) attached to SLG, as well as the interfacial electronic interactions between Cu NPs and SLG both in air and in a mildly alkaline aqueous solution under electrochemical control. A combination of techniques, including in-situ Kelvin probe force microscopy (KPFM) and ex-situ electron microscopy, are applied. When Cu NPs are metallic at cathodic potentials, there is a relatively bias-independent offset in the SLG work function due to charge transfer at the Cu-SLG contact. When Cu NPs are oxidized at anodic potentials, on the other hand, the work function of SLG also depends on the applied bias in a quasi-linear fashion due to electrochemical gating, in addition to charge transfer at the CuOx-SLG contact. Furthermore, Cu NPs were found to oxidize and detach from SLG when kept under anodic potentials for a few hours, whereas they remain adhered to SLG at cathodic potentials. This is attributed to water intercalation at the CuO-SLG interface associated with the enhanced hydrophilicity of positively polarized graphene, as supported by the absence of Cu detachment following oxidation by galvanic corrosion in air.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extremized nonlinear and linearized responses in soft metamaterials enabled by gradient-based design and grayscale digital light processing

In this study, we develop a gradient-based design approach that exploits grayscale digital light processing (DLP) 3D printing for extremizing the nonlinear and linearized response of soft metamaterials — materials that harness engineered geometric instabilities to undergo large and programmable changes in configuration. Grayscale DLP approaches modulate local mechanical properties at the pixel scale by tuning the light intensity within a single grayscale image, unlocking an exceptionally large design space. To effectively navigate this space, we develop smooth mappings between local light intensity values and global quantities of interest that characterize the behavior of soft metamaterials. Enabling these smooth mappings are robust and differentiable nonlinear finite element simulations powered by a trust region solver. A PDE-constrained optimization problem is then solved to invert these mappings and produce light intensity distributions that endow the printed part with varying stiffness and flexibility in distinctive regions. It is shown that optimizing the distribution of soft and stiff phases throughout a metamaterial structure results in markedly different buckling and self-contact configurations to drive extremized nonlinear compression and linearized vibration responses. Optimized light intensity distributions are translated to grayscale images and directly used to print soft metamaterial samples, showing remarkable agreement between the buckling and self-contact response in simulated and measured deformed configurations.

Additive manufacturing↗

Unveiling stimulation fluid-driven alterations in shale pore architecture through combined interpretation of TD-NMR and multi-component gas adsorption

The hydraulic fracturing stimulation process can alter petrophysical and flow properties of the shale matrix in unconventional reservoirs. These alterations can result from interactions with native and exogenous fluids, and can be detrimental to reservoir productivity. Different experimental procedures have been proposed to study reactive processes and thereby understand the aforementioned interactions. Most published works do not consider the rock initial saturation with formation fluid, or the effect of organic matter (OM). To gain insights into these effects, two reactive experimental sets, differing in experimental conditions, were conducted. In the first set, at 125°C and 45 MPa, a stimulation fluid surrogate was put in contact with Baxter Shale samples, previously saturated and equilibrated with synthetic formation fluid. A baseline consisted of exposure to formation fluid exclusively. In the second set, at 0.08 MPa and 25°C, contact with stimulation fluid took place as in the first set. Two additional assays in this second set either skipped the initial formation fluid saturation or removed the OM, respectively. Petrophysical properties analyses, focused on the pore architecture, relied on the combined interpretation of gas adsorption and Time-Domain Nuclear Magnetic Resonance (TD-NMR). Results show that the presence of the formation fluid significantly restricts pore space accessibility to the stimulation fluid and mitigates the pore-architecture alteration to the macro-pore region. Mineral dissolution and precipitation dominate over surface or wettability alteration. Next, results suggest that the OM restricts accessibility to some pores. Ultimately, the equilibration process with synthetic formation fluid drives mainly surface alteration, unhindered by OM.

58 GEOSCIENCES↗

Radiological analysis and transmutation calculation of representative castable nanostructured alloys

In this work, radiological analysis regarding the waste management was performed on two current reduced-activation ferritic-martensitic (RAFM) steels Eurofer 97 and F82H and two castable nanostructured alloys (CNAs) CNA1 and CNA3 using the European DEMO first wall spectrum. CNA1 and CNA3 are a new generation of RAFM steels, which have increased amounts of nanoscale carbonitrides and carbides precipitates, respectively, via small composition modifications. Furthermore, the increased nitrogen in CNA1 leads to higher activity at ≥50 years, and thus to the disposal concern under the current NRC regulation. However, the calculation using Fetter’s evaluation leads to reduced waste disposal rating (WDR) of down to 0.25 depending on the damage dose from ~0.05 wt percent nitrogen in CNA1, and thus less concern on the waste disposal. Contact dose rate after 10 years from the first wall replacement primarily depends on the initial Co and Nb impurity levels. Other added elements, such as V and Si in CNA1 and Ti and Mn in CNA3, do not have significant effect on the specific activity, contact dose rate, or decay heat. Composition limits of CNAs are given based on the shallow burial disposal and materials recycling. In addition, composition evolution of typical precipitates in the CNAs was calculated using the DEMO first wall and the HFIR spectra. Larger transmutation-induced composition evolution occurs in the V/Ta-rich carbonitrides in CNA1 than that in the Ti-rich carbides in CNA3.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extraction of lithium from battery recycling wastewater using synergistic D2EHPA and TBP

The recycling of spent lithium-ion batteries (LIBs) poses significant challenges, including the generation of large volumes of chemically complex wastewater. The composition of this wastewater is influenced by both the intrinsic chemistry of the batteries and the specific recycling processes employed. Notably, this wastewater contains economically valuable components, such as lithium, which can be recovered. Here, in this study, a solvent extraction (SX) process was investigated as a method to recover lithium from battery recycling wastewater, especially from anode-washing stage. Initially, various commercial extractants were evaluated, including di(2-ethylhexyl)phosphoric acid (D2EHPA), mono-2-ethylhexyl (2-ethylhexyl)phosphonate (PC88A), bis(2,4,4-trimethylpentyl)phosphinic acid (Cyanex 272), 2-hydroxy-5-nonylacetophenone oxime (LIX 84-I), tri-butyl phosphate (TBP), and their combinations. Among these, D2EHPA + TBP demonstrated a synergism to advance and maximize the lithium extraction. Subsequently, the effects of key parameters, including D2EHPA concentration, TBP concentration, contact time, initial pH, and aqueous-to-organic (A/O) phase ratio, were systematically investigated and optimized. A two-stage SX approach was employed to enhance lithium recovery. Under the optimized conditions of 30.0 vol% D2EHPA, 10.0 vol% TBP, 10 min of contact time, and a 1:1 A/O phase ratio, a lithium extraction efficiency of more than 88% in a two-stage solvent extraction was achieved. Lithium was subsequently stripped from the loaded organic solution using sulfuric acid (H 2 SO 4 ). Using 2.0 M H 2 SO 4 , lithium stripping was achieved after two counter-current stripping stages at an organic-to-aqueous (O/A) phase ratio of 6:1. This stripping process enriched the lithium concentration by a factor of four compared to the original lithium concentration in the anode-washing wastewater. The recyclability of the synergistic D2EHPA + TBP system was also evaluated over four extraction-stripping cycles. The results demonstrated that the system maintained high extraction and stripping efficiencies.

D2EHPA↗

Creative absorber design and optimization for CO 2 capture with aqueous piperazine

The Texas Carbon Management Program is preparing a Front-End Engineering Design (FEED) of the PZASTM (Piperazine with the Advanced Stripper TM ) for CO 2 capture at a natural gas combined cycle power plant in West Texas. This paper develops the absorber design for the capture system and studies the operational strategies of the absorber in response to ambient temperature swing. The design uses 2 absorber trains, and each train integrates absorption and a water-wash section in one column. Each column has a cross-section area of 174 m 2 and consists of 3 sections of MellapakPlus® 252Y structured packing. Packing heights for each section are 3, 4.9, and 2.7 m from top to bottom. A pump-around intercooler applied at the bottom bed circulates solvent and eliminates the direct contact cooler. The baseline CO 2 removal rate is 90% and will vary between 78% and 90% with ambient temperature. The system uses air cooling to eliminate cooling water usage. The water balance in the system will be managed by operating the water wash at a high circulation rate (1300 kg/h) to condense water during the night and make up for water loss during the day. Finally, the proposed design reduces the capital cost of the absorber column by minimizing the packing height and eliminating the direct contact cooler (DCC), and improves the energy cost by lowering the pump-around intercooling temperature and increasing the pump-around rate.

03 NATURAL GAS↗

Unified descriptor for enhanced critical heat flux during pool boiling of hemi-wicking surfaces

Boiling heat transfer is dictated by interfacial phenomena at the three-phase contact line where vapor bubbles form on the surface. Structured surfaces have shown significant enhancement in critical heat flux (CHF) during pool boiling by tailoring interfacial phenomena. This CHF enhancement has been primarily explained by two structural effects: roughness, which extends the contact line length at the bubble base, and wickability, the ability to imbibe liquid through surface structures by capillary pumping. Here, in this work, we show that CHF enhancement on structured surfaces cannot be described by roughness or wickability alone. This result was confirmed using systematically designed micropillar surfaces with controlled roughness and wickability. Further, we performed a scaling analysis and derived a unified descriptor, which represents the combined effects of thin film density and volumetric wicking rate. This unified descriptor shows a reasonable correlation with CHF values with our experiments and literature data. This work provides important insights in understanding the role of surface structures on CHF enhancement, thereby providing guidelines for the systematic design of surface structures for enhanced pool boiling heat transfer.

42 ENGINEERING↗

Droplet Evaporation on Hot Micro-Structured Superhydrophobic Surfaces: Analysis of Evaporation from Droplet Cap and Base Surfaces

In this study, evaporation of sessile water droplets on hot micro-structured superhydrophobic surfaces is experimentally and theoretically investigated. Water droplets of 4 µL are placed on micro-pillared silicon substrates with the substrate temperature heated up to 120°C. A comprehensive thermal circuit model is developed to analyze the effects of substrate roughness and substrate temperature on the sessile droplet evaporation. For the first time, two components of heat and mass transfer, i.e., one from the droplet cap surface and the other from the droplet base surface, during droplet evaporation are distinguished and systematically studied. As such, the evaporation heat transfer rates from both the droplet cap surface and the interstitial liquid-vapor interface between micropillars at the droplet base are calculated in various conditions. For droplet evaporation on the heated substrates in the range of 40°C – 80°C, the predicted droplet cap temperature matches well with the experimental results. Furthermore, during the constant contact radius mode of droplet evaporation, the decrease of evaporation rate from the droplet base contributes most to the continuously decreasing overall evaporation heat transfer rate, whereas the decrease of evaporation rate from the droplet cap surface is dominant in the constant contact angle mode. The influence of internal fluid flow is considered for droplet evaporation on substrates heated above 100°C, and an effective thermal conductivity is adopted as a correction factor to account for the effect of convection heat transfer inside the droplet. Temperature differences between the droplet base and the substrate base are estimated to be about 2°C, 5°C, 8°C, 13°C and 18°C for droplet evaporation on substrates heated at 40°C, 60°C, 80°C, 100°C, and 120°C, respectively, elucidating the delayed or depressed boiling of water droplets on a heated rough surface due to evaporative cooling.

42 ENGINEERING↗

Robust silane self-assembled monolayer coatings on plasma-engineered copper surfaces promoting dropwise condensation

Dropwise condensation is well known to result in better heat transfer performance owing to efficient condensate/droplet removal, which can be harnessed in various industrial heat/mass transfer applications such as power generation and conversion, water harvesting/desalination, and electronics thermal management. The key to enhancing condensation via the dropwise mode is thin low surface energy coatings (<100 nm) with low contact angle hysteresis. Ultrathin (<5 nm) silane self-assembled monolayers (or SAMs) have been widely studied to promote dropwise condensation due to their minimal thermal resistance and scalable integration processes. Such thin coatings typically degrade within an hour during condensation of water vapor. After coating failure, water vapor condensation transitions to the inefficient filmwise mode with poor heat transfer performance. We enhance silane SAM quality and durability during water vapor condensation on copper compared to state-of-the-art silane coatings on metal surfaces. We achieve this via (i) surface polishing to sub-10 nm levels, (ii) pure oxygen plasma surface treatment, and (iii) silane coating integration with the copper substrate in an anhydrous/moisture-free environment. The resulting silane SAM has low contact angle hysteresis (≈20°) and promotes efficient dropwise condensation of water for >360 hours without any visible sign of coating failure/degradation in the absence of non-condensable gases. We further demonstrate enhanced heat transfer performance (≈5-7× increase over filmwise condensation) over an extended period of time. Surface characterization data post-condensation leads us to propose that in the absence of non-condensable gases in the vapor environment, the silane SAM degrades due to reduction and subsequent dissolution of copper oxide at the oligomer-substrate interface. The experiments also indicate that the magnitude of surface subcooling (or condensation rate) affects the rate of coating degradation. Finally, this work identifies a pathway to durable dropwise promoter coatings that will enable efficient heat transfer in industrial applications.

42 ENGINEERING↗

Thermal conductivity measurement using modulated photothermal radiometry for nitrate and chloride molten salts

Molten salts are being used or explored for thermal energy storage and conversion systems in concentrating solar power and nuclear power plants. Thermal conductivity of molten salts is an important thermophysical property dictating the performance and cost of these systems, but its accurate measurement has been challenging, as evidenced by wide scattering of existing data in literature. The corrosive and conducting nature of these fluids also leads to time consuming sample preparation processes of many contact-based measurements. Here, in this work, we report the measurement of thermal conductivity of molten salts using a modulated photothermal radiometry (MPR) technique, which is a laser-based, non-contact, frequency-domain method adopted for molten salts for the first time. By unitizing the advantages of front side sensing of frequency-domain measurements and the vertical holder orientation, the technique can minimize the natural convection and salt creeping effects, thus yielding accurate molten salt thermal conductivity . The MPR technique is first calibrated using standard molten materials including paraffin wax and sulfur. It is then applied on measuring pure nitrate salts (NaNO 3 and KNO 3 ), solar salt (NaNO 3 –KNO 3 mixture), and chloride salt (NaCl–KCl–MgCl 2 ). The measurement results are compared with data from literature, especially those obtained from laser flash analysis (LFA). Our results demonstrate that the MPR is a convenient and reliable technique of measuring thermal conductivity of molten salts. Accurate thermal conductivity data of molten salts will be valuable in developing the next-generation high-temperature thermal energy storage and conversion systems.

14 SOLAR ENERGY↗

Alternating conduction and convection drying of paper – an experimental analysis with a continuous data acquisition approach

In conventional multi-cylinder drying of paper and board, both conductive drying from steam-heated dryer cylinders and convective drying by flowing air over the paper surface in the pockets are used. Conductive drying from steam-heated drying cylinders is a critical component in providing the necessary thermal energy to paper and board as they dry. Steam temperature and internal and external resistances at the contacting surface are critical process parameters influencing the conductive drying process. An experimental setup was developed to study the alternating conductive and convective drying of paper and board. Paper sheet moisture, temperature, and temperature distribution within the heated platen and the instantaneous heat flux as the sheet was being dried were measured. The instantaneous heat flux, contact heat transfer coefficient, and drying rates were determined as drying proceeds. Experimental results, as well as comparisons to literature and commercial data, are presented. The conductive heat transfer coefficients determined were compared to traditional correlations normally used in the modeling of paper drying. Similarly, the convective heat and mass transfer coefficients are also determined and compared to literature data. In addition to the evaluation of alternating conductive and convective drying characteristics of paper and board, the potential inclusion of auxiliary energy components will also be included. Experimental results from the conduction and convection drying system are presented. Furthermore, this data will be useful in process development, intensification of manufacturing processes, and modeling and simulation of paper drying processes.

42 ENGINEERING↗