Abstract: Computations with higher than the IEEE 754 standard double-precision (about 16 significant digits) are required recently. Although there are available software routines in Fortran and C for high-precision computation, users are required to implement such routines in their own computers with detailed knowledges about them. We have constructed an user-friendly online system for octupleprecision computation. In our Web system users with no knowledges about high-precision computation can easily perform octupleprecision computations, by choosing mathematical functions with argument(s) inputted, by writing simple mathematical expression(s) or by uploading C program(s). In this paper we enhance the Web system above by adding the facility of uploading Fortran programs, which have been widely used in scientific computing. To this end we construct converter routines in two stages.
Abstract: In this paper, we present an efficient numerical algorithm, namely block homotopy perturbation method, for solving fuzzy linear systems based on homotopy perturbation method. Some numerical examples are given to show the efficiency of the algorithm.
Abstract: This paper deals with a numerical analysis of the
transient response of composite beams with strain rate dependent
mechanical properties by use of a finite difference method. The
equations of motion based on Timoshenko beam theory are derived.
The geometric nonlinearity effects are taken into account with von
Kármán large deflection theory. The finite difference method in
conjunction with Newmark average acceleration method is applied to
solve the differential equations. A modified progressive damage
model which accounts for strain rate effects is developed based on
the material property degradation rules and modified Hashin-type
failure criteria and added to the finite difference model. The
components of the model are implemented into a computer code in
Mathematica 6. Glass/epoxy laminated composite beams with
constant and strain rate dependent mechanical properties under
dynamic load are analyzed. Effects of strain rate on dynamic
response of the beam for various stacking sequences, load and
boundary conditions are investigated.
Abstract: This paper presents a novel two-phase hybrid optimization algorithm with hybrid genetic operators to solve the optimal control problem of a single stage hybrid manufacturing system. The proposed hybrid real coded genetic algorithm (HRCGA) is developed in such a way that a simple real coded GA acts as a base level search, which makes a quick decision to direct the search towards the optimal region, and a local search method is next employed to do fine tuning. The hybrid genetic operators involved in the proposed algorithm improve both the quality of the solution and convergence speed. The phase–1 uses conventional real coded genetic algorithm (RCGA), while optimisation by direct search and systematic reduction of the size of search region is employed in the phase – 2. A typical numerical example of an optimal control problem with the number of jobs varying from 10 to 50 is included to illustrate the efficacy of the proposed algorithm. Several statistical analyses are done to compare the validity of the proposed algorithm with the conventional RCGA and PSO techniques. Hypothesis t – test and analysis of variance (ANOVA) test are also carried out to validate the effectiveness of the proposed algorithm. The results clearly demonstrate that the proposed algorithm not only improves the quality but also is more efficient in converging to the optimal value faster. They can outperform the conventional real coded GA (RCGA) and the efficient particle swarm optimisation (PSO) algorithm in quality of the optimal solution and also in terms of convergence to the actual optimum value.
Abstract: Control chart pattern recognition is one of the most important tools to identify the process state in statistical process control. The abnormal process state could be classified by the recognition of unnatural patterns that arise from assignable causes. In this study, a wavelet based neural network approach is proposed for the recognition of control chart patterns that have various characteristics. The procedure of proposed control chart pattern recognizer comprises three stages. First, multi-resolution wavelet analysis is used to generate time-shape and time-frequency coefficients that have detail information about the patterns. Second, distance based features are extracted by a bi-directional Kohonen network to make reduced and robust information. Third, a back-propagation network classifier is trained by these features. The accuracy of the proposed method is shown by the performance evaluation with numerical results.
Abstract: The wind resource in the Italian site of Lendinara
(RO) is analyzed through a systematic anemometric campaign
performed on the top of the bell tower, at an altitude of over 100 m
above the ground. Both the average wind speed and the Weibull
distribution are computed. The resulting average wind velocity is in
accordance with the numerical predictions of the Italian Wind Atlas,
confirming the accuracy of the extrapolation of wind data adopted for
the evaluation of wind potential at higher altitudes with respect to the
commonly placed measurement stations.
Abstract: Dengue fever is an important human arboviral disease. Outbreaks are now reported quite often from many parts of the world. The number of cases involving pregnant women and infant cases are increasing every year. The illness is often severe and complications may occur. Deaths often occur because of the difficulties in early diagnosis and in the improper management of the diseases. Dengue antibodies from pregnant women are passed on to infants and this protects the infants from dengue infections. Antibodies from the mother are transferred to the fetus when it is still in the womb. In this study, we formulate a mathematical model to describe the transmission of this disease in pregnant women. The model is formulated by dividing the human population into pregnant women and non-pregnant human (men and non-pregnant women). Each class is subdivided into susceptible (S), infectious (I) and recovered (R) subclasses. We apply standard dynamical analysis to our model. Conditions for the local stability of the equilibrium points are given. The numerical simulations are shown. The bifurcation diagrams of our model are discussed. The control of this disease in pregnant women is discussed in terms of the threshold conditions.
Abstract: This paper describes simple implementation of
homotopy (also called continuation) algorithm for determining the proper resistance of the resistor to dissipate energy at a specified rate of an electric circuit. Homotopy algorithm can be considered as a developing of the classical methods in numerical computing such as Newton-Raphson and fixed
point methods. In homoptopy methods, an embedding
parameter is used to control the convergence. The method purposed in this work utilizes a special homotopy called Newton homotopy. Numerical example solved in MATLAB is given to show the effectiveness of the purposed method
Abstract: Thermoelastic temperature, displacement, and
stress in heat transfer during laser surface hardening are solved
in Eulerian formulation. In Eulerian formulations the heat flux
is fixed in space and the workpiece is moved through a control
volume. In the case of uniform velocity and uniform heat flux
distribution, the Eulerian formulations leads to a steady-state
problem, while the Lagrangian formulations remains transient.
In Eulerian formulations the reduction to a steady-state
problem increases the computational efficiency. In this study
also an analytical solution is developed for an uncoupled
transient heat conduction equation in which a plane slab is
heated by a laser beam. The thermal result of the numerical
model is compared with the result of this analytical model.
Comparing the results shows numerical solution for uncoupled
equations are in good agreement with the analytical solution.
Abstract: Phase-Contrast MR imaging methods are widely used
for measurement of blood flow velocity components. Also there are
some other tools such as CT and Ultrasound for velocity map
detection in intravascular studies. These data are used in deriving
flow characteristics. Some clinical applications are investigated
which use pressure distribution in diagnosis of intravascular disorders
such as vascular stenosis. In this paper an approach to the problem of
measurement of intravascular pressure field by using velocity field
obtained from flow images is proposed. The method presented in this
paper uses an algorithm to calculate nonlinear equations of Navier-
Stokes, assuming blood as an incompressible and Newtonian fluid.
Flow images usually suffer the lack of spatial resolution. Our
attempt is to consider the effect of spatial resolution on the pressure
distribution estimated from this method. In order to achieve this aim,
velocity map of a numerical phantom is derived at six different
spatial resolutions. To determine the effects of vascular stenoses on
pressure distribution, a stenotic phantom geometry is considered. A
comparison between the pressure distribution obtained from the
phantom and the pressure resulted from the algorithm is presented. In
this regard we also compared the effects of collocated and staggered
computational grids on the pressure distribution resulted from this
algorithm.
Abstract: This article illustrates a model selection management approach for virtual prototypes in interactive simulations. In those numerical simulations, the virtual prototype and its environment are modelled as a multiagent system, where every entity (prototype,human, etc.) is modelled as an agent. In particular, virtual prototyp ingagents that provide mathematical models of mechanical behaviour inform of computational methods are considered. This work argues that selection of an appropriate model in a changing environment,supported by models? characteristics, can be managed by the deter-mination a priori of specific exploitation and performance measures of virtual prototype models. As different models exist to represent a single phenomenon, it is not always possible to select the best one under all possible circumstances of the environment. Instead the most appropriate shall be selecting according to the use case. The proposed approach consists in identifying relevant metrics or indicators for each group of models (e.g. entity models, global model), formulate their qualification, analyse the performance, and apply the qualification criteria. Then, a model can be selected based on the performance prediction obtained from its qualification. The authors hope that this approach will not only help to inform engineers and researchers about another approach for selecting virtual prototype models, but also assist virtual prototype engineers in the systematic or automatic model selection.
Abstract: Design for Disassembly (DfD) aims to reuse the
structural components instead of demolition followed by recycling of
the demolition debris. This concept preserves the invested embodied
energy of materials, thus reducing inputs of new embodied energy
during materials reprocessing or remanufacturing. Both analytical and
experimental research on a proposed DfD beam-column connection
for use in residential apartments is currently investigated at the
National University of Singapore in collaboration with the Housing
and Development Board of Singapore. The present study reports on
the results of a numerical analysis of the proposed connection utilizing
finite element analysis. The numerical model was calibrated and
validated by comparison against experimental results. Results of a
parametric study will also be presented and discussed.
Abstract: Recently, as information industry and mobile
communication technology are developing, this study is conducted on
the new concept of intelligent structures and maintenance techniques
that applied wireless sensor network, USN (Ubiquitous Sensor
Network), to social infrastructures such as civil and architectural
structures on the basis of the concept of Ubiquitous Computing that
invisibly provides human life with computing, along with mutually
cooperating, compromising and connecting networks each other by
having computers within all objects around us.
Therefore, the purpose of this study is to investigate the capability
of wireless communication of sensor node embedded in reinforced
concrete structure with a basic experiment on an electric wave
permeability of sensor node by fabricating molding with variables of
concrete thickness and steel bars that are mostly used in constructing
structures to determine the feasibility of application to constructing
structures with USN.
At this time, with putting the pitches of steel bars, the thickness of
concrete placed, and the intensity of RF signal of a
transmitter-receiver as variables and when wireless communication
module was installed inside, the possible communication distance of
plain concrete and the possible communication distance by the pitches
of steel bars was measured in the horizontal and vertical direction
respectively. Besides, for the precise measurement of diminution of an
electric wave, the magnitude of an electric wave in the range of used
frequencies was measured by using Spectrum Analyzer. The
phenomenon of diminution of an electric wave was numerically
analyzed and the effect of the length of wavelength of frequencies was
analyzed by the properties of a frequency band area.
As a result of studying the feasibility of an application to
constructing structures with wireless sensor, in case of plain concrete,
it shows 45cm for the depth of permeability and in case of reinforced
concrete with the pitches of 5cm, it shows 37cm and 45cm for the
pitches of 15cm.
Abstract: The ultimate goal of this article is to develop a robust and accurate numerical method for solving hyperbolic conservation laws in one and two dimensions. A hybrid numerical method, coupling a cheap fourth order total variation diminishing (TVD) scheme [1] for smooth region and a Robust seventh-order weighted non-oscillatory (WENO) scheme [2] near discontinuities, is considered. High order multi-resolution analysis is used to detect the high gradients regions of the numerical solution in order to capture the shocks with the WENO scheme, while the smooth regions are computed with fourth order total variation diminishing (TVD). For time integration, we use the third order TVD Runge-Kutta scheme. The accuracy of the resulting hybrid high order scheme is comparable with these of WENO, but with significant decrease of the CPU cost. Numerical demonstrates that the proposed scheme is comparable to the high order WENO scheme and superior to the fourth order TVD scheme. Our scheme has the added advantage of simplicity and computational efficiency. Numerical tests are presented which show the robustness and effectiveness of the proposed scheme.
Abstract: In this work, the plate bending formulation of the boundary element method - BEM, based on the Reissner?s hypothesis, is extended to the analysis of plates reinforced by beams taking into account the membrane effects. The formulation is derived by assuming a zoned body where each sub-region defines a beam or a slab and all of them are represented by a chosen reference surface. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. In order to reduce the number of degrees of freedom, the problem values defined on the interfaces are written in terms of their values on the beam axis. Initially are derived separated equations for the bending and stretching problems, but in the final system of equations the two problems are coupled and can not be treated separately. Finally are presented some numerical examples whose analytical results are known to show the accuracy of the proposed model.
Abstract: In this paper, a delayed predator–prey system with stage
structure is investigated. Sufficient conditions for the system to have
multiple periodic solutions are obtained when the delay is sufficiently
large by applying Bendixson-s criterion. Further, some numerical
examples are given.
Abstract: The time dependent progress of a chemical reaction over a flat horizontal plate is here considered. The problem is solved through the group similarity transformation method which reduces the number of independent by one and leads to a set of nonlinear ordinary differential equation. The problem shows a singularity at the chemical reaction order n=1 and is analytically solved through the perturbation method. The behavior of the process is then numerically investigated for n≠1 and different Schmidt numbers. Graphical results for the velocity and concentration of chemicals based on the analytical and numerical solutions are presented and discussed.
Abstract: Contour filter strips planted with perennial vegetation
can be used to improve surface and ground water quality by reducing
pollutant, such as NO3-N, and sediment outflow from cropland to a
river or lake. Meanwhile, the filter strips of perennial grass with biofuel
potentials also have economic benefits of producing ethanol. In
this study, The Soil and Water Assessment Tool (SWAT) model was
applied to the Walnut Creek Watershed to examine the effectiveness
of contour strips in reducing NO3-N outflows from crop fields to the
river or lake. Required input data include watershed topography,
slope, soil type, land-use, management practices in the watershed and
climate parameters (precipitation, maximum/minimum air
temperature, solar radiation, wind speed and relative humidity).
Numerical experiments were conducted to identify potential
subbasins in the watershed that have high water quality impact, and
to examine the effects of strip size and location on NO3-N reduction
in the subbasins under various meteorological conditions (dry,
average and wet). Variable sizes of contour strips (10%, 20%, 30%
and 50%, respectively, of a subbasin area) planted with perennial
switchgrass were selected for simulating the effects of strip size and
location on stream water quality. Simulation results showed that a
filter strip having 10%-50% of the subbasin area could lead to 55%-
90% NO3-N reduction in the subbasin during an average rainfall
year. Strips occupying 10-20% of the subbasin area were found to be
more efficient in reducing NO3-N when placed along the contour
than that when placed along the river. The results of this study can
assist in cost-benefit analysis and decision-making in best water
resources management practices for environmental protection.
Abstract: This paper proposes the numerical simulation of the
investment casting of gold jewelry. It aims to study the behavior of
fluid flow during mould filling and solidification and to optimize the
process parameters, which lead to predict and control casting defects
such as gas porosity and shrinkage porosity. A finite difference
method, computer simulation software FLOW-3D was used to
simulate the jewelry casting process. The simplified model was
designed for both numerical simulation and real casting production.
A set of sensor acquisitions were allocated on the different positions
of the wax tree of the model to detect filling times, while a set of
thermocouples were allocated to detect the temperature during
casting and cooling. Those detected data were applied to validate the
results of the numerical simulation to the results of the real casting.
The resulting comparisons signify that the numerical simulation can
be used as an effective tool in investment-casting-process
optimization and casting-defect prediction.
Abstract: In this article, the phenomenon of nonlinear
consolidation in saturated and homogeneous clay layer is studied.
Considering time-varied drainage model, the excess pore water
pressure in the layer depth is calculated. The Generalized Differential
Quadrature (GDQ) method is used for the modeling and numerical
analysis. For the purpose of analysis, first the domain of independent
variables (i.e., time and clay layer depth) is discretized by the
Chebyshev-Gauss-Lobatto series and then the nonlinear system of
equations obtained from the GDQ method is solved by means of the
Newton-Raphson approach. The obtained results indicate that the
Generalized Differential Quadrature method, in addition to being
simple to apply, enjoys a very high accuracy in the calculation of
excess pore water pressure.