1,720 research outputs found
Towards a new generation of multi-dimensional stellar evolution models: development of an implicit hydrodynamic code
This paper describes the first steps of development of a new multidimensional
time implicit code devoted to the study of hydrodynamical processes in stellar
interiors. The code solves the hydrodynamical equations in spherical geometry
and is based on the finite volume method. Radiation transport is taken into
account within the diffusion approximation. Realistic equation of state and
opacities are implemented, allowing the study of a wide range of problems
characteristic of stellar interiors. We describe in details the numerical
method and various standard tests performed to validate the method. We present
preliminary results devoted to the description of stellar convection. We first
perform a local simulation of convection in the surface layers of a A-type star
model. This simulation is used to test the ability of the code to address
stellar conditions and to validate our results, since they can be compared to
similar previous simulations based on explicit codes. We then present a global
simulation of turbulent convective motions in a cold giant envelope, covering
80% in radius of the stellar structure. Although our implicit scheme is
unconditionally stable, we show that in practice there is a limitation on the
time step which prevent the flow to move over several cells during a time step.
Nevertheless, in the cold giant model we reach a hydro CFL number of 100. We
also show that we are able to address flows with a wide range of Mach numbers
(10^-3 < Ms< 0.5), which is impossible with an anelastic approach. Our first
developments are meant to demonstrate that the use of an implicit scheme
applied to a stellar evolution context is perfectly thinkable and to provide
useful guidelines to optimise the development of an implicit multi-D
hydrodynamical code.Comment: 21 pages, 18 figures, accepted for publication in A&
Monitoring a PGD solver for parametric power flow problems with goal-oriented error assessment
This is the peer reviewed version of the following article: [García-Blanco, R., Borzacchiello, D., Chinesta, F., and Diez, P. (2017) Monitoring a PGD solver for parametric power flow problems with goal-oriented error assessment. Int. J. Numer. Meth. Engng, 111: 529–552. doi: 10.1002/nme.5470], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nme.5470/full. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.The parametric analysis of electric grids requires carrying out a large number of Power Flow computations. The different parameters describe loading conditions and grid properties. In this framework, the Proper Generalized Decomposition (PGD) provides a numerical solution explicitly accounting for the parametric dependence. Once the PGD solution is available, exploring the multidimensional parametric space is computationally inexpensive. The aim of this paper is to provide tools to monitor the error associated with this significant computational gain and to guarantee the quality of the PGD solution. In this case, the PGD algorithm consists in three nested loops that correspond to 1) iterating algebraic solver, 2) number of terms in the separable greedy expansion and 3) the alternated directions for each term. In the proposed approach, the three loops are controlled by stopping criteria based on residual goal-oriented error estimates. This allows one for using only the computational resources necessary to achieve the accuracy prescribed by the end- user. The paper discusses how to compute the goal-oriented error estimates. This requires linearizing the error equation and the Quantity of Interest to derive an efficient error representation based on an adjoint problem. The efficiency of the proposed approach is demonstrated on benchmark problems.Peer ReviewedPostprint (author's final draft
CFD investigation of a complete floating offshore wind turbine
This chapter presents numerical computations for floating offshore wind turbines for a machine of 10-MW rated power. The rotors were computed using the Helicopter Multi-Block flow solver of the University of Glasgow that solves the Navier-Stokes equations in integral form using the arbitrary Lagrangian-Eulerian formulation for time-dependent domains with moving boundaries. Hydrodynamic loads on the support platform were computed using the Smoothed Particle Hydrodynamics method. This method is mesh-free, and represents the fluid by a set of discrete particles. The motion of the floating offshore wind turbine is computed using a Multi-Body Dynamic Model of rigid bodies and frictionless joints. Mooring cables are modelled as a set of springs and dampers. All solvers were validated separately before coupling, and the loosely coupled algorithm used is described in detail alongside the obtained results
- …