Search results for " Numerical simulation"

Multi-Phase epidemic model by a Markov chain

2008

Abstract In this paper we propose a continuous-time Markov chain to describe the spread of an infective and non-mortal disease into a community numerically limited and subjected to an external infection. We make a numerical simulation that shows tendencies for recurring epidemic outbreaks and for fade-out or extinction of the infection.

Newton algorithm for Hamiltonian characterization in quantum control

2014

We propose a Newton algorithm to characterize the Hamiltonian of a quantum system interacting with a given laser field. The algorithm is based on the assumption that the evolution operator of the system is perfectly known at a fixed time. The computational scheme uses the Crank-Nicholson approximation to explicitly determine the derivatives of the propagator with respect to the Hamiltonians of the system. In order to globalize this algorithm, we use a continuation method that improves its convergence properties. This technique is applied to a two-level quantum system and to a molecular one with a double-well potential. The numerical tests show that accurate estimates of the unknown paramete…

Surface waves on cylindrical solids: numerical and experimental study.

2013

The use of Rayleigh waves enables the solution of several important inspection problems. Propagation of surface waves along straight boundaries has been properly studied but investigations about their propagation on cylindrical surfaces are not sufficient, despite they can be still of interest for NDE applications. It has been proved experimentally that a surface wave pulse suffers a phase shift during its propagation along a cylindrical surface. A numerical approach has been developed to efficiently study these effects for different materials, curvatures and frequencies. This study can help the scientific community to better understand the phenomenon, quite complex and not yet fully explor…

Steady, periodic, quasi-periodic and chaotic flow regimes in toroidal pipes

2012

Incompressible flow in a toroidal pipe was investigated by direct numerical simulation. The curvature a/c (radius of the cross section / radius of the torus) was 0.3 or 0.1 and the bulk Reynolds number ranged between 3500 and 14 700. The study revealed a rich scenario of transition to turbulence. For the higher curvature a/c = 0.3, a supercritical transition from stationary to periodic flow (Hopf bifurcation) was observed at Re=4600. The periodic flow was characterized by a travelling wave which, in the whole periodic Re range, took the form of a varicose modulation of the twin Dean vortex rings, included 8 wavelengths along the axis of the torus, and exhibited instantaneous anti-symmetry a…

Boundary-layer Flows Past an Hemispherical Roughness Element: DNS, Global Stability and Sensitivity Analysis

2015

Abstract We investigate the full three-dimensional instability mechanism arising in the wake of an hemispherical roughness element immersed in a laminar Blasius boundary layer. The inherent three-dimensional flow pattern beyond the critical Reynolds number is characterized by coherent vortical structures called hairpin vortices. Direct numerical simulation is used to analyze the formation and the shedding of hairpin packets inside the shear layer. The first bifurcation characteristics are investigated by global stability tools. We show the spatial structure of the linear direct and adjoint global eigenmodes of the linearized Navier-Stokes operator and use structural sensitivity analysis to …

A dynamic subgrid-scale tensorial eddy viscosity model

1999

In the Navier-Stokes equations the removal of the turbulent fluctuating velocities with a frequency above a certain fixed threshold, employed in the Large Eddy Simulation (LES), causes the appearance of a turbulent stress tensor that requires a number of closure assumptions. In this paper insufficiencies are demonstrated for those closure models which are based on a scalar eddy viscosity coefficient. A new model, based on a tensorial eddy viscosity, is therefore proposed; it employs the Germano identity [1] and allows dynamical evaluation of the single required input coefficient. The tensorial expression for the eddy viscosity is deduced by removing the widely used scalar assumption of the …

Direct Numerical Simulation of Pulsatile Turbulent Channel Flow

2008

Simulation of laser generated ultrasound with application to defect detection

2008

Laser generated ultrasound holds substantial promise for use as a tool for defect detection in remote inspection thanks to its ability to produce frequencies in the MHz range, enabling fine spatial resolution of defects. Despite the potential impact of laser generated ultrasound in many areas of science and industry, robust tools for studying the phenomenon are lacking and thus limit the design and optimization of non-destructive testing and evaluation techniques. The laser generated ultrasound propagation in complex structures is an intricate phenomenon and is extremely hard to analyze. Only simple geometries can be studied analytically. Numerical techniques found in the literature have pr…

Optimal valves regulation for calibration of pipe networks models

1998

The calibration of a pipe network model is faced by means of a two-steps iterative procedure. The first step of the procedure is the estimation of the parameters in the network simulation model, obtained from a first set of water head measures. The second step of the procedure is the change of the resistances provided by the valves in the network and the acquisition of a new set of measures. The inverse problem is then solved again and the procedure is repeated, until the set of computed parameters converge to stable values. The new resistances are chosen to minimize the sum of the variation coefficients, computed from the Fisher information matrix. The effect of different correlation funct…

Numerical model for composite material with polymer matrix reinforced by carbon nanotubes

2008

Due to the high stiffness and strength, as well as their ability to act as conductors, carbon nanotubes are under intense investigation as fillers in polymeric materials. The nature of the carbon nanotube/polymer bonding and the curvature of the carbon nanotubes within the polymer have arisen as particular factors in the efficacy of the carbon nanotubes to actually provide any enhanced stiffness or strength to the nanocomposite. Here the effects of carbon nanotube curvature and interface interaction with the matrix on the nanocomposite stiffness are investigated using nanomechanical analysis. In particular, the effects of poor bonding and thus poor shear lag load transfer to the carbon nano…