Comparison of natural and forced amplification regimes in plasma-based soft-x-ray lasers seeded by high-order harmonic

The amplification of high-order harmonics (HOH) in a plasma-based amplifier is a multiscale, temporal phenomenon that couples plasma hydrodynamics, atomic processes, and HOH electromagnetic fields. We use a one-dimensional, time-dependent Maxwell-Bloch code to compare the natural amplification regime and another regime where plasma polarization is constantly forced by the HOH. In this regime, a 10-MW (i.e., 100 times higher than current seeded soft x-ray laser power), 1.5-μJ, 140-fs pulse free from the parasitic temporal structures appearing on the natural amplification regime can be obtained.

Equation of state for hot dense matter using a relativistic screened hydrogenic model

The study of matter under conditions of high density, pressure, and temperature is a valuable subject for inertial confinement fusion (ICF), astrophysical phenomena, high-power laser interaction with matter, etc. In all these cases, matter is heated and compressed by strong shocks to high pressures and temperatures, becomes partially or completely ionized via thermal or pressure ionization, and is in the form of dense plasma. The thermodynamics and the hydrodynamics of hot dense plasmas cannot be predicted without the knowledge of the equation of state (EOS) that describes how a material reacts to pressure and how much energy is involved. Therefore, the equation of state often takes the form of pressure and energy as functions of density and temperature. Furthermore, EOS data must be obtained in a timely manner in order to be useful as input in hydrodynamic codes. By this reason, the use of fast, robust and reasonably accurate atomic models, is necessary for computing the EOS of a material.

WCSPH viscosity diffusion processes in vortex flows

The aim of this paper is to clarify the role played by the most commonly used viscous terms in simulating viscous laminar flows using the weakly compressible approach in the context of smooth particle hydrodynamics (WCSPH). To achieve this, Takeda et al. (Prog. Theor. Phys. 1994; 92(5):939–960), Morris et al. (J. Comput. Phys. 1997; 136:214–226) and Monaghan–Cleary–Gingold's (Appl. Math. Model. 1998; 22(12):981–993; Monthly Notices of the Royal Astronomical Society 2005; 365:199–213) viscous terms will be analysed, discussing their origins, structures and conservation properties. Their performance will be monitored with canonical flows of which related viscosity phenomena are well understood, and in which boundary effects are not relevant. Following the validation process of three previously published examples, two vortex flows of engineering importance have been studied. First, an isolated Lamb–Oseen vortex evolution where viscous effects are dominant and second, a pair of co-rotating vortices in which viscous effects are combined with transport phenomena. The corresponding SPH solutions have been compared to finite-element numerical solutions. The SPH viscosity model's behaviour in modelling the viscosity related effects for these canonical flows is adequate

Comparison of experimental and numerical sloshing loads in partially filled tanks

Sloshing describes the movement of liquids inside partially filled tanks, generating dynamic loads on the tank structure. The resulting impact pressures are of great importance in assessing structural strength, and their correct evaluation still represents a challenge for the designer due to the high level of nonlinearities involved, with complex free surface deformations, violent impact phenomena and influence of air trapping. In the present paper, a set of two-dimensional cases, for which experimental results are available, is considered to assess the merits and shortcomings of different numerical methods for sloshing evaluation, namely two commercial RANS solvers (FLOW-3D and LS-DYNA), and two academic software (Smoothed Particle Hydrodynamics and RANS). Impact pressures at various critical locations and global moment induced by water motion in a partially filled rectangular tank, subject to a simple harmonic rolling motion, are evaluated and predictions are compared with experimental measurements. 2012 Copyright Taylor and Francis Group, LLC.

Systematic Analysis of Direct-Drive Baseline Designs for Shock-Ignition with the Laser Megajoule

Direct-drive inertial confinement thermonuclear fusion consists in illuminating a shell of cryogenic Deuterium and Tritium (DT) mixture with many intense beams of laser light. Capsule is composed of DT gassurrounded by cryogenic DT as combustible fuel. Basic rules are used to define shell geometry from aspect ratio, fuel mass and layers densities. We define baseline designs using two aspect ratio (A=3 and A=5) who complete HiPER baseline design (A=7.7). Aspect ratio is defined as the ratio of ice DT shell inner radius over DT shell thickness. Low aspect ratio improves hydrodynamics stabilities of imploding shell. Laser impulsion shape and ablator thickness are initially defined by using Lindl (1995) pressure ablation and mass ablation formulae for direct-drive using CH layer as ablator. In flight adiabat parameter is close to one during implosion. Velocitie simplosions chosen are between 260 km/s and 365 km/s. More than thousand calculations are realized for each aspect ratio in order to optimize the laser pulse shape. Calculations are performed using the one-dimensional version of the Lagrangian radiation hydrodynamics FCI2. We choose implosion velocities for each initial aspect ratio, and we compute scaled-target family curves for each one to find self-ignition threshold. Then, we pick points on each curves that potentially product high thermonuclear gain and compute shock ignition in the context of Laser MegaJoule. This systematic analyze reveals many working points which complete previous studies ´allowing to highlight baseline designs, according to laser intensity and energy, combustible mass and initial aspect ratio to be relevant for Laser MegaJoule.

Particle methods parallel implementations by GP-GPU strategies

This paper outlines the problems found in the parallelization of SPH (Smoothed Particle Hydrodynamics) algorithms using Graphics Processing Units. Different results of some parallel GPU implementations in terms of the speed-up and the scalability compared to the CPU sequential codes are shown. The most problematic stage in the GPU-SPH algorithms is the one responsible for locating neighboring particles and building the vectors where this information is stored, since these specific algorithms raise many dificulties for a data-level parallelization. Because of the fact that the neighbor location using linked lists does not show enough data-level parallelism, two new approaches have been pro- posed to minimize bank conflicts in the writing and subsequent reading of the neighbor lists. The first strategy proposes an efficient coordination between CPU-GPU, using GPU algorithms for those stages that allow a straight forward parallelization, and sequential CPU algorithms for those instructions that involve some kind of vector reduction. This coordination provides a relatively orderly reading of the neighbor lists in the interactions stage, achieving a speed-up factor of x47 in this stage. However, since the construction of the neighbor lists is quite expensive, it is achieved an overall speed-up of x41. The second strategy seeks to maximize the use of the GPU in the neighbor's location process by executing a specific vector sorting algorithm that allows some data-level parallelism. Al- though this strategy has succeeded in improving the speed-up on the stage of neighboring location, the global speed-up on the interactions stage falls, due to inefficient reading of the neighbor vectors. Some changes to these strategies are proposed, aimed at maximizing the computational load of the GPU and using the GPU texture-units, in order to reach the maximum speed-up for such codes. Different practical applications have been added to the mentioned GPU codes. First, the classical dam-break problem is studied. Second, the wave impact of the sloshing fluid contained in LNG vessel tanks is also simulated as a practical example of particle methods

Impulsive motions of the floating zone

A set of problems concerning the behaviour of a suddenly disturbed ideal floating zone is considered. Mathematical techniques of asymptotic expansions arc used to solve these problems. It is seen that many already available solutions, most of them concerning liquids enclosed in cavities, will be regarded as starting approximations which are valid except in the proximity of the free surface which laterally bounds the floating zone. In particular, the problem of the linear spin-up of an initially cylindrical floating zone is considered in some detail. The presence of a recircuiating fluid pattern near the free surface is detected. This configuration is attributed to the interplay between Coriolis forces and the azimuthal component of the viscous forces.

Conceptual design of the liquid metal laboratory of the TECHNOFUSION facility

The application of liquid metal technology in fusion devices requires R&D related to many phenomena: interaction between liquid metals and structural material as corrosion, erosion and passivation techniques; magneto-hydrodynamics; free surface fluid-dynamics and any other physical aspect that will be needed for their safe reliable operation. In particular, there is a significant shortage of experimental facilities dedicated to the development of the lithium technology. In the framework of the TECHNOFUSION project, an experimental laboratory devoted to the lithium technology development is proposed, in order to shed some light in the path to IFMIF and the design of chamber's first wall and divertors. The conceptual design foresee a development in two stages, the first one consisting on a material testing loop. The second stage proposes the construction of a mock-up of the IFMIF target that will allow to assess the behaviour of a free-surface lithium target under vacuum conditions. In this paper, such conceptual design is addressed.

Self-similar motion of laser fusion plasmas. Absorption in an unbounded plasma

The one-dimensional motion generated in a cold, infinite, uniform plasma of density na by the absorption, in a certain plane, of a linear pulse of energy per unit time and area = 4>0t/r, 0< t< r, is considered, the analysis allows for thermal conduction and viscosity of ions and electrons, their energy exchange, and an electron heat flux limiter The resulting motion is self-similar and governed by a single nondimensional parameter a«(n0 2T/0)2/3 Detailed asymptotic results are obtained for both a < l and a > l , the general behavior of the solution for arbitrary a is discussed The analysis can be extended to the case of a plasma initially occupying a half-space, and throws light on how to optimize the hydrodynamics of laser fusion plasmas Known approximate results corresponding to motion of a plasma submitted to constant irradiation (<()) are recovered in the present work under appropriate limiting processes

New steps in the development of the second generation TEC GESMEY

The paper shows the results of the new steps that have been done in the development of the tidal energy converter GESMEY. These are the design, construction and trials into the sea of a 1/10 scale prototype and also the construction with the same scale of the buoy BOSCEM, that anchors the device and lets it in the correct work position and depth, along the two directions of the flow that the daily tidal cycle have. Inside the paper is described the objectives and the methodology of the experimental trials that were ca rry out the last summer with the scale prototype. GESMEY is a new type of tidal energy converter (TEC) that has the capability to exploit currents in waters over forty meters by itself and it gets only using its internal ballast system the necessary equilibrium between hy drostatics and hydrodynamics forces to make the emersion and the immersion procedures without any other help. Finally the paper shows the description of the results obtained over the performance of the devices along the immersion, emersion and floating transport manoeuvres and afterwards the results, that were obtained along the generation power tests that were carried out, are shown.

Finite volume modelling of an icefield with subglacial lake

Momentum, mass and energy balance laws provide the tools for the study of the evolution of an icefield covering a subglacial lake. The ice is described as a non-Newtonian fluid with a power-law constitutive relationship with temperature- and stress-dependent viscosity (Glen?s law) [1]. The phase transition mechanisms at the air/ice and ice/water interfaces yield moving boundary formulations, and lake hydrodynamics requires equation reduction for treating the turbulence.

Modelling of landslides propagation with SPH : effects of rheology and pore water pressure

El estudio desarrollado en este trabajo de tesis se centra en la modelización numérica de la fase de propagación de los deslizamientos rápidos de ladera a través del método sin malla Smoothed Particle Hydrodynamics (SPH). Este método tiene la gran ventaja de permitir el análisis de problemas de grandes deformaciones evitando operaciones costosas de remallado como en el caso de métodos numéricos con mallas tal como el método de los Elementos Finitos. En esta tesis, particular atención viene dada al rol que la reología y la presión de poros desempeñan durante estos eventos. El modelo matemático utilizado se basa en la formulación de Biot-Zienkiewicz v - pw, que representa el comportamiento, expresado en términos de velocidad del esqueleto sólido y presiones de poros, de la mezcla de partículas sólidas en un medio saturado. Las ecuaciones que gobiernan el problema son: • la ecuación de balance de masa de la fase del fluido intersticial, • la ecuación de balance de momento de la fase del fluido intersticial y de la mezcla, • la ecuación constitutiva y • una ecuación cinemática. Debido a sus propiedades geométricas, los deslizamientos de ladera se caracterizan por tener una profundidad muy pequeña frente a su longitud y a su anchura, y, consecuentemente, el modelo matemático mencionado anteriormente se puede simplificar integrando en profundidad las ecuaciones, pasando de un modelo 3D a 2D, el cual presenta una combinación excelente de precisión, sencillez y costes computacionales. El modelo propuesto en este trabajo se diferencia de los modelos integrados en profundidad existentes por incorporar un ulterior modelo capaz de proveer información sobre la presión del fluido intersticial a cada paso computacional de la propagación del deslizamiento. En una manera muy eficaz, la evolución de los perfiles de la presión de poros está numéricamente resuelta a través de un esquema explicito de Diferencias Finitas a cada nodo SPH. Este nuevo enfoque es capaz de tener en cuenta la variación de presión de poros debida a cambios de altura, de consolidación vertical o de cambios en las tensiones totales. Con respecto al comportamiento constitutivo, uno de los problemas principales al modelizar numéricamente deslizamientos rápidos de ladera está en la dificultad de simular con la misma ley constitutiva o reológica la transición de la fase de iniciación, donde el material se comporta como un sólido, a la fase de propagación donde el material se comporta como un fluido. En este trabajo de tesis, se propone un nuevo modelo reológico basado en el modelo viscoplástico de Perzyna, pensando a la viscoplasticidad como a la llave para poder simular tanto la fase de iniciación como la de propagación con el mismo modelo constitutivo. Con el fin de validar el modelo matemático y numérico se reproducen tanto ejemplos de referencia con solución analítica como experimentos de laboratorio. Finalmente, el modelo se aplica a casos reales, con especial atención al caso del deslizamiento de 1966 en Aberfan, mostrando como los resultados obtenidos simulan con éxito estos tipos de riesgos naturales. The study developed in this thesis focuses on the modelling of landslides propagation with the Smoothed Particle Hydrodynamics (SPH) meshless method which has the great advantage of allowing to deal with large deformation problems by avoiding expensive remeshing operations as happens for mesh methods such as, for example, the Finite Element Method. In this thesis, special attention is given to the role played by rheology and pore water pressure during these natural hazards. The mathematical framework used is based on the v - pw Biot-Zienkiewicz formulation, which represents the behaviour, formulated in terms of soil skeleton velocity and pore water pressure, of the mixture of solid particles and pore water in a saturated media. The governing equations are: • the mass balance equation for the pore water phase, • the momentum balance equation for the pore water phase and the mixture, • the constitutive equation and • a kinematic equation. Landslides, due to their shape and geometrical properties, have small depths in comparison with their length or width, therefore, the mathematical model aforementioned can then be simplified by depth integrating the equations, switching from a 3D to a 2D model, which presents an excellent combination of accuracy, computational costs and simplicity. The proposed model differs from previous depth integrated models by including a sub-model able to provide information on pore water pressure profiles at each computational step of the landslide's propagation. In an effective way, the evolution of the pore water pressure profiles is numerically solved through a set of 1D Finite Differences explicit scheme at each SPH node. This new approach is able to take into account the variation of the pore water pressure due to changes of height, vertical consolidation or changes of total stress. Concerning the constitutive behaviour, one of the main issues when modelling fast landslides is the difficulty to simulate with the same constitutive or rheological model the transition from the triggering phase, where the landslide behaves like a solid, to the propagation phase, where the landslide behaves in a fluid-like manner. In this work thesis, a new rheological model is proposed, based on the Perzyna viscoplastic model, thinking of viscoplasticity as the key to close the gap between the triggering and the propagation phase. In order to validate the mathematical model and the numerical approach, benchmarks and laboratory experiments are reproduced and compared to analytical solutions when possible. Finally, applications to real cases are studied, with particular attention paid to the Aberfan flowslide of 1966, showing how the mathematical model accurately and successfully simulate these kind of natural hazards.

Modelling of short runout propagation landslides and debris flows

This paper proposes an extension of methods used to predict the propagation of landslides having a long runout to smaller landslides with much shorter propagation distances. The method is based on: (1) a depth-integrated mathematical model including the coupling between the soil skeleton and the pore fluids, (2) suitable rheological models describing the relation between the stress and the rate of deformation tensors for fluidised soils and (3) a meshless numerical method, Smooth Particle Hydrodynamics, which separates the computational mesh (or set of computational nodes) from the mesh describing the terrain topography, which is of structured type – thus accelerating search operations. The proposed model is validated using two examples for which there are analytical solutions, and then it is applied to two short runout landslides which happened in Hong Kong in 1995, for which there is available information.