Abstract: The authors present an algorithm for order reduction of linear time invariant dynamic systems using the combined advantages of the eigen spectrum analysis and the error minimization by particle swarm optimization technique. Pole centroid and system stiffness of both original and reduced order systems remain same in this method to determine the poles, whereas zeros are synthesized by minimizing the integral square error in between the transient responses of original and reduced order models using particle swarm optimization technique, pertaining to a unit step input. It is shown that the algorithm has several advantages, e.g. the reduced order models retain the steady-state value and stability of the original system. The algorithm is illustrated with the help of two numerical examples and the results are compared with the other existing techniques.
Abstract: to simulate the phenomenon of electronic transport in semiconductors, we try to adapt a numerical method, often and most frequently it’s that of Monte Carlo. In our work, we applied this method in the case of a ternary alloy semiconductor GaInP in its cubic form; The Calculations are made using a non-parabolic effective-mass energy band model. We consider a band of conduction to three valleys (ΓLX), major of the scattering mechanisms are taken into account in this modeling, as the interactions with the acoustic phonons (elastic collisions) and optics (inelastic collisions). The polar optical phonons cause anisotropic collisions, intra-valleys, very probable in the III-V semiconductors. Other optical phonons, no polar, allow transitions inter-valleys. Initially, we present the full results obtained by the simulation of Monte Carlo in GaInP in stationary regime. We consider thereafter the effects related to the application of an electric field varying according to time, we thus study the transient phenomenon which make their appearance in ternary material
Abstract: In contrast to existing methods which do not take into account multiconnectivity in a broad sense of this term, we develop mathematical models and highly effective combination (BIEM and FDM) numerical methods of calculation of stationary and quasistationary temperature field of a profile part of a blade with convective cooling (from the point of view of realization on PC). The theoretical substantiation of these methods is proved by appropriate theorems. For it, converging quadrature processes have been developed and the estimations of errors in the terms of A.Ziqmound continuity modules have been received. For visualization of profiles are used: the method of the least squares with automatic conjecture, device spline, smooth replenishment and neural nets. Boundary conditions of heat exchange are determined from the solution of the corresponding integral equations and empirical relationships. The reliability of designed methods is proved by calculation and experimental investigations heat and hydraulic characteristics of the gas turbine first stage nozzle blade.
Abstract: The paper presents a numerical investigation on the
rapid gas decompression in pure nitrogen which is made by using the
one-dimensional (1D) and three-dimensional (3D) mathematical
models of transient compressible non-isothermal fluid flow in pipes.
A 1D transient mathematical model of compressible thermal multicomponent
fluid mixture flow in pipes is presented. The set of the
mass, momentum and enthalpy conservation equations for gas phase
is solved in the model. Thermo-physical properties of multicomponent
gas mixture are calculated by solving the Equation of
State (EOS) model. The Soave-Redlich-Kwong (SRK-EOS) model is
chosen. This model is successfully validated on the experimental data
[1] and shows a good agreement with measurements. A 3D transient
mathematical model of compressible thermal single-component gas
flow in pipes, which is built by using the CFD Fluent code (ANSYS),
is presented in the paper. The set of unsteady Reynolds-averaged
conservation equations for gas phase is solved. Thermo-physical
properties of single-component gas are calculated by solving the Real
Gas Equation of State (EOS) model. The simplest case of gas
decompression in pure nitrogen is simulated using both 1D and 3D
models. The ability of both models to simulate the process of rapid
decompression with a high order of agreement with each other is
tested. Both, 1D and 3D numerical results show a good agreement
between each other. The numerical investigation shows that 3D CFD
model is very helpful in order to validate 1D simulation results if the
experimental data is absent or limited.
Abstract: We discuss the convergence property of the minimum residual (MINRES) method for the solution of complex shifted Hermitian system (αI + H)x = f. Our convergence analysis shows that the method has a faster convergence than that for real shifted Hermitian system (Re(α)I + H)x = f under the condition Re(α) + λmin(H) > 0, and a larger imaginary part of the shift α has a better convergence property. Numerical experiments show such convergence properties.
Abstract: A nonlinear model of two-beam free-electron laser
(FEL) in the absence of slippage is presented. The two beams are
assumed to be cold with different energies and the fundamental
resonance of the higher energy beam is at the third harmonic of lower
energy beam. By using Maxwell-s equations and full Lorentz force
equations of motion for the electron beams, coupled differential
equations are derived and solved numerically by the fourth order
Runge–Kutta method. In this method a considerable growth of third
harmonic electromagnetic field in the XUV and X-ray regions is
predicted.
Abstract: The Constraints imposed by non-thermal
leptogenesis on the survival of the neutrino mass models describing
the presently available neutrino mass patterns, are studied
numerically. We consider the Majorana CP violating phases coming
from right-handed Majorana mass matrices to estimate the baryon
asymmetry of the universe, for different neutrino mass models
namely quasi-degenerate, inverted hierarchical and normal
hierarchical models, with tribimaximal mixings. Considering two
possible diagonal forms of Dirac neutrino mass matrix as either
charged lepton or up-quark mass matrix, the heavy right-handed
mass matrices are constructed from the light neutrino mass matrix.
Only the normal hierarchical model leads to the best predictions of
baryon asymmetry of the universe, consistent with observations in
non-thermal leptogenesis scenario.
Abstract: We present a numerical study of the sensitivity of the so called time relaxation family of models of fluid motion with respect to the time relaxation parameter χ on the two dimensional cavity problem. The goal of the study is to compute and compare the sensitivity of the model using finite difference method (FFD) and sensitivity equation method (SEM).
Abstract: In order to enhance the aircraft survivability, the
infrared signatures emitted by hot engine parts should be determined
exactly. For its reduction it is necessary for the rear fuselage
temperature to be decreased. In this study, numerical modeling of flow
fields and heat transfer characteristics of an aircraft nozzle is
performed and its temperature distribution along each component wall
is predicted. The radiation shield is expected to reduce the skin
temperature of rear fuselage. The effect of material characteristic of
radiation shield on the heat transfer is also investigated. Through this
numerical analysis, design parameters related to the susceptibility of
aircraft are examined.
Abstract: This paper addresses the problem of trajectory
tracking control of an underactuated autonomous underwater vehicle
(AUV) in the horizontal plane. The underwater vehicle under
consideration is not actuated in the sway direction, and the system
matrices are not assumed to be diagonal and linear, as often found in
the literature. In addition, the effect of constant bias of environmental
disturbances is considered. Using backstepping techniques and the
tracking error dynamics, the system states are stabilized by forcing
the tracking errors to an arbitrarily small neighborhood of zero. The
effectiveness of the proposed control method is demonstrated through
numerical simulations. Simulations are carried out for an
experimental vehicle for smooth, inertial, two dimensional (2D)
reference trajectories such as constant velocity trajectory (a circle
maneuver – constant yaw rate), and time varying velocity trajectory
(a sinusoidal path – sinusoidal yaw rate).
Abstract: This purpose of this paper is to present the acceptance single sampling plan when the fraction of nonconforming items is a fuzzy number and being modeled based on the fuzzy Poisson distribution. We have shown that the operating characteristic (oc) curves of the plan is like a band having a high and low bounds whose width depends on the ambiguity proportion parameter in the lot when that sample size and acceptance numbers is fixed. Finally we completed discuss opinion by a numerical example. And then we compared the oc bands of using of binomial with the oc bands of using of Poisson distribution.
Abstract: In this paper, free vibration analysis of carbon nanotube (CNT) reinforced laminated composite panels is presented. Three types of panels such as flat, concave and convex are considered for study. Numerical simulation is carried out using commercially available finite element analysis software ANSYS. Numerical homogenization is employed to calculate the effective elastic properties of randomly distributed carbon nanotube reinforced composites. To verify the accuracy of the finite element method, comparisons are made with existing results available in the literature for conventional laminated composite panels and good agreements are obtained. The results of the CNT reinforced composite materials are compared with conventional composite materials under different boundary conditions.
Abstract: A numerical study of flow in a horizontally channel
partially filled with a porous screen with non-uniform inlet has been
performed by lattice Boltzmann method (LBM). The flow in porous
layer has been simulated by the Brinkman-Forchheimer model.
Numerical solutions have been obtained for variable porosity models
and the effects of Darcy number and porosity have been studied in
detail. It is found that the flow stabilization is reliant on the Darcy
number. Also the results show that the stabilization of flow field and
heat transfer is depended to Darcy number. Distribution of stream
field becomes more stable by decreasing Darcy number. Results
illustrate that the effect of variable porosity is significant just in the
region of the solid boundary. In addition, difference between constant
and variable porosity models is decreased by decreasing the Darcy
number.
Abstract: Numerical parametric study is conducted to study the effects of ampoule rotation on the flows and the dopant segregation in vertical bridgman (vb) crystal growth. Calculations were performed in unsteady state. The extended darcy model, which includes the time derivative and coriolis terms, has been employed in the momentum equation. It’s found that the convection, and dopant segregation can be affected significantly by ampoule rotation, and the effect is similar to that by an axial magnetic field. Ampoule rotation decreases the intensity of convection and stretches the flow cell axially. When the convection is weak, the flow can be suppressed almost completely by moderate ampoule rotation and the dopant segregation becomes diffusion-controlled. For stronger convection, the elongated flow cell by ampoule rotation may bring dopant mixing into the bulk melt reducing axial segregation at the early stage of the growth. However, if the cellular flow cannot be suppressed completely, ampoule rotation may induce larger radial segregation due to poor mixing.
Abstract: We study the semiconvergence of Gauss-Seidel iterative
methods for the least squares solution of minimal norm of rank
deficient linear systems of equations. Necessary and sufficient conditions
for the semiconvergence of the Gauss-Seidel iterative method
are given. We also show that if the linear system of equations is
consistent, then the proposed methods with a zero vector as an initial
guess converge in one iteration. Some numerical results are given to
illustrate the theoretical results.
Abstract: In this paper a new method is suggested for risk
management by the numerical patterns in data-mining. These patterns
are designed using probability rules in decision trees and are cared to
be valid, novel, useful and understandable. Considering a set of
functions, the system reaches to a good pattern or better objectives.
The patterns are analyzed through the produced matrices and some
results are pointed out. By using the suggested method the direction
of the functionality route in the systems can be controlled and best
planning for special objectives be done.
Abstract: Internal combustion engines rejects 30-40% of the
energy supplied by fuel to the environment through exhaust gas. thus, there is a possibility for further significant improvement of efficiency with the utilization of exhaust gas energy and its conversion to mechanical energy or electrical energy. The Thermo-Electric
Generator (TEG) will be located in the exhaust system and will make use of an energy flow between the warmer exhaust gas and the external environment. Predict to th optimum position of temperature
distribution and the performance of TEG through numerical analysis.
The experimental results obtained show that the power output significantly increases with the temperature difference between cold
and hot sides of a thermoelectric generator.
Abstract: In this paper, solution of fuzzy differential equation
under general differentiability is obtained by simulink. The simulink
solution is equivalent or very close to the exact solution of the
problem. Accuracy of the simulink solution to this problem is
qualitatively better. An illustrative numerical example is presented
for the proposed method.
Abstract: The effect of different combinations of response
feedback on the performance of active control system on nonlinear
frames has been studied in this paper. To this end different feedback
combinations including displacement, velocity, acceleration and full
response feedback have been utilized in controlling the response of
an eight story bilinear hysteretic frame which has been subjected to a
white noise excitation and controlled by eight actuators which could
fully control the frame. For active control of nonlinear frame
Newmark nonlinear instantaneous optimal control algorithm has been
used which a diagonal matrix has been selected for weighting
matrices in performance index. For optimal design of active control
system while the objective has been to reduce the maximum drift to
below the yielding level, Distributed Genetic Algorithm (DGA) has
been used to determine the proper set of weighting matrices. The
criteria to assess the effect of each combination of response feedback
have been the minimum required control force to reduce the
maximum drift to below the yielding drift. The results of numerical
simulation show that the performance of active control system is
dependent on the type of response feedback where the velocity
feedback is more effective in designing optimal control system in
comparison with displacement and acceleration feedback. Also using
full feedback of response in controller design leads to minimum
control force amongst other combinations. Also the distributed
genetic algorithm shows acceptable convergence speed in solving the
optimization problem of designing active control systems.
Abstract: The indoor airflow with a mixed natural/forced convection
was numerically calculated using the laminar and turbulent
approach. The Boussinesq approximation was considered for a simplification
of the mathematical model and calculations. The results
obtained, such as mean velocity fields, were successfully compared
with experimental PIV flow visualizations. The effect of the distance
between the cooled wall and the heat exchanger on the temperature
and velocity distributions was calculated. In a room with a simple
shape, the computational code OpenFOAM demonstrated an ability to
numerically predict flow patterns. Furthermore, numerical techniques,
boundary type conditions and the computational grid quality were
examined. Calculations using the turbulence model k-omega had a
significant effect on the results influencing temperature and velocity
distributions.