Abstract: The influence of injector attitude on wall heat flux plays an important role in predicting the start-up transient and also determining the combustion chamber wall durability of liquid rockets. In this paper comprehensive numerical studies have been carried out on an idealized liquid rocket combustion chamber to examine the transient wall heat flux during its start-up transient at different injector attitude. Numerical simulations have been carried out with the help of a validated 2d axisymmetric, double precision, pressure-based, transient, species transport, SST k-omega model with laminar finite rate model for governing turbulent-chemistry interaction for four cases with different jet intersection angles, viz., 0o, 30o, 45o, and 60o. We concluded that the jets intersection angle is having a bearing on the time and location of the maximum wall-heat flux zone of the liquid rocket combustion chamber during the start-up transient. We also concluded that the wall heat flux mapping in liquid rocket combustion chamber during the start-up transient is a meaningful objective for the chamber wall material selection and the lucrative design optimization of the combustion chamber for improving the payload capability of the rocket.
Abstract: We present in this paper a fully implicit finite element
method tailored for the numerical modeling of inextensible fluidic
membranes in a surrounding Newtonian fluid. We consider a highly
simplified version of the Canham-Helfrich model for phospholipid
membranes, in which the bending force and spontaneous curvature
are disregarded. The coupled problem is formulated in a fully
Eulerian framework and the membrane motion is tracked using
the level set method. The resulting nonlinear problem is solved
by a Newton-Raphson strategy, featuring a quadratic convergence
behavior. A monolithic solver is implemented, and we report several
numerical experiments aimed at model validation and illustrating
the accuracy of the proposed method. We show that stability is
maintained for significantly larger time steps with respect to an
explicit decoupling method.
Abstract: The importance of fibre reinforced plastics continually
increases due to the excellent mechanical properties, low material
and manufacturing costs combined with significant weight reduction.
Today, components are usually designed and calculated numerically
by using finite element methods (FEM) to avoid expensive laboratory
tests. These programs are based on material models including
material specific deformation characteristics. In this research project,
material models for short glass fibre reinforced plastics are presented
to simulate the visco-elasto-plastic deformation behaviour. Prior
to modelling specimens of the material EMS Grivory HTV-5H1,
consisting of a Polyphthalamide matrix reinforced by 50wt.-% of
short glass fibres, are characterized experimentally in terms of
the highly time dependent deformation behaviour of the matrix
material. To minimize the experimental effort, the cyclic deformation
behaviour under tensile and compressive loading (R = −1) is
characterized by isothermal complex low cycle fatigue (CLCF)
tests. Combining cycles under two strain amplitudes and strain
rates within three orders of magnitude and relaxation intervals
into one experiment the visco-elastic deformation is characterized.
To identify visco-plastic deformation monotonous tensile tests
either displacement controlled or strain controlled (CERT) are
compared. All relevant modelling parameters for this complex
superposition of simultaneously varying mechanical loadings are
quantified by these experiments. Subsequently, two different material
models are compared with respect to their accuracy describing the
visco-elasto-plastic deformation behaviour. First, based on Chaboche
an extended 12 parameter model (EVP-KV2) is used to model cyclic
visco-elasto-plasticity at two time scales. The parameters of the
model including a total separation of elastic and plastic deformation
are obtained by computational optimization using an evolutionary
algorithm based on a fitness function called genetic algorithm.
Second, the 12 parameter visco-elasto-plastic material model by
Launay is used. In detail, the model contains a different type of a
flow function based on the definition of the visco-plastic deformation
as a part of the overall deformation. The accuracy of the models is
verified by corresponding experimental LCF testing.
Abstract: There exist a wide range of failure modes in composite
structures due to the increased usage of the structures especially in
aerospace industry. Moreover, temperature dependent wave response
of composite and layered structures have been continuously studied,
though still limited, in the last decade mainly due to the broad
operating temperature range of aerospace structures. A wave finite
element (WFE) and finite element (FE) based computational method
is presented by which the temperature dependent wave dispersion
characteristics and interaction phenomenon in composite structures
can be predicted. Initially, the temperature dependent mechanical
properties of the panel in the range of -100 ◦C to 150 ◦C are
measured experimentally using the Thermal Mechanical Analysis
(TMA). Temperature dependent wave dispersion characteristics of
each waveguide of the structural system, which is discretized as a
system of a number of waveguides coupled by a coupling element, is
calculated using the WFE approach. The wave scattering properties,
as a function of temperature, is determined by coupling the WFE
wave characteristics models of the waveguides with the full FE
modelling of the coupling element on which defect is included.
Numerical case studies are exhibited for two waveguides coupled
through a coupling element.
Abstract: A wave finite element (WFE) and finite element
(FE) based computational method is presented by which the
dispersion properties as well as the wave interaction coefficients for
one-dimensional structural system can be predicted. The structural
system is discretized as a system comprising a number of waveguides
connected by a coupling joint. Uniform nodes are ensured at the
interfaces of the coupling element with each waveguide. Then,
equilibrium and continuity conditions are enforced at the interfaces.
Wave propagation properties of each waveguide are calculated using
the WFE method and the coupling element is modelled using the
FE method. The scattering of waves through the coupling element,
on which damage is modelled, is determined by coupling the FE and
WFE models. Furthermore, the central aim is to evaluate the effect of
pressurization on the wave dispersion and scattering characteristics
of the prestressed structural system compared to that which is not
prestressed. Numerical case studies are exhibited for two waveguides
coupled through a coupling joint.
Abstract: This paper is concerned with the development of a
fully implicit and purely Eulerian fluid-structure interaction method
tailored for the modeling of the large deformations of elastic
membranes in a surrounding Newtonian fluid. We consider a
simplified model for the mechanical properties of the membrane, in
which the surface strain energy depends on the membrane stretching.
The fully Eulerian description is based on the advection of a modified
surface tension tensor, and the deformations of the membrane are
tracked using a level set strategy. The resulting nonlinear problem
is solved by a Newton-Raphson method, featuring a quadratic
convergence behavior. A monolithic solver is implemented, and we
report several numerical experiments aimed at model validation and
illustrating the accuracy of the presented method. We show that
stability is maintained for significantly larger time steps.
Abstract: Out of all methods for ground improvement, stone column became more popular these days due to its simple construction and economic consideration. Installation of stone column especially in loose fine graded soil causes increasing in load bearing capacity and settlement reduction. Encased granular stone columns (EGCs) are commonly subjected to vertical load. However, they may also be subjected to significant amount of shear loading. In this study, three-dimensional finite element (FE) analyses were conducted to estimate the shear load capacity of EGCs in sandy soil. Two types of different cases, stone column and geosynthetic encased stone column were studied at different normal pressures varying from 15 kPa to 75 kPa. Also, the effect of diameter in two cases was considered. A close agreement between the experimental and numerical curves of shear stress - horizontal displacement trend line is observed. The obtained result showed that, by increasing the normal pressure and diameter of stone column, higher shear strength is mobilized by soil; however, in the case of encased stone column, increasing the diameter had more dominated effect in mobilized shear strength.
Abstract: This article is an attempt to present a numerically study of the behaviour of an eccentrically loaded circular footing resting on sand to determine its ultimate bearing capacity. A surface circular footing of diameter 12 cm (D) was used as shallow foundation. For this purpose, three dimensional models consist of foundation, and medium sandy soil was modelled by ABAQUS software. Bearing capacity of footing was evaluated and the effects of the load eccentricity on bearing capacity, its settlement, and modulus of subgrade reaction were studied. Three different values of load eccentricity with equal space from inside the core on the core boundary and outside the core boundary, which were respectively e=0.75, 1.5, and 2.25 cm, were considered. The results show that by increasing the load eccentricity, the ultimate load and the modulus of subgrade reaction decreased.
Abstract: In this work, numerical simulations were carried out using a specific CFD code in order to study the performance of an innovative Scraped Surface Heat Exchanger (SSHE) with helical ribbons for Bingham fluids (threshold fluids). The resolution of three-dimensional form of the conservation equations (continuity, momentum and energy equations) was carried out basing on the finite volume method (FVM). After studying the effect of dimensionless numbers (axial Reynolds, rotational Reynolds and Oldroyd numbers) on the hydrodynamic and thermal behaviors within SSHE, a parametric study was developed, by varying the width of the helical ribbon, the clearance between the stator wall and the tip of the ribbon and the number of turns of the helical ribbon, in order to improve the heat transfer inside the exchanger. The effect of these geometrical numbers on the hydrodynamic and thermal behaviors was discussed.
Abstract: Environment and surrounding plays a pivotal rule in structuring life-style of the consumers. Living standards intern effect the energy consumption of the consumers. In smart grid paradigm, climate drifts, weather parameter and green environmental directly relates to the energy profiles of the various consumers, such as residential, commercial and industrial. Considering above factors helps policy in shaping utility load curves and optimal management of demand and supply. Thus, there is a pressing need to develop correlation models of load and weather parameters and critical analysis of the factors effecting energy profiles of smart grid consumers. In this paper, we elaborated various environment and weather parameter factors effecting demand of consumers. Moreover, we developed correlation models, such as Pearson, Spearman, and Kendall, an inter-relation between dependent (load) parameter and independent (weather) parameters. Furthermore, we validated our discussion with real-time data of Texas State. The numerical simulations proved the effective relation of climatic drifts with energy consumption of smart grid consumers.
Abstract: So far, the conventional experimental and theoretical analysis in fracture mechanics have been applied to study concrete flexural- cracked beams, which are strengthened using fiber reinforced polymer (FRP) composite sheets. However, there is still little knowledge about the shear capacity of a side face FRP- strengthened shear-cracked beam. A numerical analysis is herein presented to model the fracture mechanics of a four-point RC beam, with two inclined initial notch on the supports, which is strengthened with side face FRP sheets. In the present study, the shear crack is forced to conduct by using an initial notch in supports. The ABAQUS software is used to model crack propagation by conventional cohesive elements. It is observed that the FRP sheets play important roles in preventing the propagation of shear cracks.
Abstract: In this numerical work, natural convection and entropy generation of Al2O3–water nanofluid in square cavity have been studied. A two-dimensional steady laminar natural convection in a differentially heated square cavity of length L, filled with a nanofluid is investigated numerically. The horizontal walls are considered adiabatic. Vertical walls corresponding to x=0 and x=L are respectively maintained at hot temperature, Th and cold temperature, Tc. The resolution is performed by the CFD code "FLUENT" in combination with GAMBIT as mesh generator. These simulations are performed by maintaining the Rayleigh numbers varied as 103 ≤ Ra ≤ 106, while the solid volume fraction varied from 1% to 5%, the particle size is fixed at dp=33 nm and a range of the temperature from 20 to 70 °C. We used models of thermophysical nanofluids properties based on experimental measurements for studying the effect of adding solid particle into water in natural convection heat transfer and entropy generation of nanofluid. Such as models of thermal conductivity and dynamic viscosity which are dependent on solid volume fraction, particle size and temperature. The average Nusselt number is calculated at the hot wall of the cavity in a different solid volume fraction. The most important results is that at low temperatures (less than 40 °C), the addition of nanosolids Al2O3 into water leads to a decrease in heat transfer and entropy generation instead of the expected increase, whereas at high temperature, heat transfer and entropy generation increase with the addition of nanosolids. This behavior is due to the contradictory effects of viscosity and thermal conductivity of the nanofluid. These effects are discussed in this work.
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: This paper copes with the numerical simulation for convective heat transfer in the stator disk of an axial flux permanent magnet (AFPM) electrical machine. Overheating is one of the main issues in the design of AFMPs, which mainly occurs in the stator disk, so that it needs to be prevented. A rotor-stator configuration with 16 magnets at the periphery of the rotor is considered. Air is allowed to flow through openings in the rotor disk and channels being formed between the magnets and in the gap region between the magnets and the stator surface. The rotating channels between the magnets act as a driving force for the air flow. The significant non-dimensional parameters are the rotational Reynolds number, the gap size ratio, the magnet thickness ratio, and the magnet angle ratio. The goal is to find correlations for the Nusselt number on the stator disk according to these non-dimensional numbers. Therefore, CFD simulations have been performed with the multiple reference frame (MRF) technique to model the rotary motion of the rotor and the flow around and inside the machine. A minimization method is introduced by a pattern-search algorithm to find the appropriate values of the reference temperature. It is found that the correlations are fast, robust and is capable of predicting the stator heat transfer with a good accuracy. The results reveal that the magnet angle ratio diminishes the stator heat transfer, whereas the rotational Reynolds number and the magnet thickness ratio improve the convective heat transfer. On the other hand, there a certain gap size ratio at which the stator heat transfer reaches a maximum.
Abstract: This study focused on the contribution of micro-mechanical parameters on the macro-mechanical response of short fiber composites, namely polypropylene matrix reinforced by glass fibers. In the framework of this paper, an attention has been given to the glass fibers length, as micromechanical parameter influences the overall macroscopic material’s behavior. Three dimensional numerical models were developed and analyzed through the concept of a Representative Volume Element (RVE). Results of the RVE-based approach were compared with analytical Halpin-Tsai’s model.
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: In this paper, there is an implementation, verification, and graphical demonstration of a software application, which can be used swiftly over different preliminary orbit determination methods. A passive orbit determination method is used in this study to determine the location of a satellite or a flying body. It is named a passive orbit determination because it depends on observation without the use of any aids (radio and laser) installed on satellite. In order to understand how these methods work and how their output is accurate when compared with available verification data, the built models help in knowing the different inputs used with each method. Output from the different orbit determination methods (Gibbs, Lambert, and Gauss) will be compared with each other and verified by the data obtained from Satellite Tool Kit (STK) application. A modified model including all of the orbit determination methods using the same input will be introduced to investigate different models output (orbital parameters) for the same input (azimuth, elevation, and time). Simulation software is implemented using MATLAB. A Graphical User Interface (GUI) application named OrDet is produced using the GUI of MATLAB. It includes all the available used inputs and it outputs the current Classical Orbital Elements (COE) of satellite under observation. Produced COE are then used to propagate for a complete revolution and plotted on a 3-D view. Modified model which uses an adapter to allow same input parameters, passes these parameters to the preliminary orbit determination methods under study. Result from all orbit determination methods yield exactly the same COE output, which shows the equality of concept in determination of satellite’s location, but with different numerical methods.
Abstract: In this paper, we present a model and an algorithm for
the calculation of the optimal control limit, average cost, sample size,
and the sampling interval for an optimal Bayesian chart to control
the proportion of defective items produced using a semi-Markov
decision process approach. Traditional p-chart has been widely
used for controlling the proportion of defectives in various kinds
of production processes for many years. It is well known that
traditional non-Bayesian charts are not optimal, but very few optimal
Bayesian control charts have been developed in the literature, mostly
considering finite horizon. The objective of this paper is to develop
a fast computational algorithm to obtain the optimal parameters of a
Bayesian p-chart. The decision problem is formulated in the partially
observable framework and the developed algorithm is illustrated by
a numerical example.
Abstract: With the rapid development of national modern industry, people begin to pay attention to environmental pollution and harm caused by industrial dust. Based on above, a numerical study on the dedusting technology of industrial environment was conducted. The dynamic models of multicomponent particles collision and coagulation, breakage and deposition are developed, and the interaction of water droplet and aerosol particle in 2-Dimension flow field was researched by Eulerian-Lagrangian method and Multi-Monte Carlo method. The effects of the droplet scale, movement speed of droplet and the flow field structure on scavenging efficiency were analyzed. The results show that under the certain condition, 30μm of droplet has the best scavenging efficiency. At the initial speed 1m/s of droplets, droplets and aerosol particles have more time to interact, so it has a better scavenging efficiency for the particle.