<p>Fractional calculus is the integration and differentiation to an arbitrary or fractional order. The techniques of fractional calculus are not commonly taught in engineering curricula since physical laws are expressed in integer order notation. Dr. Richard Magin (2006) notes how engineers occasionally encounter dynamic systems in which the integer order methods do not properly model the physical characteristics and lead to numerous mathematical operations. In the following study, the application of fractional order calculus to approximate the angular position of the disk oscillating in a Newtonian fluid was experimentally validated. The proposed experimental study was conducted to model the nonlinear response of an oscillating system using fractional order calculus. The integer and fractional order mathematical models solved the differential equation of motion specific to the experiment. The experimental results were compared to the integer order and the fractional order analytical solutions. The fractional order mathematical model in this study approximated the nonlinear response of the designed system by using the Bagley and Torvik fractional derivative. The analytical results of the experiment indicate that either the integer or fractional order methods can be used to approximate the angular position of the disk oscillating in the homogeneous solution. The following research was in collaboration with Dr. Richard Mark French, Dr. Garcia Bravo, and Rajarshi Choudhuri, and the experimental design was derived from the previous experiments conducted in 2018.</p>

Particle suspensions occur in many situations in nature and industry. In this master’s thesis, the motion of a single rigid spheroidal particle immersed in Stokes flow is studied numerically using a boundary integral method and a new specialized quadrature method known as quadrature by expansion (QBX). This method allows the spheroid to be massless or inertial, and placed in any kind of underlying Stokesian flow.   A parameter study of the QBX method is presented, together with validation cases for spheroids in linear shear flow and quadratic flow. The QBX method is able to compute the force and torque on the spheroid as well as the resulting rigid body motion with small errors in a short time, typically less than one second per time step on a regular desktop computer. Novel results are presented for the motion of an inertial spheroid in quadratic flow, where in contrast to linear shear flow the shear rate is not constant. It is found that particle inertia induces a translational drift towards regions in the fluid with higher shear rate. / Partikelsuspensioner förekommer i många sammanhang i naturen och industrin. I denna masteruppsats studeras rörelsen hos en enstaka stel sfäroidisk partikel i Stokesflöde numeriskt med hjälp av en randintegralmetod och en ny specialiserad kvadraturmetod som kallas quadrature by expansion (QBX). Metoden fungerar för masslösa eller tröga sfäroider, som kan placeras i ett godtyckligt underliggande Stokesflöde.   En parameterstudie av QBX-metoden presenteras, tillsammans med valideringsfall för sfäroider i linjärt skjuvflöde och kvadratiskt flöde. QBX-metoden kan beräkna kraften och momentet på sfäroiden samt den resulterande stelkroppsrörelsen med små fel på kort tid, typiskt mindre än en sekund per tidssteg på en vanlig persondator. Nya resultat presenteras för rörelsen hos en trög sfäroid i kvadratiskt flöde, där skjuvningen till skillnad från linjärt skjuvflöde inte är konstant. Det visar sig att partikeltröghet medför en drift i sidled mot områden i fluiden med högre skjuvning.

Indiana University-Purdue University Indianapolis (IUPUI) / Our ability to engineer materials is limited by our capacity to tailor the material’s microstructure morphology and predict resulting properties. The insufficient knowledge on microstructure-property relationship is due to complexity and randomness in all materials at different scales. The objective of this research is to establish a design optimization methodology for microstructured materials. The material design problem is stated as finding the optimum microstructure to maximize the desired performance satisfying material processing constrains. This problem has been solved in this thesis by means of numerical techniques through four main steps: microstructure characterization, model reconstruction, property evaluation, and optimization. Two methods of microstructure characterizations have been investigated along with the advantages and disadvantages of each method. The first microstructure characterization method is a statistical method which utilizes correlation functions to extract the microstructural information. Algorithms for calculating these correlations functions have been developed and optimized based on their computational cost using MATLAB software. The second microstructure characterization method is physical characterization which works based on evaluation of physical features in microstructured domain. These features have been measured by means of MATLAB codes. Three model reconstruction techniques are proposed based on these characterization methods and employed to generate material models for further evaluation. The first reconstructing algorithm uses statistical functions to reconstruct the statistical equivalent model through simulating annealing optimization method. The second algorithm uses cellular automaton concepts to simulate the grain growth utilizing physical descriptors, and the third one generates elliptical inclusions in a material matrix using physical characteristic of microstructure. The finite element method is used to analysis the mechanical behavior of material models. Several material samples with different microstructural characteristics have been generated to model the micro-scale design domain of AZ31 magnesium alloy and magnesium matrix composite with silicon carbide fibers. Then, surrogate models have been created based on these samples to approximate the entire design domain and demonstrate the sensitivity of the desired mechanical property to two independent microstructural features. Finally, the optimum microstructure characteristics of material samples for fracture strength maximization have been obtained.

Le contexte d’étude est celui de la Compatibilité ÉlectroMagnétique (CEM). L’objectif de la CEM est, comme son nom l’indique, d’assurer la compatibilité entre une source de perturbation électromagnétique et un système électronique victime. Or, la prédiction de ces niveaux de perturbation ne peut pas s’effectuer à l’aide d’un simple calcul analytique, en raison de la géométrie qui est généralement complexe pour le système que l’on étudie, tel que le champ à l’intérieur d’un cockpit d’avion par exemple. En conséquence, nous sommes contraints d’employer des méthodes numériques, dans le but de prédire ce niveau de couplage entre les sources et les victimes. Parmi les nombreuses méthodes numériques existantes à ce jour, les méthodes Multi-Domaines (MD) sont très prisées. En effet, elles offrent la liberté aux utilisateurs de choisir la méthode numérique la plus adaptée, en fonction de la zone géométrique à calculer. Au sein de ces méthodes MD, la « Domain Decomposition Method » (DDM) présente l’avantage supplémentaire de découpler chacun de ces domaines. En conséquence, la DDM est particulièrement intéressante, vis-à-vis des méthodes concurrentes, en particulier sur l’aspect du coût numérique. Pour preuve, l’ONERA continue de développer cette méthode qui ne cesse de montrer son efficacité depuis plusieurs années, notamment pour le domaine des Surfaces Équivalentes Radar (SER) et des antennes. L’objectif de l’étude est de tirer profit des avantages de cette méthode pour des problématiques de CEM. Jusqu’à maintenant, de nombreuses applications de CEM, traitées par le code DDM, fournissaient des résultats fortement bruités. Même pour des problématiques électromagnétiques très simples, des problèmes subsistaient, sans explication convaincante. Ceci justifie cette étude. Le but de cette thèse est de pouvoir appliquer ce formalisme DDM à des problématiques de CEM. Dans cette optique, nous avons été amenés à redéfinir un certain nombre de conventions, qui interviennent au sein de la DDM. Par ailleurs, nous avons développé un modèle spécifique pour les ouvertures, qui sont des voies de couplage privilégiées par les ondes, à l’intérieur des cavités que représentent les blindages. Comme les ouvertures sont, en pratique, de petites dimensions devant la longueur d’onde, on s’est intéressé à un modèle quasi-statique. Nous proposons alors un modèle, qui a été implémenté, puis validé. Suite à ce modèle, nous avons développé une méthode originale, basée sur un calcul en deux étapes, permettant de ne plus discrétiser le support des ouvertures dans les calculs 3D. / The context of the study is the ElectroMagnetic Compatibility (EMC). Principal aim of the EMC is to ensure the compatibility between an electromagnetic perturbance source and an electronic device victim. Unfortunately, the perturbation levels prediction can not be made using an analytic formula, because the geometry which is generally complex for the interesting system, for example the field inside an aircraft’s cockpit. Therefore, we are contrained to use numerical methods, to be able to evaluate this coupling level between sources and victims. Among several existing numerical methods, Multi-Domains (MD) methods are very interesting. They offer to users the freedom to choose the most powerfull numerical method, in terms of the geometrical zone evaluated. With the MD methods, « Domain Decomposition Method » (DDM) has the avantage of decouplingeach of theses areas. Therefore, DDM is very interesting, compared to other methods, in particular on the numerical cost. ONERA keeps on developing this method, which has not stop showing his efficiency since several years, in particular in Radar Cross Section (RCS) and antennas. The objective of this study is to take the benefits of this method for EMC problems. Up to now, several EMC applications treated by the DDM code provided results strongly noisy. Even for with very simple electromagnetic cases, some problems remained without convincing explanations. This justifies this study. The aim of this thesis is to can be able to apply DDM formalism to EMC problems. Then, we have been induced to redefine a number of conventions which are involved in the DDM. Otherwise, we have developed a specific model for the apertures which are privilegied tracts of the coupling by the penetration of waves inside cavities (shieldings). As the apertures have in practice smaller dimensions compared to the wavelength, we have been interested to a quasistatic model which was developped, implemented and validated. Following this model, we have developed an original method, based on a two step calculation, able to do not discretize the apertures support in 3D computations.

No description available.