• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.

• Coupled systems mechanics
• /
• v.9 no.2
• /
• pp.125-145
• /
• 2020

In this paper, we present a 3D thermo-hydro-mechanical coupled discrete beam lattice model of structure built of the nonisothermal saturated poro-plastic medium subjected to mechanical loads and nonstationary heat transfer conditions. The proposed model is based on Voronoi cell representation of the domain with cohesive links represented as inelastic Timoshenko beam finite elements enhanced with additional kinematics in terms of embedded strong discontinuities in axial and both transverse directions. The enhanced Timoshenko beam finite element is capable of modeling crack formation in mode I, mode II and mode III. Mode I relates to crack opening, mode II relates to in-plane crack sliding, and mode III relates to the out-of-plane shear sliding. The pore fluid flow and heat flow in the proposed model are governed by Darcy's law and Fourier's law for heat conduction, respectively. The pore pressure field and temperature field are approximated with linear tetrahedral finite elements. By exploiting nodal point quadrature rule for numerical integration on tetrahedral finite elements and duality property between Voronoi diagram and Delaunay tetrahedralization, the numerical implementation of the coupling results with additional pore pressure and temperature degrees of freedom placed at each node of a Timoshenko beam finite element. The results of several numerical simulations are presented and discussed.

• International Journal of Aeronautical and Space Sciences
• /
• v.17 no.3
• /
• pp.423-431
• /
• 2016

In this study, wing design optimization for long-endurance unmanned aerial vehicles (UAVs) is investigated. The fluid-structure integration (FSI) analysis is carried out to simulate the aeroelastic characteristics of a high-aspect ratio wing for a long-endurance UAV. High-fidelity computational codes, FLUENT and DIAMOND/IPSAP, are employed for the loose coupling FSI optimization. In addition, this optimization procedure is improved by adopting the design of experiment (DOE) and Kriging model. A design optimization tool, PIAnO, integrates with an in-house codes, CAE simulation and an optimization process for generating the wing geometry/computational mesh, transferring information, and finding the optimum solution. The goal of this optimization is to find the best high-aspect ratio wing shape that generates minimum drag at a cruise condition of $C_L=1.0$. The result shows that the optimal wing shape produced 5.95 % less drag compared to the initial wing shape.

• Wind and Structures
• /
• v.13 no.2
• /
• pp.169-189
• /
• 2010

This paper describes the use of coupled Computational Fluid Dynamics (CFD) and Rigid Body Dynamics (RBD) in modelling the aerodynamic behaviour of wind-borne plate type objects. Unsteady 2D and 3D Reynolds Averaged Navier-Stokes (RANS) CFD models are used to simulate the unsteady and non-uniform flow field surrounding static, forced rotating, auto-rotating and free-flying plates. The auto-rotation phenomenon itself is strongly influenced by vortex shedding, and the realisable k-epsilon turbulence modelling approach is used, with a second order implicit time advancement scheme and equal or higher order advection schemes for the flow variables. Sequentially coupling the CFD code with a RBD solver allows a more detailed modelling of the Fluid-Structure Interaction (FSI) behaviour of the plate and how this influences plate motion. The results are compared against wind tunnel experiments on auto-rotating plates and an existing 3D analytical model.

• Transactions of Materials Processing
• /
• v.25 no.1
• /
• pp.49-55
• /
• 2016

Electrohydraulic forming (EHF) is a high-speed forming process that uses an electric arc discharge in water. Shock waves resulting from the electric arc discharge are propagated to the blank through water and the blank moves toward the die. Advantages of EHF include improved formability due to the high-speed process and reduction of the bouncing effect. In the current study, a numerical simulation of EHF was developed using LS-DYNA. In the simulation, the model for the electric arc was assumed as an adiabatic gas expansion and an Arbitrary Lagrange-Eulerian (ALE) multi material formulation was used to describe the interaction between the electric arc and the water. In order to model the Fluid-Structure Interaction (FSI), a coupling mechanism was used. The blank of Al 1100-O was simulated using shell elements. The results of the simulation showed that the blank was deformed due to the pressure propagation of water and the bouncing effect did not affect the formability of blank.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.4
• /
• pp.1024-1040
• /
• 2014

The theoretical background and a numerical solution procedure for a time domain hydroelastic code are presented in this paper. The code combines a VOF-based free surface flow solver with a flexible body motion solver where the body linear elastic deformation is described by a modal superposition of dry mode shapes expressed in a local floating frame of reference. These mode shapes can be obtained from any finite element code. The floating frame undergoes a pseudo rigid-body motion which allows for a large rigid body translation and rotation and fully preserves the coupling with the local structural deformation. The formulation relies on the ability of the flow solver to provide the total fluid action on the body including e.g. the viscous forces, hydrostatic and hydrodynamic forces, slamming forces and the fluid damping. A numerical simulation of a flexible barge is provided and compared to experiments to show that the VOF-based flow solver has this ability and the code has the potential to predict the global hydroelastic responses accurately.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.1
• /
• pp.72-78
• /
• 2008

In this study, advanced (fluid-structure interaction (FSI)) analysis system has been developed in order to predict turbine cascade performance with blade deformation effect due to aerodynamic loads. Intereference effects due to the relative movement of the rotor cascade with respect to the stator cascade are also considered. Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation k-ω SST turbulence models are solved to accurately predict fluid dynamic loads considering flow separation effects. A fully implicit time marching scheme based on the (coupled Newmark time-integration method) with high artificial damping is efficiently used to compute the complex fluid-structure interaction problem. Predicted aerodynamic performance considering structural deformation effect of the blade shows somewhat different results compared to the case of rigid blade model. Cascade performance evaluations for different elastic axis positions are importantly presented and its aeroelastic effects are investigated.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2011.02a
• /
• pp.21-21
• /
• 2011

A series of Quasi-Elastic Neutron Scattering (QENS) experiments helps us to understand the single-particle (hydrogen atom) dynamics of a globular protein and its hydration water and strong coupling between them. We also performed Molecular Dynamics (MD) simulations on a realistic model of the hydrated hen-egg Lysozyme powder having two proteins in the periodic box. We found the existence of a Fragile-to-Strong dynamic Crossover (FSC) phenomenon in hydration water around a protein occurring at TL=$225{\pm}5K$ by analyzing Intermediate Scattering Function (ISF). On lowering of the temperature toward FSC, the structure of hydration water makes a transition from predominantly the High Density Liquid (HDL) form, a more fluid state, to predominantly the Low Density Liquid (LDL) form, a less fluid state, derived from the existence of a liquid?liquid critical point at an elevated pressure. We showed experimentally and confirmed theoretically that this sudden switch in the mobility of the hydration water around a protein triggers the dynamic transition (so-called glass transition) of the protein, at a temperature TD=220 K. Mean Square Displacement (MSD) is the important factor to show that the FSC is the key to the strong coupling between a protein and its hydration water by suggesting TL${\fallingdotseq}$TD. MD simulations with TIP4P force field for water were performed to understand hydration level dependency of the FSC temperature. We added water molecules to increase hydration level of the protein hydration water, from 0.30, 0.45, 0.60 and 1.00 (1.00 is the bulk water). These confirm the existence of the FSC and the hydration level dependence of the FSC temperature: FSC temperature is decreased upon increasing hydration level. We compared the hydration water around Lysozyme, B-DNA and RNA. Similarity among those suggests that the FSC and this coupling be universal for globular proteins, biopolymers.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.16 no.2
• /
• pp.58-65
• /
• 2008

The finite element programs have been developed for structure, collision, flow, dynamics, heat transfer, acoustics, electromagnetism, MEMS (Micro Electro Mechanical Systems), and etc. These programs can be classified as either "package" program or "single purpose" program. Single purpose programs usually have convenient and powerful functions, but these programs have limited expandability to different fields of analysis. Therefore, the method to converter the analysis results of single purpose program to other programs is needed. In the research, multi-body data conversion methods of 1) finite element model and 2) solid model were created to convert fluid analysis result of CFD-ACE+ to ANSYS data structure. Automatic boundary condition algorithms were developed for blade, and finite element model was compared with solid model. It is expected that, by sealess data transfer, the Multi-body Data Conversion Program could reduce the development period of torque converters.

• Journal of Power System Engineering
• /
• v.6 no.3
• /
• pp.57-63
• /
• 2002

This paper describes the acoustic analysis of mid duty truck. The focus of the analysis is on structure borne engine noise with major contributions of 2nd order. It has been previously recognized that the noise contribution of each transfer path of structure borne noise can be varied with the charateristics of each mounts and vibro acoustic sensitivity of car body. The structure of car body will be split up into three major sub components, which are modeled separately, the engine, the frame and the cab. The acoustic performance is evaluated on three levels: engine to frame transfer, frame to cab transfer, and panel contribution from cab to driver. In order to perform these analyses, analytical models are created for the engine, frame, cab and acoustic cavity. The models are linked through a coupled fluid structure calculation, and through FRF Based Substructuring for the structural couplings. Based on the structural coupling calculations, a transfer path analysis is performed to identify the most important transfer paths. These paths are then the focussing points for applying modifications to the structure or the mount system. Finally, a number of modification are proposed and their effect is quantified.