Computation of Mathematical Models for Complex Industrial Processes

Designed for undergraduate and postgraduate students, academic researchers and industrial practitioners, this book provides comprehensive case studies on numerical computing of industrial processes and step-by-step procedures for conducting industrial computing. It assumes minimal knowledge in numerical computing and computer programming, making it easy to read, understand and follow. Topics discussed include fundamentals of industrial computing, finite difference methods, the Wavelet-Collocation Method, the Wavelet-Galerkin Method, High Resolution Methods, and comparative studies of various methods. These are discussed using examples of carefully selected models from real processes of industrial significance. The step-by-step procedures in all these case studies can be easily applied to other industrial processes without a need for major changes and thus provide readers with useful frameworks for the applications of engineering computing in fundamental research problems and practical development scenarios.

Vibration of simply-supported trapezoidal plates I. Symmetric trapezoids

A detailed investigation of the natural frequencies and mode shapes of simply supported symmetric trapezoidal plates is undertaken in this paper. For numerical calculations, the relationship that exists between the eigenvalue problem of a polygonal simply supported plate and the eigenvalue problem of polygonal membrane of the same shape is utilized with advantage. The deflection surface is expressed in terms of a Fourier sine series in transformed coordinates and the Galerkin method is used. Results are presented in the form of tables and graphs. Several features like the crossing of frequency curves and the metamorphosis of some of the nodal patterns are observed. By a suitable interpretation of the modes of those symmetric trapezoidal plates which have the median as the nodal line, the results for some of the modes of unsymmetrical trapezoidal plates are also deduced.

Vibration of simply-supported trapezoidal plates II. Unsymmetric trapezoids

This paper deals with the investigation of the vibration characteristics of simply-supported unsymmetric trapezoidal plates. For numerical calculations, the relationship between the eigenvalue problems of a polygonal simply-supported plate and polygonal membrane is again effectively utilized. The Galerkin method is applied, with the deflection surface expressed in terms of a Fourier sine series in transformed coordinates. Numerical values for the first seven to eight frequencies for different geometries of the unsymmetric trapezoid are presented in the form of tables. Also the nodal patterns for a few representative configurations are presented.

A comparative study of dragonfly inspired flapping wings actuated by single crystal piezoceramic

A dragonfly inspired flapping wing is investigated in this paper. The flapping wing is actuated from the root by a PZT-5H and PZN-7%PT single crystal unimorph in the piezofan configuration. The nonlinear governing equations of motion of the smart flapping wing are obtained using the Hamilton's principle. These equations are then discretized using the Galerkin method and solved using the method of multiple scales. Dynamic characteristics of smart flapping wings having the same size as the actual wings of three different dragonfly species Aeshna Multicolor, Anax Parthenope Julius and Sympetrum Frequens are analyzed using numerical simulations. An unsteady aerodynamic model is used to obtain the aerodynamic forces. Finally, a comparative study of performances of three piezoelectrically actuated flapping wings is performed. The numerical results in this paper show that use of PZN-7%PT single crystal piezoceramic can lead to considerable amount of wing weight reduction and increase of lift and thrust force compared to PZT-5H material. It is also shown that dragonfly inspired smart flapping wings actuated by single crystal piezoceramic are a viable contender for insect scale flapping wing micro air vehicles.

Frequency Domain Based Solution for Certain Class of Wave Equations: An exhaustive study of Numerical Solutions

The paper discusses the frequency domain based solution for a certain class of wave equations such as: a second order partial differential equation in one variable with constant and varying coefficients (Cantilever beam) and a coupled second order partial differential equation in two variables with constant and varying coefficients (Timoshenko beam). The exact solution of the Cantilever beam with uniform and varying cross-section and the Timoshenko beam with uniform cross-section is available. However, the exact solution for Timoshenko beam with varying cross-section is not available. Laplace spectral methods are used to solve these problems exactly in frequency domain. The numerical solution in frequency domain is done by discretisation in space by approximating the unknown function using spectral functions like Chebyshev polynomials, Legendre polynomials and also Normal polynomials. Different numerical methods such as Galerkin Method, Petrov- Galerkin method, Method of moments and Collocation method or the Pseudo-spectral method in frequency domain are studied and compared with the available exact solution. An approximate solution is also obtained for the Timoshenko beam with varying cross-section using Laplace Spectral Element Method (LSEM). The group speeds are computed exactly for the Cantilever beam and Timoshenko beam with uniform cross-section and is compared with the group speeds obtained numerically. The shear mode and the bending modes of the Timoshenko beam with uniform cross-section are separated numerically by applying a modulated pulse as the shear force and the corresponding group speeds for varying taper parameter in are obtained numerically by varying the frequency of the input pulse. An approximate expression for calculating group speeds corresponding to the shear mode and the bending mode, and also the cut-off frequency is obtained. Finally, we show that the cut-off frequency disappears for large in, for epsilon > 0 and increases for large in, for epsilon < 0.

Faster computation of the Karhunen-Loeve expansion using its domain independence property

The goal of this work is to reduce the cost of computing the coefficients in the Karhunen-Loeve (KL) expansion. The KL expansion serves as a useful and efficient tool for discretizing second-order stochastic processes with known covariance function. Its applications in engineering mechanics include discretizing random field models for elastic moduli, fluid properties, and structural response. The main computational cost of finding the coefficients of this expansion arises from numerically solving an integral eigenvalue problem with the covariance function as the integration kernel. Mathematically this is a homogeneous Fredholm equation of second type. One widely used method for solving this integral eigenvalue problem is to use finite element (FE) bases for discretizing the eigenfunctions, followed by a Galerkin projection. This method is computationally expensive. In the current work it is first shown that the shape of the physical domain in a random field does not affect the realizations of the field estimated using KL expansion, although the individual KL terms are affected. Based on this domain independence property, a numerical integration based scheme accompanied by a modification of the domain, is proposed. In addition to presenting mathematical arguments to establish the domain independence, numerical studies are also conducted to demonstrate and test the proposed method. Numerically it is demonstrated that compared to the Galerkin method the computational speed gain in the proposed method is of three to four orders of magnitude for a two dimensional example, and of one to two orders of magnitude for a three dimensional example, while retaining the same level of accuracy. It is also shown that for separable covariance kernels a further cost reduction of three to four orders of magnitude can be achieved. Both normal and lognormal fields are considered in the numerical studies. (c) 2014 Elsevier B.V. All rights reserved.

A stabilization approach for mesh-free simulations of systems developing shocks or extreme strain localizations

A new stabilization scheme, based on a stochastic representation of the discretized field variables, is proposed with a view to reduce or even eliminate unphysical oscillations in the mesh-free numerical simulations of systems developing shocks or exhibiting localized bands of extreme deformation in the response. The origin of the stabilization scheme may be traced to nonlinear stochastic filtering and, consistent with a class of such filters, gain-based additive correction terms are applied to the simulated solution of the system, herein achieved through the element-free Galerkin method, in order to impose a set of constraints that help arresting the spurious oscillations. The method is numerically illustrated through its Applications to inviscid Burgers' equations, wherein shocks may develop as a result of intersections of the characteristics, and to a gradient plasticity model whose response is often characterized by a developing shear band as the external load is gradually increased. The potential of the method in stabilized yet accurate numerical simulations of such systems involving extreme gradient variations in the response is thus brought forth. (C) 2014 Elsevier Ltd. All rights reserved.

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Pull-in Instability Study of Carbon Nanotube Tweezers Under the Influence of Van der Waals Forces

The pull-in instability of two nanotubes under van der Waals force is studied. The cantilever beam with large deformation model is used. The influence of nanotube parameters such as the interior radius, the gap distance between the two nanotubes, etc, on the pull-in instability is studied. The critical nanotube length is determined for each specific set of nanotube parameters. The Galerkin method is applied to discretize the governing equations, and it shows good convergence.

Numerical and Analytical Study on the Pull-in Instability of Micro-Structure under Electrostatic Loading

The one-mode analysis method on the pull-in instability of micro-structure under electrostatic loading is presented. Taylor series are used to expand the electrostatic loading term in the one-mode analysis method, which makes analytical solution available. The one-mode analysis is the combination of Galerkin method and Cardan solution of cubic equation. The one-mode analysis offers a direct computation method on the pull-in voltage and displacement. In low axial loading range, it shows little difference with the established multi-mode analysis on predicting the pull-in voltages for three different structures (cantilever, clamped-clamped beams and the plate with four edges simply-supported) studied here. For numerical multi-mode analysis, we also show that using the structural symmetry to select the symmetric mode can greatly reduce both the computation effort and the numerical fluctuation.

Vibration of an Adhered Microbeam Under a Periodically Shaking Electrical Force

The vibration analysis of an adhered S-shaped microbeam under alternating sinusoidal voltage is presented. The shaking force is the electrical force due to the sinusoidal voltage. During vibration, both the microbeam deflection and the adhesion length keep changing. The microbeam deflection and adhesion length are numerically determined by the iteration method. As the adhesion length keeps changing, the domain of the equation of motion for the microbeam (unadhered part) changes correspondingly, which results in changes of the structure natural frequencies. For this reason, the system can never reach a steady state. The transient behaviors of the microbeam under different shaking frequencies are compared. We deliberately choose the initial conditions to compare our dynamic results with the existing static theory. The paper also analyzes the changing behavior of adhesion length during vibration and an asymmetric pattern of adhesion length change is revealed, which may be used to guide the dynamic de-adhering process. The abnormal behavior of the adhered microbeam vibrating at almost the same frequency under two quite different shaking frequencies is also shown. The Galerkin method is used to discretize the equation of motion and its convergence study is also presented. The model is only applicable in the case that the peel number is equal to 1. Some other model limitations are also discussed.

On the topological bifurcation of flows around a rotating circular cylinder

Flow fields around a rotating circular cylinder in a uniform stream are computed using a low dimensional Galerkin method. Results show that the formation of a Fopple vortex pair behind a stationary circular cylinder is caused by the structural instability in the vicinity of the saddle located at the rear of the cylinder. For rotating cylinder a bifurcation diagram with the consideration of two parameters, Reynolds number Re and rotation parameter a, is built by a kinematic analysis of the steady flow fields.

Hopf bifurcation in wakes behind a rotating and translating circular cylinder

A low-dimensional Galerkin method, initiated by Noack and Eckelmann [Physica D 56, 151 (1992)], for the prediction of the flow field around a stationary two-dimensional circular cylinder in a uniform stream at low Reynolds number is generalized to the case of a rotating and translating cylinder. The Hopf bifurcation describing the transition from steady to time-periodic solution is investigated. A curve indicating the transitional boundary is given in the two-dimensional parameter plane of Reynolds number Re and rotating parameter alpha. Our results show that rotation may delay the onset of vortex street and decrease the vortex-shedding frequency. (C) 1996 American Institute of Physics.

Métodos espectrais aplicados à relatividade numérica: determinação dos dados iniciais

Neste trabalho aplicamos métodos espectrais para a determinação da configuração inicial de três espaços-tempos contendo buracos negros. Para isto apresentamos primeiro a foliação do espaço-tempo em hipersuperfícies tridimensionais espaciais parametrizadas pela função temporal t. Este processo é chamado de decomposição 3+1 [2] [5]. O resultado deste processo são dois conjuntos de equações classificadas em equações de vínculo e evolução [4]. As equações de vínculo podem ser divididas em vínculos Hamiltoniano e dos momentos. Para a obtenção dos dados iniciais dos problemas estudados aqui, apenas a equação de vínculo Hamiltoniano será resolvida numericamente, pois as equações de vínculo dos momentos possuem solução analítica nestes casos. Uma pequena descrição dos métodos espectrais é apresentada, destacando-se os método de Galerkin, método pseudoespectral ou de colocação e método de Tau, que são empregados na resolução das equações de vínculo Hamiltoniano dos problemas estudados. Verificamos que os resultados obtidos neste trabalho superam aqueles produzidos por Kidder e Finn [15], devido a uma escolha diferente das funções de base, que aqui satisfazem uma das condições de contorno.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.