1. |
- Frassinetti, Lorenzo, et al.
(author)
-
Vektoranalys
- 2019. - 1
-
Book (peer-reviewed)abstract
- Läroböcker i vektoranalys är ofta kortfattade. Denna bok, som kan användas för såväl grundläggande som mer avancerade kurser, behandlar ämnet mer utförligt.Bokens pedagogiska idé skiljer sig markant från liknande böcker. Återkommande inslag är tydligt formulerade problem som fångar det centrala i vektoranalysen. Syftet med problemen är dels att väcka intresse för den teori och de metoder som behandlas, dels att stimulera till aktivt lärande. Boken innehåller genomarbetade och lättillgängliga teoriavsnitt - som börjar med grundläggande vektoralgebra och slutar med kartesiska tensorer och en härledning av vektoranalysens huvudsats. Dessutom ingår ett stort antal konkreta exempel och många tillämpningar. Sist i varje kapitel finns en sammanfattning av den viktigaste teorin och övningsuppgifter med svar. Ledningar och fullständiga lösningar finns på Libers webbplats. Där finns även ett Appendix med tillämpningar.
|
|
2. |
- Lindvall, Kristoffer, 1990-, et al.
(author)
-
Can the Time-Spectral Method GWRM Advance Fusion Transport Modelling?
- 2017
-
In: 59th Annual Meeting of the APS Division of Plasma Physics.
-
Conference paper (peer-reviewed)abstract
- Transport phenomena in fusion plasma pose a daunting task for both real-time experiments and numerical modelling. The transport is driven by micro-instabilities caused by a host of unstable modes, for example ion temperature gradient and trapped electron modes. These modes can be modelled using fluid or gyrokinetic equations. However, the equations are characterised by high degrees of freedom and high temporal and spatial numerical requirements. Thus, a time-spectral method GWRM has been developed in order to efficiently solve these multiple time scale equations. The GWRM assumes a multivariate Chebyshev expansion ansatz in time, space, and parameter domain. Advantages are that time constraining CFL criteria no longer apply and that the solution accurately averages over small time-scale dynamics. For benchmarking, a two-fluid 2D drift wave turbulence model has been solved in order to study toroidal ion temperature gradient growth rates and nonlinear behaviour.
|
|
3. |
- Lindvall, Kristoffer, 1990-, et al.
(author)
-
Spectral Representation of Time and Physical Parameters in Numerical Weather Prediction
- 2018
-
In: Understanding of Atmospheric Systems with Efficient Numerical Methods for Observation and Prediction. - : IntechOpen.
-
Book chapter (peer-reviewed)abstract
- Numerical weather prediction (NWP) is a difficult task in chaotic dynamical regimes because of the strong sensitivity to initial conditions and physical parameters. As a result, high numerical accuracy is usually necessary. In this chapter, an accurate and efficient alternative to the traditional time stepping solution methods is presented; the time-spectral method. The generalized weighted residual method (GWRM) solves systems of nonlinear ODEs and PDEs using a spectral representation of time. Not being subject to CFL-like criteria, the GWRM typically employs time intervals two orders of magnitude larger than those of time-stepping methods. As an example, efficient solution of the chaotic Lorenz 1984 equations is demonstrated. The results indicate that the method has strong potential for NWP. Furthermore, employing spectral representations of physical parameters and initial values, families of solutions are obtained in a single computation. Thus, the GWRM is conveniently used for studies of system parameter dependency and initial condition error growth in NWP.
|
|
4. |
- Lindvalll, Kristoffer, et al.
(author)
-
A time-spectral method for initial-value problems using a novel spatial subdomain scheme
- 2018
-
In: Cogent Mathematics & Statistics. - : Informa UK Limited. - 2331-1835 .- 2574-2558. ; 5:1
-
Journal article (peer-reviewed)abstract
- We analyse a novel subdomain scheme for time-spectral solution of initial-value partial differential equations. In numerical modelling spectral methods are commonplace for spatially dependent systems, whereas finite difference schemes are typically applied for the temporal domain. The Generalized Weighted Residual Method (GWRM) is a fully spectral method that spectrally decomposes all specified domains, including the temporal domain, using multivariate Chebyshev polynomials. The Common Boundary-Condition method (CBC) here proposed is a spatial subdomain scheme for the GWRM. It solves the physical equations independently from the global connection of subdomains in order to reduce the total number of modes. Thus, it is a condensation procedure in the spatial domain that allows for a simultaneous global temporal solution. It is here evaluated against the finite difference methods of Crank-Nicolson and Lax-Wendroff for two example linear PDEs. The CBC-GWRM is also applied to the linearised ideal magnetohydrodynamic (MHD) equations for a screw pinch equilibrium. The growth rate of the most unstable mode was efficiently computed with an error <0.1%.
|
|
5. |
- Scheffel, Jan, 1954-, et al.
(author)
-
Generalized Weighted Residual Method; Advancements and Current Studies
- 2016
-
In: 58th Annual Meeting of the APS Division of Plasma Physics, San Jose, California, USA, 31 October – 4 November 2016.
-
Conference paper (peer-reviewed)abstract
- The Generalized Weighted Residual Method (GWRM) is a time-spectral method for solving initial value partial differential equations. The GWRM treats the temporal, spatial, and parameter domains by projecting the residual to a Chebyshev polynomial space, with the variational principle being that the residual is zero. This treatment provides a global semi-analytical solution. However, straightforward global solution is not economical. One remedy is the inclusion of spatial and temporal sub-domains with coupled internal boundary conditions, which decreases memory requirements and introduces sparse matrices. Only the equations pertaining to the boundary conditions need be solved globally, making the method parallelizable in time. Efficient solution of the linearized ideal MHD stability equations of screw-pinch equilibria are proved possible. The GWRM has also been used to solve strongly nonlinear ODEs such as the Lorenz equations (1984), and is capable of competing with finite time difference schemes in terms of both accuracy and efficiency. GWRM solutions of linear and nonlinear model problems of interest for stability and turbulence modelling will be presented, including detailed comparisons with time stepping methods.
|
|
6. |
- Scheffel, Jan, 1954-, et al.
(author)
-
Optimizing Time-Spectral Solution of Initial-Value Problems
- 2018
-
In: American Journal of Computational Mathematics. - : Scientific Research Publishing, Inc.. - 2161-1203 .- 2161-1211. ; 8:1, s. 7-26
-
Journal article (peer-reviewed)abstract
- Time-spectral solution of ordinary and partial differential equations is often regarded as an inefficient approach. The associated extension of the time domain, as compared to finite difference methods, is believed to result in uncomfortably many numerical operations and high memory requirements. It is shown in this work that performance is substantially enhanced by the introduction of algorithms for temporal and spatial subdomains in combination with sparse matrix methods. The accuracy and efficiency of the recently developed time spectral, generalized weighted residual method (GWRM) are compared to that of the explicit Lax-Wendroff and implicit Crank-Nicolson methods. Three initial-value PDEs are employed as model problems; the 1D Burger equation, a forced 1D wave equation and a coupled system of 14 linearized ideal magnetohydrodynamic (MHD) equations. It is found that the GWRM is more efficient than the time-stepping methods at high accuracies. The advantageous scalings Nt**1.0*Ns**1.43 and Nt**0.0*Ns**1.08 were obtained for CPU time and memory requirements, respectively, with Nt and Ns denoting the number of temporal and spatial subdomains. For time-averaged solution of the two-time-scales forced wave equation, GWRM performance exceeds that of the finite difference methods by an order of magnitude both in terms of CPU time and memory requirement. Favorable subdomain scaling is demonstrated for the MHD equations, indicating a potential for efficient solution of advanced initial-value problems in, for example, fluid mechanics and MHD.
|
|
7. |
- Scheffel, Jan, 1954-, et al.
(author)
-
SIR—An efficient solver for systems of equations
- 2018
-
In: Software Quality Professional. - : Elsevier B.V.. - 1522-0540. ; 7, s. 59-62
-
Journal article (peer-reviewed)abstract
- The Semi-Implicit Root solver (SIR) is an iterative method for globally convergent solution of systems of nonlinear equations. We here present MATLAB and MAPLE codes for SIR, that can be easily implemented in any application where linear or nonlinear systems of equations need be solved efficiently. The codes employ recently developed efficient sparse matrix algorithms and improved numerical differentiation. SIR convergence is quasi-monotonous and approaches second order in the proximity of the real roots. Global convergence is usually superior to that of Newton's method, being a special case of the method. Furthermore the algorithm cannot land on local minima, as may be the case for Newton's method with line search.
|
|