Pricing in stochastic-local volatility models with default

In recent years is becoming increasingly important to handle credit risk. Credit risk is the risk associated with the possibility of bankruptcy. More precisely, if a derivative provides for a payment at cert time T but before that time the counterparty defaults, at maturity the payment cannot be effectively performed, so the owner of the contract loses it entirely or a part of it. It means that the payoff of the derivative, and consequently its price, depends on the underlying of the basic derivative and on the risk of bankruptcy of the counterparty. To value and to hedge credit risk in a consistent way, one needs to develop a quantitative model. We have studied analytical approximation formulas and numerical methods such as Monte Carlo method in order to calculate the price of a bond. We have illustrated how to obtain fast and accurate pricing approximations by expanding the drift and diffusion as a Taylor series and we have compared the second and third order approximation of the Bond and Call price with an accurate Monte Carlo simulation. We have analysed JDCEV model with constant or stochastic interest rate. We have provided numerical examples that illustrate the effectiveness and versatility of our methods. We have used Wolfram Mathematica and Matlab.

Testing of subgrid scale (SGS) models for large-eddy simulation (LES) of turbulent channel flow

Sub-grid scale (SGS) models are required in order to model the influence of the unresolved small scales on the resolved scales in large-eddy simulations (LES), the flow at the smallest scales of turbulence. In the following work two SGS models are presented and deeply analyzed in terms of accuracy through several LESs with different spatial resolutions, i.e. grid spacings. The first part of this thesis focuses on the basic theory of turbulence, the governing equations of fluid dynamics and their adaptation to LES. Furthermore, two important SGS models are presented: one is the Dynamic eddy-viscosity model (DEVM), developed by \cite{germano1991dynamic}, while the other is the Explicit Algebraic SGS model (EASSM), by \cite{marstorp2009explicit}. In addition, some details about the implementation of the EASSM in a Pseudo-Spectral Navier-Stokes code \cite{chevalier2007simson} are presented. The performance of the two aforementioned models will be investigated in the following chapters, by means of LES of a channel flow, with friction Reynolds numbers $Re_\tau=590$ up to $Re_\tau=5200$, with relatively coarse resolutions. Data from each simulation will be compared to baseline DNS data. Results have shown that, in contrast to the DEVM, the EASSM has promising potentials for flow predictions at high friction Reynolds numbers: the higher the friction Reynolds number is the better the EASSM will behave and the worse the performances of the DEVM will be. The better performance of the EASSM is contributed to the ability to capture flow anisotropy at the small scales through a correct formulation for the SGS stresses. Moreover, a considerable reduction in the required computational resources can be achieved using the EASSM compared to DEVM. Therefore, the EASSM combines accuracy and computational efficiency, implying that it has a clear potential for industrial CFD usage.

One dimensional extended Hubbard model: two-particle bound states and resonances.

In questo lavoro di tesi è stato svolto uno studio analitico sul modello di Hubbard esteso unidimensionale al fine di osservare la presenza di eventuali risonanze che possano dare origine alla formazione di stati legati di due particelle. L'esistenza di uno stato legato stabile ha suscitato grande interesse negli ultimi anni, sia in ambito teorico che sperimentale, poichè è alla base di molti fenomeni che vengono osservati nei sistemi a molti corpi a basse temperature, come il BCS-BEC crossover. Pertanto si è ritenuto utile studiare il problema a due corpi nel modello di Hubbard esteso, che in generale non è integrabile. Il modello considerato contiene interazioni a primi e secondi vicini, in aggiunta all'interazione di contatto presente nel modello di Hubbard. Il problema è stato indagato analiticamente attraverso il Bethe ansatz, che consente di trovare tutti gli autovalori e le autofunzioni dell'Hamiltoniana. L'ansatz di Bethe sulla funzione d'onda è stato generalizzato per poter tener conto dei termini di interazione a più lungo raggio rispetto all'interazione di contatto. Si trova che, in questo modello, nel limite termodinamico, possono avvenire delle risonanze (o quasi-risonanze) in cui la lunghezza di scattering diverge, contrariamente a quanto avviene nel modello di Hubbard. Tale fenomeno si verifica quando il livello energetico discreto degli stati legati “tocca” la banda di scattering. Inoltre, con l'aggiunta di nuovi termini di interazione emergono nuovi stati legati. Nel caso in esame, si osservano due famiglie di stati legati, se lo spin totale delle due particelle è 1, e tre famiglie di stati legati, se lo spin totale è 0.