Doctor of Philosophy / Department of Mechanical and Nuclear Engineering / Prakash Krishnaswami / Xiao J. Xin / Over the past two decades meshless methods have attracted much attention owing to their advantages in adaptivity, higher degree of solution field continuity, and capability to handle moving boundary and changing geometry. In this work, a meshless integral method based on the regularized boundary integral equation has been developed and applied to two-dimensional linear elasticity and elastoplasticity with small or large deformation. The development of the meshless integral method and its application to two-dimensional linear elasticity is described first. The governing integral equation is obtained from the weak form of elasticity over a local sub-domain, and the moving least-squares approximation is employed for meshless function approximation. This formulation incorporates: a subtraction method for singularity removal in the boundary integral equation, a special numerical integration for the calculation of integrals with weak singularity which further improves accuracy, a collocation method for the imposition of essential boundary conditions, and a method for incorporation of natural boundary conditions in the system governing equation. Next, elastoplastic material behavior with small deformation is introduced into the meshless integral method. The constitutive law is rate-independent flow theory based on von Mises yielding criterion with isotropic hardening. The method is then extended to large deformation plasticity based on Green-Naghdi’s theory using updated Lagrangian description. The Green-Lagrange strain is decomposed into the elastic and plastic part, and the elastoplastic constitutive law is employed that relates the Green-Lagrange strain to the second Piola-Kirchhoff stress. Finally, a pre- and post-processor for the meshless method using node- and pixel-based approach is presented. Numerical results from the meshless integral method agree well with available analytical solutions or finite element results, and the comparisons demonstrate that the meshless integral method is accurate and robust. This research lays the foundation for modeling and simulation of metal cutting processes.

Time-domain methods are preferred over their frequency-domain counterparts for solving acoustic and electromagnetic scattering problems since they can produce wide- band data from a single simulation. Among the time-domain methods, time-domain surface integral equation solvers have recently found widespread use because they offer several benefits over differential equation solvers. This dissertation develops several second-kind surface integral equation solvers for analyzing transient acoustic scattering from rigid and penetrable objects and transient electromagnetic scattering from perfect electrically conducting and dielectric objects. For acoustically rigid, perfect electrically conducting, and dielectric scatterers, fully explicit marching-on-in-time schemes are developed for solving time domain Kirchhoff, magnetic field, and scalar potential integral equations, respectively. The unknown quantity (e.g., velocity potential, electric current, or scalar potential) on the scatterer surface is discretized using a higher-order method in space and Lagrange interpolation in time. The resulting system is cast in the form of an ordinary differen- tial equation and integrated in time using a predictor-corrector scheme to obtain the unknown expansion coefficients. The explicit scheme can use the same time step size as its implicit counterpart without sacrificing from the stability of the solution and is much faster under low-frequency excitation (i.e., for large time step size). In addition, low-frequency behavior of vector potential integral equations for perfect electrically conducting scatterers is also investigated in this dissertation. For acoustically penetrable scatterers, presence of spurious interior resonance modes in the solutions of two forms of time domain surface integral equations is investigated. Numerical results demonstrate that the solution of the form that is widely used in the literature is corrupted by the interior resonance modes. But, the amplitude of these modes in the time domain can be suppressed by increasing the accuracy of discretization especially in time. On the other hand, the proposed one in the combined form shows a resonance-free performance verified via numerical experiments. In addition to providing detailed formulations of these solvers, the dissertation presents numerical examples, which demonstrate the solvers’ accuracy, efficiency, and applicability in real-life scenarios.

京都大学 / 0048 / 新制・課程博士 / 博士(エネルギー科学) / 甲第19092号 / エネ博第316号 / 新制||エネ||64(附属図書館) / 32043 / 京都大学大学院エネルギー科学研究科エネルギー基礎科学専攻 / (主査)教授 木下 正弘, 教授 森井 孝, 教授 片平 正人 / 学位規則第4条第1項該当 / Doctor of Energy Science / Kyoto University / DFAM

京都大学 / 0048 / 新制・課程博士 / 博士(情報学) / 甲第19121号 / 情博第567号 / 新制||情||100(附属図書館) / 32072 / 京都大学大学院情報学研究科複雑系科学専攻 / (主査)教授 木上 淳, 教授 磯 祐介, 教授 西村 直志 / 学位規則第4条第1項該当 / Doctor of Informatics / Kyoto University / DFAM

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第20632号 / 工博第4370号 / 新制||工||1679(附属図書館) / 京都大学大学院工学研究科分子工学専攻 / (主査)教授 佐藤 啓文, 教授 白川 昌宏, 教授 山本 量一 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DGAM

No description available.

Significant research is ongoing in the world to meet the needs of social and environmental crisis by harnessing wind and solar energy at high altitudes. One such approach is the use of an inflatable High Altitude Aerial Platform (HAAP). In the presented work, such periodically supported beam-like structures are analyzed using various mathematical models primarily modeling them as an equivalent beam using one-dimensional theories. The Euler-Bernoulli Theory (EBT) has been widely used for high aspect ratio beams, whereas the First Order Shear Deformation Theory (FSDT), or the Timoshenko beam theory, considers transverse shear effects and hence is superior in modeling low aspect ratio beams. First, an Isogeometric Analysis (IGA) is conducted using both FSDT and EBT to predict thermal buckling of periodically supported composite beams. Isogeometric analysis overcomes the limitations of the Gibbs phenomenon at discontinuities for a periodically supported beam using a higher order textit{k}-refinement. Next, an Integral Equation Approach (IEA) is implemented using EBT to obtain natural frequencies and buckling loads of periodically supported non-prismatic beams. Ill-conditioning errors were alleviated using admissible orthogonal Chebychev polynomials to obtain higher modes. We also present the prediction of the onset of flutter instability for metal plate and inflatable wing shaped foam test articles analyzed using finite element analysis (FEA). FEA updating based on modal testing and by conducting a geometrically nonlinear analysis resulted in a good agreement against the experiment tests. Furthermore, a nonlinear co-rotational large displacement/rotation FEA including the effects of the pressure as a follower forces was implemented to predict deformations of an inflatable structures. The developed FEA based tool namely Structural Analysis of Inflatables using FEA (SAIF) was compared with the experiments and available literature. It is concluded that the validity of the developed tool depends on the flexibility of the beam, which further depends upon the length of the beam and the bending rigidity of the beam. Inflatable structures analyzed with materials with high value of the Young's modulus and low to medium slenderness ratio tend to perform better against the experimental data. This is attributed to the presence of wrinkling and/or the Brazier effect (ovalling of the cross section) for flexible beams. The presented work has applications in programmable buckling, uncertainty quantification, and design of futuristic HAAP models to help face the upcoming environmental crises and meet the societal needs. / Doctor of Philosophy / In the future, developed countries are projected to face an increase in renewable energy demands due to environmental crises and increasing societal needs for energy due to urbanization. Wind energy, a renewable source, has received increasing attention. Wind farms require large land space and offshore wind energy harvesting is prone to unstable environments. Crosswind kite power is one of the promising and emerging fields where one can harvest energy from the wind farm inaccessible and apparently endless winds at high altitudes. In this dissertation, we present analysis and experiments on investigating complex structures, such as inflatable high altitude aerial platforms (HAAP) by using various mathematical models, primarily modeling them as an equivalent beam using one-dimensional theories. We investigate the effects of internal pressure on such structures, which unlike many other types of applied loads, follow the direction of the deflections. When supported on multiple supports, these structures are more efficient in terms of increased payload capacity due to a better distribution of loads, despite the increased weight penalty. To name a few, there are direct applications of periodic supports in design of bridges and railway sleepers. To avoid violent vibrations or failure, we also investigate the effect of multiple supports on the so-called natural frequency, vibration frequency under absence of applied loads, and buckling loads, instabilities under compression, of such beam-like structures. The presented work will aid in the design of futuristic HAAP models to help face the upcoming environmental crises and meet the energy demands of society due to urbanization.

No description available.

With the rapid increase in available compute power and memory, and bolstered by the advent of efficient formulations and algorithms, the role of 3D full-wave computational methods for accurate modelling of complex electromagnetic (EM) structures has gained in significance. The range of problems includes Radar Cross Section (RCS) computation, analysis and design of antennas and passive microwave circuits, bio-medical non-invasive detection and therapeutics, energy harvesting etc. Further, with the rapid advances in technology trends like System-in-Package (SiP) and System-on-Chip (SoC), the fidelity of chip-to-chip communication and package-board electrical performance parameters like signal integrity (SI), power integrity (PI), electromagnetic interference (EMI) are becoming increasingly critical. Rising pin-counts to satisfy functionality requirements and decreasing layer-counts to maintain cost-effectiveness necessitates 3D full wave electromagnetic solution for accurate system modelling. Method of Moments (MoM) is one such widely used computational technique to solve a 3D electromagnetic problem with full-wave accuracy. Due to lesser number of mesh elements or discretization on the geometry, MoM has an advantage of a smaller matrix size. However, due to Green's Function interactions, the MoM matrix is dense and its solution presents a time and memory challenge. The thesis focuses on formulation and development of novel techniques that aid in fast MoM based electromagnetic solutions. With the recent paradigm shift in computer hardware architectures transitioning from single-core microprocessors to multi-core systems, it is of prime importance to parallelize the serial electromagnetic formulations in order to leverage maximum computational benefits. Therefore, the thesis explores the possibilities to expedite an electromagnetic simulation by scalable parallelization of near-linear complexity algorithms like Fast Multipole Method (FMM) on a multi-core platform. Secondly, with the best of parallelization strategies in place and near-linear complexity algorithms in use, the solution time of a complex EM problem can still be exceedingly large due to over-meshing of the geometry to achieve a desired level of accuracy. Hence, the thesis focuses on judicious placement of mesh elements on the geometry to capture the physics of the problem without compromising on accuracy- a technique called Adaptive Mesh Refinement. This facilitates a reduction in the number of solution variables or degrees of freedom in the system and hence the solution time. For multi-scale structures as encountered in chip-package-board systems, the MoM formulation breaks down for parts of the geometry having dimensions much smaller as compared to the operating wavelength. This phenomenon is popularly known as low-frequency breakdown or low-frequency instability. It results in an ill-conditioned MoM system matrix, and hence higher iteration count to converge when solved using an iterative solver framework. This consequently increases the solution time of simulation. The thesis thus proposes novel formulations to improve the spectral properties of the system matrix for real-world complex conductor and dielectric structures and hence form well-conditioned systems. This reduces the iteration count considerably for convergence and thus results in faster solution. Finally, minor changes in the geometrical design layouts can adversely affect the time-to-market of a commodity or a product. This is because the intermediate design variants, in spite of having similarities between them are treated as separate entities and therefore have to follow the conventional model-mesh-solve workflow for their analysis. This is a missed opportunity especially for design variant problems involving near-identical characteristics when the information from the previous design variant could have been used to expedite the simulation of the present design iteration. A similar problem occurs in the broadband simulation of an electromagnetic structure. The solution at a particular frequency can be expedited manifold if the matrix information from a frequency in its neighbourhood is used, provided the electrical characteristics remain nearly similar. The thesis introduces methods to re-use the subspace or Eigen-space information of a matrix from a previous design or frequency to solve the next incremental problem faster.

Understanding, predicting, and controlling electromagnetic field interactions on and between complex RF platforms requires high fidelity computational electromagnetic (CEM) simulation. The primary CEM tool within NASA is GEMINI, an integral equation based method-of-moments (MoM) code for frequency domain electromagnetic modeling. However, GEMINI is currently limited in the size and complexity of problems that can be effectively handled. To extend GEMINI’S CEM capabilities beyond those currently available, primary research is devoted to integrating the MFDlib library developed at the University of Kentucky with GEMINI for efficient filling, factorization, and solution of large electromagnetic problems formulated using integral equation methods. A secondary research project involves the hybrid parallelization of GEMINI for the efficient speedup of the impedance matrix filling process. This thesis discusses the research, development, and testing of the secondary research project on the High Performance Computing DLX Linux supercomputer cluster. Initial testing of GEMINI’s existing MPI parallelization establishes the benchmark for speedup and reveals performance issues subsequently solved by the NASA CEM Lab. Implementation of hybrid parallelization incorporates GEMINI’s existing course level MPI parallelization with Open MP fine level parallel threading. Simple and nested Open MP threading are compared. Final testing documents the improvements realized by hybrid parallelization.