Search results for "finite difference"
showing 10 items of 122 documents
Combining spectral and shock-capturing methods: A new numerical approach for 3D relativistic core collapse simulations
2005
We present a new three-dimensional general relativistic hydrodynamics code which is intended for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. Contrary to the common approach followed in most existing three-dimensional numerical relativity codes which are based in Cartesian coordinates, in this code both the metric and the hydrodynamics equations are formulated and solved numerically using spherical polar coordinates. A distinctive feature of this new code is the combination of two types of accurate numerical schemes specifically designed to solve each system of equations. More precisely, the code uses spectra…
Efficient finite difference formulation of a geometrically nonlinear beam element
2021
The article is focused on a two-dimensional geometrically nonlinear formulation of a Bernoulli beam element that can accommodate arbitrarily large rotations of cross sections. The formulation is based on the integrated form of equilibrium equations, which are combined with the kinematic equations and generalized material equations, leading to a set of three first-order differential equations. These equations are then discretized by finite differences and the boundary value problem is converted into an initial value problem using a technique inspired by the shooting method. Accuracy of the numerical approximation is conveniently increased by refining the integration scheme on the element lev…
Introduction: Signals and Transforms
2015
In this chapter we outline some well known facts about periodic signals and transforms, which are needed throughout the book. For details we refer to the classical textbook Oppenheim and Schafer [2].
Coupled fluid-flow and magnetic-field simulation of the Riga dynamo experiment
2006
Magnetic fields of planets, stars, and galaxies result from self-excitation in moving electroconducting fluids, also known as the dynamo effect. This phenomenon was recently experimentally confirmed in the Riga dynamo experiment [ A. Gailitis et al., Phys. Rev. Lett. 84, 4365 (2000) ; A. Gailitis et al., Physics of Plasmas 11, 2838 (2004) ], consisting of a helical motion of sodium in a long pipe followed by a straight backflow in a surrounding annular passage, which provided adequate conditions for magnetic-field self-excitation. In this paper, a first attempt to simulate computationally the Riga experiment is reported. The velocity and turbulence fields are modeled by a finite-volume Navi…
Numerical Study of Forced MHD Convection Flow and Temperature Around Periodically Placed Cylinders
2016
In this paper we consider 2D stationary boundary value problems for the system of magnetohydrodynamic (MHD) equations and the heat transfer equation. The viscous electrically conducting incompressible liquid moves between infinite cylinders with square or round sections placed periodically. We also consider similar 2D MHD channel flow with periodically placed obstacles on the channel walls. We analyse the 2D forced and free MHD convection flow and temperature around cylinders and obstacles in homogeneous external magnetic field. The cylinders, obstacles and walls of the channel with constant temperature are heated. The distributions of electromagnetic fields, forces, velocity and temperatur…
Higher Order Sobolev-Type Spaces on the Real Line
2014
This paper gives a characterization of Sobolev functions on the real line by means of pointwise inequalities involving finite differences. This is also shown to apply to more general Orlicz-Sobolev, Lorentz-Sobolev, and Lorentz-Karamata-Sobolev spaces.
Adjoint-based inversion for porosity in shallow reservoirs using pseudo-transient solvers for non-linear hydro-mechanical processes
2020
Abstract Porous flow is of major importance in the shallow subsurface, since it directly impacts on reservoir-scale processes such as waste fluid sequestration or oil and gas exploration. Coupled and non-linear hydro-mechanical processes describe the motion of a low-viscous fluid interacting with a higher viscous porous rock matrix. This two-phase flow may trigger the initiation of solitary waves of porosity, further developing into vertical high-porosity pipes or chimneys. These preferred fluid escape features may lead to localised and fast vertical flow pathways potentially problematic in the case of for instance CO2 sequestration. Constraining the porosity and the non-linearly related pe…
An Efficient Numerical Method for Time Domain Computational Electromagnetic Simulation
2018
In this paper an efficient numerical method in approximating the electric and magnetic fields is provided. The method is based on an implicit leapfrog arrangement in time and without mesh in space. Moreover, a projection scheme is introduced in order to improve the accuracy of the proposed approach and applied into the computational electromagnetic (CEM) framework. The PDEs governing the process are solved and some numerical results are reported to validate the numerical process.
Chronoamperometry of prussian blue films on ITO electrodes: Ohmic drop and film thickness effect
1999
Abstract The chronoamperograms associated with the reduction of prussian blue films deposited onto indium tin oxide (ITO) electrodes to the Everitt’s salt form, are influenced by the ohmic drop effect. These chronoamperometric curves have been simulated by means of a numerical finite difference model which is able to explain their shape and their dependence on the thickness of the film and on the uncompensated resistance. An analytical expression which describes the dependence of current against time at initial times considering the ohmic drop effect has also been proved when applied to these chronoamperometric curves at short times.
A numerical approach to the voltammograms of the reduction of Prussian Blue films on ITO electrodes
1997
The uncompensated resistance, mainly due to the ITO electrode, modifies the shape of voltammetric curves of the system Prussian Blue ⇄ Everitt's Salt films deposited on this transparent electrode. A numerical finite difference model which is able to explain the shape of these voltammetric curves is studied in this paper. This model explains the dependence of voltammetric curves on the film thickness and uncompensated resistance.