Pub Date : 2004-09-01DOI: 10.1142/S1465876304002563
Z. Xu, C. Picu, J. Fish
The classical approach to linking lattice dynamics properties to continuum equations of motion, the "method of long waves," is extended to include higher order terms. The additional terms account for non-local and non-linear effects. In the first part of the article, the derivation is made within the harmonic approximation for the perfect lattice response. Higher order terms are included in the continuum equation of motion to account for non-linear dispersion effects. Wave propagation coefficients as well as fourth order dispersion coefficients are obtained. In the second part, the lattice anharmonicity is considered and nonlinear macroscopic equations of motion are obtained within the local approximation. Both continuum solutions are particularized to the one-dimensional case and are compared with the lattice response in order to establish the accuracy of the approximation.
{"title":"Higher order continuum wave equation calibrated on lattice dynamics","authors":"Z. Xu, C. Picu, J. Fish","doi":"10.1142/S1465876304002563","DOIUrl":"https://doi.org/10.1142/S1465876304002563","url":null,"abstract":"The classical approach to linking lattice dynamics properties to continuum equations of motion, the \"method of long waves,\" is extended to include higher order terms. The additional terms account for non-local and non-linear effects. In the first part of the article, the derivation is made within the harmonic approximation for the perfect lattice response. Higher order terms are included in the continuum equation of motion to account for non-linear dispersion effects. Wave propagation coefficients as well as fourth order dispersion coefficients are obtained. In the second part, the lattice anharmonicity is considered and nonlinear macroscopic equations of motion are obtained within the local approximation. Both continuum solutions are particularized to the one-dimensional case and are compared with the lattice response in order to establish the accuracy of the approximation.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"496 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"132479281","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2004-09-01DOI: 10.1142/S146587630400254X
Ramana Kumar Kaza, S. Saikumar, M. Wang
Most of the work in the field of topology optimization is concentrated on using sensitivity analysis and optimality criteria methods that need explicit formulation. The design systems are often hard-coded for a specific problem with specialized optimization and FEM routines. This paper presents a work that uses a system approach to solid free form design. It attempts to develop a general topology optimization system that has a wide range of applicability by making use of sophisticated optimization and FEM packages available. A computer design system is implemented with an integration of commercial codes CFSQP and NASTRAN. A pre-processor and a post-processor are developed to assist the optimal design process. The system is tested with benchmark cases for minimum mean compliance and minimum weight designs. The results for the cases are presented, demonstrating the ability of the system to handle complex cases with practical feasibility. The implementation is evaluated with a parametric study of its performance. The key factors for the common problems of topology optimization are examined, including the mesh dependency and numerical instability. The computational efficiency is further studied to indicate the direction for further improvement of the system.
{"title":"A system approach to solid free form design of optimal structures","authors":"Ramana Kumar Kaza, S. Saikumar, M. Wang","doi":"10.1142/S146587630400254X","DOIUrl":"https://doi.org/10.1142/S146587630400254X","url":null,"abstract":"Most of the work in the field of topology optimization is concentrated on using sensitivity analysis and optimality criteria methods that need explicit formulation. The design systems are often hard-coded for a specific problem with specialized optimization and FEM routines. This paper presents a work that uses a system approach to solid free form design. It attempts to develop a general topology optimization system that has a wide range of applicability by making use of sophisticated optimization and FEM packages available. A computer design system is implemented with an integration of commercial codes CFSQP and NASTRAN. A pre-processor and a post-processor are developed to assist the optimal design process. The system is tested with benchmark cases for minimum mean compliance and minimum weight designs. The results for the cases are presented, demonstrating the ability of the system to handle complex cases with practical feasibility. The implementation is evaluated with a parametric study of its performance. The key factors for the common problems of topology optimization are examined, including the mesh dependency and numerical instability. The computational efficiency is further studied to indicate the direction for further improvement of the system.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"9 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"116282448","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2004-09-01DOI: 10.1142/S1465876304002502
Qiang Yu, S. Esche
This paper introduces an object-oriented programming approach into the development of a general-purpose simulation framework for the numerical modeling of microstructure evolution at the mesoscopic length level. The design requirements of this framework are discussed, and an analysis from the standpoint of object-oriented programming is performed. Based on this framework, a Monte Carlo grain growth algorithm with improved efficiency and accuracy has been developed. This simulation framework is expected to represent a reliable software module of a multi-scale microstructure prediction system for materials processing.
{"title":"Mesoscopic Computer Modeling Of Microstructure Evolution Within An Object-Oriented Simulation Framework","authors":"Qiang Yu, S. Esche","doi":"10.1142/S1465876304002502","DOIUrl":"https://doi.org/10.1142/S1465876304002502","url":null,"abstract":"This paper introduces an object-oriented programming approach into the development of a general-purpose simulation framework for the numerical modeling of microstructure evolution at the mesoscopic length level. The design requirements of this framework are discussed, and an analysis from the standpoint of object-oriented programming is performed. Based on this framework, a Monte Carlo grain growth algorithm with improved efficiency and accuracy has been developed. This simulation framework is expected to represent a reliable software module of a multi-scale microstructure prediction system for materials processing.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"18 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125261398","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2004-09-01DOI: 10.1142/S1465876304002526
K. Elmer
Realistic mechanical wave propagation problems and acoustic problems in transient dynamics with a broad frequency range require very large systems and many timesteps, to obtain reliable numerical solutions with sufficient high accuracy. New algorithms are to be developed to make effective use of todays hardware features and the visualization of complex prozesses like intensity and energy flow in transient dynamics. The propagation of mechanical waves is characterized by the change and interaction of strain energy and kinetic energy in space and time. A fast explicit FD-algorithm for wave propagation problems and acoustic problems is developed to analyse and visualise the complex behavior of traveling waves. This inherent parallel algorithm is based on the solution of the three-dimensional wave equation as a first order formulation in terms of stresses and velocities or acoustic pressure and particle velocity, representing both forms of energy in a direct way. Because of the small storage required and the short computational time, the algorithm allows numerical investigations of large systems with online graphic simulations to analyse the complex real physical behavior of propagating waves and to make numerical results comparable to measured results. Local mesh refinement helps to minimize numerical errors of the discrete model. Examples of applications are given with dispersional effects of traveling waves, instantaneous intensity distribution and local energy flow of propagating and standing waves. The complex behavior of traveling waves in a bar with a crack is analysed as a three-dimensional system. As a result, a non-destructive testing method is described using impact hammer for the detection, localization and quantification of cracks. The size of the defects can be of some order smaller than the used wave length.
{"title":"Three-dimensional wave propagation and energy flow","authors":"K. Elmer","doi":"10.1142/S1465876304002526","DOIUrl":"https://doi.org/10.1142/S1465876304002526","url":null,"abstract":"Realistic mechanical wave propagation problems and acoustic problems in transient dynamics with a broad frequency range require very large systems and many timesteps, to obtain reliable numerical solutions with sufficient high accuracy. New algorithms are to be developed to make effective use of todays hardware features and the visualization of complex prozesses like intensity and energy flow in transient dynamics. The propagation of mechanical waves is characterized by the change and interaction of strain energy and kinetic energy in space and time. A fast explicit FD-algorithm for wave propagation problems and acoustic problems is developed to analyse and visualise the complex behavior of traveling waves. This inherent parallel algorithm is based on the solution of the three-dimensional wave equation as a first order formulation in terms of stresses and velocities or acoustic pressure and particle velocity, representing both forms of energy in a direct way. Because of the small storage required and the short computational time, the algorithm allows numerical investigations of large systems with online graphic simulations to analyse the complex real physical behavior of propagating waves and to make numerical results comparable to measured results. Local mesh refinement helps to minimize numerical errors of the discrete model. Examples of applications are given with dispersional effects of traveling waves, instantaneous intensity distribution and local energy flow of propagating and standing waves. The complex behavior of traveling waves in a bar with a crack is analysed as a three-dimensional system. As a result, a non-destructive testing method is described using impact hammer for the detection, localization and quantification of cracks. The size of the defects can be of some order smaller than the used wave length.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"146 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"116117248","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2004-09-01DOI: 10.1142/S1465876304002575
H. Jin, E. Li, W. Yuan, Lewei Li
The paper presents and discusses the parallelized 3D finite-difference time-domain algorithm based on single-program multiple-data architecture using the MPI protocol for electromagnetic compatibility and signal integrity analysis and evaluation in high-speed designs. The efficiency for parallelization is verified through numerical experiments at IBM P640 and P690 multiprocessor high performance computers with up to 20 processors at different sizes of the problems.
{"title":"Parallel FDTD computing for emc simulation in high-speed electronics","authors":"H. Jin, E. Li, W. Yuan, Lewei Li","doi":"10.1142/S1465876304002575","DOIUrl":"https://doi.org/10.1142/S1465876304002575","url":null,"abstract":"The paper presents and discusses the parallelized 3D finite-difference time-domain algorithm based on single-program multiple-data architecture using the MPI protocol for electromagnetic compatibility and signal integrity analysis and evaluation in high-speed designs. The efficiency for parallelization is verified through numerical experiments at IBM P640 and P690 multiprocessor high performance computers with up to 20 processors at different sizes of the problems.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"46 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"124309796","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2004-09-01DOI: 10.1142/S1465876304002642
Y. Ogata, T. Yabe
We propose a multi-dimensional scheme to solve the shallow water equations by using the CIP(Constrained Interpolation Profile/Cubic Interpolated Pseudoparticle) method which has the accuracy by one-order of magnitude better than the cubic Lagrange. This third-order accuracy is kept even in non-uniform grid, while the cubic Lagrange becomes the first order. The semi-Lagrangian approach to the characteristic equations has been used for solving gravity waves because it can make time step be much longer than the Courant-Friedrichs-Lewy (CFL) condition. Even in two dimensions, the method can give a symmetrical wave propagation under the large CFL condition.
{"title":"Multi-dimensional semi-lagrangian characteristic approach to the shallow water equations by the CIP method","authors":"Y. Ogata, T. Yabe","doi":"10.1142/S1465876304002642","DOIUrl":"https://doi.org/10.1142/S1465876304002642","url":null,"abstract":"We propose a multi-dimensional scheme to solve the shallow water equations by using the CIP(Constrained Interpolation Profile/Cubic Interpolated Pseudoparticle) method which has the accuracy by one-order of magnitude better than the cubic Lagrange. This third-order accuracy is kept even in non-uniform grid, while the cubic Lagrange becomes the first order. The semi-Lagrangian approach to the characteristic equations has been used for solving gravity waves because it can make time step be much longer than the Courant-Friedrichs-Lewy (CFL) condition. Even in two dimensions, the method can give a symmetrical wave propagation under the large CFL condition.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"70 6 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"129639586","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2004-09-01DOI: 10.1142/S1465876304002599
G. Giese
In this paper we present an efficient numerical one-step method of high order in space and time for solving the elastic-plastic wave equation in three space dimensions. The basic idea is to decompose the hyperbolic PDE into advection equations, which can be solved numerically, Furthermore, the occurrence of plasticity makes it necessary to solve an ODE for the stress-strain relationship at every point.
{"title":"High-resolution simulation of the elastic-plastic wave equation in three space dimensions","authors":"G. Giese","doi":"10.1142/S1465876304002599","DOIUrl":"https://doi.org/10.1142/S1465876304002599","url":null,"abstract":"In this paper we present an efficient numerical one-step method of high order in space and time for solving the elastic-plastic wave equation in three space dimensions. The basic idea is to decompose the hyperbolic PDE into advection equations, which can be solved numerically, Furthermore, the occurrence of plasticity makes it necessary to solve an ODE for the stress-strain relationship at every point.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"10 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125335406","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2004-09-01DOI: 10.1142/S1465876304002629
A. A. M. Dayem, H. Shalaby
The salt-gradient solar ponds are considered in this work in order to study its performance numerically. Governing equations of the upper convective zone, nonconvective zone, and lower convective zone of a salt gradient solar pond are deduced as a set of non-linear partial differential equations. The equations are solved numerically to predict the thermal performance of the solar pond over a long time. The meteorological data of Egypt such as incident solar radiation, ambient temperature, air velocity, and relative humidity are taken into considerations in the model. Heat transfer modes considered between the upper convective zone and the ambient are convection, evaporation, and radiation. The present model is used to study the development of temperature inside the three zones of salt gradient solar pond. The optimum thickness of each layer is obtained with close agreement of previous results.
{"title":"Numerical simulation of salt gradient solar ponds","authors":"A. A. M. Dayem, H. Shalaby","doi":"10.1142/S1465876304002629","DOIUrl":"https://doi.org/10.1142/S1465876304002629","url":null,"abstract":"The salt-gradient solar ponds are considered in this work in order to study its performance numerically. Governing equations of the upper convective zone, nonconvective zone, and lower convective zone of a salt gradient solar pond are deduced as a set of non-linear partial differential equations. The equations are solved numerically to predict the thermal performance of the solar pond over a long time. The meteorological data of Egypt such as incident solar radiation, ambient temperature, air velocity, and relative humidity are taken into considerations in the model. Heat transfer modes considered between the upper convective zone and the ambient are convection, evaporation, and radiation. The present model is used to study the development of temperature inside the three zones of salt gradient solar pond. The optimum thickness of each layer is obtained with close agreement of previous results.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"45 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"124142910","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2004-09-01DOI: 10.1142/S1465876304002538
P. Jafarali, Lalitha Chattopadhyay, G. Prathap, S. Rajendran
We review critically the performance of an ingeniously designed hybrid beam element that uses a stiffness matrix based on Timoshenko theory but retains the mass matrix from classical beam theory. This clever engineering trick gives seemingly very accurate results in thin beam situations. However, the physics of thick beam behavior is consequently misrepresented. A careful study reveals that cancellation of errors is responsible for the apparent "accurate" performance.
{"title":"Error analysis of a hybrid beam element with Timoshenko stiffness and classical mass","authors":"P. Jafarali, Lalitha Chattopadhyay, G. Prathap, S. Rajendran","doi":"10.1142/S1465876304002538","DOIUrl":"https://doi.org/10.1142/S1465876304002538","url":null,"abstract":"We review critically the performance of an ingeniously designed hybrid beam element that uses a stiffness matrix based on Timoshenko theory but retains the mass matrix from classical beam theory. This clever engineering trick gives seemingly very accurate results in thin beam situations. However, the physics of thick beam behavior is consequently misrepresented. A careful study reveals that cancellation of errors is responsible for the apparent \"accurate\" performance.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"38 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"129038802","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2004-09-01DOI: 10.1142/S1465876304002551
Zheng H. Zhu, S. Meguid
Most existing curved beam elements suffer from poor convergence difficulties and a heavy computational burden while limit themselves to 2D problems. In this paper, we address and overcome these difficulties by developing a new three-noded locking-free 3D curved beam element. The element formulations, which employ coupled consistent polynomial displacement fields, satisfy the membrane locking-free requirement of being able to recover the inextensible bending mode of the curved beam. Quintic transverse displacement interpolation functions are used to represent the bending deformation of the beam, while the axial and torsional displacement fields are derived by integration of the presumably linear membrane and torsional shear strain fields, which are coupled with the transverse displacement fields. Numerical results of two- and three-dimensional applications are presented to demonstrate the superior accuracy and high convergence rate of the newly developed curved beam element compared with existing ones.
{"title":"Analysis of three-dimensional locking-free curved beam element","authors":"Zheng H. Zhu, S. Meguid","doi":"10.1142/S1465876304002551","DOIUrl":"https://doi.org/10.1142/S1465876304002551","url":null,"abstract":"Most existing curved beam elements suffer from poor convergence difficulties and a heavy computational burden while limit themselves to 2D problems. In this paper, we address and overcome these difficulties by developing a new three-noded locking-free 3D curved beam element. The element formulations, which employ coupled consistent polynomial displacement fields, satisfy the membrane locking-free requirement of being able to recover the inextensible bending mode of the curved beam. Quintic transverse displacement interpolation functions are used to represent the bending deformation of the beam, while the axial and torsional displacement fields are derived by integration of the presumably linear membrane and torsional shear strain fields, which are coupled with the transverse displacement fields. Numerical results of two- and three-dimensional applications are presented to demonstrate the superior accuracy and high convergence rate of the newly developed curved beam element compared with existing ones.","PeriodicalId":331001,"journal":{"name":"Int. J. Comput. Eng. Sci.","volume":"8 6 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2004-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125002129","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：