Kantorovich-Rubinstein misfit for inverting gravity-gradient data by the level-set method

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. G55-G73
Author(s):  
Guanghui Huang ◽  
Xinming Zhang ◽  
Jianliang Qian
Keyword(s):  
Level Set Method ◽  
Level Set ◽  
Random Noise ◽  
Magnetic Data ◽  
Gravity Gradient ◽  
Local Minima ◽  
Target Model ◽  
Level Set Function ◽  
Level Set Evolution ◽  
Set Function

We have developed a novel Kantorovich-Rubinstein (KR) norm-based misfit function to measure the mismatch between gravity-gradient data for the inverse gradiometry problem. Under the assumption that an anomalous mass body has an unknown compact support with a prescribed constant value of density contrast, we implicitly parameterize the unknown mass body by a level-set function. Because the geometry of an underlying anomalous mass body may experience various changes during inversion in terms of level-set evolution, the classic least-squares ([Formula: see text]-norm-based) and the [Formula: see text]-norm-based misfit functions for governing the level-set evolution may potentially induce local minima if an initial guess of the level-set function is far from that of the target model. The KR norm from the optimal transport theory computes the data misfit by comparing the modeled data and the measured data in a global manner, leading to better resolution of the differences between the inverted model and the target model. Combining the KR norm with the level-set method yields a new effective methodology that is not only able to mitigate local minima but is also robust against random noise for the inverse gradiometry problem. Numerical experiments further demonstrate that the new KR norm-based misfit function is able to recover deep dipping flanks of SEG/EAGE salt models even at extremely low signal-to-noise ratios. The new methodology can be readily applied to gravity and magnetic data as well.

A local level-set method for 3D inversion of gravity-gradient data

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. G35-G51 ◽  
Author(s):  
Wangtao Lu ◽  
Jianliang Qian
Keyword(s):  
Inverse Problem ◽  
Level Set Method ◽  
Level Set ◽  
Local Level ◽  
Density Contrast ◽  
Gravity Gradient ◽  
Gravity Center ◽  
Level Set Function ◽  
Local Level Set ◽  
Set Function

We have developed a local level-set method for inverting 3D gravity-gradient data. To alleviate the inherent nonuniqueness of the inverse gradiometry problem, we assumed that a homogeneous density contrast distribution with the value of the density contrast specified a priori was supported on an unknown bounded domain [Formula: see text] so that we may convert the original inverse problem into a domain inverse problem. Because the unknown domain [Formula: see text] may take a variety of shapes, we parametrized the domain [Formula: see text] by a level-set function implicitly so that the domain inverse problem was reduced to a nonlinear optimization problem for the level-set function. Because the convergence of the level-set algorithm relied heavily on initializing the level-set function to enclose the gravity center of a source body, we applied a weighted [Formula: see text]-regularization method to locate such a gravity center so that the level-set function can be properly initialized. To rapidly compute the gradient of the nonlinear functional arising in the level-set formulation, we made use of the fact that the Laplacian kernel in the gravity force relation decayed rapidly off the diagonal so that matrix-vector multiplications for evaluating the gradient can be accelerated significantly. We conducted extensive numerical experiments to test the performance and effectiveness of the new method.

Numerical Simulation of Breaking Waves Using a Level Set Method

2002 ◽  
Author(s):  
Pablo Go´mez ◽  
Julio Herna´ndez ◽  
Joaqui´n Lo´pez ◽  
Fe´lix Faura
Keyword(s):  
Free Surface ◽  
Level Set Method ◽  
Level Set ◽  
Stokes Equations ◽  
Numerical Study ◽  
Operating Conditions ◽  
Breaking Wave ◽  
Level Set Function ◽  
Set Function ◽  
The Impact

A numerical study of the initial stages of wave breaking processes in shallow water is presented. The waves considered are assumed to be generated by moving a piston in a two-dimensional channel, and may appear, for example, in the injection chamber of a high-pressure die casting machine under operating conditions far from the optimal. A numerical model based on a finite-difference discretization of the Navier-Stokes equations in a Cartesian grid and a second-order approximate projection method has been developed and used to carry out the simulations. The evolution of the free surface is described using a level set method, with a reinitialization procedure of the level set function which uses a local grid refinement near the free surface. The ability of different algorithms to improve mass conservation in the reinitialization step of the level set function has been tested in a time-reversed single vortex flow. The results for the breaking wave profiles show the flow characteristics after the impact of the first plunging jet onto the wave’s forward face and during the subsequent splash-up.

Parametric Shape and Topology Optimization: A New Level Set Approach Based on Cardinal Kernel Functions

2017 ◽  
Author(s):  
Long Jiang ◽  
Shikui Chen ◽  
Xiangmin Jiao
Keyword(s):  
Topology Optimization ◽  
Level Set Method ◽  
Level Set ◽  
Material Property ◽  
Kernel Function ◽  
Level Set Methods ◽  
Level Set Function ◽  
Set Function ◽  
Distance Regularization ◽  
Parametric Level Set Method

The parametric level set method is an extension of the conventional level set methods for topology optimization. By parameterizing the level set function, conventional levels let methods can be easily coupled with mathematical programming to achieve better numerical robustness and computational efficiency. Furthermore, the parametric level set scheme not only can inherit the original advantages of the conventional level set methods, such as clear boundary representation and high topological changes handling flexibility but also can alleviate some un-preferred features from the conventional level set methods, such as needing re-initialization. However, in the RBF-based parametric level set method, it was difficult to determine the range of the design variables. Moreover, with the mathematically driven optimization process, the level set function often results in significant fluctuations during the optimization process. This brings difficulties in both numerical stability control and material property interpolation. In this paper, an RBF partition of unity collocation method is implemented to create a new type of kernel function termed as the Cardinal Basis Function (CBF), which employed as the kernel function to parameterize the level set function. The advantage of using the CBF is that the range of the design variable, which was the weight factor in conventional RBF, can be explicitly specified. Additionally, a distance regularization energy functional is introduced to maintain a desired distance regularized level set function evolution. With this desired distance regularization feature, the level set evolution is stabilized against significant fluctuations. Besides, the material property interpolation from the level set function to the finite element model can be more accurate.

Design of Compliant Thermal Actuators Using Structural Optimization Based on the Level Set Method

2011 ◽  
Vol 11 (1) ◽  
Author(s):  
Takayuki Yamada ◽  
Shintaro Yamasaki ◽  
Shinji Nishiwaki ◽  
Kazuhiro Izui ◽  
Masataka Yoshimura
Keyword(s):  
Structural Optimization ◽  
Level Set Method ◽  
Level Set ◽  
Compliant Mechanisms ◽  
Optimization Method ◽  
Equilibrium Equations ◽  
Thermal Actuators ◽  
Level Set Function ◽  
Initialization Scheme ◽  
Set Function

Compliant mechanisms are designed to be flexible to achieve a specified motion as a mechanism. Such mechanisms can function as compliant thermal actuators in micro-electromechanical systems by intentionally designing configurations that exploit thermal expansion effects in elastic material when appropriate portions of the mechanism structure are heated or are subjected to an electric potential. This paper presents a new structural optimization method for the design of compliant thermal actuators based on the level set method and the finite element method (FEM). First, an optimization problem is formulated that addresses the design of compliant thermal actuators considering the magnitude of the displacement at the output location. Next, the topological derivatives that are used when introducing holes during the optimization process are derived. Based on the optimization formulation, a new structural optimization algorithm is constructed that employs the FEM when solving the equilibrium equations and updating the level set function. The re-initialization of the level set function is performed using a newly developed geometry-based re-initialization scheme. Finally, several design examples are provided to confirm the usefulness of the proposed structural optimization method.

A Meshless Level Set Method for Shape and Topology Optimization

2011 ◽  
Vol 308-310 ◽  
pp. 1046-1049 ◽  
Author(s):  
Yu Wang ◽  
Zhen Luo
Keyword(s):  
Topology Optimization ◽  
Free Boundary ◽  
Level Set Method ◽  
Level Set ◽  
Level Set Methods ◽  
Discrete Level ◽  
Shape And Topology Optimization ◽  
Level Set Function ◽  
And Topology ◽  
Set Function

This paper proposes a meshless Galerkin level set method for structural shape and topology optimization of continua. To taking advantage of the implicit free boundary representation scheme, structural design boundary is represented through the introduction of a scalar level set function as its zero level set, to flexibly handle complex shape fidelity and topology changes by maintaining concise and smooth interface. Compactly supported radial basis functions (CSRBFs) are used to parameterize the level set function and also to construct the shape functions for mesh free function approximation. The meshless Galerkin global weak formulation is employed to implement the discretization of the state equations. This provides a pathway to simplify two numerical procedures involved in most conventional level set methods in propagating the discrete level set functions and in approximating the discrete equations, by unifying the two different stages at two sets of grids just in terms of one set of scattered nodes. The proposed level set method has the capability of describing the implicit moving boundaries without remeshing for discontinuities. The motion of the free boundary is just a question of advancing the discrete level set function by finding the design variables of the size optimization in time. One benchmark example is used to demonstrate the effectiveness of the proposed method. The numerical results showcase that this method has the ability to simplify numerical procedures and to avoid numerical difficulties happened in most conventional level set methods. It is straightforward to apply the present method to more advanced shape and topology optimization problems.

scholarly journals Using the level set method in geodynamical modeling of multi-material flows and Earth's free surface

2014 ◽  
Vol 6 (2) ◽  
pp. 1523-1554 ◽  
Author(s):  
B. Hillebrand ◽  
C. Thieulot ◽  
T. Geenen ◽  
A. P. van den Berg ◽  
W. Spakman
Keyword(s):  
Free Surface ◽  
Level Set Method ◽  
Level Set ◽  
Material Flow ◽  
Flow Modeling ◽  
Material Interfaces ◽  
Level Set Function ◽  
Geophysical Processes ◽  
Rayleigh Taylor Instability ◽  
Set Function

Abstract. The level set method allows for tracking material surfaces in 2-D and 3-D flow modeling and is well suited for applications of multi-material flow modeling. The level set method utilizes smooth level set functions to define material interfaces, which makes the method stable and free of oscillations that are typically observed in case step-like functions parameterize interfaces. By design the level set function is a signed distance function and gives for each point in the domain the exact distance to the interface and on which side it is located. In this paper we present four benchmarks which show the validity, accuracy and simplicity of using the level set method for multi-material flow modeling. The benchmarks are simplified setups of dynamical geophysical processes such as a Rayleigh–Taylor instability, post glacial rebound, subduction and slab detachment. We also demonstrate the benefit of using the level set method for modeling a free surface with the sticky air approach. Our results show that the level set method allows for accurate material flow modeling and that the combination with the sticky air approach works well in mimicking Earth's free surface. Since the level set method tracks material interfaces instead of materials themselves, it has the advantage that the location of these interfaces is accurately known and that it represents a viable alternative to the more commonly used tracer method.

A three-dimensional one-layer particle level set method

2019 ◽  
Vol 30 (7) ◽  
pp. 3653-3684
Author(s):  
LanHao Zhao ◽  
Kailong Mu ◽  
Jia Mao ◽  
Khuc Hongvan ◽  
Dawei Peng
Keyword(s):  
Level Set Method ◽  
Level Set ◽  
Three Dimensional ◽  
Single Vortex ◽  
Content Type ◽  
Moving Interface ◽  
Level Set Function ◽  
Particle Level ◽  
Moving Interface Problems ◽  
Set Function

Purpose Moving interface problems exist commonly in nature and industry, and the main difficulty is to represent the interface. The purpose of this paper is to capture the accurate interface, a novel three-dimensional one-layer particle level set (OPLS) method is presented by introducing Lagrangian particles to reconstruct the seriously distorted level set function. Design/methodology/approach First, the interface is captured by the level set method. Then, the interface is corrected with only one-layer particles advected with the flow to ensure that the level set function value of the particle is equal to 0. When interfaces are merged, all particles in merged regions are deleted, while the added particles near the generated interface are used to determine the interface as the interface is separated. Findings The OPLS method is validated with well-known benchmark examples, such as the long-term advection of a sphere, the rotation of a three-dimensional slotted disk and sphere, single vortex in a box, sphere merging and separation, deformation of a sphere. The simulation results indicate that the proposed method is found to be highly reliable and accurate. Originality/value This method exhibits excellent conservation of the area bounded by the interface. The extraordinary performance is also shown in dealing with complex interface topological changes.

scholarly journals A multi-discretization scheme for topology optimization based on the parameterized level set method

2020 ◽  
Vol 11 ◽  
pp. 3
Author(s):  
Peng Wei ◽  
Yang Liu ◽  
Zuyu Li
Keyword(s):  
Structural Analysis ◽  
Level Set Method ◽  
Level Set ◽  
Three Dimension ◽  
Theoretical Level ◽  
Structure Topology ◽  
Topological Design ◽  
Level Set Function ◽  
Set Function ◽  
Design Freedom

In the framework of the parameterized level set method, the structural analysis and topology representation can be implemented in a decoupling way. A parameterized level set function, typically, using radial basis functions (RBFs), is a linear combination of a set of prescribed RBFs and coefficients. Once the coefficients are determined, the theoretical level set function is determined. Exploiting this inherent property, we propose a multi-discretization method based on the parameterized level set method. In this approach, a coarse discretization is applied to do the structural analysis whereas another dense discretization is employed to represent the structure topology. As a result, both efficient analysis and high-resolution topological design are available. Note that the dense discretization only accounts for a more precise and smooth description of the theoretical level set function rather than introduce extra design freedom or incur interference to structural analysis or the optimization process. In other words, this decoupling way will not add to the computational burden of structural analysis or result in non-uniqueness of converged results for a particular analysis setting. Numerical examples in both two-dimension and three-dimension show effectiveness and applicability of the proposed method.

Simulation of Wave-Body Interaction Using a Single-Phase Level Set Function in the SWENSE Method

2013 ◽  
Author(s):  
Gabriel Reliquet ◽  
Aurélien Drouet ◽  
Pierre-Emmanuel Guillerm ◽  
Erwan Jacquin ◽  
Lionel Gentaz ◽  
...  
Keyword(s):  
Free Surface ◽  
Level Set Method ◽  
Level Set ◽  
Single Phase ◽  
Stokes Equations ◽  
Flow Theory ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Level Set Function ◽  
Set Function

The purpose of this paper is to present combination of the SWENSE (Spectral Wave Explicit Navier-Stokes Equations – [1]) method — an original method to treat fully nonlinear wave-body interactions — and a free surface RANSE (Reynolds Averaged Navier-Stokes Equations) solver using a single-phase Level Set method to capture the interface. The idea is to be able to simulate wave-body interactions under viscous flow theory with strong deformations of the interface (wave breaking in the vicinity of the body, green water on ship decks…), while keeping the advantages of the SWENSE scheme. The SWENSE approach is based on a physical decomposition by combining incident waves described by a nonlinear spectral scheme based on potential flow theory and an adapted Navier-Stokes solver where only the diffracted part of the flow is solved, incident flow parameters seen as forcing terms. In the single-phase Level Set method [2, 3], the air phase is neglected. Thus, only the liquid phase is solved considering a fluid with uniform properties. The location of the free surface is determined by a Level Set function initialised as the signed distance. The accuracy of simulation depends essentially on the pressure scheme used to impose free surface dynamic boundary condition. Comparisons of numerical results with experimental and numerical data for US navy combatant DTMB 5415 in calm water and in head waves are presented.

scholarly journals Adaptive Regularized Level Set Method for Weak Boundary Object Segmentation

2012 ◽  
Vol 2012 ◽  
pp. 1-16 ◽  
Author(s):  
Meng Li ◽  
Chuanjiang He ◽  
Yi Zhan
Keyword(s):  
Level Set Method ◽  
Level Set ◽  
Active Contours ◽  
Boundary Object ◽  
Level Curve ◽  
External Energy ◽  
Zero Level ◽  
Level Set Function ◽  
Wide Range ◽  
Level Set Evolution

An adaptive regularized level set method for image segmentation is proposed. A weightedp(x)-Dirichlet integral is presented as a geometric regularization on zero level curve, which is used to diminish the influence of image noise on level set evolution while ensuring the active contours not to pass through weak object boundaries. The idea behind the new energy integral is that the amount of regularization on the zero level curve can be adjusted automatically by the variable exponentp(x)to fit the image data. This energy is then incorporated into a level set formulation with an external energy term that drives the motion of the zero level set toward the desired objects boundaries, and a level set function regularization term that is necessary for maintaining stable level set evolution. The proposed model has been applied to a wide range of both real and synthetic images with promising results.

Sign in / Sign up

Export Citation Format

Share Document