京都大学 / 0048 / 新制・論文博士 / 博士(農学) / 乙第13021号 / 論農博第2831号 / 新制||農||1042(附属図書館) / 学位論文||H28||N4967(農学部図書室) / 32949 / 京都大学大学院農学研究科地域環境科学専攻 / (主査)教授 藤原 正幸, 教授 村上 章, 准教授 宇波 耕一 / 学位規則第4条第2項該当 / Doctor of Agricultural Science / Kyoto University / DFAM

No description available.

With the algorithm’s suitability for exploiting current petascale and next-generation exascale supercomputers, stable and structure-preserving properties are necessary to develop predictive computational tools. In this dissertation, summation-by-parts (SBP) operators and a new relaxation Runge–Kutta (RRK) scheme are used to construct mimetic and structure-preserving full discretization for non-reactive compressible computational fluid dynamics (CFD) and reaction-diffusion models. In the first chapter, we provide the necessary background and a literature survey that forms the basis of this dissertation. Next, we provide a short overview of entropy stability for general conservation laws. The second chapter covers the analysis of the Eulerian model for compressible and heat-conducting flows. We provide the necessary background of the new system of parabolic partial differential equation (PDE). Then, we present the entropy stability analysis of the model at the continuous level. Subsequently, using the SBP, we construct an entropy-stable discretization of any order for unstructured grids with tensor-product elements. The third chapter discusses the implementation of RRK methods. We start by reviewing the RRK scheme constructed to guarantee conservation or stability with respect to any inner-product norm. Then, we present the extension and generalization of RRK schemes to general convex functionals and their application to compressible fluid flow problems. The final chapter demonstrates the far-reaching capabilities of the SBP operators and RRK schemes presenting the development of a novel fully discrete Lyapunov stable discretization for reaction models with spatial diffusion. Finally, we conclude this dissertation with an overview of our achievements and future research directions.

No description available.

No description available.

No description available.

There is an increasing need for accurate models describing the electrical behaviour of individual biological cells exposed to electromagnetic fields. In this area of solving linear problem, the most frequently used technique for computing the EM field is the Finite-Difference Time-Domain (FDTD) method. When modelling objects that are small compared with the wavelength, for example biological cells at radio frequencies, the standard Finite-Difference Time-Domain (FDTD) method requires extremely small time-step sizes, which may lead to excessive computation times. The problem can be overcome by implementing a quasi-static approximate version of FDTD, based on transferring the working frequency to a higher frequency and scaling back to the frequency of interest after the field has been computed. An approach to modeling and analysis of biological cells, incorporating the Hodgkin and Huxley membrane model, is presented here. Since the external medium of the biological cell is lossy material, a modified Berenger absorbing boundary condition is used to truncate the computation grid. Linear assemblages of cells are investigated and then Floquet periodic boundary conditions are imposed to imitate the effect of periodic replication of the assemblages. Thus, the analysis of a large structure of cells is made more computationally efficient than the modeling of the entire structure. The total fields of the simulated structures are shown to give reasonable and stable results at 900MHz, 1800MHz and 2450MHz. This method will facilitate deeper investigation of the phenomena in the interaction between EM fields and biological systems. Moreover, the nonlinear response of biological cell exposed to a 0.9GHz signal was discussed on observing the second harmonic at 1.8GHz. In this, an electrical circuit model has been proposed to calibrate the performance of nonlinear RF energy conversion inside a high quality factor resonant cavity with known nonlinear device. Meanwhile, the first and second harmonic responses of the cavity due to the loading of the cavity with the lossy material will also be demonstrated. The results from proposed mathematical model, give good indication of the input power required to detect the weakly effects of the second harmonic signal prior to perform the measurement. Hence, this proposed mathematical model will assist to determine how sensitivity of the second harmonic signal can be detected by placing the required specific input power.

Aeolian dune fields form some of the most striking patterns on Earth and Mars. These patterns reflect the internal dune dynamics of self-organization within boundary conditions, which are the unique set of environmental variables within which each dune field evolves. Dune-field pattern self-organization occurs because of interactions between the dunes themselves and the rich diversity of dune-field patterns arises because boundary conditions alter the type and frequency of dune interactions. These hypotheses are explored in three parts. First, source-area geometry and areal limits are two newly recognized boundary conditions. Measurements of crest length and spacing from satellite images of dune patterns with point and line source-area geometries show an increase in crest length and spacing over distance, whereas crest length and spacing in plane-sourced patterns emerge equally across the dune field. The areal limit boundary condition is the size and shape of the dune field itself. Empirical measurements from ten dune fields ranging over four orders of magnitude in area show that spacing increases and defect density decreases as the area of the dune field increases. A simple analytical model indicates that dune fields that are five times longer in the dune migration direction can achieve the greatest spacing for a given area. Second, time-series aerial photographs and airborne LiDAR show that fully developed, crescentic aeolian dunes at White Sands, New Mexico, interact and the dune pattern organizes in systematically similar ways as wind ripples and subaqueous dunes and ripples. Interaction type, classified as constructive, regenerative or neutral in terms of pattern development, changes spatially with the pattern because of the imposition of the line-source area and sediment availability boundary conditions. Upwind dominance by constructive interactions at the field line-source yields to neutral and regenerative interactions in the sediment availability-limited field center. Third, the dune-field pattern in the Olympia Undae Dune Field on Mars is comprised of two generations of dunes. This scenario of pattern reformation with a new wind regime shows that the emergence of the younger pattern is controlled by the boundary condition of the antecedent dune topography imposed upon the interaction between the younger and older patterns. / text

Piezoelectric materials in the forms of both bulk and thin-film have been widely used as actuators and sensors due to their electromechanical coupling. The characterization of piezoelectric materials plays an important role in determining device performance and reliability. Instrumented indentation is a promising method for probing mechanical as well as electrical properties of piezoelectric materials. The use of instrumented indentation to characterize the properties of piezoelectric materials requires analytical relations. Finite element methods are used to analyze the indentation of piezoelectric materials under different mechanical and electrical boundary conditions. For indentation of a piezoelectric half space, a three-dimensional finite element model is used due to the anisotropy and geometric nonlinearity. The analysis is focused on the effect of angle between poling direction and indentation-loading direction on indentation responses. For the indentation by a flat-ended cylindrical indenter, both insulating indenter and conducting indenter without a prescribed electric potential are considered. The results reveal that both the indentation load and the magnitude of the indentation-induced potential at the contact center increase linearly with the indentation depth. For the indentation by an insulating Berkovich indenter, both frictionless and frictional contact between the indenter and indented surface are considered. The results show the indentation load is proportional to the square of the indentation depth, while the indentation-induced potential at the contact center is proportional to the indentation depth. Spherical indentation of piezoelectric thin films is analyzed in an axisymmetric finite element model, in which the poling direction is anti-parallel to the indentation-loading direction. Six different combinations of electrical boundary conditions are considered for a thin film perfectly bonded to a rigid substrate under the condition of the contact radius being much larger than the film thickness. The indentation load is found to be proportional to the square of the indentation depth. To analyze the decohesion problem between a piezoelectric film and an elastic substrate, a traction-separation law is used to control the interfacial behavior between a thin film and an electrically grounded elastic substrate. The discontinuous responses at the initiation of interfacial decohesion are found to depend on interface and substrate properties.

Cette thèse traite de la diffusion par des surfaces rugueuses monodimensionnelles. Les surfaces présentant des petites échelles de variations nécessitent une discrétisation fine pour représenter les effets de diffusion sur le champ diffracté, ce qui augmente les coûts numériques. Deux aspects sont considérés : la réduction de la taille du problème en construisant une condition aux limiteséquivalente traduisant les effets des variations rapides et la réduction du nombre d’itérations nécessaires pour résoudre le système linéaire issu de la méthode des moments par une méthode basée sur les sous-espaces de Krylov. En ce qui concerne la réduction de la taille du problème, une technique d’homogénéisation est utilisée pour transformer la condition aux limites posée sur lasurface rugueuse par des paramètres effectifs. Ces paramètres sont déterminés par des problèmes auxiliaires qui tiennent compte des échelles fines de la surface. Dans le cas de surfaces parfaitement métalliques, la procédure est appliquée en polarisation Transverse Magnétique (TM) et Transverse Électrique (TE). Une impédance équivalente de Léontovich d’ordre 1 est déduite.Le procédure est automatique et les ordres supérieurs sont dérivés pour la polarisation TM. La procédure d’homogénéisation est aussi appliquée pour des interfaces rugueuses séparant deux milieux diélectriques. En ce qui concerne la réduction du nombre d’itérations, un préconditionneur, basé sur des considérations physiques, est construit à partir des modes de Floquet. Bien que le préconditionneur soit initialement élaboré pour des surfaces périodiques, nous montrons qu’il est aussi efficace pour des surfaces tronquées éclairées par une onde plane. L’efficacité des deux aspects présentés dans cette thèse est numériquement illustrée pour des configurations d’intérêt. / This work is about the scattering by monodimensional rough surfaces. Surfaces presenting small scales of variations need a very refined mesh to finally capture the scattering field behaviour what increases the computational cost. Two aspects are considered : the reduction of the problemsize through an effective boundary condition incorporating the effect of rapid variations and the reduction of the number of iterations to solve the linear system arising from method of moments by a method based on Krylov subspace. Firstly, an homogenization process is used to convert the boundary condition on the rough interface into effective parameters. These parameters are determined by the solutions of auxiliary problems which involve the detailed profile of the interface. In the case of perfectly metallic surfaces, the process is applied to the E- and H-polarization and an Leontovich impedance of order 1 is deduced. The process is automatic and higher orders are derived for E-polarization. The homogenization process is also applied to dielectric rough interfaces. Secondly, a physically-based preconditioner is built with Floquet’s modes. Although the preconditioner has been designed for periodical surfaces, it was shown to be efficient in the case of truncated surfaces illuminated by a plane wave. The efficiency of both aspects is numerically illustrated for some configurations of interest.