Search results for "Solver"
showing 10 items of 157 documents
A marching in space and time (MAST) solver of the shallow water equations. Part II: The 2D model
2007
Abstract A novel methodology for the solution of the 2D shallow water equations is proposed. The algorithm is based on a fractional step decomposition of the original system in (1) a convective prediction, (2) a convective correction, and (3) a diffusive correction step. The convective components are solved using a Marching in Space and Time (MAST) procedure, that solves a sequence of small ODEs systems, one for each computational cell, ordered according to the cell value of a scalar approximated potential. The scalar potential is sought after computing first the minimum of a functional via the solution of a large linear system and then refining locally the optimum search. Model results are…
High-quality computational tools for linear-algebra problems in FEM electromagnetic simulation [EM Programmer's Notebook]
2004
A key ingredient of finite-element analysis programs is the linear-algebra solver, typically either a linear-system solver or an eigensolver. The first part of This work tries to justify why it is important to have recourse to publicly available software for addressing this part of the computation. A number of libraries are mentioned as successful examples that exhibit a series of desirable qualities. Although some of these libraries force the programmer to somewhat change the programming style and may be difficult to learn, the benefits usually pay off the extra effort. The second part of the paper describes one of these libraries in some detail, namely SLEPc, the Scalable Library for Eige…
A Column Generation Approach to Scheduling of Periodic Tasks
2011
We present an algorithm based on column generation for a real time scheduling problem, in which all tasks appear regularly after a given period. Furthermore, the tasks exchange messages, which have to be transferred over a bus, if the tasks involved are executed on different ECUs. Experiments show that for large instances our preliminary implementation is faster than the previous approach based on an integer linear programming formulation using a state-of-the-art solver.
Detection of local tourism systems by threshold accepting
2015
Despite the importance of tourism as a leading industry in the development of a country’s economy, there is a lack of criteria and methodologies for the detection, promotion, and governance of local tourism systems. We propose a quantitative approach for the detection of local tourism systems the size of which is optimal with respect to geographical, economic, and demographical criteria: we formulate the problem as an optimisation problem and we solve it by a metaheuristic approach; then we compare the obtained results with standard clustering approaches and with an exact optimisation solver. Results show that our approach requires low computational times to provide results that are better …
Multicopter UAV design optimization
2014
Designing and selecting hardware for a multirotor can be challenging in order to get the best flight performance out of the system. In addition to selecting the hardware, the number of actuators can also be altered. For a 4 actuator (quadrotor) setup, one set of hardware can give the optimal design, while for a 6 actuator setup (hexarotor) the same hardware may not necessarily give the same response. In this paper we present a design optimization process of a multirotor, where the hardware is selected from a set of low-cost off-the-shelf standard RC hobby parts. Constraining the problem to a given hardware ensures existence of the selected hardware, and the design can be implemented. Also t…
Efficient numerical integration of neutrino oscillations in matter
2016
A special purpose solver, based on the Magnus expansion, well suited for the integration of the linear three neutrino oscillations equations in matter is proposed. The computations are speeded up to two orders of magnitude with respect to a general numerical integrator, a fact that could smooth the way for massive numerical integration concomitant with experimental data analyses. Detailed illustrations about numerical procedure and computer time costs are provided.
High-order methods for the simulation of hydromagnetic instabilities in core-collapse supernovae
2011
AbstractWe present an assessment of the accuracy of a recently developed MHD code used to study hydromagnetic flows in supernovae and related events. The code, based on the constrained transport formulation, incorporates unprecedented ultra-high-order methods (up to 9th order) for the reconstruction and the most accurate approximate Riemann solvers. We estimate the numerical resistivity of these schemes in tearing instability simulations.
IBSIMU: a three-dimensional simulation software for charged particle optics.
2010
A general-purpose three-dimensional (3D) simulation code IBSIMU for charged particle optics with space charge is under development at JYFL. The code was originally developed for designing a slit-beam plasma extraction and nanosecond scale chopping for pulsed neutron generator, but has been developed further and has been used for many applications. The code features a nonlinear FDM Poisson's equation solver based on fast stabilized biconjugate gradient method with ILU0 preconditioner for solving electrostatic fields. A generally accepted nonlinear plasma model is used for plasma extraction. Magnetic fields can be imported to the simulations from other programs. The particle trajectories are …
Numerical Relativistic Hydrodynamics
2008
High Resolution Shock Capturing (HRSC) techniques achieve highly accurate numerical approximations (formally second order or better) in smooth regions of the flow, and capture the motion of unresolved steep gradients without creating spurious oscillations. I will show how these techniques have been extended to relativistic hydrodynamics, making it possible to explore some challenging astrophysical scenarios. I will review recent literature concerning the main properties of different special relativistic Riemann solvers, and discuss several test problems which are commonly used to evaluate the performance of numerical methods in relativistic hydrodynamics. In the second part, I will illustra…
Correlation effects in the total energy, the bulk modulus, and the lattice constant of a transition metal: Combined local-density approximation and d…
2009
We present an accurate implementation of total-energy calculations into the local-density approximation plus dynamical mean-field theory $(\text{LDA}+\text{DMFT})$ method. The electronic structure problem is solved through the full-potential linear muffin-tin orbital and Korringa-Kohn-Rostoker methods with a perturbative solver for the effective impurity suitable for moderately correlated systems. We have tested the method in detail for the case of Ni, and investigated the sensitivity of the results to the computational scheme and to the complete self-consistency. It is demonstrated that the $\text{LDA}+\text{DMFT}$ method can resolve a long-standing controversy between the LDA/generalized …