| 1. |
- Almquist, Martin, et al.
(författare)
-
Atmospheric Sound Propagation Over Large-Scale Irregular Terrain
- 2014
-
Ingår i: Journal of Scientific Computing. - : Springer Science and Business Media LLC. - 0885-7474 .- 1573-7691. ; 61:2, s. 369-397
-
Tidskriftsartikel (refereegranskat)abstract
- A benchmark problem on atmospheric sound propagation over irregular terrain has been solved using a stable fourth-order accurate finite difference approximation of a high-fidelity acoustic model. A comparison with the parabolic equation method and ray tracing methods is made. The results show that ray tracing methods can potentially be unreliable in the presence of irregular terrain.
|
|
| 2. |
|
|
| 3. |
- Almquist, Martin, 1988-
(författare)
-
Efficient Simulation of Wave Phenomena
- 2017
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Wave phenomena appear in many fields of science such as acoustics, geophysics, and quantum mechanics. They can often be described by partial differential equations (PDEs). As PDEs typically are too difficult to solve by hand, the only option is to compute approximate solutions by implementing numerical methods on computers. Ideally, the numerical methods should produce accurate solutions at low computational cost. For wave propagation problems, high-order finite difference methods are known to be computationally cheap, but historically it has been difficult to construct stable methods. Thus, they have not been guaranteed to produce reasonable results.In this thesis we consider finite difference methods on summation-by-parts (SBP) form. To impose boundary and interface conditions we use the simultaneous approximation term (SAT) method. The SBP-SAT technique is designed such that the numerical solution mimics the energy estimates satisfied by the true solution. Hence, SBP-SAT schemes are energy-stable by construction and guaranteed to converge to the true solution of well-posed linear PDE. The SBP-SAT framework provides a means to derive high-order methods without jeopardizing stability. Thus, they overcome most of the drawbacks historically associated with finite difference methods.This thesis consists of three parts. The first part is devoted to improving existing SBP-SAT methods. In Papers I and II, we derive schemes with improved accuracy compared to standard schemes. In Paper III, we present an embedded boundary method that makes it easier to cope with complex geometries. The second part of the thesis shows how to apply the SBP-SAT method to wave propagation problems in acoustics (Paper IV) and quantum mechanics (Papers V and VI). The third part of the thesis, consisting of Paper VII, presents an efficient, fully explicit time-integration scheme well suited for locally refined meshes.
|
|
| 4. |
|
|
| 5. |
|
|
| 6. |
- Dao, Tuan Anh, et al.
(författare)
-
Energy stable and accurate coupling of finite element methods and finite difference methods
- 2022
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 449
-
Tidskriftsartikel (refereegranskat)abstract
- We introduce a hybrid method to couple continuous Galerkin finite element methods and high-order finite difference methods in a nonconforming multiblock fashion. The aim is to optimize computational efficiency when complex geometries are present. The proposed coupling technique requires minimal changes in the existing schemes while maintaining strict stability, accuracy, and energy conservation. Results are demonstrated on linear and nonlinear scalar conservation laws in two spatial dimensions.
|
|
| 7. |
- Duru, Kenneth, 1975-
(författare)
-
Perfectly Matched Layers and High Order Difference Methods for Wave Equations
- 2012
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- The perfectly matched layer (PML) is a novel technique to simulate the absorption of waves in unbounded domains. The underlying equations are often a system of second order hyperbolic partial differential equations. In the numerical treatment, second order systems are often rewritten and solved as first order systems. There are several benefits with solving the equations in second order formulation, though. However, while the theory and numerical methods for first order hyperbolic systems are well developed, numerical techniques to solve second order hyperbolic systems are less complete.We construct a strongly well-posed PML for second order systems in two space dimensions, focusing on the equations of linear elasto-dynamics. In the continuous setting, the stability of both first order and second order formulations are linearly equivalent. We have found that if the so-called geometric stability condition is violated, approximating the first order PML with standard central differences leads to a high frequency instability at most resolutions. In the second order setting growth occurs only if growing modes are well resolved. We determine the number of grid points that can be used in the PML to ensure a discretely stable PML, for several anisotropic elastic materials.We study the stability of the PML for problems where physical boundaries are important. First, we consider the PML in a waveguide governed by the scalar wave equation. To ensure the accuracy and the stability of the discrete PML, we derived a set of equivalent boundary conditions. Second, we consider the PML for second order symmetric hyperbolic systems on a half-plane. For a class of stable boundary conditions, we derive transformed boundary conditions and prove the stability of the corresponding half-plane problem. Third, we extend the stability analysis to rectangular elastic waveguides, and demonstrate the stability of the discrete PML.Building on high order summation-by-parts operators, we derive high order accurate and strictly stable finite difference approximations for second order time-dependent hyperbolic systems on bounded domains. Natural and mixed boundary conditions are imposed weakly using the simultaneous approximation term method. Dirichlet boundary conditions are imposed strongly by injection. By constructing continuous strict energy estimates and analogous discrete strict energy estimates, we prove strict stability.
|
|
| 8. |
|
|
| 9. |
|
|
| 10. |
- Eriksson, Gustav, et al.
(författare)
-
Boundary and interface methods for energy stable finite difference discretizations of the dynamic beam equation
- 2023
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 476
-
Tidskriftsartikel (refereegranskat)abstract
- We consider energy stable summation by parts finite difference methods (SBP-FD) for the homogeneous and piecewise homogeneous dynamic beam equation (DBE). Previously the constant coefficient problem has been solved with SBP-FD together with penalty terms (SBP-SAT) to impose boundary conditions. In this work, we revisit this problem and compare SBP-SAT to the projection method (SBP-P). We also consider the DBE with discontinuous coefficients and present novel SBP-SAT, SBP-P, and hybrid SBP-SAT-P discretizations for imposing interface conditions. To demonstrate the methodology for two-dimensional problems, we also present a discretization of the piecewise homogeneous dynamic Kirchoff-Love plate equation based on the hybrid SBP-SAT-P method. Numerical experiments show that all methods considered are similar in terms of accuracy, but that SBP-P can be more computationally efficient (less restrictive time step requirement for explicit time integration methods) for both the constant and piecewise constant coefficient problems.
|
|
