As a result of the increase in accessibility of computational resources and the increase in the power of the computers during the last two decades, designers are able to create computer models to simulate the behavior of a complex products. To address global competitiveness, companies are forced to optimize their designs and products. Optimizing the design needs several runs of computationally expensive simulation models. Therefore, using metamodels as an efficient and sufficiently accurate approximate of the simulation model is necessary. Radial basis functions (RBF) is one of the several metamodeling methods that can be found in the literature. The established approach is to add a bias to RBF in order to obtain a robust performance. The a posteriori bias is considered to be unknown at the beginning and it is defined by imposing extra orthogonality constraints. In this thesis, a new approach in constructing RBF with the bias to be set a priori by using the normal equation is proposed. The performance of the suggested approach is compared to the classic RBF with a posteriori bias. Another comprehensive comparison study by including several modeling criteria, such as problem dimension, sampling technique and size of samples is conducted. The studies demonstrate that the suggested approach with a priori bias is in general as good as the performance of RBF with a posteriori bias. Using the a priori RBF, it is clear that the global response is modeled with the bias and that the details are captured with radial basis functions. Multi-objective optimization and the approaches used in solving such problems are briefly described in this thesis. One of the methods that proved to be efficient in solving multi-objective optimization problems (MOOP) is the strength Pareto evolutionary algorithm (SPEA2). Multi-objective optimization of a disc brake system of a heavy truck by using SPEA2 and RBF with a priori bias is performed. As a result, the possibility to reduce the weight of the system without extensive compromise in other objectives is found. Multi-objective optimization of material model parameters of an adhesive layer with the aim of improving the results of a previous study is implemented. The result of the original study is improved and a clear insight into the nature of the problem is revealed.

Philosophiae Doctor - PhD / Many problems in financial world are being modeled by means of differential equation. These problems are time dependent, highly nonlinear, stochastic and heavily depend on the previous history of time. A variety of financial products exists in the market, such as forwards, futures, swaps and options. Our main focus in this thesis is to use the numerical analysis tools to solve some option pricing problems. Depending upon the inter-relationship of the financial derivatives, the dimension of the associated problem increases drastically and hence conventional methods (for example, the finite difference methods or finite element methods) for solving them do not provide satisfactory results. To resolve this issue, we use a special class of numerical methods, namely, the mesh free methods. These methods are often better suited to cope with changes in the geometry of the domain of interest than classical discretization techniques. In this thesis, we apply these methods to solve problems that price standard and non-standard options. We then extend the proposed approach to solve Heston' volatility model. The methods in each of these cases are analyzed for stability and thorough comparative numerical results are provided. / South Africa

Au cours de cette thèse, un modèle potentiel résolvant les équations d’Euler-Zakharov a été développé dans le but de simuler la propagation de vagues et d’états de mer irréguliers et multi-directionnels, du large jusqu’à la côte, sur des bathymétries variables. L’objectif est de représenter les effets non-linéaires et dispersifs le plus précisément possible pour des domainescôtiers bidimensionnels (dans le plan horizontal) de l’ordre de quelques kilomètres.La version 1DH initiale du modèle, résolvant le problème aux limites de Laplace à l’aide de schémas aux différences finies d’ordre élevé dans la direction horizontale combinés à une approche spectrale sur la verticale, a été améliorée et validée. L’implémentation de conditions aux limites de type Dirichlet et Neumann pour générer des vagues dans le domaine a été étudiée en détail. Dans la pratique, une zone de relaxation a été utilisée en complément deces conditions pour améliorer la stabilité du modèle.L’expression analytique de la relation de dispersion a été établie dans le cas d’un fond plat. Son analyse a montré que la représentation des effets dispersifs s’améliorait significativement avec l’augmentation de la résolution sur la direction verticale (i.e. avec le degré maximal de la basede polynômes de Tchebyshev utilisée pour projeter le potentiel des vitesses sur la verticale).Une étude de convergence menée pour des ondes solitaires modérément à fortement non-linéaires a confirmé la convergence exponentielle avec la résolution verticale grâce à l’approche spectrale, ainsi que les convergences algébriques en temps et en espace sur l’horizontale avec des ordres d’environ 4 (ou plus) en accord avec les schémas numériques utilisés.La comparaison des résultats du modèle à plusieurs jeux de données expérimentales a démontré les capacités du modèle à représenter les effets non-linéaires induits par les variations de bathymétrie, notamment les transferts d’énergie entre les composantes harmoniques, ainsi que la représentation précise des propriétés dispersives. Une formulation visco-potentielle a également été implémentée afin de prendre en compte les effets visqueux induits par la dissipation interne et le frottement sur le fond. Cette formulation a été validée dans le cas d’une faible viscosité avec un fond plat ou présentant une faible pente.Dans le but de représenter des champs de vagues 2DH, le modèle a été étendu en utilisant une discrétisation non-structurée (par nuage de points) du plan horizontal. Les dérivées horizontales ont été estimées à l’aide de la méthode RBF-FD (Radial Basis Function-Finite Difference), en conservant l’approche spectrale sur la verticale. Une étude numérique de sensibilité a été menée afin d’évaluer la robustesse de la méthode RBF-FD, en comparant différents types de RBFs, avec ou sans paramètre de forme et l’ajout éventuel d’un polynôme. La version 2DH du modèle a été utilisée pour simuler deux expériences en bassin, validant ainsi l’approche choisie et démontrant son applicabilité pour simuler la propagation 3D des vagues faisant intervenir des effets non-linéaires. Dans le but de réduire le temps de calcul et de pouvoir appliquer le code à des simulations sur de grands domaines, le code a été modifié pour utiliser le solveur linéaire direct en mode parallèle / In this work, a potential flow model based on the Euler-Zakharov equations was developed with the objective of simulating the propagation of irregular and multidirectional sea states from deep water conditions to the coast over variable bathymetry. A highly accurate representation of nonlinear and dispersive effects for bidimensional (2DH) nearshore and coastal domains on the order of kilometers is targeted.The preexisting 1DH version of the model, resolving the Laplace Boundary Value problem using a combination of high-order finite difference schemes in the horizontal direction and a spectral approach in the vertical direction, was improved and validated. The generation of incident waves through the implementation of specific Dirichlet and Neumann boundary conditions was studied in detail. In practice, these conditions were used in combination witha relaxation zone to improve the stability of the model.The linear dispersion relation of the model was derived analytically for the flat bottom case. Its analysis showed that the accuracy of the representation of dispersive effects improves significantly by increasing the vertical resolution (i.e. the maximum degree of the Chebyshev polynomial basis used to project the potential in the vertical). A convergence study conducted for moderate to highly nonlinear solitary waves confirmed the exponential convergence in the vertical dimension owing to the spectral approach, and the algebraic convergence in time and in space (horizontal dimension) with orders of about 4 (or higher) in agreement with the numerical schemes used.The capability of the model to represent nonlinear effects induced by variable bathymetry, such as the transfer of energy between harmonic components, as well as the accurate representation of dispersive properties, were demonstrated with comparisons to several experimental data sets. A visco-potential flow formulation was also implemented to take into account viscous effects induced by bulk viscosity and bottom friction. This formulation was validated inthe limit of small viscosity for mild slope bathymetries.To represent 2DH wave fields in complex nearshore domains, the model was extended using an unstructured discretization (scattered nodes) in the horizontal plane. The horizontal derivatives were estimated using the RBF-FD (Radial Basis Function - Finite Difference) method, while the spectral approach in the vertical remained unchanged. A series of sensitivity tests were conducted to evaluate numerically the robustness of the RBF-FD method, including a comparison of a variety of RBFs with or without shape factors and augmented polynomials. The 2DH version of the model was used to simulate two wave basin experiments, validating the approach and demonstrating the applicability of this method for 3D wave propagation, including nonlinear effects. As an initial attempt to improve the computational efficiency ofthe model for running simulations of large spatial domains, the code was adapted to use a parallelized direct linear solver

This thesis deals with neural modeling of electromagnetic field inside small aircrafts, witch can contain composite materials in their construction. Introduction to neural networks and its application in EMC of small airplanes is discussed in the first part of the text. In the second part of this thesis we design a simple EM model of small airplane. The airplane is simulated by two parallel dielectric layers (the left-hand side wall and the right hand side wall of the airplane). The layers are put into a rectangular metallic waveguide terminated by the absorber in order to simulate the illumination of the airplane by the external wave (both of the harmonic nature and pulse one). Numerical analyses are performed to search the relations between the distribution of an electromagnetic field inside the aircraft and electric parameters of model walls. The results of numerical analyses are used to train two types of neural network. In this way we can obtain accurate continuous model of electromagnetic field inside the aircraft. For the comparison with neural networks a multi-dimensional cubic spline interpolation is provided also. Neural classifiers are also investigated. We use them for classification of imaginary composite materials in terms of EMC. The nearest neighbour algorithm is applied as a classic approach to problem of classification.

Philosophiae Doctor - PhD / Many problems in financial world are being modeled by means of differential equation. These problems are time dependent, highly nonlinear, stochastic and heavily depend on the previous history of time. A variety of financial products exists in the market, such as forwards, futures, swaps and options. Our main focus in this thesis is to use the numerical analysis tools to solve some option pricing problems. Depending upon the inter-relationship of the financial derivatives, the dimension of the associated problem increases drastically and hence conventional methods (for example, the finite difference methods or finite element methods) for solving them do not provide satisfactory results. To resolve this issue, we use a special class of numerical methods, namely, the mesh free methods. These methods are often better suited to cope with changes in the geometry of the domain of interest than classical discretization techniques. In this thesis, we apply these methods to solve problems that price standard and non-standard options. We then extend the proposed approach to solve Heston's volatility model. The methods in each of these cases are analyzed for stability and thorough comparative numerical results are provided.

No description available.

In this study, the effect of a passive trapped-vortex cell on lift to drag (L/D) ratio of an FFA-W3-301 airfoil is studied. The upper surface of the airfoil was modified to incorporate a cavity defined by seven parameters. The L/D ratio of the airfoil is modeled using a radial basis function metamodel. This model is used to find the optimal design parameter values that give the highest L/D. The numerical results indicate that the L/D ratio is most sensitive to the position on an airfoil’s upper surface at which the cavity starts, the position of the end point of the cavity, and the vertical distance of the cavity end point relative to the airfoil surface. The L/D ratio can be improved by locating the cavity start point at the point of separation for a particular angle of attack. The optimal cavity shape (o19_aXX) is also tested for a NACA0024 airfoil.

Numerical methods for solving partial differential equations are commonplace in the engineering community and their popularity can be attributed to the rapid performance improvement of modern workstations and desktop computers. The ubiquity of computer technology has allowed all areas of engineering to have access to detailed thermal, stress, and fluid flow analysis packages capable of performing complex studies of current and future designs. The rapid pace of computer development, however, has begun to outstrip efforts to reduce analysis overhead. As such, most commercially available software packages are now limited by the human effort required to prepare, develop, and initialize the necessary computational models. Primarily due to the mesh-based analysis methods utilized in these software packages, the dependence on model preparation greatly limits the accessibility of these analysis tools. In response, the so-called meshless or mesh-free methods have seen considerable interest as they promise to greatly reduce the necessary human interaction during model setup. However, despite the success of these methods in areas demanding high degrees of model adaptability (such as crack growth, multi-phase flow, and solid friction), meshless methods have yet to gain notoriety as a viable alternative to more traditional solution approaches in general solution domains. Although this may be due (at least in part) to the relative youth of the techniques, another potential cause is the lack of focus on developing robust methodologies. The failure to approach development from a practical perspective has prevented researchers from obtaining commercially relevant meshless methodologies which reach the full potential of the approach. The primary goal of this research is to present a novel meshless approach called MIMS (Model Integrated Meshless Solver) which establishes the method as a generalized solution technique capable of competing with more traditional PDE methodologies (such as the finite element and finite volume methods). This was accomplished by developing a robust meshless technique as well as a comprehensive model generation procedure. By closely integrating the model generation process into the overall solution methodology, the presented techniques are able to fully exploit the strengths of the meshless approach to achieve levels of automation, stability, and accuracy currently unseen in the area of engineering analysis. Specifically, MIMS implements a blended meshless solution approach which utilizes a variety of shape functions to obtain a stable and accurate iteration process. This solution approach is then integrated with a newly developed, highly adaptive model generation process which employs a quaternary triangular surface discretization for the boundary, a binary-subdivision discretization for the interior, and a unique shadow layer discretization for near-boundary regions. Together, these discretization techniques are able to achieve directionally independent, automatic refinement of the underlying model, allowing the method to generate accurate solutions without need for intermediate human involvement. In addition, by coupling the model generation with the solution process, the presented method is able to address the issue of ill-constructed geometric input (small features, poorly formed faces, etc.) to provide an intuitive, yet powerful approach to solving modern engineering analysis problems.

Diese Dissertation beschreibt einen Fortschritt in der Vorhersage elektronischer Materialeigenschaften durch quantenchemische Methoden höherer Ordnung, mit Fokus auf der Berechnung elektronischer Bandlücken innerhalb des FHI-aims-Softwarepakets. Im Rahmen dieser Arbeit werden Coupled-Cluster-Methoden durch eine neue Schnittstelle zwischen FHI-aims und dem Coupled cluster for solids (Cc4s) Code implementiert. Dabei werden zwei zentrale Herausforderungen in Equation-of-Motion-Coupled-Cluster-Rechnungen mit Einfach- und Zweifachanregungen diskutiert: Konvergenz im Bezug auf den Basissatz und auf die Systemgröße. Für die Untersuchung der Basissatzunvollständigkeit werden numerische atomzentrierte Funktionen verwendet und Extrapolationstechniken unter Verwendung von valenz-korrelationskonsistenter Basissätze untersucht. Für Effekte endlicher Systemgröße wird der Equation-of-Motion-Coupled-Cluster-Strukturfaktor hergeleitet, um die Konvergenz der elektronischen Korrelation im langreichweitigen Limes zu analysieren, wobei Konvergenzraten für Ionisationspotentiale und Elektronenaffinitäten in verschiedenen Dimensionen hergeleitet werden. Die Arbeit identifiziert einen Zusammenhang zwischen der Konvergenz der Bandlücke der Equation-of-Motion-Coupled-Cluster- und der G0W0-Methode, was eine GW-gestützte Abschätzung verbleibender systemgrößenbedingter Fehler ermöglicht. Diese Methodik wird durch Bandlückenberechnungen für Materialien wie Diamant, Borphosphid, Lithiumhydrid und Magnesiumoxid validiert, wobei für einige Materialien eine verbesserte Genauigkeit nachgewiesen wird, während verbleibende Herausforderungen für Isolatoren mit großer Bandlücke identifiziert werden. / This thesis advances the prediction of electronic properties in materials through high-level quantum chemistry methods, focusing on electronic band gap calculations within the FHI-aims software package. The work implements coupled-cluster methods through a new interface between FHI-aims and the Coupled cluster for solids (Cc4s) code. The research addresses two primary challenges in equation-of-motion coupled-cluster calculations with single and double excitations: basis set convergence and system size convergence. For basis set incompleteness, the study evaluates numeric atom-centered functions and explores extrapolation techniques using valence-correlation consistent basis sets. For finite-size effects, the equation-of-motion coupled-cluster structure factor is derived to analyze electronic correlation convergence in the long-range limit, establishing convergence rates for ionization potentials and electron affinities across different dimensions. The work identifies a correlation between finite-size convergence in equation-of-motion coupled-cluster and G0W0 methods, enabling GW-aided estimation of remaining finite-size errors. The methodology is validated through band gap calculations for materials including diamond, boron phosphide, lithium hydride, and magnesium oxide, demonstrating improved accuracy for some materials while identifying persistent challenges for wide-gap insulators.

Engine mapping is the process of empirically modelling engine behaviour as a function of adjustable engine parameters, predicting the output of the engine. The aim is to calibrate the electronic engine controller to meet decreasing emission requirements and increasing fuel economy demands. Modern engines have an increasing number of control parameters that are having a dramatic impact on time and e ort required to obtain optimal engine calibrations. These are further complicated due to transient engine operating mode. A new model-based transient calibration method has been built on the application of hierarchical statistical modelling methods, and analysis of repeated experiments for the application of engine mapping. The methodology is based on two-stage regression approach, which organise the engine data for the mapping process in sweeps. The introduction of time-dependent covariates in the hierarchy of the modelling led to the development of a new approach for the problem of transient engine calibration. This new approach for transient engine modelling is analysed using a small designed data set for a throttle body inferred air ow phenomenon. The data collection for the model was performed on a transient engine test bed as a part of this work, with sophisticated software and hardware installed on it. Models and their associated experimental design protocols have been identi ed that permits the models capable of accurately predicting the desired response features over the whole region of operability. Further, during the course of the work, the utility of multi-layer perceptron (MLP) neural network based model for the multi-covariate case has been demonstrated. The MLP neural network performs slightly better than the radial basis function (RBF) model. The basis of this comparison is made on assessing relevant model selection criteria, as well as internal and external validation ts. Finally, the general ability of the model was demonstrated through the implementation of this methodology for use in the calibration process, for populating the electronic engine control module lookup tables.