Global ETD Search

11	High-Accuracy and Stable Finite Difference Methods for Solving the Acoustic Wave Equation Boughanmi, Aimen January 2024 (has links) This report presents a comprehensive investigation into the accuracy and stability of Finite Difference Methods (FDM) when applied to the acoustic wave equation. The analysis focuses on comparing the classical 2nd order FDM with highly-accurate computational stencil of order 2p = 2,4 and 6 with Summation-by-Parts (SBP) and Simultaneous Approximation Term (SAT) technique of the Finite Difference Method. The objective of the study is to investigate complex numerical techniques that contributes to highly-accurate and stable solutions to many hyperbolic PDEs. The report starts by introducing the governing problem and studies its well-posedness to ensure stable and unique solutions of the governing equations. It continues with basic introduction to the classic spatial discretization of the FDM and introduces the SBP-SAT implementation of the method. The governing equations are rewritten as a semi-discrete problem, such that it can be written as a system of ordinary differential equations (ODEs) only dependent on the temporal evolution. This system can be solved with classic Runge-Kutta methods to ensure robust and accurate time-stepping schemes. The results show that the implementation of the higher-order SBP-SAT Finite Difference Method provides highly accurate solutions of the acoustic wave equation compared to the classic FDM. The results also show that the method provides stable solutions with no visible oscillations (dispersion), which can be a challenge for higher order methods. Overall, this paper contributes with valuable insights into the analysis of accuracy and stability in finite difference methods for acoustic wave equation. Acoustic Wave Equation Stability Accuracy Finite Difference Methods Summation-by-Parts (SBP) Simultaneous Approximation Term (SAT) Mathematics Matematik
12	Efficient Algorithms for Future Aircraft Design: Contributions to Aerodynamic Shape Optimization Hicken, Jason 24 September 2009 (has links) Advances in numerical optimization have raised the possibility that efficient and novel aircraft configurations may be ``discovered'' by an algorithm. To begin exploring this possibility, a fast and robust set of tools for aerodynamic shape optimization is developed. Parameterization and mesh-movement are integrated to accommodate large changes in the geometry. This integrated approach uses a coarse B-spline control grid to represent the geometry and move the computational mesh; consequently, the mesh-movement algorithm is two to three orders faster than a node-based linear elasticity approach, without compromising mesh quality. Aerodynamic analysis is performed using a flow solver for the Euler equations. The governing equations are discretized using summation-by-parts finite-difference operators and simultaneous approximation terms, which permit nonsmooth mesh continuity at block interfaces. The discretization results in a set of nonlinear algebraic equations, which are solved using an efficient parallel Newton-Krylov-Schur strategy. A gradient-based optimization algorithm is adopted. The gradient is evaluated using adjoint variables for the flow and mesh equations in a sequential approach. The flow adjoint equations are solved using a novel variant of the Krylov solver GCROT. This variant of GCROT is flexible to take advantage of non-stationary preconditioners and is shown to outperform restarted flexible GMRES. The aerodynamic optimizer is applied to several studies of induced-drag minimization. An elliptical lift distribution is recovered by varying spanwise twist, thereby validating the algorithm. Planform optimization based on the Euler equations produces a nonelliptical lift distribution, in contrast with the predictions of lifting-line theory. A study of spanwise vertical shape optimization confirms that a winglet-up configuration is more efficient than a winglet-down configuration. A split-tip geometry is used to explore nonlinear wake-wing interactions: the optimized split-tip demonstrates a significant reduction in induced drag relative to a single-tip wing. Finally, the optimal spanwise loading for a box-wing configuration is investigated. induced drag shape optimization aircraft design computational fluid dynamics b-spline volume Newton Krylov GMRES GCROT adjoint summation-by-parts simultaneous approximation terms 0538
13	Coupled High-Order Finite Difference and Unstructured Finite Volume Methods for Earthquake Rupture Dynamics in Complex Geometries O'Reilly, Ossian January 2011 (has links) The linear elastodynamic two-dimensional anti-plane stress problem, where deformations occur in only one direction is considered for one sided non-planar faults. Fault dynamics are modeled using purely velocity dependent friction laws, and applied on boundaries with complex geometry. Summation-by-parts operators and energy estimates are used to couple a high-order finite difference method with an unstructured finite volume method. The unstructured finite volume method is used near the fault and the high-order finite difference method further away from the fault where no complex geometry is present. Boundary conditions are imposed weakly on characteristic form using the simultaneous approximation term technique, allowing explicit time integration to be used. Numerical computations are performed to verify the accuracy and time stability, of the method.
14	Efficient Algorithms for Future Aircraft Design: Contributions to Aerodynamic Shape Optimization Hicken, Jason 24 September 2009 (has links) Advances in numerical optimization have raised the possibility that efficient and novel aircraft configurations may be ``discovered'' by an algorithm. To begin exploring this possibility, a fast and robust set of tools for aerodynamic shape optimization is developed. Parameterization and mesh-movement are integrated to accommodate large changes in the geometry. This integrated approach uses a coarse B-spline control grid to represent the geometry and move the computational mesh; consequently, the mesh-movement algorithm is two to three orders faster than a node-based linear elasticity approach, without compromising mesh quality. Aerodynamic analysis is performed using a flow solver for the Euler equations. The governing equations are discretized using summation-by-parts finite-difference operators and simultaneous approximation terms, which permit nonsmooth mesh continuity at block interfaces. The discretization results in a set of nonlinear algebraic equations, which are solved using an efficient parallel Newton-Krylov-Schur strategy. A gradient-based optimization algorithm is adopted. The gradient is evaluated using adjoint variables for the flow and mesh equations in a sequential approach. The flow adjoint equations are solved using a novel variant of the Krylov solver GCROT. This variant of GCROT is flexible to take advantage of non-stationary preconditioners and is shown to outperform restarted flexible GMRES. The aerodynamic optimizer is applied to several studies of induced-drag minimization. An elliptical lift distribution is recovered by varying spanwise twist, thereby validating the algorithm. Planform optimization based on the Euler equations produces a nonelliptical lift distribution, in contrast with the predictions of lifting-line theory. A study of spanwise vertical shape optimization confirms that a winglet-up configuration is more efficient than a winglet-down configuration. A split-tip geometry is used to explore nonlinear wake-wing interactions: the optimized split-tip demonstrates a significant reduction in induced drag relative to a single-tip wing. Finally, the optimal spanwise loading for a box-wing configuration is investigated. induced drag shape optimization aircraft design computational fluid dynamics b-spline volume Newton Krylov GMRES GCROT adjoint summation-by-parts simultaneous approximation terms 0538
15	Weak Boundary and Interface Procedures for Wave and Flow Problems Abbas, Qaisar January 2011 (has links) In this thesis, we have analyzed the accuracy and stability aspects of weak boundary and interface conditions (WBCs) for high order finite difference methods on Summations-By-Parts (SBP) form. The numerical technique has been applied to wave propagation and flow problems. The advantage of WBCs over strong boundary conditions is that stability of the numerical scheme can be proven. The boundary procedures in the advection-diffusion equation for a boundary layer problem is analyzed. By performing Navier-Stokes calculations, it is shown that most of the conclusions from the model problem carries over to the fully nonlinear case. The work was complemented to include the new idea of using WBCs on multiple grid points in a region, where the data is known, instead of at a single point. It was shown that we can achieve high accuracy, an increased rate of convergence to steady-state and non-reflecting boundary conditions by using this approach. Using the SBP technique and WBCs, we have worked out how to construct conservative and energy stable hybrid schemes for shocks using two different approaches. In the first method, we combine a high order finite difference scheme with a second order MUSCL scheme. In the second method, a procedure to locally change the order of accuracy of the finite difference schemes is developed. The main purpose is to obtain a higher order accurate scheme in smooth regions and a low order non-oscillatory scheme in the vicinity of shocks. Furthermore, we have analyzed the energy stability of the MUSCL scheme, by reformulating the scheme in the framework of SBP and artificial dissipation operators. It was found that many of the standard slope limiters in the MUSCL scheme do not lead to a negative semi-definite dissipation matrix, as required to get pointwise stability. Finally, high order simulations of shock diffracting over a convex wall with two facets were performed. The numerical study is done for a range of Reynolds numbers. By monitoring the velocities at the solid wall, it was shown that the computations were resolved in the boundary layer. Schlieren images from the computational results were obtained which displayed new interesting flow features. weak boundary conditions multiple penalty finite difference methods summation-by-parts high order scheme hybrid methods MUSCL scheme shocks stability energy estimate steady-state non-reflecting
16	Numerical Simulation of Soliton Tunneling Tiberg, Matilda, Estensen, Elias, Seger, Amanda January 2020 (has links) This project studied two different ways of imposing boundary conditions weakly with the finite difference summation-by-parts (SBP) operators. These operators were combined with the boundary handling methods of simultaneous-approximation-terms (SAT) and the Projection to impose homogeneous Neumann and Dirichlet boundary conditions. The convergence rate of both methods was analyzed for different boundary conditions for the one-dimensional (1D) Schrödinger equation, without potential, which resulted in both methods performing similarly. A multi-block discretization was then implemented and different combinations of SBP-SAT and SBP-Projection were applied to impose inner boundary conditions of continuity between the blocks. A convergence study of the different methods of imposing the inner BC:s was conducted for the 1D Schrödinger equation without potential. The resulting convergence was the same for all methods and it was concluded that they performed similarly. Methods involving SBP-Projection had the slight advantage of faster computation time. Finally, the 1D Gross-Pitaevskii equation (GPE) and the 1D Schrödinger equation were analyzed with a step potential. The waves propagating towards the potential barrier were in both cases partially transmitted and partially reflected. The waves simulated with the Schrödinger equation dispersed, while the solitons simulated with the GPE kept their shape due to the equations reinforcing non-linear term. The bright soliton was partly transmitted and partly reflected. The dark soliton was either totally reflected or totally transmitted. soliton GPE gross pitaevskii Schrödinger step potential summation by parts SBP-SAT Simultaneous approximation term SAT Projection SBP-Projection finite differences boundary treatment interface treatment Computational Mathematics Beräkningsmatematik
17	Numerical Simulation of the Generalized Modified Benjamin-Bona-Mahony Equation Using SBP-SAT in Time Kjelldahl, Vilma January 2023 (has links) This paper describes simulations of the generalized modified Benjamin-Bona-Mahony (BBM) equation, using finite difference methods (FDM). Well-posed boundary conditions (BCs) as well as stable semi-discrete approximations are derived using summations-by-parts (SBP) operators combined with the projection method. For time integration, explicit Runge-Kutta 4 (RK4) is used, as well as SBP-SAT, which discretizes the temporal domain using SBP operators and imposes initial conditions using simultaneous approximation term (SAT). These time-marching methods are evaluated and compared in terms of accuracy and computing times, and soliton-boundary interaction is studied. It is shown that SBP-SAT time-marching perform well and is more suitable than RK4 for this type of non-linear, dispersive problem. Generalized summation-by-parts (GSBP) time-marching perform particularly well, due to high accuracy with few solution points. Benjamin-Bona-Mahony equation non-linear initial boundary value problem summation-by-parts simultaneous-approximation-term time integration simulation Computational Mathematics Beräkningsmatematik
18	Numerical simulation of acoustic wave propagation with a focus on modeling sediment layers and large domains Estensen, Elias January 2022 (has links) In this report, we study how finite differences can be used to simulate acoustic wave propagation originating from a point source in the ocean using the Helmholtz equation. How to model sediment layers and the vast size of the ocean is studied in particular. The finite differences are implemented with summation by parts operators with boundary conditions enforced with simultaneous approximation terms and projection. The numerical solver is combined with the WaveHoltz method to improve the performance. Sediment layers are handled with interface conditions and the domain is artificially expanded using absorbing layers. The absorbing layer is implemented with an alternative approach to the super-grid method where the domain expansion is accomplished by altering the wave speed rather than with coordinate transformations. To isolate these issues, other parameters such as variations in the ocean floor are neglected. With this simplification, cylindrical coordinates are used and the angular variation is assumed to be zero. This reduces the problem to a quasi-three-dimensional system. We study how the parameters of the alternative absorbing layer approach affect its quality. The numerical solver is verified on several test cases and appears to work according to theory. Finally, a semi-realistic simulation is carried out and the solution seems correct in this setting. Helmholtz wave equation acoustics ocean acoustics numerical analysis computational science finite differences summation by parts sbp waveholtz cylindrical coordinates sediment layers interface super-grid super grid absorbing layer Helmholtz vågekvationen akustik undervattensakustik numerisk analys beräkningsvetenskap finita differenser summation by parts sbp cylindriska koordinater waveholtz sedimentskikt interface super-grid absorberande lager Computational Mathematics Beräkningsmatematik
19	High order summation-by-parts methods in time and space Lundquist, Tomas January 2016 (has links) This thesis develops the methodology for solving initial boundary value problems with the use of summation-by-parts discretizations. The combination of high orders of accuracy and a systematic approach to construct provably stable boundary and interface procedures makes this methodology especially suitable for scientific computations with high demands on efficiency and robustness. Most classes of high order methods can be applied in a way that satisfies a summation-by-parts rule. These include, but are not limited to, finite difference, spectral and nodal discontinuous Galerkin methods. In the first part of this thesis, the summation-by-parts methodology is extended to the time domain, enabling fully discrete formulations with superior stability properties. The resulting time discretization technique is closely related to fully implicit Runge-Kutta methods, and may alternatively be formulated as either a global method or as a family of multi-stage methods. Both first and second order derivatives in time are considered. In the latter case also including mixed initial and boundary conditions (i.e. conditions involving derivatives in both space and time). The second part of the thesis deals with summation-by-parts discretizations on multi-block and hybrid meshes. A new formulation of general multi-block couplings in several dimensions is presented and analyzed. It collects all multi-block, multi-element and hybrid summation-by-parts schemes into a single compact framework. The new framework includes a generalized description of non-conforming interfaces based on so called summation-by-parts preserving interpolation operators, for which a new theoretical accuracy result is presented. summation-by-parts time integration stiff problems weak initial conditions high order methods simultaneous-approximation-term finite difference discontinuous Galerkin spectral methods conservation energy stability complex geometries non-conforming grid interfaces interpolation
20	Summation By Parts Finite Difference Methods with Simultaneous Approximation Terms for the Heat Equation with Discontinuous Coefficients Kåhlman, Niklas January 2019 (has links) In this thesis we will investigate how the SBP-SAT finite difference method behave with and without an interface. As model problem, we consider the heat equation with piecewise constant coefficients. The thesis is split in two main parts. In the first part we look at the heat equation in one-dimension, and in the second part we expand the problem to a two-dimensional domain. We show how the SAT-parameters are chosen such that the scheme is dual consistent and stable. Then, we perform numerical experiments, now looking at the static case. In the one-dimensional case we see that the second order SBP-SAT method with an interface converge with an order of two, while the second order SBP-SAT method without an interface converge with an order of one. Finite Difference Method FDM Summation By Parts Simultaneous Approximation Terms SBP-SAT Heat Equation Poisson's Equation Discontinuous Coefficients Piecewise Constant Coefficients Convergence Dual Consistency Stability Computational Mathematics Beräkningsmatematik Discrete Mathematics Diskret matematik Mathematical Analysis Matematisk analys Other Mathematics Annan matematik

Search results