scholarly journals A Newton-Krylov solver for fully implicit 3D extended MHD

2007 ◽  
Author(s):  
Luis Chacon
Keyword(s):  
2010 ◽  
Vol 229 (21) ◽  
pp. 8130-8143 ◽  
Author(s):  
Hinrich Lütjens ◽  
Jean-François Luciani
Keyword(s):  

Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1639
Author(s):  
Abdelkrim Aharmouch ◽  
Brahim Amaziane ◽  
Mustapha El Ossmani ◽  
Khadija Talali

We present a numerical framework for efficiently simulating seawater flow in coastal aquifers using a finite volume method. The mathematical model consists of coupled and nonlinear partial differential equations. Difficulties arise from the nonlinear structure of the system and the complexity of natural fields, which results in complex aquifer geometries and heterogeneity in the hydraulic parameters. When numerically solving such a model, due to the mentioned feature, attempts to explicitly perform the time integration result in an excessively restricted stability condition on time step. An implicit method, which calculates the flow dynamics at each time step, is needed to overcome the stability problem of the time integration and mass conservation. A fully implicit finite volume scheme is developed to discretize the coupled system that allows the use of much longer time steps than explicit schemes. We have developed and implemented this scheme in a new module in the context of the open source platform DuMu X . The accuracy and effectiveness of this new module are demonstrated through numerical investigation for simulating the displacement of the sharp interface between saltwater and freshwater in groundwater flow. Lastly, numerical results of a realistic test case are presented to prove the efficiency and the performance of the method.


2021 ◽  
Vol 11 (12) ◽  
pp. 5509
Author(s):  
Hongjin Choi ◽  
Seonghwan Choi ◽  
Soo-Chang Kang ◽  
Myoung-Gyu Lee

A fully implicit stress integration algorithm is developed for the distortional hardening model, namely the e−HAH model, capable of simulating cross−hardening/softening under orthogonal loading path changes. The implicit algorithm solves a complete set of residuals as nonlinear functions of stress, a microstructure deviator, and plastic state variables of the constitutive model, and provides a consistent tangent modulus. The number of residuals is set to be 20 or 14 for the continuum or shell elements, respectively. Comprehensive comparison programs are presented regarding the predictive accuracy and stability with different numerical algorithms, strain increments, material properties, and loading conditions. The flow stress and r−value evolutions under reverse/cross−loading conditions prove that the algorithm is robust and accurate, even with large strain increments. By contrast, the cutting−plane method and partially implicit Euler backward method, which are characterized by a reduced number of residuals, result in unstable responses under abrupt loading path changes. Finally, the algorithm is implemented into the finite element modeling of large−size, S−rail forming and the springback for two automotive steel sheets, which is often solved by a hybrid dynamic explicit–implicit scheme. The fully implicit algorithm performs well for the whole simulation with the solely static implicit scheme.


2021 ◽  
Vol 33 (7) ◽  
pp. 071909
Author(s):  
Pei-Chuan Chao ◽  
Ali Gürbüz ◽  
Frederick Sachs ◽  
M. V. Sivaselvan

Sign in / Sign up

Export Citation Format

Share Document