Abstract: Numerical computation of wave propagation in a large
domain usually requires significant computational effort. Hence, the
considered domain must be truncated to a smaller domain of interest.
In addition, special boundary conditions, which absorb the outward
travelling waves, need to be implemented in order to describe the
system domains correctly. In this work, the linear one dimensional
wave equation is approximated by utilizing the Fourier Galerkin
approach. Furthermore, the artificial boundaries are realized with
absorbing boundary conditions. Within this work, a systematic work
flow for setting up the wave problem, including the absorbing
boundary conditions, is proposed. As a result, a convenient modal
system description with an effective absorbing boundary formulation
is established. Moreover, the truncated model shows high accuracy
compared to the global domain.
Abstract: Conditions corresponding to the unconditional stability
of convection in a mechanically anisotropic fluid saturated porous
medium of infinite horizontal extent are determined. The medium
is heated from below and its bounding surfaces are subjected to
temperature modulation which consists of a steady part and a
time periodic oscillating part. The Brinkman model is employed
in the momentum equation with the Bousinessq approximation.
The stability region is found for arbitrary values of modulational
frequency and amplitude using the energy method. Higher order
numerical computations are carried out to find critical boundaries
and subcritical instability regions more accurately.
Abstract: Migration of a core-shell soft particle under the
influence of an external electric field in an electrolyte solution is
studied numerically. The soft particle is coated with a positively
charged polyelectrolyte layer (PEL) and the rigid core is having
a uniform surface charge density. The Darcy-Brinkman extended
Navier-Stokes equations are solved for the motion of the ionized
fluid, the non-linear Nernst-Planck equations for the ion transport and
the Poisson equation for the electric potential. A pressure correction
based iterative algorithm is adopted for numerical computations. The
effects of convection on double layer polarization (DLP) and diffusion
dominated counter ions penetration are investigated for a wide range
of Debye layer thickness, PEL fixed surface charge density, and
permeability of the PEL. Our results show that when the Debye
layer is in order of the particle size, the DLP effect is significant
and produces a reduction in electrophoretic mobility. However, the
double layer polarization effect is negligible for a thin Debye layer
or low permeable cases. The point of zero mobility and the existence
of mobility reversal depending on the electrolyte concentration are
also presented.
Abstract: The nonlinear time history analysis of seismically base-isolated structures can require a significant computational effort when the behavior of each seismic isolator is predicted by adopting the widely used differential equation Bouc-Wen model. In this paper, a nonlinear exponential model, able to simulate the response of seismic isolation bearings within a relatively large displacements range, is described and adopted in order to reduce the numerical computations and speed up the nonlinear dynamic analysis. Compared to the Bouc-Wen model, the proposed one does not require the numerical solution of a nonlinear differential equation for each time step of the analysis. The seismic response of a 3d base-isolated structure with a lead rubber bearing system subjected to harmonic earthquake excitation is simulated by modeling each isolator using the proposed analytical model. The comparison of the numerical results and computational time with those obtained by modeling the lead rubber bearings using the Bouc-Wen model demonstrates the good accuracy of the proposed model and its capability to reduce significantly the computational effort of the analysis.
Abstract: In order to reduce numerical computations in the
nonlinear dynamic analysis of seismically base-isolated structures, a
Mixed Explicit-Implicit time integration Method (MEIM) has been
proposed. Adopting the explicit conditionally stable central
difference method to compute the nonlinear response of the base
isolation system, and the implicit unconditionally stable Newmark’s
constant average acceleration method to determine the superstructure
linear response, the proposed MEIM, which is conditionally stable
due to the use of the central difference method, allows to avoid the
iterative procedure generally required by conventional monolithic
solution approaches within each time step of the analysis. The main
aim of this paper is to investigate the stability and computational
efficiency of the MEIM when employed to perform the nonlinear
time history analysis of base-isolated structures with sliding bearings.
Indeed, in this case, the critical time step could become smaller than
the one used to define accurately the earthquake excitation due to the
very high initial stiffness values of such devices. The numerical
results obtained from nonlinear dynamic analyses of a base-isolated
structure with a friction pendulum bearing system, performed by
using the proposed MEIM, are compared to those obtained adopting a
conventional monolithic solution approach, i.e. the implicit
unconditionally stable Newmark’s constant acceleration method
employed in conjunction with the iterative pseudo-force procedure.
According to the numerical results, in the presented numerical
application, the MEIM does not have stability problems being the
critical time step larger than the ground acceleration one despite of
the high initial stiffness of the friction pendulum bearings. In
addition, compared to the conventional monolithic solution approach,
the proposed algorithm preserves its computational efficiency even
when it is adopted to perform the nonlinear dynamic analysis using a
smaller time step.
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: 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: In this work we use the Discrete Proper Orthogonal Decomposition transform to characterize the properties of coupled dynamics in thin-walled beams by exploiting numerical simulations obtained from finite element simulations. The outcomes of the will improve our understanding of the linear and nonlinear coupled behavior of thin-walled beams structures. Thin-walled beams have widespread usage in modern engineering application in both large scale structures (aeronautical structures), as well as in nano-structures (nano-tubes). Therefore, detailed knowledge in regard to the properties of coupled vibrations and buckling in these structures are of great interest in the research community. Due to the geometric complexity in the overall structure and in particular in the cross-sections it is necessary to involve computational mechanics to numerically simulate the dynamics. In using numerical computational techniques, it is not necessary to over simplify a model in order to solve the equations of motions. Computational dynamics methods produce databases of controlled resolution in time and space. These numerical databases contain information on the properties of the coupled dynamics. In order to extract the system dynamic properties and strength of coupling among the various fields of the motion, processing techniques are required. Time- Proper Orthogonal Decomposition transform is a powerful tool for processing databases for the dynamics. It will be used to study the coupled dynamics of thin-walled basic structures. These structures are ideal to form a basis for a systematic study of coupled dynamics in structures of complex geometry.
Abstract: In this paper, the effect of the intake temperature (IT) and intake pressure (IP) on ignition timing and pollutants emission of Homogeneous Charge Compression Ignition (HCCI) engine is investigated. Numerical computations are performed using the CHEMKIN computer code. The numerical temperature obtained using different boundary conditions is compared to published data and a good agreement is assigned. Results show that the HCCI combustion engine is significantly improved by increasing the IT. With a value of IT lower than 390 K, combustion cannot occur. However, with an IT greater than 420 K, the cylinder pressure decreases. An optimum crank rotation angle is achieved by using IT of 420 K. So, we can conclude that the variation of the IT and IP influence notably the emission concentration.
Abstract: In recent decades, flapping wing aerodynamics has attracted great interest. Understanding the physics of biological flyers such as birds and insects can help improve the performance of micro air vehicles. The present research focuses on the aerodynamics of insect-like flapping wing flight with the approach of numerical computation. Insect model of hawkmoth is adopted in the numerical study with rigid wing assumption currently. The numerical model integrates the computational fluid dynamics of the flow and active control of wing kinematics to achieve stable flight. The computation grid is a hybrid consisting of background Cartesian nodes and clouds of mesh-free grids around immersed boundaries. The generalized finite difference method is used in conjunction with single value decomposition (SVD-GFD) in computational fluid dynamics solver to study the dynamics of a free hovering hummingbird hawkmoth. The longitudinal dynamics of the hovering flight is governed by three control parameters, i.e., wing plane angle, mean positional angle and wing beating frequency. In present work, a PID controller works out the appropriate control parameters with the insect motion as input. The controller is adjusted to acquire desired maneuvering of the insect flight. The numerical scheme in present study is proven to be accurate and stable to simulate the flight of the hummingbird hawkmoth, which has relatively high Reynolds number. The PID controller is responsive to provide feedback to the wing kinematics during the hovering flight. The simulated hovering flight agrees well with the real insect flight. The present numerical study offers a promising route to investigate the free flight aerodynamics of insects, which could overcome some of the limitations of experiments.
Abstract: On this research, experimental analyses have been performed in order to determine the oil pump mechanism dynamics and stability from an oil unit mechanical structure. The experimental tests were focused on the vibrations which occur inside of the rod element during functionality of the oil pump unit. The oil pump mechanism dynamic parameters were measured and also determined through numerical computations. Entire research is based on the oil pump unit mechanical system virtual prototyping. For a complete analysis of the mechanism, the frequency dynamic response was identified, mainly for the mechanism driven element, based on two methods: processing and virtual simulations with MSC Adams aid and experimental analysis. In fact, through this research, a complete methodology is presented where numerical simulations of a mechanism with deformed elements are developed on a dynamic mode and these can be correlated with experimental tests.
Abstract: Three-dimensional incompressible turbulent fluid flow and heat transfer of pin fin heat sinks using air as a cooling fluid are numerically studied in this study. Two different kinds of pin fins are compared in the thermal performance, including circular and square cross sections, both are in-line and staggered arrangements. The turbulent governing equations are solved using a control-volume- based finite-difference method. Subsequently, numerical computations are performed with the realizable k - ԑ turbulence for the parameters studied, the fin height H, fin diameter D, and Reynolds number (Re) in the range of 7 ≤ H ≤ 10, 0.75 ≤ D ≤ 2, 2000 ≤ Re ≤ 126000 respectively. The numerical results are validated with available experimental data in the literature and good agreement has been found. It indicates that circular pin fins are streamlined in comparing with the square pin fins, the pressure drop is small than that of square pin fins, and heat transfer is not as good as the square pin fins. The thermal performance of the staggered pin fins is better than that of in-line pin fins because the staggered arrangements produce large disturbance. Both in-line and staggered arrangements show the same behavior for thermal resistance, pressure drop, and the entropy generation.
Abstract: This paper develops a meshless approach, called Element Free Galerkin (EFG) method, which is based on the weak form Moving Least Squares (MLS) of the partial differential governing equations and employs the interpolation to construct the meshless shape functions. The variation weak form is used in the EFG where the trial and test functions are approximated bye the MLS approximation. Since the shape functions constructed by this discretization have the weight function property based on the randomly distributed points, the essential boundary conditions can be implemented easily. The local weak form of the partial differential governing equations is obtained by the weighted residual method within the simple local quadrature domain. The spline function with high continuity is used as the weight function. The presently developed EFG method is a truly meshless method, as it does not require the mesh, either for the construction of the shape functions, or for the integration of the local weak form. Several numerical examples of two-dimensional static structural analysis are presented to illustrate the performance of the present EFG method. They show that the EFG method is highly efficient for the implementation and highly accurate for the computation. The present method is used to analyze the static deflection of beams and plate hole
Abstract: In designing a low-energy-consuming buildings, the heat transfer through a large glass or wall becomes critical. Multiple layers of the window glasses and walls are employed for the high insulation. The gravity driven air flow between window glasses or wall layers is a natural heat convection phenomenon being a key of the heat transfer. For the first step of the natural heat transfer analysis, in this study the development and application of a finite volume method for the numerical computation of viscous incompressible flows is presented. It will become a part of the natural convection analysis with high-order scheme, multi-grid method, and dual-time step in the future. A finite volume method based on a fully-implicit second-order is used to discretize and solve the fluid flow on unstructured grids composed of arbitrary-shaped cells. The integrations of the governing equation are discretised in the finite volume manner using a collocated arrangement of variables. The convergence of the SIMPLE segregated algorithm for the solution of the coupled nonlinear algebraic equations is accelerated by using a sparse matrix solver such as BiCGSTAB. The method used in the present study is verified by applying it to some flows for which either the numerical solution is known or the solution can be obtained using another numerical technique available in the other researches. The accuracy of the method is assessed through the grid refinement.
Abstract: To ensure the gas transmittal GCU's efficient operation, leakages through the labyrinth packings (LP) should be minimized. Leakages can be minimized by decreasing the LP gap, which in turn depends on thermal processes and possible rotor vibrations and is designed to ensure absence of mechanical contact. Vibration mitigation allows to minimize the LP gap. It is advantageous to research influence of processes in the dynamic gas-structure system on LP vibrations. This paper considers influence of rotor vibrations on LP gas dynamics and influence of the latter on the rotor structure within the FSI unidirectional dynamical coupled problem. Dependences of nonstationary parameters of gas-dynamic process in LP on rotor vibrations under various gas speeds and pressures, shaft rotation speeds and vibration amplitudes, and working medium features were studied. The programmed multi-processor ANSYS CFX was chosen as a numerical computation tool. The problem was solved using PNRPU high-capacity computer complex. Deformed shaft vibrations are replaced with an unyielding profile that moves in the fixed annulus "up-and-down" according to set harmonic rule. This solves a nonstationary gas-dynamic problem and determines time dependence of total gas-dynamic force value influencing the shaft. Pressure increase from 0.1 to 10 MPa causes growth of gas-dynamic force oscillation amplitude and frequency. The phase shift angle between gas-dynamic force oscillations and those of shaft displacement decreases from 3π/4 to π/2. Damping constant has maximum value under 1 MPa pressure in the gap. Increase of shaft oscillation frequency from 50 to 150 Hz under P=10 MPa causes growth of gas-dynamic force oscillation amplitude. Damping constant has maximum value at 50 Hz equaling 1.012. Increase of shaft vibration amplitude from 20 to 80 µm under P=10 MPa causes the rise of gas-dynamic force amplitude up to 20 times. Damping constant increases from 0.092 to 0.251. Calculations for various working substances (methane, perfect gas, air at 25 ˚С) prove the minimum gas-dynamic force persistent oscillating amplitude under P=0.1 MPa being observed in methane, and maximum in the air. Frequency remains almost unchanged and the phase shift in the air changes from 3π/4 to π/2. Calculations for various working substances (methane, perfect gas, air at 25 ˚С) prove the maximum gas-dynamic force oscillating amplitude under P=10 MPa being observed in methane, and minimum in the air. Air demonstrates surging. Increase of leakage speed from 0 to 20 m/s through LP under P=0.1 MPa causes the gas-dynamic force oscillating amplitude to decrease by 3 orders and oscillation frequency and the phase shift to increase 2 times and stabilize. Increase of leakage speed from 0 to 20 m/s in LP under P=1 MPa causes gas-dynamic force oscillating amplitude to decrease by almost 4 orders. The phase shift angle increases from π/72 to π/2. Oscillations become persistent. Flow rate proved to influence greatly on pressure oscillations amplitude and a phase shift angle. Work medium influence depends on operation conditions. At pressure growth, vibrations are mostly affected in methane (of working substances list considered), and at pressure decrease, in the air at 25 ˚С.
Abstract: The modelling of physical phenomena, such as the
earth’s free oscillations, the vibration of strings, the interaction of
atomic particles, or the steady state flow in a bar give rise to Sturm-
Liouville (SL) eigenvalue problems. The boundary applications of
some systems like the convection-diffusion equation, electromagnetic
and heat transfer problems requires the combination of Dirichlet and
Neumann boundary conditions. Hence, the incorporation of Robin
boundary condition in the analyses of Sturm-Liouville problem. This
paper deals with the computation of the eigenvalues and
eigenfunction of generalized Sturm-Liouville problems with Robin
boundary condition using the finite element method. Numerical
solution of classical Sturm–Liouville problem is presented. The
results show an agreement with the exact solution. High results
precision is achieved with higher number of elements.
Abstract: In this paper numerical studies have been carried out
to examine the pre-ignition flow features of high-performance solid
propellant rocket motors with two different port geometries but with
same propellant loading density. Numerical computations have been
carried out using a validated 3D, unsteady, 2nd-order implicit, SST k-
ω turbulence model. In the numerical study, a fully implicit finite
volume scheme of the compressible, Reynolds-Averaged, Navier-
Stokes equations is employed. We have observed from the numerical
results that in solid rocket motors with highly loaded propellants
having divergent port geometry the hot igniter gases can create preignition
pressure oscillations leading to thrust oscillations due to the
flow unsteadiness and recirculation. We have also observed that the
igniter temperature fluctuations are diminished rapidly thereby
reaching the steady state value faster in the case of solid propellant
rocket motors with convergent port than the divergent port
irrespective of the igniter total pressure. We have concluded that the
prudent selection of the port geometry, without altering the propellant
loading density, for damping the total temperature fluctuations within
the motor is a meaningful objective for the suppression and control of
instability and/or thrust oscillations often observed in solid propellant
rocket motors with non-uniform port geometry.
Abstract: In this paper comprehensive studies have been carried
out for the design optimization of a waste heat recovery system for
effectively utilizing the domestic air conditioner heat energy for
producing hot water. Numerical studies have been carried for the
geometry optimization of a waste heat recovery system for domestic
air conditioners. Numerical computations have been carried out using
a validated 2d pressure based, unsteady, 2nd-order implicit, SST k-ω
turbulence model. In the numerical study, a fully implicit finite
volume scheme of the compressible, Reynolds-Averaged, Navier-
Stokes equations is employed. At identical inflow and boundary
conditions various geometries were tried and effort has been taken for
proposing the best design criteria. Several combinations of pipe line
shapes viz., straight and spiral with different number of coils for the
radiator have been attempted and accordingly the design criteria has
been proposed for the waste heat recovery system design. We have
concluded that, within the given envelope, the geometry optimization
is a meaningful objective for getting better performance of waste heat
recovery system for air conditioners.
Abstract: This paper deals with the study of reflection and
transmission characteristics of acoustic waves at the interface of a
semiconductor half-space and elastic solid. The amplitude ratios
(reflection and transmission coefficients) of reflected and transmitted
waves to that of incident wave varying with the incident angles have
been examined for the case of quasi-longitudinal wave. The special
cases of normal and grazing incidence have also been derived with
the help of Gauss elimination method. The mathematical model
consisting of governing partial differential equations of motion and
charge carriers’ diffusion of n-type semiconductors and elastic solid
has been solved both analytically and numerically in the study. The
numerical computations of reflection and transmission coefficients
has been carried out by using MATLAB programming software for
silicon (Si) semiconductor and copper elastic solid. The computer
simulated results have been plotted graphically for Si
semiconductors. The study may be useful in semiconductors,
geology, and seismology in addition to surface acoustic wave (SAW)
devices.
Abstract: In this report we have discussed the theoretical aspects
of the flow transformation, occurring through a series of bifurcations.
The parameters and their continuous diversion, the intermittent bursts
in the transition zone, variation of velocity and pressure with time,
effect of roughness in turbulent zone, and changes in friction factor
and head loss coefficient as a function of Reynolds number for a
transverse flow across a cylinder have been discussed. An analysis of
the variation in the wake length with Reynolds number was done in
FORTRAN.