- Geological and Geochemical Analysis
- Advanced Numerical Methods in Computational Mathematics
- High-pressure geophysics and materials
- earthquake and tectonic studies
- Advanced Mathematical Modeling in Engineering
- Seismic Imaging and Inversion Techniques
- Hydraulic Fracturing and Reservoir Analysis
- Hydrocarbon exploration and reservoir analysis
- Electromagnetic Simulation and Numerical Methods
- Numerical methods in engineering
- Methane Hydrates and Related Phenomena
- Geochemistry and Geologic Mapping
- Parallel Computing and Optimization Techniques
- Lattice Boltzmann Simulation Studies
- Enhanced Oil Recovery Techniques
- Advanced Numerical Analysis Techniques
- Computational Geometry and Mesh Generation
- Distributed and Parallel Computing Systems
- Big Data and Business Intelligence
- COVID-19 diagnosis using AI
- Radioactive contamination and transfer
- Solidification and crystal growth phenomena
- Advanced Bandit Algorithms Research
- Neural Networks and Applications
- Groundwater and Watershed Analysis
University of Cambridge
2007-2023
Bridge University
2019
Nagoya University
1997-1998
Ifremer
1997
Observatoire Midi-Pyrénées
1997
Centre National de la Recherche Scientifique
1997
Université de Bretagne Occidentale
1997
We develop a method for generating degree-of-freedom maps arbitrary order finite element spaces any cell shape. The approach is based on the composition of permutations and transformations by sub-entity. Current approaches to problems typically rely consistent orientation entities that permits definition common local coordinate system shared edges faces. However, while mesh straightforward simplex cells operation, it not strictly operation quadrilateral in case hexahedral all meshes are...
The finite element method (FEM) (Ciarlet, 1978) is a widely used numerical for approximating the solution of partial differential equations (PDEs).Solving problem using FEM involves discretising and searching in dimensional space: these spaces are created by defining on each cell mesh.Following Ciarlet (1978), commonly defined triple (R, V, L), where:• R reference cell, example triangle with vertices at (0,0), (1,0) (0,1); • V polynomial space, span{1, x, y, x 2 , xy, y }; L basis dual space...
Longitudinal river profiles, where elevation of a bed is plotted as function distance along the bed, contain information about uplift rate. When region adjacent to reference level (e.g., sea level) uplifted, rapid change in gradient occurs near mouth. The erosional process causes this migrate upstream. Thus profile effectively ‘tape recording’ rate history, provided that can be adequately parameterized. Here, we use non‐linear equation relate shape profile, z ( x ), U t ). If erosion assumed...
We consider the onset of thermosolutal (double-diffusive) convection a binary fluid in horizontal porous layer subject to fixed temperatures and chemical equilibrium on bounding surfaces, case when solubility dissolved component depends temperature. use linear stability analysis investigate how dissolution or precipitation this affects selection an unstable wavenumber; we extend using Galerkin method predict structure initial bifurcation compare our analytical results with numerical...
Research Article| September 01, 2012 Disequilibrium dihedral angles in dolerite sills: A new proxy for cooling rate Marian B. Holness; Holness 1Department of Earth Sciences, University Cambridge, Downing Street, Cambridge CB2 3EQ, UK Search other works by this author on: GSW Google Scholar Chris Richardson; Richardson 2BP Institute, Madingley Road, CB3 0EZ, Rosalind T. Helz 3U.S. Geological Survey National Center, 12201 Sunrise Valley Drive, Reston, Virginia 20192, USA Author and Article...
In this paper we investigate salt transport during the evaporation and upflow of saline groundwater. We describe a model in which sharp evaporation–precipitation front separates regions soil saturated with an air–vapour mixture water. then consider two idealized problems. first equilibrium configurations freshwater system when depth layer is finite, obtaining results for location water induced by evaporation. develop solution propagating infinite gravitational stability salinity profile...
It is often assumed that the molten regions of planetary interiors can be modeled as a compacting viscous matrix permeated by less melt phase. In order to assess true behaviour, it important test many different aspects these models possible. calculations presented here, viscosity made decrease proportion increases. New results in two dimensions show if only gravity drives flow, pockets travel upwards solitary waves, and behaviour similar well known constant case. However, when an external...
Because of the proximity Euler poles rotation Pacific and Antarctic plates, small variations in plate kinematics are fully recorded axial morphology geometry Pacific-Antarctic Ridge south Udintsev fracture zone. Swath bathymetry magnetic data show that clockwise rotations relative motion between plates over last 6 million years resulted rift propagation or linkage ridge segments, with transitions from transform faults to giant overlapping spreading centers. This bimodal rearrangement has...
Mixed dimensional partial differential equations (PDEs) are coupling unknown fields defined over domains of differing topological dimension. Such naturally arise in a wide range scientific including geology, physiology, biology, and fracture mechanics. PDEs also commonly encountered when imposing non-standard conditions subspace lower dimension, e.g., through Lagrange multiplier. In this article, we present general abstractions algorithms for finite element discretizations mixed domain...
Abstract Contact metamorphism caused by the Glenmore plug in Ardnamurchan, a magma conduit active for 1 month, resulted partial melting, with melt now preserved as glass. The pristine nature of much aureole provides natural laboratory which to investigate distribution melt. A simple thermal model, based on first appearance quartz–feldspar grain boundaries, quartz paramorphs after tridymite and plausible intrusion temperature, time‐scale melting. onset melting boundaries was initially rapid,...
This paper introduces LEoPart, an add-on for the open source finite element software library FEniCS to seamlessly integrate Lagrangian particle functionality with (Eulerian) mesh-based (FE) approaches. LEoPart - which is so much as say: `Lagrangian-Eulerian on Particles' contains tools efficient, accurate and scalable advection of particles arbitrary polyhedral meshes. In addition, comes several projection operators exchanging information between scattered mesh \textit{vice versa}. These are...
Although melting is inferred to start at depths of 120 km or more in the mantle, little known about style melt transport from source region surface. Geochemical observations require depth be fast, and that not chemical equilibrium with matrix through which it passes. We present a model allows fast by an open conduit viscous permeable discuss possibility isolation such conduit. show disturbances on surface propagate slowly upwards as solitary waves cause flow between surrounding porous matrix.
Abstract The crystalline monolayer of 1,2-bis(4-pyridyl)ethylene physisorbed on a graphite surface at 0.44 monolayers coverage has been observed and characterized by synchrotron X-ray diffraction differential scanning calorimetry. experimentally determined structure p2 symmetry with lattice parameters = 17.77 Å, b 13.69 Å ν 39.7°. unit cell contains two molecules, which are oriented in plane parallel to the surface. It is proposed that molecules arranged such they able form weak C–H ··· N...
Breeding blankets are designed to ensure tritium self-sufficiency in deuterium–tritium fusion power plants. In addition this, breeder play a vital role shielding key components of the reactor, and provide main source heat which will ultimately be used generate electricity. Blanket design is critical success reactors integral process. Neutronic simulations regularly performed ascertain performance particular design. An iterative process improvements parametric studies required optimize meet...
The Skaergaard Intrusion of East Greenland is cut by several generations dykes, the earliest which thought to have intruded shortly after solidification Skaergaard. Two ~ 6 m wide doleritic dykes from generation are exposed in campsite area near Homestead Bay Peninsula. One (the Campsite Dyke) locally contains abundant xenoliths troctolitic cumulate. other Plagioclase-phyric large plagioclase phenocrysts. Cross-cutting relationships between two not exposed. median...
We examine the performance of matrix-free SYCL implementations conjugate gradient method for solving sparse linear systems equations. Performance is tested on an NVIDIA A100-80GB device and a dual socket Intel Ice Lake CPU node using different implementations, compared to CUDA BLAS (cuBLAS) A100 GPU MKL node. All considered kernels in implementation are memory bandwidth limited, simple model applied estimate asymptotic latency. Our experiments show that most cases match reference...