Abstract: The paper presents a one-dimensional transient
mathematical model of thermal oil-water two-phase emulsion flows
in pipes. The set of the mass, momentum and enthalpy conservation
equations for the continuous fluid and droplet phases are solved. Two
friction correlations for the continuous fluid phase to wall friction are
accounted for in the model and tested. The aerodynamic drag force
between the continuous fluid phase and droplets is modeled, too. The
density and viscosity of both phases are assumed to be constant due
to adiabatic experimental conditions. The proposed mathematical
model is validated on the experimental measurements of oil-water
emulsion flows in horizontal pipe [1,2]. Numerical analysis on
single- and two-phase oil-water flows in a pipe is presented in the
paper. The continuous oil flow having water droplets is simulated.
Predictions, which are performed by using the presented model, show
excellent agreement with the experimental data if the water fraction is
equal or less than 10%. Disagreement between simulations and
measurements is increased if the water fraction is larger than 10%.
Abstract: This paper investigates experimental and numerical study of the airflow characteristics for vortex, round and square ceiling diffusers and its effect on the thermal comfort in a ventilated room. Three different thermal comfort criteria namely; Mean Age of the Air (MAA), ventilation effectiveness (E), and Effective Draft Temperature (EDT) have been used to predict the thermal comfort zone inside the room. In experimental work, a sub-scale room is set-up to measure the temperature field in the room. In numerical analysis, unstructured grids have been used to discretize the numerical domain. Conservation equations are solved using FLUENT commercial flow solver. The code is validated by comparing the numerical results obtained from three different turbulence models with the available experimental data. The comparison between the various numerical models shows that the standard k-ε turbulence model can be used to simulate these cases successfully. After validation of the code, effect of supply air velocity on the flow and thermal field could be investigated and hence the thermal comfort. The results show that the pressure coefficient created by the square diffuser is 1.5 times greater than that created by the vortex diffuser. The velocity decay coefficient is nearly the same for square and round diffusers and is 2.6 times greater than that for the vortex diffuser.
Abstract: The modified Arcan fixture was used in order to
investigate the mixed mode fracture properties of high strength steel
butt weld through experimental and numerical analysis. The fixture
consisted of a central section with "butterfly-shaped" specimen that
had central crack. The specimens were under pure mode I (opening),
pure mode II (shearing) and all in plane mixed mode loading angles
starting from 0 to 90 degrees. The geometric calibration factors were
calculated with the aid of finite element analysis for various loading
mode and different crack length (0.45≤ a/w ≤0.55) and the critical
fracture loads obtained experimentally. The critical fracture
toughness (KIC & KIIC) estimated with experimental and numerical
analysis under mixed mode loading conditions.
Abstract: This paper focuses on a technique for identifying the geological boundary of the ground strata in front of a tunnel excavation site using the first order adjoint method based on the optimal control theory. The geological boundary is defined as the boundary which is different layers of elastic modulus. At tunnel excavations, it is important to presume the ground situation ahead of the cutting face beforehand. Excavating into weak strata or fault fracture zones may cause extension of the construction work and human suffering. A theory for determining the geological boundary of the ground in a numerical manner is investigated, employing excavating blasts and its vibration waves as the observation references. According to the optimal control theory, the performance function described by the square sum of the residuals between computed and observed velocities is minimized. The boundary layer is determined by minimizing the performance function. The elastic analysis governed by the Navier equation is carried out, assuming the ground as an elastic body with linear viscous damping. To identify the boundary, the gradient of the performance function with respect to the geological boundary can be calculated using the adjoint equation. The weighed gradient method is effectively applied to the minimization algorithm. To solve the governing and adjoint equations, the Galerkin finite element method and the average acceleration method are employed for the spatial and temporal discretizations, respectively. Based on the method presented in this paper, the different boundary of three strata can be identified. For the numerical studies, the Suemune tunnel excavation site is employed. At first, the blasting force is identified in order to perform the accuracy improvement of analysis. We identify the geological boundary after the estimation of blasting force. With this identification procedure, the numerical analysis results which almost correspond with the observation data were provided.
Abstract: In this study, the precision heading process of
spur gears has been investigated by means of numerical
analysis. The effect of some parameters such as teeth number
and module on the forming force and material flow were
presented. The simulation works were performed rigid-plastic
finite element method using DEFORM 3D software. In order
to validate the estimated numerical results, they were
compared with those obtained experimentally during heading
of spur gear using lead as a model material. Results showed
that the optimum number of gear teeth is between 10 to 20,
that is because of being the specific pressure in its minimum
value.
Abstract: The paper presents a one-dimensional transient
mathematical model of compressible thermal multi-component gas
mixture flows in pipes. The set of the mass, momentum and enthalpy
conservation equations for gas phase is solved. Thermo-physical
properties of multi-component gas mixture are calculated by solving
the Equation of State (EOS) model. The Soave-Redlich-Kwong
(SRK-EOS) model is chosen. Gas mixture viscosity is calculated on
the basis of the Lee-Gonzales-Eakin (LGE) correlation. Numerical
analysis on rapid decompression in conventional dry gases is
performed by using the proposed mathematical model. The model is
validated on measured values of the decompression wave speed in
dry natural gas mixtures. All predictions show excellent agreement
with the experimental data at high and low pressure. The presented
model predicts the decompression in dry natural gas mixtures much
better than GASDECOM and OLGA codes, which are the most
frequently-used codes in oil and gas pipeline transport service.
Abstract: Nowadays, there is little information, concerning the
heat shield systems, and this information is not completely reliable to
use in so many cases. for example, the precise calculation cannot be
done for various materials. In addition, the real scale test has two
disadvantages: high cost and low flexibility, and for each case we
must perform a new test. Hence, using numerical modeling program
that calculates the surface recession rate and interior temperature
distribution is necessary. Also, numerical solution of governing
equation for non-charring material ablation is presented in order to
anticipate the recession rate and the heat response of non-charring
heat shields. the governing equation is nonlinear and the Newton-
Rafson method along with TDMA algorithm is used to solve this
nonlinear equation system. Using Newton- Rafson method for
solving the governing equation is one of the advantages of the
solving method because this method is simple and it can be easily
generalized to more difficult problems. The obtained results
compared with reliable sources in order to examine the accuracy of
compiling code.
Abstract: In this research work, investigations are carried out on
Continuous Wave (CW) Nd:YAG laser welding system after
preliminary experimentation to understand the influencing parameters
associated with laser welding of AISI 304. The experimental
procedure involves a series of laser welding trials on AISI 304
stainless steel sheets with various combinations of process parameters
like beam power, beam incident angle and beam incident angle. An
industrial 2 kW CW Nd:YAG laser system, available at Welding
Research Institute (WRI), BHEL Tiruchirappalli, is used for
conducting the welding trials for this research. After proper tuning of
laser beam, laser welding experiments are conducted on AISI 304
grade sheets to evaluate the influence of various input parameters on
weld bead geometry i.e. bead width (BW) and depth of penetration
(DOP). From the laser welding results, it is noticed that the beam
power and welding speed are the two influencing parameters on
depth and width of the bead. Three dimensional finite element
simulation of high density heat source have been performed for laser
welding technique using finite element code ANSYS for predicting
the temperature profile of laser beam heat source on AISI 304
stainless steel sheets. The temperature dependent material properties
for AISI 304 stainless steel are taken into account in the simulation,
which has a great influence in computing the temperature profiles.
The latent heat of fusion is considered by the thermal enthalpy of
material for calculation of phase transition problem. A Gaussian
distribution of heat flux using a moving heat source with a conical
shape is used for analyzing the temperature profiles. Experimental
and simulated values for weld bead profiles are analyzed for stainless
steel material for different beam power, welding speed and beam
incident angle. The results obtained from the simulation are
compared with those from the experimental data and it is observed
that the results of numerical analysis (FEM) are in good agreement
with experimental results, with an overall percentage of error
estimated to be within ±6%.
Abstract: The objective of this paper is to study the analysis and testing for determining the torsional stiffness of the student formula-s space frame. From past study, the space frame for Chulalongkorn University Student Formula team used in 2011 TSAE Auto Challenge Student Formula in Thailand was designed by considering required mass and torsional stiffness based on the numerical method and experimental method. The numerical result was compared with the experimental results to verify the torsional stiffness of the space frame. It can be seen from the large error of torsional stiffness of 2011 frame that the experimental result can not verify by the numerical analysis due to the different between the numerical model and experimental setting. In this paper, the numerical analysis and experiment of the same 2011 frame model is performed by improving the model setting. The improvement of both numerical analysis and experiment are discussed to confirm that the models from both methods are same. After the frame was analyzed and tested, the results are compared to verify the torsional stiffness of the frame. It can be concluded that the improved analysis and experiments can used to verify the torsional stiffness of the space frame.
Abstract: Real options theory suggests that managerial flexibility embedded within irreversible investments can account for a significant value in project valuation. Although the argument has become the dominant focus of capital investment theory over decades, yet recent survey literature in capital budgeting indicates that corporate practitioners still do not explicitly apply real options in investment decisions. In this paper, we explore how real options decision criteria can be transformed into equivalent capital budgeting criteria under the consideration of uncertainty, assuming that underlying stochastic process follows a geometric Brownian motion (GBM), a mixed diffusion-jump (MX), or a mean-reverting process (MR). These equivalent valuation techniques can be readily decomposed into conventional investment rules and “option impacts", the latter of which describe the impacts on optimal investment rules with the option value considered. Based on numerical analysis and Monte Carlo simulation, three major findings are derived. First, it is shown that real options could be successfully integrated into the mindset of conventional capital budgeting. Second, the inclusion of option impacts tends to delay investment. It is indicated that the delay effect is the most significant under a GBM process and the least significant under a MR process. Third, it is optimal to adopt the new capital budgeting criteria in investment decision-making and adopting a suboptimal investment rule without considering real options could lead to a substantial loss in value.
Abstract: In this paper, numerical simulation is used to
investigate the thermal performance of liquid cooling heatsink with
microchannels due to geometric arrangement. Commercial software
ICEPAK is utilized for the analysis. The considered parameters
include aspect ratio, porosity and the length and height of
microchannel. The aspect ratio varies from 3 to 16 and the length of
microchannel is 10mm, 14mm, and 18mm. The height of
microchannel is 2mm, 3mm and 4mm. It is found short channel have
better thermal efficiency than long channel at 490Pa. No matter the
length of channel the best aspect ratio is 4. It is also noted that pressure
difference at 2940Pa the best aspect ratio from 4 to 8, it means pressure
difference affect aspect ratio, effective thermal resistance at low
pressure difference but lower effective thermal resistance at high
pressure difference.
Abstract: This paper presents ageing experiments controlled by the evolution of junction parameters. The deterioration of the device is related to high injection effects which modified the transport mechanisms in the space charge region of the junction. Physical phenomena linked to the degradation of junction parameters that affect the devices reliability are reported and discussed. We have used the method based on numerical analysis of experimental current-voltage characteristic of the junction, in order to extract the electrical parameters. The simultaneous follow-up of the evolutions of the series resistance and of the transition voltage allow us to introduce a new parameter for reliability evaluation.
Abstract: In recent years, environment regulation forcing
manufactures to consider recovery activity of end-of- life products
and/or return products for refurbishing, recycling,
remanufacturing/repair and disposal in supply chain management. In
this paper, a mathematical model is formulated for single product
production-inventory system considering remanufacturing/reuse of
return products and rate of return products follows a demand like
function, dependent on purchasing price and acceptance quality level.
It is useful in decision making to determine whether to go for
remanufacturing or disposal of returned products along with newly
produced products to satisfy a stationary demand. In addition, a
modified genetic algorithm approach is proposed, inspired by particle
swarm optimization method. Numerical analysis of the case study is
carried out to validate the model.
Abstract: In this paper, experimental testing and numerical analysis were used to investigate the effect of tube thickness on the face bending for concrete filled hollow sections connected to other structural members using Extended Hollobolts. Six samples were tested experimentally by applying pull-out load on the bolts. These samples were designed to fail by column face bending. The main variable in all tests is the column face thickness. Finite element analyses were also performed using ABAQUS 6.11 to extend the experimental results and to quantify the effect of column face thickness. Results show that, the column face thickness has a clear impact on the connection strength and stiffness. However, the amount of improvement in the connection stiffness by changing the column face thickness from 5mm to 6.3mm seems to be higher than that when increasing it from 6.3mm to 8mm. The displacement at which the bolts start pulling-out from their holes increased with the use of thinner column face due to the high flexibility of the section. At the ultimate strength, the yielding of the column face propagated to the column corner and there was no yielding in its walls. After the ultimate resistance is reached, the propagation of the yielding was mainly in the column face with a miner yielding in the walls.
Abstract: In the present paper some recommendations for the
use of software package “Mathematica" in a basic numerical analysis
course are presented. The methods which are covered in the course
include solution of systems of linear equations, nonlinear equations
and systems of nonlinear equations, numerical integration,
interpolation and solution of ordinary differential equations. A set of
individual assignments developed for the course covering all the
topics is discussed in detail.
Abstract: The optimal grid spacing and turbulence model for the
2D numerical analysis of a vertical-axis water turbine (VAWaterT)
operating in a 2 m/s freestream current has been investigated. The
results of five different spatial domain discretizations and two
turbulence models (k-ω SST and k-ε RNG) have been compared, in
order to gain the optimal y+ parameter distribution along the blade
walls during a full rotor revolution. The resulting optimal mesh has
appeared to be quite similar to that obtained for the numerical
analysis of a vertical-axis wind turbine.
Abstract: In present work the problem of the ITER fusion
plasma neutron source parameter reconstruction using only the
Vertical Neutron Camera data was solved. The possibility of neutron
source parameter reconstruction was estimated by the numerical
simulations and the analysis of adequateness of mathematic model
was performed. The neutron source was specified in a parametric
form. The numerical analysis of solution stability with respect to data
distortion was done. The influence of the data errors on the
reconstructed parameters is shown:
• is reconstructed with errors less than 4% at all examined values
of δ (until 60%);
• is determined with errors less than 10% when δ do not overcome
5%;
• is reconstructed with relative error more than 10 %;
• integral intensity of the neutron source is determined with error
10% while δ error is less than 15%;
where -error of signal measurements, (R0,Z0), the plasma center
position,- /parameter of neutron source profile.
Abstract: This paper presents a numerical analysis of the
seismic behaviour of earth dams. Analysis is conducted for the solid
phase. It may correspond to the response of the dam before water
filling. Analysis is conducted for a simple case which concerns the
elastic response of the dam. Numerical analyses are conducted using
the FLAC3D program. The behaviour of the Shell and core of the
dam and the foundation behaviour is assumed to be elastic. Result
shows the influence of the variation of the shear modulus of the core
and shell on the seismic amplification of the dam. It can be observed
that the variation of the shearing modulus of the core leads to a
moderate increase in the dynamic amplification and the increase in
the shell shearing modulus leads to a significant increase in the
dynamic amplification.
Abstract: This paper describes dynamic analysis using proposed
fast finite element method for a shock absorbing structure including a
sponge. The structure is supported by nonlinear concentrated springs.
The restoring force of the spring has cubic nonlinearity and linear
hysteresis damping. To calculate damping properties for the structures
including elastic body and porous body, displacement vectors as
common unknown variable are solved under coupled condition. Under
small amplitude, we apply asymptotic method to complex eigenvalue
problem of this system to obtain modal parameters. And then
expressions of modal loss factor are derived approximately. This
approach was proposed by one of the authors previously. We call this
method as Modal Strain and Kinetic Energy Method (MSKE method).
Further, using the modal loss factors, the discretized equations in
physical coordinate are transformed into the nonlinear ordinary
coupled equations using normal coordinate corresponding to linear
natural modes. This transformation yields computation efficiency. As
a numerical example of a shock absorbing structures, we adopt double
skins with a sponge. The double skins are supported by nonlinear
concentrated springs. We clarify influences of amplitude of the input
force on nonlinear and chaotic responses.
Abstract: In this paper, all-optical signal processors that perform
both microwave mixing and bandpass filtering in a radio-over-fiber
(RoF) link are presented. The key device is a Mach-Zehnder
modulator (MZM) which performs all-optical microwave mixing. An
up-converted microwave signal is obtained and other unwanted
frequency components are suppressed at the end of the fiber span.