Abstract: Passive resonant vibration absorbers are among the most widely used dynamic control systems in civil engineering. They typically consist in a single-degree-of-freedom mechanical appendage of the main structure, tuned to one structural target mode through frequency and damping optimization. One classical scheme is the pendulum absorber, whose mass is constrained to move along a curved trajectory and is damped by viscous dashpots. Even though the principle is well known, the search for improved arrangements is still under way. In recent years this investigation inspired a type of bidirectional pendulum absorber (BPA), consisting of a mass constrained to move along an optimal three-dimensional (3D) concave surface. For such a BPA, the surface principal curvatures are designed to ensure a bidirectional tuning of the absorber to both principal modes of the main structure, while damping is produced either by horizontal viscous dashpots or by vertical friction dashpots, connecting the BPA to the main structure. In this paper, a variant of BPA is proposed, where damping originates from the variable tangential friction force which develops between the pendulum mass and the 3D surface as a result of a spatially-varying friction coefficient pattern. Namely, a friction coefficient is proposed that varies along the pendulum surface in proportion to the modulus of the 3D surface gradient. With such an assumption, the dissipative model of the absorber can be proven to be nonlinear homogeneous in the small displacement domain. The resulting homogeneous BPA (HBPA) has a fundamental advantage over conventional friction-type absorbers, because its equivalent damping ratio results independent on the amplitude of oscillations, and therefore its optimal performance does not depend on the excitation level. On the other hand, the HBPA is more compact than viscously damped BPAs because it does not need the installation of dampers. This paper presents the analytical model of the HBPA and an optimal methodology for its design. Numerical simulations of single- and multi-story building structures under wind and earthquake loads are presented to compare the HBPA with classical viscously damped BPAs. It is shown that the HBPA is a promising alternative to existing BPA types and that homogeneous tangential friction is an effective means to realize systems provided with amplitude-independent damping.
Abstract: The control problem of underactuated spacecrafts has
attracted a considerable amount of interest. The control method for
a spacecraft equipped with less than three control torques is useful
when one of the three control torques had failed. On the other hand,
the quantized control of systems is one of the important research
topics in recent years. The random dither quantization method that
transforms a given continuous signal to a discrete signal by adding
artificial random noise to the continuous signal before quantization
has also attracted a considerable amount of interest. The objective of
this study is to develop the control method based on random dither
quantization method for stabilizing the rotational motion of a rigid
spacecraft with two control inputs. In this paper, the effectiveness of
random dither quantization control method for the stabilization of
rotational motion of spacecrafts with two torque inputs is verified
by numerical simulations.
Abstract: Generation of high DC voltages is necessary for testing the insulation material of high voltage AC transmission lines with long lengths. The harmonic and ripple contents of the output DC voltage supplied by high voltage DC circuits require the use of costly capacitors to smooth the output voltage after rectification. This paper proposes a new modular multiplier high voltage DC generator with embedded Cockcroft-Walton circuits that achieve a negligible harmonic and ripple contents of the output DC voltage without the need for costly filters to produce a nearly constant output voltage. In this new topology, Cockcroft-Walton modules are connected in series to produce a high DC output voltage. The modules are supplied by low input AC voltage sources that have the same magnitude and frequency and shifted from each other by a certain angle to eliminate the harmonics from the output voltage. The small ripple factor is provided by the smoothing column capacitors and the phase shifted input voltages of the cascaded modules. The constituent harmonics within each module are determined using Fourier analysis. The viability of the proposed DC generator for testing purposes and the effectiveness of the cascaded connection are confirmed by numerical simulations using MATLAB/Simulink.
Abstract: For any country in the world, it has become a priority to protect the critical infrastructure from looming risks of terrorism. In any infrastructure system, the structural elements like lower floors, exterior columns, walls etc. are key elements which are the most susceptible to damage due to blast load. The present study revisits the state of art review of the design and analysis of reinforced concrete panels subjected to blast loading. Various aspects in association with blast loading on structure, i.e. estimation of blast load, experimental works carried out previously, the numerical simulation tools, various material models, etc. are considered for exploring the current practices adopted worldwide. Discussion on various parametric studies to investigate the effect of reinforcement ratios, thickness of slab, different charge weight and standoff distance is also made. It was observed that for the simulation of blast load, CONWEP blast function or equivalent numerical equations were successfully employed by many researchers. The study of literature indicates that the researches were carried out using experimental works and numerical simulation using well known generalized finite element methods, i.e. LS-DYNA, ABAQUS, AUTODYN. Many researchers recommended to use concrete damage model to represent concrete and plastic kinematic material model to represent steel under action of blast loads for most of the numerical simulations. Most of the studies reveal that the increase reinforcement ratio, thickness of slab, standoff distance was resulted in better blast resistance performance of reinforced concrete panel. The study summarizes the various research results and appends the present state of knowledge for the structures exposed to blast loading.
Abstract: In response to widely used wearable medical devices equipped with a continuous glucose monitor (CGM) and insulin pump, the advanced control methods are still demanding to get the full benefit of these devices. Unlike costly clinical trials, implementing effective insulin-glucose control strategies can provide significant contributions to the patients suffering from chronic diseases such as diabetes. This study deals with a key role of two-layer insulin-glucose regulator based on model-predictive-control (MPC) scheme so that the patient’s predicted glucose profile is in compliance with the insulin level injected through insulin pump automatically. It is achieved by iterative optimization algorithm which is called an integrated perturbation analysis and sequential quadratic programming (IPA-SQP) solver for handling uncertainties due to unexpected variations in glucose-insulin values and body’s characteristics. The feasibility evaluation of the discussed control approach is also studied by means of numerical simulations of two case scenarios via measured data. The obtained results are presented to verify the superior and reliable performance of the proposed control scheme with no negative impact on patient safety.
Abstract: Zone Plates (ZPs) are used in many areas of physics where planar fabrication is advantageous in comparison with conventional curved lenses. There are several types of ZPs, such as the well-known Fresnel ZPs or the more recent Fractal ZPs and Fibonacci ZPs. The material selection of the lens plays a very important role in the beam modulation control. This work presents a comparison between two Fresnel ZP made from different materials in the ultrasound domain: Polylactic Acid (PLA) and brass. PLA is the most common material used in commercial 3D-printers due to its high design flexibility and low cost. Numerical simulations based on Finite Element Method (FEM) and experimental results are shown, and they prove that the focusing capabilities of brass ZPs and PLA ZPs are similar. For this reason, PLA is proposed as a Magnetic Resonance Imaging (MRI) compatible material with great potential for therapeutic ultrasound focusing applications.
Abstract: We present a new class of numerical techniques to
solve shallow water flows over dry areas including run-up. Many
recent investigations on wave run-up in coastal areas are based on
the well-known shallow water equations. Numerical simulations have
also performed to understand the effects of several factors on tsunami
wave impact and run-up in the presence of coastal areas. In all these
simulations the shallow water equations are solved in entire domain
including dry areas and special treatments are used for numerical
solution of singularities at these dry regions. In the present study we
propose a new method to deal with these difficulties by reformulating
the shallow water equations into a new system to be solved only in the
wetted domain. The system is obtained by a change in the coordinates
leading to a set of equations in a moving domain for which the
wet/dry interface is the reconstructed using the wave speed. To solve
the new system we present a finite volume method of Lax-Friedrich
type along with a modified method of characteristics. The method is
well-balanced and accurately resolves dam-break problems over dry
areas.
Abstract: We develop a fast, user-friendly implementation of
a potential flow solver based on the unsteady vortex lattice
method (UVLM). The computational framework uses the Python
programming language which has easy integration with the scripts
requiring computationally-expensive operations written in Fortran.
The mixed-language approach enables high performance in terms
of solution time and high flexibility in terms of easiness of code
adaptation to different system configurations and applications. This
computational tool is intended to predict the unsteady aerodynamic
behavior of multiple moving bodies (e.g., flapping wings, rotating
blades, suspension bridges...) subject to an incoming air. We
simulate different aerodynamic problems to validate and illustrate
the usefulness and effectiveness of the developed computational tool.
Abstract: The physical effects of upstream flow obstructions such
as vegetation on cross-ventilation phenomena of a building are
important for issues such as indoor thermal comfort. Modelling such
effects in Computational Fluid Dynamics simulations may also be
challenging. The aim of this work is to establish the cross-ventilation
jet behaviour in such complex terrain conditions as well as to provide
guidelines on the implementation of CFD numerical simulations in
order to model complex terrain features such as vegetation in an
efficient manner. The methodology consists of onsite measurements
on a test cell coupled with numerical simulations. It was found
that the cross-ventilation flow is highly turbulent despite the very
low velocities encountered internally within the test cells. While no
direct measurement of the jet direction was made, the measurements
indicate that flow tends to be reversed from the leeward to the
windward side. Modelling such a phenomenon proves challenging
and is strongly influenced by how vegetation is modelled. A solid
vegetation tends to predict better the direction and magnitude of the
flow than a porous vegetation approach. A simplified terrain model
was also shown to provide good comparisons with observation. The
findings have important implications on the study of cross-ventilation
in complex terrain conditions since the flow direction does not remain
trivial, as with the traditional isolated building case.
Abstract: In this paper, we studied the optimal portfolio selection in a defined contribution (DC) pension scheme with multiple contributors under constant elasticity of variance (CEV) model and the impact of stochastic additional voluntary contribution on the investment strategies. We assume that the voluntary contributions are stochastic and also consider investments in a risk free asset and a risky asset to increase the expected returns of the contributing members. We derived a stochastic differential equation which consists of the members’ monthly contributions and the invested fund and obtained an optimized problem with the help of Hamilton Jacobi Bellman equation. Furthermore, we find an explicit solution for the optimal investment strategy with stochastic voluntary contribution using power transformation and change of variables method and the corresponding optimal fund size was obtained. We discussed the impact of the voluntary contribution on the optimal investment strategy with numerical simulations and observed that the voluntary contribution reduces the optimal investment strategy of the risky asset.
Abstract: The influence of the geometric parameters of trapezoidal labyrinth channel on the emitter discharge is investigated in this work. The impact of the dentate angle, the dentate spacing, and the dentate height are studied among the geometric parameters of the labyrinth channel. Numerical simulations of the water flow movement are performed according to central cubic composite design using Commercial codes GAMBIT and FLUENT. Inlet pressure of the dripper is set up to be 1 bar. The objective of this paper is to derive a mathematical model of the emitter discharge depending on the dentate angle, the dentate spacing, the dentate height of the labyrinth channel. As a result, the obtained mathematical model is a second-order polynomial reporting 2-way interactions among the geometric parameters. The dentate spacing has the most important and positive influence on the emitter discharge, followed by the simultaneous impact of the dentate spacing and the dentate height. The dentate angle in the observed interval has no significant effect on the emitter discharge. The obtained model can be used as a basis for a future emitter design.
Abstract: The ultimate load analysis of RC pile groups has assumed a lot of significance under liquefying soil conditions, especially due to post-earthquake studies of 1964 Niigata, 1995 Kobe and 2001 Bhuj earthquakes. The present study reports the results of numerical simulations on pile groups subjected to monotonically increasing lateral loads under design amounts of pile axial loading. The soil liquefaction has been considered through the non-linear p-y relationship of the soil springs, which can vary along the depth/length of the pile. This variation again is related to the liquefaction potential of the site and the magnitude of the seismic shaking. As the piles in the group can reach their extreme deflections and rotations during increased amounts of lateral loading, a precise modeling of the inelastic behavior of the pile cross-section is done, considering the complete stress-strain behavior of concrete, with and without confinement, and reinforcing steel, including the strain-hardening portion. The possibility of the inelastic buckling of the individual piles is considered in the overall collapse modes. The model is analysed using Riks analysis in finite element software to check the post buckling behavior and plastic collapse of piles. The results confirm the kinds of failure modes predicted by centrifuge test results reported by researchers on pile group, although the pile material used is significantly different from that of the simulation model. The extension of the present work promises an important contribution to the design codes for pile groups in liquefying soils.
Abstract: Gas release from the pipelines is one of the main factors in the gas industry accidents. Released gas ejects from the pipeline as a free jet and in the growth process, the fuel gets mixed with the ambient air. Accordingly, an accidental spark will release the chemical energy of the mixture with an explosion. Gas explosion damages the equipment and endangers the life of staffs. So due to importance of safety in gas industries, prevision of accident can reduce the number of the casualties. In this paper, natural gas leakages from the low pressure pipelines are studied in two steps: 1) the simulation of mixing process and identification of flammable zones and 2) the simulation of wind effects on the mixing process. The numerical simulations were performed by using the finite volume method and the pressure-based algorithm. Also, for the grid generation the structured method was used. The results show that, in just 6.4 s after accident, released natural gas could penetrate to 40 m in vertical and 20 m in horizontal direction. Moreover, the results show that the wind speed is a key factor in dispersion process. In fact, the wind transports the flammable zones into the downstream. Hence, to improve the safety of the people and human property, it is preferable to construct gas facilities and buildings in the opposite side of prevailing wind direction.
Abstract: Energy transmission pipelines are one of the most vital parts of each country which several strict laws have been conducted to enhance the safety of these lines and their vicinity. One of these laws is the safety distance around high pressure gas pipelines. Safety distance refers to the minimum distance from the pipeline where people and equipment do not confront with serious damages. In the present study, safety distance around high pressure gas transmission pipelines were determined by using numerical methods. For this purpose, gas leakages from cracked pipeline and created jet fires were simulated as continuous ignition, three dimensional, unsteady and turbulent cases. Numerical simulations were based on finite volume method and turbulence of flow was considered using k-ω SST model. Also, the combustion of natural gas and air mixture was applied using the eddy dissipation method. The results show that, due to the high pressure difference between pipeline and environment, flow chocks in the cracked area and velocity of the exhausted gas reaches to sound speed. Also, analysis of the incident radiation results shows that safety distances around 42 inches high pressure natural gas pipeline based on 5 and 15 kW/m2 criteria are 205 and 272 meters, respectively.
Abstract: Traditional techniques for analyzing time series are based on the notion of stationarity of phenomena under study, but in reality most economic and financial series do not verify this hypothesis, which implies the implementation of specific tools for the detection of such behavior. In this paper, we study nonstationary non-seasonal time series tests in a non-exhaustive manner. We formalize the problem of nonstationary processes with numerical simulations and take stock of their statistical characteristics. The theoretical aspects of some of the most common unit root tests will be discussed. We detail the specification of the tests, showing the advantages and disadvantages of each. The empirical study focuses on the application of these tests to the exchange rate (USD/TND) and the Consumer Price Index (CPI) in Tunisia, in order to compare the Power of these tests with the characteristics of the series.
Abstract: A mechanical wave or vibration propagating through
granular media exhibits a specific signature in time. A coherent
pulse or wavefront arrives first with multiply scattered waves (coda)
arriving later. The coherent pulse is micro-structure independent i.e.
it depends only on the bulk properties of the disordered granular
sample, the sound wave velocity of the granular sample and hence
bulk and shear moduli. The coherent wavefront attenuates (decreases
in amplitude) and broadens with distance from its source. The
pulse attenuation and broadening effects are affected by disorder
(polydispersity; contrast in size of the granules) and have often been
attributed to dispersion and scattering. To study the effect of disorder
and initial amplitude (non-linearity) of the pulse imparted to the
system on the coherent wavefront, numerical simulations have been
carried out on one-dimensional sets of particles (granular chains).
The interaction force between the particles is given by a Hertzian
contact model. The sizes of particles have been selected randomly
from a Gaussian distribution, where the standard deviation of this
distribution is the relevant parameter that quantifies the effect of
disorder on the coherent wavefront. Since, the coherent wavefront is
system configuration independent, ensemble averaging has been used
for improving the signal quality of the coherent pulse and removing
the multiply scattered waves. The results concerning the width of the
coherent wavefront have been formulated in terms of scaling laws. An
experimental set-up of photoelastic particles constituting a granular
chain is proposed to validate the numerical results.
Abstract: To date, high-performance structural steel has been widely used for columns in construction practices due to its significant advantages over conventional steel. However, the same design approach with conventional steel columns is still adopted in the design of high-performance steel columns. As a result, its superior properties cannot be fully considered in design. This paper conducts a test and finite element analysis on the overall stability behaviour of welded Q460GJ steel box columns. In the test, four steel columns with different slenderness and width-to-thickness ratio were compressed under an axial compression testing machine. And finite element models were established in which material nonlinearity and residual stress distributions of test columns were included. Then, comparisons were made between test results and finite element result, it showed that finite element analysis results are agree well with the test result. It means that the test and finite element model are reliable. Then, we compared the test result with the design value calculated by current code, the result showed that Q460GJ steel box columns have the higher overall buckling capacity than the design value. It is necessary to update the design curves for Q460GJ steel columns so that the overall stability capacity of Q460GJ box columns can be designed appropriately.
Abstract: Computations for two-dimensional flow past a stationary and harmonically pitching wind turbine airfoil at a moderate value of Reynolds number (400000) are carried out by progressively increasing the angle of attack for stationary airfoil and at fixed pitching frequencies for rotary one. The incompressible Navier-Stokes equations in conjunction with Unsteady Reynolds Average Navier-Stokes (URANS) equations for turbulence modeling are solved by OpenFOAM package to investigate the aerodynamic phenomena occurred at stationary and pitching conditions on a NACA 6-series wind turbine airfoil. The aim of this study is to enhance the accuracy of numerical simulation in predicting the aerodynamic behavior of an oscillating airfoil in OpenFOAM. Hence, for turbulence modelling, k-ω-SST with low-Reynolds correction is employed to capture the unsteady phenomena occurred in stationary and oscillating motion of the airfoil. Using aerodynamic and pressure coefficients along with flow patterns, the unsteady aerodynamics at pre-, near-, and post-static stall regions are analyzed in harmonically pitching airfoil, and the results are validated with the corresponding experimental data possessed by the authors. The results indicate that implementing the mentioned turbulence model leads to accurate prediction of the angle of static stall for stationary airfoil and flow separation, dynamic stall phenomenon, and reattachment of the flow on the surface of airfoil for pitching one. Due to the geometry of the studied 6-series airfoil, the vortex on the upper surface of the airfoil during upstrokes is formed at the trailing edge. Therefore, the pattern flow obtained by our numerical simulations represents the formation and change of the trailing-edge vortex at near- and post-stall regions where this process determines the dynamic stall phenomenon.
Abstract: The random dither quantization method enables us
to achieve much better performance than the simple uniform
quantization method for the design of quantized control systems.
Motivated by this fact, the stochastic model predictive control
method in which a performance index is minimized subject to
probabilistic constraints imposed on the state variables of systems
has been proposed for linear feedback control systems with random
dither quantization. In other words, a method for solving optimal
control problems subject to probabilistic state constraints for linear
discrete-time control systems with random dither quantization has
been already established. To our best knowledge, however, the
feasibility of such a kind of optimal control problems has not
yet been studied. Our objective in this paper is to investigate the
feasibility of stochastic model predictive control problems for linear
discrete-time control systems with random dither quantization. To
this end, we provide the results of numerical simulations that verify
the feasibility of stochastic model predictive control problems for
linear discrete-time control systems with random dither quantization.
Abstract: This paper presents the Colzate-Vertova landslide, a Deep Seated Gravitational Slope Deformation (DSGSD) located in the Seriana Valley, Northern Italy. The paper aims at describing the development as well as evaluating the factors that influence the evolution of the landslide. After defining the conceptual model of the landslide, numerical simulations were developed using a finite element numerical model, first with a two-dimensional domain, and later with a three-dimensional one. The results of the 2-D model showed a displacement field typical of a sackung, as a consequence of the erosion along the Seriana Valley. The analysis also showed that the groundwater flow could locally affect the slope stability, bringing about a reduction in the safety factor, but without reaching failure conditions. The sensitivity analysis carried out on the strength parameters pointed out that slope failures could be reached only for relevant reduction of the geotechnical characteristics. Such a result does not fit the real conditions observed on site, where a number of small failures often develop all along the hillslope. The 3-D model gave a more comprehensive analysis of the evolution of the DSGSD, also considering the border effects. The results showed that the convex profile of the slope favors the development of displacements along the lateral valley, with a relevant reduction in the safety factor, justifying the existing landslides.