Numerical modeling of multiple length scales in thermal transport processes

Author(s):  
Yogesh Jaluria

Purpose – Multiple length and time scales arise in a wide variety of practical and fundamental problems. It is important to obtain accurate and validated numerical simulation results, considering the different scales that exist, in order to predict, design and optimize the behavior of practical thermal processes and systems. The purpose of this paper is to present modeling at the different length scales and then addresses the question of coupling the different models to obtain the overall model for the system or process. Design/methodology/approach – Both numerical and experimental methods to obtain results at the different length scales, particularly at micro and nanoscales, are considered. Even though the paper focusses on length scales, multiple time scales lead to similar concerns and are also considered. The two circumstances considered in detail are multiple length scales in different domains and those in the same domain. These two cases have to be modeled quite differently in order to obtain a model for the overall process or system. The basic considerations involved in such a modeling are discussed. A wide range of thermal processes are considered and the methods that may be used are presented. The models employed must be validated and the accuracy of the simulation results established if the simulation results are to be used for prediction, control and design. Findings – Of particular interest are concerns like verification and validation, imposition of appropriate boundary conditions, and modeling of complex, multimode transport phenomena in multiple scales. Additional effects such as viscous dissipation, surface tension, buoyancy and rarefaction that could arise and complicate the modeling are discussed. Uncertainties that arise in material properties and in boundary conditions are also important in design and optimization. Large variations in the geometry and coupled multiple regions are also discussed. Research limitations/implications – The paper is largely focussed on multiple-scale considerations in thermal processes. Both numerical modeling/simulation and experimentation are considered, with the latter being used for validation and physical insight. Practical implications – Several examples from materials processing, environmental flows and electronic systems, including data centers, are given to present the different techniques that may be used to achieve the desired level of accuracy and predictability. Originality/value – Present state of the art and future needs in this interesting and challenging area are discussed, providing the impetus for further work. Different methods for treating multiscale problems are presented.

Author(s):  
G. Rossini ◽  
A. Caimi ◽  
A. Redaelli ◽  
E. Votta

AbstractA Finite Element workflow for the multiscale analysis of the aortic valve biomechanics was developed and applied to three physiological anatomies with the aim of describing the aortic valve interstitial cells biomechanical milieu in physiological conditions, capturing the effect of subject-specific and leaflet-specific anatomical features from the organ down to the cell scale. A mixed approach was used to transfer organ-scale information down to the cell-scale. Displacement data from the organ model were used to impose kinematic boundary conditions to the tissue model, while stress data from the latter were used to impose loading boundary conditions to the cell level. Peak of radial leaflet strains was correlated with leaflet extent variability at the organ scale, while circumferential leaflet strains varied over a narrow range of values regardless of leaflet extent. The dependency of leaflet biomechanics on the leaflet-specific anatomy observed at the organ length-scale is reflected, and to some extent emphasized, into the results obtained at the lower length-scales. At the tissue length-scale, the peak diastolic circumferential and radial stresses computed in the fibrosa correlated with the leaflet surface area. At the cell length-scale, the difference between the strains in two main directions, and between the respective relationships with the specific leaflet anatomy, was even more evident; cell strains in the radial direction varied over a relatively wide range ($$0.36-0.87$$ 0.36 - 0.87 ) with a strong correlation with the organ length-scale radial strain ($$R^{2}= 0.95$$ R 2 = 0.95 ); conversely, circumferential cell strains spanned a very narrow range ($$0.75-0.88$$ 0.75 - 0.88 ) showing no correlation with the circumferential strain at the organ level ($$R^{2}= 0.02$$ R 2 = 0.02 ). Within the proposed simulation framework, being able to account for the actual anatomical features of the aortic valve leaflets allowed to gain insight into their effect on the structural mechanics of the leaflets at all length-scales, down to the cell scale.


Author(s):  
David O. Kazmer ◽  
Stephen P. Johnston ◽  
Mary E. Moriarty ◽  
Christopher Santeufemio

Methods are presented for self-alignment and assembly of objects with micron and nanometer-level features. The approach is a combination of kinematic coupling and elastic averaging in which mating alignment features spanning multiple length scales are successively brought into contact. When the objects are pressed together, the larger alignment features cause necessary deformation to ensure adequate alignment at the smaller length scales. Analytical and numerical modeling indicate that the largest alignment features can be designed to generally resolve global systematic errors while the smaller alignment features can correct local errors to achieve sub-micron alignment. Physical realization with ion beam etching, deposition, and thermal imprint lithography are also discussed.


2020 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Yogesh Jaluria

Purpose This paper aims to discuss inverse problems that arise in a variety of practical thermal processes and systems. It presents some of the approaches that may be used to obtain results that lie within a small region of uncertainty. Therefore, the non-uniqueness of the solution is reduced so that the final design and boundary conditions may be determined. Optimization methods that may be used to reduce the uncertainty and to select locations for experimental data and for minimizing the error are presented. A few examples of thermal systems are given to illustrate the applicability of these methods and the challenges that must be addressed in solving inverse problems. Design/methodology/approach In most analytical and numerical solutions, the basic equations that describe the process, as well as the relevant and appropriate boundary conditions, are known. The interest lies in obtaining a unique solution that satisfies the equations and boundary conditions. This may be termed as a direct or forward solution. However, there are many problems, particularly in practical systems, where the desired result is known but the conditions needed for achieving it are not known. These are generally known as inverse problems. In manufacturing, for instance, the temperature variation to which a component must be subjected to obtain desired characteristics is prescribed, but the means to achieve this variation are not known. An example of this circumstance is the annealing, tempering or hardening of steel. In such cases, the boundary and initial conditions are not known and must be determined by solving the inverse problem to obtain the desired temperature variation in the component. The solutions, thus, obtained are generally not unique. This is a review paper, which discusses inverse problems that arise in a variety of practical thermal processes and systems. It presents some of the approaches or strategies that may be used to obtain results that lie within a small region of uncertainty. It is important to realize that the solution is not unique, and this non-uniqueness must be reduced so that the final design and boundary conditions may be determined with acceptable accuracy and repeatability. Optimization techniques are often used for minimizing the error. This review presents several methods that may be applied to reduce the uncertainty and to select locations for experimental data for the best results. A few examples of thermal systems are given to illustrate the applicability of these methods and the challenges that must be addressed in solving inverse problems. By considering a variety of systems, the paper also shows the importance of solving inverse problems to obtain results that may be used to model and design thermal processes and systems. Findings The solution of inverse problems, which frequently arise in thermal processes, is discussed. Different strategies to obtain the conditions that lead to the desired result are given. The goal of these approaches is to reduce uncertainty and obtain essentially unique solutions for different circumstances. The error of the method can be checked against known conditions to see if it is acceptable for the given problem. Several examples are given to illustrate the use of these methods. Originality/value The basic strategies presented here for solving inverse problems that arise in thermal processes and systems, as well as the optimization techniques used to reduce the domain of uncertainty, are fairly original. They are used for certain challenging problems that have not been considered in detail earlier. Several methods are outlined for considering different types of problems.


Author(s):  
Iman Rashidi ◽  
Lioua Kolsi ◽  
Goodarz Ahmadi ◽  
Omid Mahian ◽  
Somchai Wongwises ◽  
...  

Purpose This study aims to investigate a three-dimensional computational modelling of free convection of Al2O3 water-based nanofluid in a cylindrical cavity under heterogeneous heat fluxes that can be used as a thermal storage tank. Design/methodology/approach Effects of different heat flux boundary conditions on heat transfer and entropy generation were examined and the optimal configuration was identified. The simulation results for nanoparticle (NP) volume fractions up to 4 per cent, and Rayleigh numbers of 104, 105 and 106 were presented. Findings The results showed that for low Ra (104) the heat transfer and entropy generation patterns were symmetric, whereas with increasing the Rayleigh number these patterns became asymmetric and more complex. Therefore, despite the symmetric boundary conditions imposed on the periphery of the enclosure (uniform in Ɵ), it was necessary to simulate the problem as three-dimensional instead of two-dimensional. The simulation results showed that by selecting the optimal values of heat flux distribution and NP volume fraction for these systems the energy consumption can be reduced, and consequently, the energy efficiency can be ameliorated. Originality/value The results of the present study can be used for the design of energy devices such as thermal storage tanks, as both first and second laws of thermodynamics have been considered. Using the optimal design will reduce energy consumption.


2012 ◽  
Vol 134 (3) ◽  
Author(s):  
Yogendra Joshi

Thermal systems often involve multiple spatial and temporal scales, where transport information from one scale is relevant at others. Optimized thermal design of such systems and their control require approaches for their rapid simulation. These activities are of increasing significance due to the need for energy efficiency in the operation of these systems. Traditional full-field simulation methodologies are typically unable to resolve these scales in a computationally efficient manner. We summarize recent work on simulations of conjugate transport processes over multiple length scales via reduced order modeling through approaches such as compact finite elements and proper orthogonal decomposition. In order to incorporate the influence of length scales beyond those explicitly considered, lumped models are invoked, with appropriate handshaking between the two frameworks. We illustrate the methodology through selected examples, with a focus on information technology systems.


2021 ◽  
Author(s):  
Dietmar Hutmacher ◽  
Benjamin Gorissen ◽  
Simon Liponsky ◽  
Katia Bertoldi ◽  
Tara Shabab ◽  
...  

Abstract The quest for an advanced soft robotic actuator technology that is fast and can execute a wide range of application-specific tasks at multiple length scales is still ongoing. Here, we demonstrate a new design strategy leveraging the concepts of miniaturisation and fibre-reinforcement to realize high-speed inflatable actuators exhibiting diverse movements. To fabricate the designs, we employ a class of additive manufacturing technology called melt electrowriting. We demonstrate 3D printing of microfibre architectures on soft elastomers with precision at unprecedently small length scales, leading to miniaturised composite actuators with highly controlled deformation characteristics. We show that owing to their small dimensions and deterministically designed fibrous networks, our actuators require extremely low amounts of fluid to inflate. We demonstrate that actuators with a length of 10 to 15mm and an inner diameter 1mm can reach their full range of motion within ~ 20ms without exploiting snapping instabilities or material non-linearities. We display the speed of our actuators by building an ultrafast, soft flycatcher.


Geosphere ◽  
2021 ◽  
Author(s):  
Samuel Angiboust ◽  
Armel Menant ◽  
Taras Gerya ◽  
Onno Oncken

Several decades of field, geophysical, analogue, and numerical modeling investigations have enabled documentation of the wide range of tectonic transport processes in accretionary wedges, which constitute some of the most dynamic plate boundary environments on Earth. Active convergent margins can exhibit basal accretion (via underplating) leading to the formation of variably thick duplex structures or tectonic erosion, the latter known to lead to the consumption of the previously accreted material and eventually the forearc continental crust. We herein review natural examples of actively underplating systems (with a focus on circum-Pacific settings) as well as field examples highlighting internal wedge dynamics recorded by fossil accretionary systems. Duplex formation in deep paleo–accretionary systems is known to leave in the rock record (1) diagnostic macro- and microscopic deformation patterns as well as (2) large-scale geochronological characteristics such as the downstepping of deformation and metamorphic ages. Zircon detrital ages have also proved to be a powerful approach to deciphering tectonic transport in ancient active margins. Yet, fundamental questions remain in order to understand the interplay of forces at the origin of mass transfer and crustal recycling in deep accretionary systems. We address these questions by presenting a suite of two-dimensional thermo-mechanical experiments that enable unravelling the mass-flow pathways and the long-term distribution of stresses along and above the subduction interface as well as investigating the importance of parameters such as fluids and slab roughness. These results suggest the dynamical instability of fluid-bearing accretionary systems causes either an episodic or a periodic character of subduction erosion and accretion processes as well as their topographic expression. The instability can be partly deciphered through metamorphic and strain records, thus explaining the relative scarcity of paleo–accretionary systems worldwide despite the tremendous amounts of material buried by the subduction process over time scales of tens or hundreds of millions of years. We finally stress that the understanding of the physical processes at the origin of underplating processes as well as the forearc topographic response paves the way for refining our vision of long-term plate-interface coupling as well as the rheological behavior of the seismogenic zone in active subduction settings.


SPE Journal ◽  
2011 ◽  
Vol 16 (03) ◽  
pp. 713-724 ◽  
Author(s):  
Juliana Y. Leung ◽  
Sanjay Srinivasan

Summary Reservoir heterogeneities occur over a wide range of length scales, and their interaction with various transport mechanisms controls the performance of subsurface flow and transport processes. Modeling these processes at large scales requires proper scaleup of petrophysical properties that are autocorrelated or heterogeneously distributed in space, and analyzing their interaction with underlying transport mechanisms. A method is proposed to investigate and quantify the uncertainty in reservoir models introduced by scaleup. It is demonstrated that when the volume support of the measurement is smaller than the representative elementary volume (REV) scale of the attribute to be modeled, there is uncertainty in the conditioning data because of scaleup and that uncertainty has to be propagated to spatial models for the attribute. This important consideration is demonstrated for mapping total porosity for a carbonate reservoir in the Gulf of Mexico. The results demonstrate that in most cases, the uncertainty distributions obtained by accounting for the scaleup procedure successfully characterize the variability in the actual core and log data observed along new wells. Conventional reservoir models considering the well data as "hard" conditioning data fail to predict the "true" values. Following this discussion on scaling of reservoir attributes, a conceptual understanding of the scaling characteristics of flow responses such as recovery factor (RF) is provided, in terms of the mean and variance of RF at different length scales. Finally, a new technique is presented to systematically quantify the scaling characteristics of transport processes by accounting for subscale heterogeneities and their interaction with various transport mechanisms based on the volume averaging approach. The objective is to provide a tool for understanding the scaling relationships for RF using detailed fine-scale compositional reservoir simulations over a subdomain of the reservoir.


2020 ◽  
Vol 37 (8) ◽  
pp. 2761-2783
Author(s):  
Łukasz Łach ◽  
Dmytro Svyetlichnyy

Purpose Some functional properties of engineering materials, i.e. physical, mechanical and thermal ones, depend directly on the microstructure, which is a result of processes occurring in the material during the forming and thermomechanical processing. The proper microstructure can be obtained in many cases by the phase transformation. This phenomenon is one of the most important processes during hot forming and heat treatment. The purpose of this paper is to develop a new comprehensive hybrid model for modeling diffusion phase transformations. A problem has been divided into several tasks and is carried out on several stages. The purpose of this stage is a development of the structure of a hybrid model, development of an algorithm used in the diffusion module and one-dimensional heat flow and diffusion modeling. Generally, the processes of phase transformations are studied well enough but there are not many tools for their complex simulations. The problems of phase transformation simulation are related to the proper consideration of diffusion, movement of phase boundaries and kinetics of transformation. The proposed new model at the final stage of development will take into account the varying grain growth rate, different shape of growing grains and will allow for proper modeling of heat flow and carbon diffusion during the transformation in many processes, where heating, annealing and cooling can be considered (e.g. homogenizing and normalizing). Design/methodology/approach One of the most suitable methods for modeling of microstructure evolution during the phase transformation is cellular automata (CA), while lattice Boltzmann method (LBM) suits for modeling of diffusion and heat flow. Then, the proposed new hybrid model is based on CA and LBM methods and uses high performing parallel computations. Findings The first simulation results obtained for one-dimensional modeling confirm the correctness of interaction between LBM and CA in common numerical solution and the possibility of using these methods for modeling of phase transformations. The advantages of the LBM method can be used for the simulation of heat flow and diffusion during the transformation taking into account the results obtained from the simulations. LBM creates completely new possibilities for modeling of phase transformations in combination with CA. Practical implications The studies are focused on diffusion phase transformations in solid state in condition of low cooling rate (e.g. transformation of austenite into ferrite and pearlite) and during the heating and annealing (e.g. transformation of the ferrite-pearlite structure into austenite, the alignment of carbon concentration in austenite and growth of austenite grains) in carbon steels within a wide range of carbon content. The paper presents the comprehensive modeling system, which can operate with the technological processes with phase transformation during heating, annealing or cooling. Originality/value A brief review of the modeling of phase transformations and a description of the structure of a new CA and LBM hybrid model and its modules are presented in the paper. In the first stage of model implementation, the one-dimensional LBM model of diffusion and heat flow was developed. The examples of simulation results for several variants of modeling with different boundary conditions are shown.


Sign in / Sign up

Export Citation Format

Share Document