• International Journal of Air-Conditioning and Refrigeration
• /
• v.11 no.1
• /
• pp.1-9
• /
• 2003

The examinations of the operating mechanism of an oscillating capillary tube heat pipe (OCHP) using the visualization method revealed that the working fluid in the OCHP oscillated to the axial direction by the contraction and expansion of vapor plugs. The contraction and expansion were due to the formation and extinction of bubbles in the evaporating and condensing part, respectively The actual physical mechanism, whereby the heat which was transferred in such an OCHP was complex and not well understood. In this study, a theoretical model of the OCHP was developed to model the oscillating motion of working fluid in the OCHP. The differential equations of two-phase flow were applied and simultaneous non-linear partial differential equations were solved. From the analysis of the numerical results, it was found that the oscillating motion Of working fluid in the OCHP was affected by the operation and design conditions such as the heat flux, the charging ratio of working fluid and the hydraulic diameter of flow channel. The simulation results showed that the proposed model and solution could be used for estimating the operating mechanism in the OCHP.

• Structural Engineering and Mechanics
• /
• v.83 no.4
• /
• pp.537-549
• /
• 2022

In this paper, a novel multi-feature model predictive control (MPC) framework with real-time and adaptive performances is proposed for intelligent structural control in which some drawbacks of the algorithm including, complex control rule and non-optimality, are alleviated. Hence, Linear Programming (LP) is utilized to simplify the resulted control rule. Afterward, the Whale Optimization Algorithm (WOA) is applied to the optimal and adaptive tuning of the LP weights independently at each time step. The stochastic control rule is also achieved using Kalman Filter (KF) to handle noisy measurements. The Extreme Learning Machine (ELM) is then adopted to develop a data-driven and real-time control algorithm. The efficiency of the developed algorithm is then demonstrated by numerical simulation of a twenty-story high-rise benchmark building subjected to earthquake excitations. The competency of the proposed method is proven from the aspects of optimality, stochasticity, and adaptivity compared to the KF-based MPC (KMPC) and constrained MPC (CMPC) algorithms in vibration suppression of building structures. The average value for performance indices in the near-field and far-field (El earthquakes demonstrates a reduction up to 38.3% and 32.5% compared with KMPC and CMPC, respectively.

• Journal of ILASS-Korea
• /
• v.28 no.4
• /
• pp.169-175
• /
• 2023

Multi-injection scheme is being applied to GDI combustion to reduce PM and PN emission to meet the EU7 regulation. However, very short injection duration encounters the ballistic injection region, which injection quantity does not increase linearly with injection duration when applying multi-injection. In this study, numerical studies were conducted to reveal the cause of ballistic injection and the effect of design parameters on ballistic region using 1-D simulation, AMESim. Injection rate and injection quantity were compared with experiment to validate the established model, which showed the accuracy with 10% error. The model revealed that the tendency of ballistic region coincides with the needle motion behavior, which means that parameters at the upper part of needle such as electro-magnetic force, needle spring force and needle friction force have dominant effect on ballistic injection. To figure out the effect of electro-magnetic and needle friction force on ballistic, those parameters were varied to plus and minus 10% with model. The result showed that those parameters clearly changed the ballistic region characteristics, however, the impact became insignificant for outside of ballistic region, which means that the ballistic injection is mainly influenced by initial motion of injector needle.

• Earthquakes and Structures
• /
• v.26 no.5
• /
• pp.327-336
• /
• 2024

Evaluation of tunnel performance in seismic-prone areas demands efficient means of estimating performance at different hazard levels. The present study introduces an innovative push-over analysis approach which employs the standard earthquake spectrum to simulate the performance of a tunnel. The numerical simulation has taken into account the lining and surrounding rock to calculate the rock-tunnel interaction subjected to a static push-over displacement regime. Elastic perfectly plastic models for the lining and hardening strain rock medium were used to portray the development of plastic hinges, nonlinear deformation, and performance of the tunnel structure. Separately using a computational algorithm, the non-linear response spectrum was approximated from the average shear strain of the rock model. A NATM tunnel in Turkey was chosen for parametric study. A seismic performance curve and two performance thresholds are introduced that are based on the proposed nonlinear seismic static loading approach and the formation of plastic hinges. The tunnel model was also subjected to a harmonic excitation with a smooth response spectrum and different amplitudes in the fully-dynamic phase to assess the accuracy of the approach. The parametric study investigated the effects of the lining stiffness and capacity and soil stiffness on the seismic performance of the tunnel.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Structural Engineering and Mechanics
• /
• v.44 no.1
• /
• pp.109-125
• /
• 2012

Post-buckling behavior of Timoshenko beams subjected to uniform temperature rising with temperature dependent physical properties are studied in this paper by using the total Lagrangian Timoshenko beam element approximation. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces end therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. The beams considered in numerical examples are made of Austenitic Stainless Steel (316). The convergence studies are made. In this study, the difference between temperature dependent and independent physical properties are investigated in detail in post-buckling case. The relationships between deflections, thermal post-buckling configuration, critical buckling temperature, maximum stresses of the beams and temperature rising are illustrated in detail in post-buckling case.

• Computers and Concrete
• /
• v.4 no.5
• /
• pp.377-402
• /
• 2007

The paper presents quasi-static plane strain FE-simulations of strain localization in reinforced concrete beams without stirrups. The material was modeled with two different isotropic continuum crack models: an elasto-plastic and a damage one. In case of elasto-plasticity, linear Drucker-Prager criterion with a non-associated flow rule was defined in the compressive regime and a Rankine criterion with an associated flow rule was adopted in the tensile regime. In the case of a damage model, the degradation of the material due to micro-cracking was described with a single scalar damage parameter. To ensure the mesh-independence and to capture size effects, both criteria were enhanced in a softening regime by nonlocal terms. Thus, a characteristic length of micro-structure was included. The effect of a characteristic length, reinforcement ratio, bond-slip stiffness, fracture energy and beam size on strain localization was investigated. The numerical results with reinforced concrete beams were quantitatively compared with corresponding laboratory tests by Walraven (1978).

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.94-103
• /
• 2015

The turbulent flow characteristics in the channel flow are investigated using large eddy simulation(LES) of FDS code, built in NIST(USA), in which the near-wall flow is solved by Werner-Wengle wall function. The periodic flow condition is applied in streamwise direction to get the fully developed turbulent flow and symmetric condition is applied in lateral direction. The height of the channel is H=1m, and the length of the channel is 6H, and the lateral length is H. The total grid is $32{\times}32{\times}32$ and $y^+$ is kept above 11 to fulfill the near-wall flow requirement. The Smagorinsky model is used to solve the sub-grid scale stress. Smagorinsky constant $C_s$ is 0.2(default in FDS). Three cases of Reynolds number(10,700, 26,000, 49,000.), based on the channel height, are analyzed. The simulated results are compared with direct numerical simulation(DNS) and particle image velocimetry(PIV) experimental data. The linear low-Re eddy viscosity model of Launder & Sharma and non-linear low-Re eddy viscosity model of Abe-Jang-Leschziner are utilized to compare the results with LES of FDS. Reynolds normal stresses, Reynolds shear stresses, turbulent kinetic energys and mean velocity flows are well compared with DNS and PIV data.

• Advances in Computational Design
• /
• v.3 no.3
• /
• pp.233-267
• /
• 2018

Conventional seismic design of concentrically braced frame (CBF) structures suggests that the gusset plate connecting a steel brace to beams and/or columns should be designed as non-dissipative in earthquakes, while the steel brace members should be designed as dissipative elements. These design intentions lead to thicker and larger gusset plates in design on one hand and a potentially under-rated contribution of gusset plates in design, on the other hand. In contrast, research has shown that compact and thinner gusset plates designed in accordance with the elliptical clearance method rather than the conventional standard linear clearance method can enhance system ductility and energy dissipation capacity in concentrically braced steel frames. In order to assess the two design methods, six cyclic push-over tests on full scale models of concentric braced steel frame structures were conducted. Furthermore, a 3D finite element (FE) shell model, incorporating state-of-the-art tools and techniques in numerical simulation, was developed that successfully replicates the response of gusset plate and bracing members under fully reversed cyclic axial loading. Direct measurements from strain gauges applied to the physical models were used primarily to validate FE models, while comparisons of hysteresis load-displacement loops from physical and numerical models were used to highlight the overall performance of the FE models. The study shows the two design methods attain structural response as per the design intentions; however, the elliptical clearance method has a superiority over the standard linear method as a fact of improving detailing of the gusset plates, enhancing resisting capacity and improving deformability of a CBF structure. Considerations were proposed for improvement of guidelines for detailing gusset plates and bracing members in CBF structures.

• Proceedings of the Korean Powder Metallurgy Institute Conference
• /
• 1997.10a
• /
• pp.7-7
• /
• 1997

During sintering of very porous green bodies, as obtained by compaction of hard powders - such as tungsten carbide or ceramics - or by injection moulding, important shrinkage occurs. Due to heterogeneous green density field, gravity effects, friction on the support, thermal gradients, etc., this shrinkage is often non-uniform, which' may induce significant shape changes. As the ratio of compact dimension to powder size is very high, the mechanics of continuum is relevant to model such phenomena. Thus numerical techniques, such as the finite element method can be used to simulate the sintering process and predict the final shape of the sintered part. Such type of simulation has much been developed in the last decade firstly for hot isostatic pressing and next for die compaction. Finite element modelling has been recently applied to free sintering. The simulation of sintering should be based on constitutive equations describing the thermo-mechanical behaviour of the material under any state of stress and any temperature which may arise within the sintering body. These equations can be drawn either from experimental data or from micromechanical models. The experiments usually consist in free sintering and sinter-forging tests. Indeed applying more complex loading conditions at high temperature under controlled atmosphere is delicate. Micromechanical models describe the constitutive behaviour of aggregates of spheres from the deformation of two-sphere contact either by viscous flow or grain boundary diffusion. Such models are not able to describe complex microstructure and mechanisms as observed in real materials but they can give some basic information on the formulation of constitutive equations. Practically both experimental and theoretical approaches can be coupled to identify the constitutive equations. Such procedure has been performed for modelling the sintering of compacts obtained by die pressing of a mixture of tungsten carbide and cobalt powders. The constitutive behaviour of this material during sintering has been described by a linear viscous constitutive model, whose functions have been fitted from results of free sintering and sinter-forging experiments. This model has next been introduced in ABAQUS finite element code to simulate the sintering of heterogeneous green compacts of various geometries at constant temperature. Examples of simulations are shown and compared with experiments.