scholarly journals Efficient Numerical Treatment of Ambipolar and Hall Drift as Hyperbolic System

2021 ◽  
Vol 923 (1) ◽  
pp. 79
Author(s):  
M. Rempel ◽  
D. Przybylski
Keyword(s):  
Time Derivative ◽  
Transport Processes ◽  
Solar Chromosphere ◽  
Explicit Integration ◽  
Time Step ◽  
Time Stepping ◽  
Computational Overhead ◽  
Partially Ionized Plasmas ◽  
Ambipolar Drift ◽  
Diffusive Time

Abstract Partially ionized plasmas, such as the solar chromosphere, require a generalized Ohm’s law including the effects of ambipolar and Hall drift. While both describe transport processes that arise from the multifluid equations and are therefore of hyperbolic nature, they are often incorporated in models as a diffusive, i.e., parabolic process. While the formulation as such is easy to include in standard MHD models, the resulting diffusive time-step constraints do require often a computationally more expensive implicit treatment or super-time-stepping approaches. In this paper we discuss an implementation that retains the hyperbolic nature and allows for an explicit integration with small computational overhead. In the case of ambipolar drift, this formulation arises naturally by simply retaining a time derivative of the drift velocity that is typically omitted. This alone leads to time-step constraints that are comparable to the native MHD time-step constraint for a solar setup including the region from photosphere to lower solar corona. We discuss an accelerated treatment that can further reduce time-step constraints if necessary. In the case of Hall drift we propose a hyperbolic formulation that is numerically similar to that for the ambipolar drift and we show that the combination of both can be applied to simulations of the solar chromosphere at minimal computational expense.

Fast Conjugate Heat Transfer Simulation of Long Transient Flexible Operations Using Adaptive Time Stepping

2018 ◽  
Vol 140 (9) ◽  
Author(s):  
R. Maffulli ◽  
L. He ◽  
P. Stein ◽  
G. Marinescu
Keyword(s):  
Heat Transfer ◽  
Conjugate Heat Transfer ◽  
Time Derivative ◽  
Computational Cost ◽  
Strong Impact ◽  
Time Step ◽  
Time Stepping ◽  
Adaptive Time ◽  
Fully Coupled ◽  
Adaptive Time Stepping

The emerging renewable energy market calls for more advanced prediction tools for turbine transient operations in fast startup/shutdown cycles. Reliable numerical analysis of such transient cycles is complicated by the disparity in time scales of the thermal responses in fluid and solid domains. Obtaining fully coupled time-accurate unsteady conjugate heat transfer (CHT) results under these conditions would require to march in both domains using the time-step dictated by the fluid domain: typically, several orders of magnitude smaller than the one required by the solid. This requirement has strong impact on the computational cost of the simulation as well as being potentially detrimental to the accuracy of the solution due to accumulation of round-off errors in the solid. A novel loosely coupled CHT methodology has been recently proposed, and successfully applied to both natural and forced convection cases that remove these requirements through a source-term based modeling (STM) approach of the physical time derivative terms in the relevant equations. The method has been shown to be numerically stable for very large time steps with adequate accuracy. The present effort is aimed at further exploiting the potential of the methodology through a new adaptive time stepping approach. The proposed method allows for automatic time-step adjustment based on estimating the magnitude of the truncation error of the time discretization. The developed automatic time stepping strategy is applied to natural convection cases under long (2000 s) transients: relevant to the prediction of turbine thermal loads during fast startups/shutdowns. The results of the method are compared with fully coupled unsteady simulations showing comparable accuracy with a significant reduction of the computational costs.

Variable Timestepping in EHD Problems

2005 ◽  
Author(s):  
Vasilios Bakolas ◽  
Wolfgang Borchers
Keyword(s):  
Control Method ◽  
Time Derivative ◽  
Second Order ◽  
Operating Conditions ◽  
Time Interval ◽  
Time Step ◽  
Simulation Code ◽  
Automotive Engine ◽  
Time Stepping

Transient EHD calculations have moved to the center of the attention in the past years as more machine elements that operate under variable operating conditions, such as automotive engine components, are examined. Almost all of the transient calculations today are carried out with a constant time step using first or second order discretisation for the determination of the time derivatives. Recently, a method was proposed that could control the time stepping of a transient calculation using first order time derivatives. In this paper, a method is presented for the control of the time step in transient calculations. The method employs a second order backward differentiation scheme with variable time stepping for the determination of the time derivatives. The time interval between successive calculations is determined depending on an estimate of the local truncation error without solving additional large linear systems, thus allowing for larger time steps when the time derivative is almost constant. This time stepping control method was implemented in an EHD simulation code that employed multi-grid and multi-level multi-summation techniques. Results presented in this paper prove that the proposed method can significantly reduce the number of necessary time steps of a transient calculation without any loss of accuracy of the results.

Fast Conjugate Heat Transfer Simulation of Long Transient Flexible Operations Using Adaptive Time Stepping

2018 ◽  
Author(s):  
R. Maffulli ◽  
L. He ◽  
P. Stein ◽  
G. Marinescu
Keyword(s):  
Heat Transfer ◽  
Conjugate Heat Transfer ◽  
Time Derivative ◽  
Computational Cost ◽  
Strong Impact ◽  
Time Step ◽  
Time Stepping ◽  
Adaptive Time ◽  
Fully Coupled ◽  
Adaptive Time Stepping

The emerging renewable energy market calls for more advanced prediction tools for turbine transient operations in fast startup/shutdown cycles. Reliable numerical analysis of such transient cycles is complicated by the disparity in time scales of the thermal responses in fluid and solid domains. Obtaining fully-coupled time-accurate unsteady Conjugate Heat Transfer (CHT) results under these conditions would require to march in both domains using the time step dictated by the fluid domain: typically, several orders of magnitude smaller than the one required by the solid. This requirement has strong impact on the computational cost of the simulation as well as being potentially detrimental to the accuracy of the solution due to accumulation of round-off errors in the solid. A novel loosely coupled CHT methodology has been recently proposed, and successfully applied to both natural and forced convection cases that removes these requirements through a source-term based modelling approach of the physical time derivative terms in the relevant equations. The method has been shown to be numerically stable for very large time steps with adequate accuracy. The present effort is aimed at further exploiting the potential of the methodology through a new adaptive time stepping approach. The proposed method allows for automatic time step adjustment based on estimating the magnitude of the truncation error of the time discretization. The developed automatic time stepping strategy is applied to natural convection cases under long (2000 s) transients: relevant to the prediction of turbine thermal loads during fast startups/shutdowns. The results of the method are compared with fully coupled unsteady simulations showing comparable accuracy with a significant reduction of the computational costs.

Using the rapid expansion method for accurate time-stepping in modeling and reverse-time migration

Geophysics ◽  
2011 ◽  
Vol 76 (4) ◽  
pp. S177-S185 ◽  
Author(s):  
Ekkehart Tessmer
Keyword(s):  
Time Derivative ◽  
Expansion Method ◽  
Rapid Expansion ◽  
Numerical Dispersion ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Step ◽  
Time Stepping ◽  
Time Migration ◽  
Spatial Derivatives

Reverse-time migration is based on seismic forward modeling algorithms, where spatial derivatives usually are calculated by finite differences or by the Fourier method. Time integration in general is done by finite-difference time stepping of low orders. If the spatial derivatives are calculated by high-order methods and time stepping is based on low-order methods, there is an imbalance that might require that the time-step size needs to be very small to avoid numerical dispersion. As a result, computing times increase. Using the rapid expansion method (REM) avoids numerical dispersion if the number of expansion terms is chosen properly. Comparisons with analytical solutions show that the REM is preferable, especially at larger propagation times. For reverse-time migration, the REM needs to be applied in a time-stepping manner. This is necessary because the original implementation based on very large time spans requires that the source term is separable in space and time. This is not appropriate for reverse-time migration where the sources have different time histories. In reverse-time migration, it might be desirable to use the Poynting vector information to estimate opening angles to improve the quality of the image. In the solution of the wave equation, this requires that one calculates not only the pressure wavefield but also its time derivative. The rapid expansion method can be extended easily to provide this time derivative with negligible extra cost.

scholarly journals Analysis of a Decoupled Time-Stepping Scheme for Evolutionary Micropolar Fluid Flows

2016 ◽  
Vol 2016 ◽  
pp. 1-13 ◽  
Author(s):  
S. S. Ravindran
Keyword(s):  
Micropolar Fluid ◽  
Stokes Equations ◽  
Fluid Model ◽  
Navier Stokes ◽  
Optimal Order ◽  
Step Size ◽  
Time Step ◽  
Order Error ◽  
Time Step Size ◽  
Time Stepping

Micropolar fluid model consists of Navier-Stokes equations and microrotational velocity equations describing the dynamics of flows in which microstructure of fluid is important. In this paper, we propose and analyze a decoupled time-stepping algorithm for the evolutionary micropolar flow. The proposed method requires solving only one uncoupled Navier-Stokes and one microrotation subphysics problem per time step. We derive optimal order error estimates in suitable norms without assuming any stability condition or time step size restriction.

A Sequential Integration Method

1988 ◽  
Vol 110 (4) ◽  
pp. 382-388
Author(s):  
Liang-Wey Chang ◽  
James F. Hamilton
Keyword(s):  
Integration Method ◽  
Equations Of Motion ◽  
Slow Motion ◽  
Explicit Integration ◽  
Time Step ◽  
Lumped Model ◽  
Guaranteed Stability ◽  
Fast Motion ◽  
The Stability ◽  
Secular Terms

This paper presents a method for simulating systems with two inertially coupled motions, i.e., a slow motion and a fast motion. The equations of motion are separated into two sets of coupled nonlinear ordinary differential equations. For each time step, the two sets of equations are integrated sequentially rather than simultaneously. Explicit integration methods are used for integrating the slow motion since the stability of the integration is not a problem and the explicit methods are very convenient for nonlinear equations. For the fast motion, the equations are linear and the implicit integrations can be used with guaranteed stability. The size of time step only needs to be chosen to provide accuracy of the solution for the modes that are excited. The interaction between the two types of motion must be treated such that secular terms do not appear due to the sequential integration method. A lumped model of a flexible pendulum will be presented in this paper to illustrate the application of the method. Numerical results for both simultaneous and sequential integration are presented for comparison.

Semi-Elimination Methodology for Simulating High Flow Features in a Reservoir

2021 ◽  
Author(s):  
Yahan Yang ◽  
Ali Samii ◽  
Zhenlong Zhao ◽  
Guotong Ren
Keyword(s):  
High Flow ◽  
Small Cells ◽  
Computational Techniques ◽  
Computational Technique ◽  
System Of Equations ◽  
Time Step ◽  
Linear System Of Equations ◽  
Time Stepping ◽  
Model Complex ◽  
Upstream Weighting

Abstract Despite the rapid rise of computing power and advances in computational techniques in past decades, it is still challenging in reservoir simulation to model complex and detailed features that are represented by small cells with large permeability values, for example, fractures, multi-segment wells, etc. While those features may carry a large amount of flow and thus have a significant impact on the performance prediction, the combination of small volume and large permeability unfortunately leads to well-known time stepping and convergence difficulties during Newton iteration. We address this issue of high flow through small cells by developing a new semi-elimination computational technique. At the beginning of simulation, we construct a set of pressure basis which is a mapping from pressures at surrounding cells in the bulk of reservoir to pressures at those small cells. Next, we start the time-stepping scheme. For each time step or iteration within a time step, small cells are first employed to provide an accurate computation of flow rates and derivatives using upstream weighting and a flow partitioning scheme. Afterwards, small cells are eliminated and a linear system of equations is assembled and solved involving only bulk cells. This semi-elimination technique allows us to fundamentally avoid the drawbacks caused by including small cells in the global system of equations, while capturing their effect on the flow of hydrocarbon in the reservoir. One of the advantages of the proposed techniques over other existing methods is that it is fully implicit and preserves upstream weighting and compositions of the flow field even after small cells are eliminated, which enhances numerical stability and accuracy of simulation results. Application of this technique to several synthetic and field models demonstrates significant performance and accuracy improvement over standard approaches. This method thus offers a practical way to model complex and dynamic flow behaviors in important features without incurring penalties in speed and robustness of the simulation.

scholarly journals A general-purpose time-step criterion for simulations with gravity

2020 ◽  
Vol 495 (4) ◽  
pp. 4306-4313 ◽  
Author(s):  
Michael Y Grudić ◽  
Philip F Hopkins
Keyword(s):  
Time Scale ◽  
General Purpose ◽  
Gravitational Acceleration ◽  
Time Step ◽  
Simulation Code ◽  
Time Stepping ◽  
Dynamical Time ◽  
Adaptive Time Step ◽  
Order Of Magnitude ◽  
True Time

ABSTRACT We describe a new adaptive time-step criterion for integrating gravitational motion, which uses the tidal tensor to estimate the local dynamical time-scale and scales the time-step proportionally. This provides a better candidate for a truly general-purpose gravitational time-step criterion than the usual prescription derived from the gravitational acceleration, which does not respect the equivalence principle, breaks down when $\boldsymbol {a}=0$, and does not obey the same dimensional scaling as the true time-scale of orbital motion. We implement the tidal time-step criterion in the simulation code gizmo, and examine controlled tests of collisionless galaxy and star cluster models, as well as galaxy merger simulations. The tidal criterion estimates the dynamical time faithfully, and generally provides a more efficient time-stepping scheme compared to an acceleration criterion. Specifically, the tidal criterion achieves order-of-magnitude smaller energy errors for the same number of force evaluations in potentials with inner profiles shallower than ρ ∝ r−1 (i.e. where $\boldsymbol {a}\rightarrow 0$), such as star clusters and cored galaxies. For a given problem these advantages must be weighed against the additional overhead of computing the tidal tensor on-the-fly, but in many cases this overhead is small.

Numerical simulation of higher-order diffusion-wave equations of variable coefficients using the meshless spectral method

2020 ◽  
pp. 2150010
Author(s):  
Manzoor Hussain ◽  
Sirajul Haq
Keyword(s):  
Wave Equations ◽  
Interpolation Method ◽  
Variable Coefficients ◽  
Higher Order ◽  
Implicit Time ◽  
Test Problems ◽  
Time Step ◽  
Diffusion Wave ◽  
Time Step Size ◽  
Time Stepping

In this paper, meshless spectral interpolation technique using implicit time stepping scheme is proposed for the numerical simulations of time-fractional higher-order diffusion wave equations (TFHODWEs) of variable coefficients. Meshless shape functions, obtained from radial basis functions (RBFs) and point interpolation method (PIM), are used for spatial approximation. Central differences coupled with quadrature rule of [Formula: see text] are employed for fractional temporal approximation. For advancement of solution, an implicit time stepping scheme is then invoked. Simulations performed for different benchmark test problems feature good agreement with exact solutions. Stability analysis of the proposed method is theoretically discussed and computationally validated to support the analysis. Accuracy and efficiency of the proposed method are assessed via [Formula: see text], [Formula: see text] and [Formula: see text] error norms as well as number of nodes [Formula: see text] and time step-size [Formula: see text].

A computational method for earthquake cycles within anisotropic media

2019 ◽  
Vol 219 (2) ◽  
pp. 816-833 ◽  
Author(s):  
Maricela Best Mckay ◽  
Brittany A Erickson ◽  
Jeremy E Kozdon
Keyword(s):  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Elastic Solid ◽  
Elliptic Partial Differential Equation ◽  
Computational Method ◽  
Parameter Study ◽  
Time Step ◽  
Time Stepping ◽  
Static Elasticity ◽  
Earthquake Cycles

SUMMARY We present a numerical method for the simulation of earthquake cycles on a 1-D fault interface embedded in a 2-D homogeneous, anisotropic elastic solid. The fault is governed by an experimentally motivated friction law known as rate-and-state friction which furnishes a set of ordinary differential equations which couple the interface to the surrounding volume. Time enters the problem through the evolution of the ordinary differential equations along the fault and provides boundary conditions for the volume, which is governed by quasi-static elasticity. We develop a time-stepping method which accounts for the interface/volume coupling and requires solving an elliptic partial differential equation for the volume response at each time step. The 2-D volume is discretized with a second-order accurate finite difference method satisfying the summation-by-parts property, with boundary and fault interface conditions enforced weakly. This framework leads to a provably stable semi-discretization. To mimic slow tectonic loading, the remote side-boundaries are displaced at a slow rate, which eventually leads to earthquake nucleation at the fault. Time stepping is based on an adaptive, fourth-order Runge–Kutta method and captures the highly varying timescales present. The method is verified with convergence tests for both the orthotropic and fully anisotropic cases. An initial parameter study reveals regions of parameter space where the systems experience a bifurcation from period one to period two behaviour. Additionally, we find that anisotropy influences the recurrence interval between earthquakes, as well as the emergence of aseismic transients and the nucleation zone size and depth of earthquakes.

Sign in / Sign up

Export Citation Format

Share Document