Abstract: In this paper, an advanced Nonlinear Exponential
Model (NEM), able to simulate the uniaxial dynamic behavior of
seismic isolators having a continuously decreasing tangent stiffness
with increasing displacement in the relatively large displacements
range and a hardening or softening behavior at large displacements, is
presented. The mathematical model is validated by comparing the
experimental force-displacement hysteresis loops obtained during
cyclic tests, conducted on a helical wire rope isolator and a recycled
rubber-fiber reinforced bearing, with those predicted analytically.
Good agreement between the experimental and simulated results
shows that the proposed model can be an effective numerical tool to
predict the force-displacement relationship of seismic isolation
devices within the large displacements range. Compared to the
widely used Bouc-Wen model, unable to simulate the response of
seismic isolators at large displacements, the proposed one allows to
avoid the numerical solution of a first order nonlinear ordinary
differential equation for each time step of a nonlinear time history
analysis, thus reducing the computation effort. Furthermore, the
proposed model can simulate the smooth transition of the hysteresis
loops from small to large displacements by adopting only one set of
five parameters determined from the experimental hysteresis loops
having the largest amplitude.
Abstract: In this paper, a one-dimensional (1d) Parallel Elasto-
Plastic Model (PEPM), able to simulate the uniaxial dynamic
behavior of seismic isolators having a continuously decreasing
tangent stiffness with increasing displacement, is presented. The
parallel modeling concept is applied to discretize the continuously
decreasing tangent stiffness function, thus allowing to simulate the
dynamic behavior of seismic isolation bearings by putting linear
elastic and nonlinear elastic-perfectly plastic elements in parallel. The
mathematical model has been validated by comparing the
experimental force-displacement hysteresis loops, obtained testing a
helical wire rope isolator and a recycled rubber-fiber reinforced
bearing, with those predicted numerically. Good agreement between
the simulated and experimental results shows that the proposed
model can be an effective numerical tool to predict the forcedisplacement
relationship of seismic isolators within relatively large
displacements. Compared to the widely used Bouc-Wen model, the
proposed one allows to avoid the numerical solution of a first order
ordinary nonlinear differential equation for each time step of a
nonlinear time history analysis, thus reducing the computation effort,
and requires the evaluation of only three model parameters from
experimental tests, namely the initial tangent stiffness, the asymptotic
tangent stiffness, and a parameter defining the transition from the
initial to the asymptotic tangent stiffness.
Abstract: The solution of the nonlinear dynamic equilibrium equations of base-isolated structures adopting a conventional monolithic solution approach, i.e. an implicit single-step time integration method employed with an iteration procedure, and the use of existing nonlinear analytical models, such as differential equation models, to simulate the dynamic behavior of seismic isolators can require a significant computational effort. In order to reduce numerical computations, a partitioned solution method and a one dimensional nonlinear analytical model are presented in this paper. A partitioned solution approach can be easily applied to base-isolated structures in which the base isolation system is much more flexible than the superstructure. Thus, in this work, the explicit conditionally stable central difference method is used to evaluate the base isolation system nonlinear response and the implicit unconditionally stable Newmark’s constant average acceleration method is adopted to predict the superstructure linear response with the benefit in avoiding iterations in each time step of a nonlinear dynamic analysis. The proposed mathematical model is able to simulate the dynamic behavior of seismic isolators without requiring the solution of a nonlinear differential equation, as in the case of widely used differential equation model. The proposed mixed explicit-implicit time integration method and nonlinear exponential model are adopted to analyze a three dimensional seismically isolated structure with a lead rubber bearing system subjected to earthquake excitation. The numerical results show the good accuracy and the significant computational efficiency of the proposed solution approach and analytical model compared to the conventional solution method and mathematical model adopted in this work. Furthermore, the low stiffness value of the base isolation system with lead rubber bearings allows to have a critical time step considerably larger than the imposed ground acceleration time step, thus avoiding stability problems in the proposed mixed method.
Abstract: The objective of this article is to improve the passive vibration damping of solar array (SA) used in space structures, by the effective application of numerical optimization. A case study of a SA is used for demonstration. A finite element (FE) model was created and verified by experimental testing. Optimization was then conducted by implementing the FE model with the genetic algorithm, to find the optimal placement of aluminum circular patches, to suppress the first two bending mode shapes. The results were verified using experimental testing. Finally, a parametric study was conducted using the FE model where patch locations, material type, and shape were varied one at a time, and the results were compared with the optimal ones. The results clearly show that through the proper application of FE modeling and numerical optimization, passive vibration damping of space structures has been successfully achieved.
Abstract: Treatment for the non-matching interface is an important computational issue. To handle this problem, the method of Lagrange multipliers including classical and localized versions are the most popular technique. It essentially imposes the interface compatibility conditions by introducing Lagrange multipliers. However, the numerical system becomes unstable and inefficient due to the Lagrange multipliers. The interface element-independent formulation that does not include the Lagrange multipliers can be obtained by modifying the independent variables mathematically. Through this modification, more efficient and stable system can be achieved while involving equivalent accuracy comparing with the conventional method. A numerical example is conducted to verify the validity of the presented method.
Abstract: The harmonic distortion of voltage is important in relation to power quality due to the interaction between the large diffusion of non-linear and time-varying single-phase and three-phase loads with power supply systems. However, harmonic distortion levels can be reduced by improving the design of polluting loads or by applying arrangements and adding filters. The application of passive filters is an effective solution that can be used to achieve harmonic mitigation mainly because filters offer high efficiency, simplicity, and are economical. Additionally, possible different frequency response characteristics can work to achieve certain required harmonic filtering targets. With these ideas in mind, the objective of this paper is to determine what size single tuned passive filters work in distribution networks best, in order to economically limit violations caused at a given point of common coupling (PCC). This article suggests that a single tuned passive filter could be employed in typical industrial power systems. Furthermore, constrained optimization can be used to find the optimal sizing of the passive filter in order to reduce both harmonic voltage and harmonic currents in the power system to an acceptable level, and, thus, improve the load power factor. The optimization technique works to minimize voltage total harmonic distortions (VTHD) and current total harmonic distortions (ITHD), where maintaining a given power factor at a specified range is desired. According to the IEEE Standard 519, both indices are viewed as constraints for the optimal passive filter design problem. The performance of this technique will be discussed using numerical examples taken from previous publications.
Abstract: The noise requirements for naval and research vessels
have seen an increasing demand for quieter ships in order to fulfil
current regulations and to reduce the effects on marine life. Hence,
new methods dedicated to the characterization of propeller noise,
which is the main source of noise in the far-field, are needed. The
study of cavitating propellers in closed-section is interesting for
analyzing hydrodynamic performance but could involve significant
difficulties for hydroacoustic study, especially due to reverberation
and boundary layer noise in the tunnel. The aim of this paper
is to present a numerical methodology for the identification of
hydroacoustic sources on marine propellers using hydrophone arrays
in a large hydrodynamic tunnel. The main difficulties are linked to the
reverberation of the tunnel and the boundary layer noise that strongly
reduce the signal-to-noise ratio. In this paper it is proposed to estimate
the reflection coefficients using an inverse method and some reference
transfer functions measured in the tunnel. This approach allows to
reduce the uncertainties of the propagation model used in the inverse
problem. In order to reduce the boundary layer noise, a cleaning
algorithm taking advantage of the low rank and sparse structure of the
cross-spectrum matrices of the acoustic and the boundary layer noise
is presented. This approach allows to recover the acoustic signal even
well under the boundary layer noise. The improvement brought by
this method is visible on acoustic maps resulting from beamforming
and DAMAS algorithms.
Abstract: In this work, the hemodynamics in the sinuses of
Valsalva after Transcatheter Aortic Valve Implantation is numerically
examined. We focus on the physical results in the two-dimensional
case. We use a finite element methodology based on a Lagrange
multiplier technique that enables to couple the dynamics of blood
flow and the leaflets’ movement. A massively parallel implementation
of a monolithic and fully implicit solver allows more accuracy and
significant computational savings. The elastic properties of the aortic
valve are disregarded, and the numerical computations are performed
under physiologically correct pressure loads. Computational results
depict that blood flow may be subject to stagnation in the lower
domain of the sinuses of Valsalva after Transcatheter Aortic Valve
Implantation.
Abstract: The heat and mass transfer characteristics of
contaminants in groundwater subjected to a biodegradation reaction
is analyzed by taking into account the thermal diffusion (Soret)
effects. This phenomenon is modulated mathematically by a
system of partial differential equations which govern the motion
of fluid (groundwater) and solid (contaminants) particles. The
numerical results are presented graphically for different values of
the parameters entering into the problem on the velocity profiles of
fluid, contaminants, temperature and concentration profile.
Abstract: The modernization of computer technology and commercial computational fluid dynamic (CFD) simulation has given better detailed results as compared to experimental investigation techniques. CFD techniques are widely used in different field due to its flexibility and performance. Evaluation of pipeline erosion is complex phenomenon to solve by numerical arithmetic technique, whereas CFD simulation is an easy tool to resolve that type of problem. Erosion wear behaviour due to solid–liquid mixture in the slurry pipeline has been investigated using commercial CFD code in FLUENT. Multi-phase Euler-Lagrange model was adopted to predict the solid particle erosion wear in 22.5° pipe bend for the flow of bottom ash-water suspension. The present study addresses erosion prediction in three dimensional 22.5° pipe bend for two-phase (solid and liquid) flow using finite volume method with standard k-ε turbulence, discrete phase model and evaluation of erosion wear rate with varying velocity 2-4 m/s. The result shows that velocity of solid-liquid mixture found to be highly dominating parameter as compared to solid concentration, density, and particle size. At low velocity, settling takes place in the pipe bend due to low inertia and gravitational effect on solid particulate which leads to high erosion at bottom side of pipeline.
Abstract: In this paper, we consider a geometric inverse source
problem for the heat equation with Dirichlet and Neumann boundary
data. We will reconstruct the exact form of the unknown source
term from additional boundary conditions. Our motivation is to
detect the location, the size and the shape of source support.
We present a one-shot algorithm based on the Kohn-Vogelius
formulation and the topological gradient method. The geometric
inverse source problem is formulated as a topology optimization
one. A topological sensitivity analysis is derived from a source
function. Then, we present a non-iterative numerical method for the
geometric reconstruction of the source term with unknown support
using a level curve of the topological gradient. Finally, we give
several examples to show the viability of our presented method.
Abstract: This paper presents a numerical investigation of regime transition of flow of ellipsoidal particles and a comparison with that of spherical particle assembly. Particle assemblies constituting spherical and ellipsoidal particle of 2.5:1 aspect ratio are examined at separate instances in similar flow conditions in a shear cell model that is numerically developed based on the discrete element method. Correlations among elastically scaled stress, kinetically scaled stress, coordination number and volume fraction are investigated, and show important similarities and differences for the spherical and ellipsoidal particle assemblies. In particular, volume fractions at points of regime transition are identified for both types of particles. It is found that compared with spherical particle assembly, ellipsoidal particle assembly has higher volume fraction for the quasistatic to intermediate regime transition and lower volume fraction for the intermediate to inertial regime transition. Finally, the relationship between coordination number and volume fraction shows strikingly distinct features for the two cases, suggesting that different from spherical particles, the effect of the shear rate on the coordination number is not significant for ellipsoidal particles. This work provides a glimpse of currently running work on one of the most attractive scopes of research in this field and has a wide prospect in understanding rheology of more complex shaped particles in light of the strong basis of simpler spherical particle rheology.
Abstract: Numerical Methods is a course that can be conducted using workshops and group discussion. This study has been implemented on undergraduate students of level two at the Faculty of Engineering, International Islamic University Malaysia. The Numerical Method course has been delivered to two Sections 1 and 2 with 44 and 22 students in each section, respectively. Systematic steps have been followed to apply the student centered learning approach in teaching Numerical Method course. Initially, the instructor has chosen the topic which was Euler’s Method to solve Ordinary Differential Equations (ODE) to be learned. The students were then divided into groups with five members in each group. Initial instructions have been given to the group members to prepare their subtopics before meeting members from other groups to discuss the subtopics in an expert group inside the classroom. For the time assigned for the classroom discussion, the setting of the classroom was rearranged to accommodate the student centered learning approach. Teacher strength was by monitoring the process of learning inside and outside the class. The students have been assessed during the migrating to the expert groups, recording of a video explanation outside the classroom and during the final examination. Euler’s Method to solve the ODE was set as part of Question 3(b) in the final exam. It is observed that none of the students from both sections obtained a zero grade in Q3(b), compared to Q3(a) and Q3(c). Also, for Section 1(44 students), 29 students obtained the full mark of 7/7, while only 10 obtained 7/7 for Q3(a) and no students obtained 6/6 for Q3(c). Finally, we can recommend that the Numerical Method course be moved toward more student-centered Learning classrooms where the students will be engaged in group discussion rather than having a teacher one man show.
Abstract: In this paper, an electromagnetic analysis is presented
for describing the influence of shielding in a rectangular waveguide.
A hybridization based on the method of moments combined to
the generalized equivalent circuit MoM-GEC is used to model the
problem. This is validated by applying the MoM-GEC hybridization
to investigate a diffraction structure. It consists of electromagnetic
diffraction by an iris in a rectangular waveguide. Numerical results
are shown and discussed and a comparison with FEM and Marcuvitz
methods is achieved.
Abstract: We examine two-dimensional oil displacement by water in a petroleum reservoir. The pore fluid is immiscible, and the porous media is homogenous and isotropic in the horizontal direction. Buckley-Leverett theory and a combination of Laplacian and Darcy’s law are used to study the fluid flow through porous media, and the Laplacian that defines the dispersion and diffusion of fluid in the sand using heavy oil is discussed. The reservoir is homogenous in the horizontal direction, as expressed by the partial differential equation. Two main factors which are observed are the water saturation and pressure distribution in the reservoir, and they are evaluated for predicting oil recovery in two dimensions by a physical and mathematical simulation model. We review the numerical simulation that solves difficult partial differential reservoir equations. Based on the numerical simulations, the saturation and pressure equations are calculated by the iterative alternating direction implicit method and the iterative alternating direction explicit method, respectively, according to the finite difference assumption. However, to understand the displacement of oil by water and the amount of water dispersion in the reservoir better, an interpolated contour line of the water distribution of the five-spot pattern, that provides an approximate solution which agrees well with the experimental results, is also presented. Finally, a computer program is developed to calculate the equation for pressure and water saturation and to draw the pressure contour line and water distribution contour line for the reservoir.
Abstract: The effects of thermal radiation, Soret and Dufour
parameters on mixed convection and nanofluid flow over a stretching
sheet in the presence of a magnetic field are investigated. The flow is
subject to temperature dependent viscosity and a chemical reaction
parameter. It is assumed that the nanoparticle volume fraction at the
wall may be actively controlled. The physical problem is modelled
using systems of nonlinear differential equations which have been
solved numerically using a spectral relaxation method. In addition
to the discussion on heat and mass transfer processes, the velocity,
nanoparticles volume fraction profiles as well as the skin friction
coefficient are determined for different important physical parameters.
A comparison of current findings with previously published results
for some special cases of the problem shows an excellent agreement.
Abstract: Hydraulic structures such as gravity dams are classified as essential structures, and have the vital role in providing strong and safe water resource management. Three major aspects must be considered to achieve an effective design of such a structure: 1) The building cost, 2) safety, and 3) accurate analysis of seepage characteristics. Due to the complexity and non-linearity relationships of the seepage process, many approximation theories have been developed; however, the application of these theories results in noticeable errors. The analytical solution, which includes the difficult conformal mapping procedure, could be applied for a simple and symmetrical problem only. Therefore, the objectives of this paper are to: 1) develop a surrogate model based on numerical simulated data using SEEPW software to approximately simulate seepage process related to a hydraulic structure, 2) develop and solve a linked simulation-optimization model based on the developed surrogate model to describe the seepage occurring under a concrete gravity dam, in order to obtain optimum and safe design at minimum cost. The result shows that the linked simulation-optimization model provides an efficient and optimum design of concrete gravity dams.
Abstract: For the stability and control demand of offshore small floating platform, a 2-HUS/U parallel mechanism was presented as offshore platform. Inverse kinematics was obtained by institutional constraint equation, and the dynamic model of offshore 2-HUS/U parallel platform was derived based on rigid body’s Lagrangian method. The equivalent moment of inertia, damping and driving force/torque variation of offshore 2-HUS/U parallel platform were analyzed. A numerical example shows that, for parallel platform of given motion, system’s equivalent inertia changes 1.25 times maximally. During the movement of platform, they change dramatically with the system configuration and have coupling characteristics. The maximum equivalent drive torque is 800 N. At the same time, the curve of platform’s driving force/torque is smooth and has good sine features. The control system needs to be adjusted according to kinetic equation during stability and control and it provides a basis for the optimization of control system.
Abstract: Thermally Activated Building Systems (TABS) have proven to be an energy-efficient solution to provide buildings with an optimal indoor thermal environment. This solution uses the structure of the building to store heat, reduce the peak loads, and decrease the primary energy demand. TABS require the heated or cooled surfaces to be as exposed as possible to the indoor space, but exposing the bare concrete surfaces has a diminishing effect on the acoustic qualities of the spaces in a building. Acoustic solutions capable of providing optimal acoustic comfort and allowing the heat exchange between the TABS and the room are desirable. In this study, the effects of free-hanging units on the cooling performance of TABS and the occupants’ thermal comfort was measured in a full-scale TABS laboratory. Investigations demonstrate that the use of free-hanging sound absorbers are compatible with the performance of TABS and the occupant’s thermal comfort, but an appropriate acoustic design is needed to find the most suitable solution for each case. The results show a reduction of 11% of the cooling performance of the TABS when 43% of the ceiling area is covered with free-hanging horizontal sound absorbers, of 23% for 60% ceiling coverage ratio and of 36% for 80% coverage. Measurements in actual buildings showed an increase of the room operative temperature of 0.3 K when 50% of the ceiling surface is covered with horizontal panels and of 0.8 to 1 K for a 70% coverage ratio. According to numerical simulations using a new TRNSYS Type, the use of comfort ventilation has a considerable influence on the thermal conditions in the room; if the ventilation is removed, then the operative temperature increases by 1.8 K for a 60%-covered ceiling.
Abstract: This research is aimed to study a two-step iteration
process defined over a finite family of σ-asymptotically
quasi-nonexpansive nonself-mappings. The strong convergence
is guaranteed under the framework of Banach spaces with some
additional structural properties including strict and uniform
convexity, reflexivity, and smoothness assumptions. With similar
projection technique for nonself-mapping in Hilbert spaces, we
hereby use the generalized projection to construct a point within
the corresponding domain. Moreover, we have to introduce the use
of duality mapping and its inverse to overcome the unavailability
of duality representation that is exploit by Hilbert space theorists.
We then apply our results for σ-asymptotically quasi-nonexpansive
nonself-mappings to solve for ideal efficiency of vector optimization
problems composed of finitely many objective functions. We also
showed that the obtained solution from our process is the closest to
the origin. Moreover, we also give an illustrative numerical example
to support our results.