A stabilized finite element method for the Stokes–Temperature coupled problem

Araya, Rodolfo; Cárcamo, Cristian; Poza, Abner H.

doi:10.1016/j.apnum.2023.02.002

Cited by 5 publications

(4 citation statements)

References 42 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Finally, we solve equations ( 8), ( 9) and ( 11) using a stabilized finite element 45 which allows us to use equal order of polynomial interpolation at each unknown point.…”

Section: Methods 1 Fault Friction and Mantle Viscositymentioning

confidence: 99%

“…To obtain the pressuretemperature path along the fault (Extended Data Fig. 1b), we build a finite-element thermal model based on near-lithostatic pore-fluid pressure at megathrust depth 45 , where the temperature satisfies the steady-state heat advection-diffusion equation (Methods: Section 2). The resulting temperature distribution (Extended Data Fig.…”

Section: Modeling the Subduction Seismic Cycle Of The Southern Andesmentioning

confidence: 99%

See 1 more Smart Citation

Recurrence time and size of Chilean earthquakes influenced by geological structure

Julve,

Barbot,

Moreno

et al. 2023

Nat. Geosci.

Self Cite

View full text Add to dashboard Cite

In 1960, the giant Valdivia earthquake (Mw 9.5), the largest earthquake ever recorded, struck the Chilean subduction zone, rupturing the entire depth of the seismogenic zone and extending for 1,000 km along strike. The first sign of new seismic energy release since 1960 occurred in 2017 with the Melinka earthquake (Mw 7.6), which affected only a portion of the deepest part of the seismogenic zone. Despite the recognition that rupture characteristics and rheology vary with depth, the mechanical controls behind such variations of earthquake size remain elusive. Here, we build quasi-dynamic simulations of the seismic cycle in Southern Chile including frictional and viscoelastic properties, drawing upon a compilation of geological and geophysical insights, to explain the recurrence times of recent, historic, and paleoseismic earthquakes and the distribution of fault slip and crustal deformation associated with the Melinka and Valdivia earthquakes. The frictional and rheological properties of the forearc, primarily controlled by the geological structure and distribution of fluids at the megathrust, govern earthquake size and recurrence patterns in Chile. Article.Subduction zones produce 90% of the world's natural seismicity and host Earth's largest earthquakes accompanied by devastating tsunamis, outlying the key importance to understand subduction dynamics during the seismic cycle. This process is controlled by the structure of the accretionary prism and the fabric of the megathrust -the seismogenic interface that separates the down-going oceanic lithosphere from the upper plate. Earthquakes originate from a frictional instability at the plate interface, primarily affected by the lithology, temperature, and fluid content of the fault zone, leading to a stratification of source properties 1-3 . Numerous studies have shown variability in rupture width and slip of earthquakes in different depth domains, including great ruptures (moment magnitude Mw > 8.5) with a trench-breaking feature 4 , followed by blind and large ruptures (7 < Mw < 8.5) nucleated within the interseismically locked region as illustrated by paleoseismic and instrumental record at the Japan 5 , Sunda 6 , and Aleutian 7 trenches. The succession of large and great ruptures at the subduction megathrust modulates the size and frequency of earthquakes, but the underlying mechanics is uncertain.

show abstract

“…Finally, we solve equations ( 8), ( 9) and ( 11) using a stabilized finite element 45 which allows us to use equal order of polynomial interpolation at each unknown point.…”

Section: Methods 1 Fault Friction and Mantle Viscositymentioning

confidence: 99%

Section: Modeling the Subduction Seismic Cycle Of The Southern Andesmentioning

confidence: 99%

Recurrence time and size of Chilean earthquakes influenced by geological structure

Julve,

Barbot,

Moreno

et al. 2023

Nat. Geosci.

Self Cite

View full text Add to dashboard Cite

show abstract

“…Finally, we solve equations (8), ( 9) and ( 11) using a stabilized nite element 84 which allows us to use equal order of polynomial interpolation at each unknown point.…”

Section: Temperature Modelmentioning

confidence: 99%

Size and timing of giant Chilean earthquake controlled by rock composition and forearc structure

Julve

Barbot

Moreno

et al. 2023

Preprint

Self Cite

View full text Add to dashboard Cite

In 1960, the Southern Chile subduction zone produced the Mw9.5 Valdivia mega-quake, the largest earthquake on record, breaking up the entire seismogenic zone for 1,000 km along the Chile Trench. Crustal deformation in the aftermath involved a combination of fault motion and viscoelastic flow. The region has recovered its seismic potential, as evidenced by the 2016 Mw7.6 Melinka earthquake that only partially affected the 1960 rupture zone below the continental shelf. The succession of such partial and full ruptures at subduction megathrusts controls seismic hazards but remains poorly understood. Here we show how the distribution of frictional and rheological properties within the forearc and the surrounding lithosphere, controls the size and timing of great and giant earthquakes at the Chile Trench. We build numerical simulations of the seismic cycle in Southern Chile using frictional and viscoelastic properties calibrated to geologic, geodetic, and geophysical information. The model explains the recurrence times of great and giant earthquakes from paleoseismic data, but also the fault slip distribution and crustal deformation associated with the Melinka and Valdivia earthquakes. Our synoptic model illuminates the structural and geological controls on seismicity, explaining the typical sequence of blind and trench-breaking ruptures at subduction zones.

show abstract

“…In conjunction with the zone method, a heat conduction model for billets affected by radiation and convection was constructed, adopting a suitable convective heat transfer empirical coefficient. This method significantly improves the uniform temperature rise of billet in a large reheating furnace [ 5 ]. Peck [ 6 ] presented a novel reduced-order transfer-function-based mathematical model for temperature uniformity of a gas powered industrial box furnace, demonstrating that fine adjustments to the controller can reduce temperature rise time while maintaining compliance with Temperature Uniformity Survey (TUS) standards.…”

Section: Introductionmentioning

confidence: 99%

Soft Sensor Modeling for 3D Transient Temperature Field of Large-Scale Aluminum Alloy Workpieces Based on Multi-Loss Consistency Optimization PINN

Shen

Chen

Wang

et al. 2023

Sensors

View full text Add to dashboard Cite

Uniform temperature distribution during quenching thermal treatment is crucial for achieving exceptional mechanical and physical properties of alloy materials. Accurate and rapid prediction of the 3D transient temperature field model of large-scale aluminum alloy workpieces is key to realizing effective thermal treatment. This paper establishes a 3D transient temperature field model of large aluminum alloy workpieces and proposes a multi-loss consistency optimization-based physics-informed neural network (MCO-PINN) to realize soft sensing of the 3D temperature field model. The method is based on a MLP structure and adopts Gaussian activation functions. A surrogate model of the partial differential equation (PDE) is first constructed, and the residuals of the PDE, initial and boundary conditions, and observed data are encoded into the loss functions of the network. By establishing a Gaussian probability distribution model of each loss function and combining it with maximum likelihood estimation, the weight consistency optimization method of each loss function is then proposed to further improve the approximation ability of the model. To optimize the training speed of the network, an adaptive initial-value-eigenvector coding clustering (AIV-ECC) algorithm is finally proposed, which quickly determines the parameters of the Gaussian activation function, reduces the dependence on the initial value and improves the generalization performance of the network. Simulation and industrial experiments demonstrate that the proposed MCO-PINN can solve the 3D transient temperature field model with high precision and high time efficiency based on sparse measurements.

show abstract

A stabilized finite element method for the Stokes–Temperature coupled problem

Cited by 5 publications

References 42 publications

Recurrence time and size of Chilean earthquakes influenced by geological structure

Recurrence time and size of Chilean earthquakes influenced by geological structure

Size and timing of giant Chilean earthquake controlled by rock composition and forearc structure

Soft Sensor Modeling for 3D Transient Temperature Field of Large-Scale Aluminum Alloy Workpieces Based on Multi-Loss Consistency Optimization PINN

Contact Info

Product

Resources

About