Abstract: Laminar natural convection in a cylindrical annular cavity filled with air and provided with two fins is studied numerically using the discretization of the governing equations with the Centered Finite Difference method based on the Alternating Direction Implicit (ADI) scheme. The fins are attached to the inner cylinder of radius ri (hot wall of temperature Ti). The outer cylinder of radius ro is maintained at a temperature To (To < Ti). Two values of the dimensionless thickness of the fins are considered: 0.015 and 0.203. We consider a low fin height equal to 0.078 and medium fin heights equal to 0.093 and 0.203. The position of the fin is 0.82π and the radius ratio is equal to 2. The effect of Rayleigh number, Ra, on the flow structure and heat transfer is analyzed for a range of Ra from 103 to 104. The results for established flow structures and heat transfer at low height indicate that the flow regime that occurs is unicellular for all Ra and fin thickness; in addition, the heat transfer rate increases with increasing Rayleigh number and is the same for both thicknesses. At median fin heights 0.093 and 0.203, the increase of Rayleigh number leads to transitions of flow structure which correspond to significant variations of the heat transfer. The critical Rayleigh numbers, Rac.app and Rac.disp corresponding to the appearance of the bicellular flow regime and its disappearance, are determined and their influence on the change of heat transfer rate is analyzed.
Abstract: Nowadays, data center industry faces strong challenges for increasing the speed and data processing capacities while at the same time is trying to keep their devices a suitable working temperature without penalizing that capacity. Consequently, the cooling systems of this kind of facilities use a large amount of energy to dissipate the heat generated inside the servers, and developing new cooling techniques or perfecting those already existing would be a great advance in this type of industry. The installation of a temperature sensor matrix distributed in the structure of each server would provide the necessary information for collecting the required data for obtaining a temperature profile instantly inside them. However, the number of temperature probes required to obtain the temperature profiles with sufficient accuracy is very high and expensive. Therefore, other less intrusive techniques are employed where each point that characterizes the server temperature profile is obtained by solving differential equations through simulation methods, simplifying data collection techniques but increasing the time to obtain results. In order to reduce these calculation times, complicated and slow computational fluid dynamics simulations are replaced by simpler and faster finite element method simulations which solve the Burgers‘ equations by backward, forward and central discretization techniques after simplifying the energy and enthalpy conservation differential equations. The discretization methods employed for solving the first and second order derivatives of the obtained Burgers‘ equation after these simplifications are the key for obtaining results with greater or lesser accuracy regardless of the characteristic truncation error.
Abstract: The neural network quantization is highly desired
procedure to perform before running neural networks on mobile
devices. Quantization without fine-tuning leads to accuracy drop of
the model, whereas commonly used training with quantization is done
on the full set of the labeled data and therefore is both time- and
resource-consuming. Real life applications require simplification and
acceleration of quantization procedure that will maintain accuracy of
full-precision neural network, especially for modern mobile neural
network architectures like Mobilenet-v1, MobileNet-v2 and MNAS. Here we present a method to significantly optimize training with
quantization procedure by introducing the trained scale factors for
discretization thresholds that are separate for each filter. Using the
proposed technique, we quantize the modern mobile architectures of
neural networks with the set of train data of only ∼ 10% of the
total ImageNet 2012 sample. Such reduction of train dataset size and
small number of trainable parameters allow to fine-tune the network
for several hours while maintaining the high accuracy of quantized
model (accuracy drop was less than 0.5%). Ready-for-use models and
code are available in the GitHub repository.
Abstract: This work is the modeling and simulation of fluid flow (liquid) through porous media. This type of flow occurs in many situations of interest in applied sciences and engineering, fluid (oil) consists of several individual substances in pure, single-phase flow is incompressible and isothermal. The porous medium is isotropic, homogeneous optionally, with the rectangular format and the flow is two-dimensional. Modeling of hydrodynamic phenomena incorporates Darcy's law and the equation of mass conservation. Correlations are used to model the density and viscosity of the fluid. A finite volume code is used in the discretization of differential equations. The nonlinearity is treated by Newton's method with relaxation coefficient. The results of the simulation of the pressure and the mobility of liquid flowing through porous media are presented, analyzed, and illustrated.
Abstract: Bioassay is the measurement of the potency of a chemical substance by its effect on a living animal or plant tissue. Bioassay data and chemical structures from pharmacokinetic and drug metabolism screening are mined from and housed in multiple databases. Bioassay prediction is calculated accordingly to determine further advancement. This paper proposes a four-step preprocessing of datasets for improving the bioassay predictions. The first step is instance selection in which dataset is categorized into training, testing, and validation sets. The second step is discretization that partitions the data in consideration of accuracy vs. precision. The third step is normalization where data are normalized between 0 and 1 for subsequent machine learning processing. The fourth step is feature selection where key chemical properties and attributes are generated. The streamlined results are then analyzed for the prediction of effectiveness by various machine learning algorithms including Pipeline Pilot, R, Weka, and Excel. Experiments and evaluations reveal the effectiveness of various combination of preprocessing steps and machine learning algorithms in more consistent and accurate prediction.
Abstract: Turbulence modelling is still evolving, and efforts are on to improve and develop numerical methods to simulate the real turbulence structures by using the empirical and experimental information. The monotonically integrated large eddy simulation (MILES) is an attractive approach for modelling turbulence in high Re flows, which is based on the solving of the unfiltered flow equations with no explicit sub-grid scale (SGS) model. In the current work, this approach has been used, and the action of the SGS model has been included implicitly by intrinsic nonlinear high-frequency filters built into the convection discretization schemes. The MILES solver is developed using the opensource CFD OpenFOAM libraries. The role of flux limiters schemes namely, Gamma, superBee, van-Albada and van-Leer, is studied in predicting turbulent statistical quantities for a fully developed channel flow with a friction Reynolds number, ReT = 180, and compared the numerical predictions with the well-established Direct Numerical Simulation (DNS) results for studying the wall generated turbulence. It is inferred from the numerical predictions that Gamma, van-Leer and van-Albada limiters produced more diffusion and overpredicted the velocity profiles, while superBee scheme reproduced velocity profiles and turbulence statistical quantities in good agreement with the reference DNS data in the streamwise direction although it deviated slightly in the spanwise and normal to the wall directions. The simulation results are further discussed in terms of the turbulence intensities and Reynolds stresses averaged in time and space to draw conclusion on the flux limiter schemes performance in OpenFOAM context.
Abstract: Dynamic analysis of composite doubly curved panels with variable thickness subjected to different pulse types using Generalized Differential Quadrature method (GDQ) is presented in this study. Panels with variable thickness are used in the construction of aerospace and marine industry. Giving variable thickness to panels can allow the designer to get optimum structural efficiency. For this reason, estimating the response of variable thickness panels is very important to design more reliable structures under dynamic loads. Dynamic equations for composite panels with variable thickness are obtained using virtual work principle. Partial derivatives in the equation of motion are expressed with GDQ and Newmark average acceleration scheme is used for temporal discretization. Several examples are used to highlight the effectiveness of the proposed method. Results are compared with finite element method. Effects of taper ratios, boundary conditions and loading type on the response of composite panel are investigated.
Abstract: The global coverage of broadband multimedia and
internet-based services in terrestrial-satellite networks demand
particular interests for satellite providers in order to enhance services
with low latencies and high signal quality to diverse users. In
particular, the delay of on-board processing is an inherent source
of latency in a satellite communication that sometimes is discarded
for the end-to-end delay of the satellite link. The frame work for this
paper includes modelling of an on-orbit satellite payload using an
agent model that can reproduce the properties of processing delays.
In essence, a comparison of different spatial interpolation methods is
carried out to evaluate physical data obtained by an GEO satellite
in order to define a discretization function for determining that
delay. Furthermore, the performance of the proposed agent and the
development of a delay discretization function are together validated
by simulating an hybrid satellite and terrestrial network. Simulation
results show high accuracy according to the characteristics of initial
data points of processing delay for Ku bands.
Abstract: The present work aims to investigate numerically the thermal and flow characteristics of a rectangular latent heat storage unit (LHSU) during the melting process of a phase change material (PCM). The LHSU consists of a number of vertical and identical plates of PCM separated by rectangular channels. The melting process is initiated when the LHSU is heated by a heat transfer fluid (HTF: water) flowing in channels in a downward or upward direction. The proposed study is motivated by the need to optimize the thermal performance of the LHSU by accelerating the charging process. A mathematical model is developed and a fixed-grid enthalpy formulation is adopted for modeling the melting process coupling with convection-conduction heat transfer. The finite volume method was used for discretization. The obtained numerical results are compared with experimental, analytical and numerical ones found in the literature and reasonable agreement is obtained. Thereafter, the numerical investigations were carried out to highlight the effects of the HTF flow direction and the aspect ratio of the PCM slabs on the heat transfer characteristics and thermal performance enhancement of the LHSU.
Abstract: The article presents an application of Fractional Model Predictive Control (FMPC) to a fractional order thermal system using Controlled Auto Regressive Integrated Moving Average (CARIMA) model obtained by discretization of a continuous fractional differential equation. Moreover, the output deviation approach is exploited to design the K -step ahead output predictor, and the corresponding control law is obtained by solving a quadratic cost function. Experiment results onto a thermal system are presented to emphasize the performances and the effectiveness of the proposed predictive controller.
Abstract: The performance of the tightening equipment will decline with the working process in manufacturing system. The main manifestations are the randomness and discretization degree increasing of the tightening performance. To evaluate the degradation tendency of the tightening performance accurately, a complexity measurement approach based on Kolmogorov entropy is presented. At first, the states of performance index are divided for calibrating the discrete degree. Then the complexity measurement model based on Kolmogorov entropy is built. The model describes the performance degradation tendency of tightening equipment quantitatively. At last, a study case is applied for verifying the efficiency and validity of the approach. The research achievement shows that the presented complexity measurement can effectively evaluate the degradation tendency of the tightening equipment. It can provide theoretical basis for preventive maintenance and life prediction of equipment.
Abstract: In this work, we present a Bayesian non-parametric
approach to model the motion control of ATVs. The motion control
model is based on a Dirichlet Process-Gaussian Process (DP-GP)
mixture model. The DP-GP mixture model provides a flexible
representation of patterns of control manoeuvres along trajectories
of different lengths and discretizations. The model also estimates the
number of patterns, sufficient for modeling the dynamics of the ATV.
Abstract: This paper develops a meshless approach, called Element Free Galerkin (EFG) method, which is based on the weak form Moving Least Squares (MLS) of the partial differential governing equations and employs the interpolation to construct the meshless shape functions. The variation weak form is used in the EFG where the trial and test functions are approximated bye the MLS approximation. Since the shape functions constructed by this discretization have the weight function property based on the randomly distributed points, the essential boundary conditions can be implemented easily. The local weak form of the partial differential governing equations is obtained by the weighted residual method within the simple local quadrature domain. The spline function with high continuity is used as the weight function. The presently developed EFG method is a truly meshless method, as it does not require the mesh, either for the construction of the shape functions, or for the integration of the local weak form. Several numerical examples of two-dimensional static structural analysis are presented to illustrate the performance of the present EFG method. They show that the EFG method is highly efficient for the implementation and highly accurate for the computation. The present method is used to analyze the static deflection of beams and plate hole
Abstract: In this paper, the influence of van der Waals, as well as electrostatic forces on the structural behavior of MEMS and NEMS actuators, has been investigated using of a Euler-Bernoulli beam continuous model. In the proposed nonlinear model, the electrostatic fringing-fields and the mid-plane stretching (geometric nonlinearity) effects have been considered. The nonlinear integro-differential equation governing the static structural behavior of the actuator has been derived. An original Galerkin-based reduced-order model has been developed to avoid problems arising from the nonlinearities in the differential equation. The obtained reduced-order model equations have been solved numerically using the Newton-Raphson method. The basic design parameters such as the pull-in parameters (voltage and deflection at pull-in), as well as the detachment length due to the van der Waals force of some investigated micro- and nano-actuators have been calculated. The obtained numerical results have been compared with some other existing methods (finite-elements method and finite-difference method) and the comparison showed good agreement among all assumed numerical techniques.
Abstract: The objective of this research work is to investigate
for one dimensional transient radiative transfer equations with
conduction using finite volume method. Within the infrastructure of
finite-volume, we obtain the conservative discretization of the terms
in order to preserve the overall conservative property of finitevolume
schemes. Coupling of conductive and radiative equation
resulting in fluxes is governed by the magnitude of emissivity,
extinction coefficient, and temperature of the medium as well as
geometry of the problem.
The problem under consideration has been solved, for a slab
dominating radiation coupled with transient conduction based on
finite volume method. The boundary conditions are also chosen so as
to give a good model of the discretized form of radiation transfer
equation. The important feature of the present method is flexibility in
specifying the control angles in the FVM, while keeping the
simplicity in the solution procedure.
Effects of various model parameters are examined on the
distributions of temperature, radiative and conductive heat fluxes and
incident radiation energy etc. The finite volume method is considered
to effectively evaluate the propagation of radiation intensity through
a participating medium.
Abstract: In this paper, we present a binary cat swarm
optimization for solving the Set covering problem. The set covering
problem is a well-known NP-hard problem with many practical
applications, including those involving scheduling, production
planning and location problems. Binary cat swarm optimization
is a recent swarm metaheuristic technique based on the behavior
of discrete cats. Domestic cats show the ability to hunt and are
curious about moving objects. The cats have two modes of behavior:
seeking mode and tracing mode. We illustrate this approach with
65 instances of the problem from the OR-Library. Moreover, we
solve this problem with 40 new binarization techniques and we select
the technical with the best results obtained. Finally, we make a
comparison between results obtained in previous studies and the new
binarization technique, that is, with roulette wheel as transfer function
and V3 as discretization technique.
Abstract: Steady three-dimensional and two free surface waves
generated by moving bodies are presented, the flow problem to be
simulated is rich in complexity and poses many modeling challenges
because of the existence of breaking waves around the ship hull, and
because of the interaction of the two-phase flow with the turbulent
boundary layer. The results of several simulations are reported. The
first study was performed for NACA0012 of hydrofoil with different
meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second
simulation a mathematically defined Wigley hull form is used to
investigate the application of a commercial CFD code in prediction of
the total resistance and its components from tangential and normal
forces on the hull wetted surface. The computed resistance and wave
profiles are used to estimate the coefficient of the total resistance for
Wigley hull advancing in calm water under steady conditions. The
commercial CFD software FLUENT version 12 is used for the
computations in the present study. The calculated grid is established
using the code computer GAMBIT 2.3.26. The shear stress k-ωSST
model is used for turbulence modeling and the volume of fluid
technique is employed to simulate the free-surface motion. The
second order upwind scheme is used for discretizing the convection
terms in the momentum transport equations, the Modified HRIC
scheme for VOF discretization. The results obtained compare well
with the experimental data.
Abstract: Steady three-dimensional and two free surface waves
generated by moving bodies are presented, the flow problem to be
simulated is rich in complexity and poses many modeling challenges
because of the existence of breaking waves around the ship hull, and
because of the interaction of the two-phase flow with the turbulent
boundary layer. The results of several simulations are reported. The
first study was performed for NACA0012 of hydrofoil with different
meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second
simulation a mathematically defined Wigley hull form is used to
investigate the application of a commercial CFD code in prediction of
the total resistance and its components from tangential and normal
forces on the hull wetted surface. The computed resistance and wave
profiles are used to estimate the coefficient of the total resistance for
Wigley hull advancing in calm water under steady conditions. The
commercial CFD software FLUENT version 12 is used for the
computations in the present study. The calculated grid is established
using the code computer GAMBIT 2.3.26. The shear stress k-ωSST
model is used for turbulence modeling and the volume of fluid
technique is employed to simulate the free-surface motion. The
second order upwind scheme is used for discretizing the convection
terms in the momentum transport equations, the Modified HRIC
scheme for VOF discretization. The results obtained compare well
with the experimental data.
Abstract: This paper reports the numerical simulation of doublediffusive
natural convection flows within a horizontal annular filled
with a saturated porous medium. The analysis concerns the influence
of the different parameters governing the problem, namely, the
Rayleigh number Ra, the Lewis number Le and the buoyancy ratio N,
on the heat and mass transfer and on the flow structure, in the case of
a fixed radius ratio R = 2. The numerical model used for the
discretization of the dimensionless equations governing the problem
is based on the finite difference method, using the ADI scheme. The
study is focused on steady-state solutions in the cooperation situation.
Abstract: In this paper, a new trend for improvement in semianalytical
method based on scale boundaries in order to solve the 2D
elastodynamic problems is provided. In this regard, only the
boundaries of the problem domain discretization are by specific subparametric
elements. Mapping functions are uses as a class of higherorder
Lagrange polynomials, special shape functions, Gauss-Lobatto-
Legendre numerical integration, and the integral form of the weighted
residual method, the matrix is diagonal coefficients in the equations
of elastodynamic issues. Differences between study conducted and
prior research in this paper is in geometry production procedure of
the interpolation function and integration of the different is selected.
Validity and accuracy of the present method are fully demonstrated
through two benchmark problems which are successfully modeled
using a few numbers of DOFs. The numerical results agree very well
with the analytical solutions and the results from other numerical
methods.