Transboundary extraction of groundwater in the presence of hydraulic fracturing

We studied transboundary groundwater management problems in the presence of hydraulic fracturing. We found that the presence of risk suggests there should be caution when considering hydraulic fracturing. Our results from the cooperative solution show a decrease in hydraulic fracturing and increase in the steady state survival rate of groundwater. We also provide a Pigouvian type tax that could be imposed on natural gas developers.     

A Phase Field Model for Martensitic Transformations

The martensitic transformation is described using a phase field model which is in mathematical terms the regularization of a sharp interface approach. In this work, up to two martensitic orientation variants are considered. The evolution of microstructure is assumed to follow a time dependent Ginzburg-Landau equation. The coupled problem of the mechanical balance equation and the evolution equations is solved using finite elements and an implicit time integration scheme. In this work, the global energy evolution during the martensitic transformation and the influence of external loads on the formation of the different martensitic phases are studied in 2d. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Phase Field Approach for Dynamic Fracture

In phase field fracture models cracks are indicated by the value of a scalar field variable which interpolates smoothly between broken and undamaged material. The evolution equation for this crack field is coupled to the mechanical field equations in order to model the mutual interaction between the crack evolution and mechanical quantities. In finite element simulations of crack growth at comparatively slow loading velocities, a quasi-static phase field model yields reasonable results. However, the simulation of fast loading or the nucleation of new cracks challenges the limits of such a formulation. Here, the quasi-static phase field model predicts brutal crack extension with an artificially high crack speed. In this work, we analyze to which extend a dynamic formulation of the mechanical part of the phase field model can overcome this paradox created by the quasi-static formulation. In finite element simulations, the impact of the dynamic effects is studied, and differences between the crack propagation behavior of the quasi-static model and the dynamic formulation are highlighted. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Robust Error Estimates for Adaptive Phase Field Simulations

Phase field equations are commonly used as a regularized model, where bulk phases are separated by interface regions that have a thickness of the order γ. Their numerical analysis is well established for a fixed parameter size γ, but conventional error estimates depend exponentially on γ⁻¹ and thus become useless in the relevant case if γ→0. Technical applications include e.g. the simulation of Sn–Cu alloys for the production of lead free solder or Ni–Al alloys used for rotor blade surfaces. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Well-Posedness and Regularity for a Parabolic-Hyperbolic Penrose-Fife Phase Field System

This work is concerned with the study of an initial boundary value problem for a non-conserved phase field system arising from the Penrose-Fife approach to the kinetics of phase transitions. The system couples a nonlinear parabolic equation for the absolute temperature with a nonlinear hyperbolic equation for the phase variable χ, which is characterized by the presence of an inertial term multiplied by a small positive coefficient μ. This feature is the main consequence of supposing that the response of χ to the generalized force (which is the functional derivative of a free energy potential and arises as a consequence of the tendency of the free energy to decay towards a minimum) is subject to delay. We first obtain well-posedness for the resulting initial-boundary value problem in which the heat flux law contains a special function of the absolute temperature ϑ, i.e. α(ϑ) ∼ ϑ − 1/ϑ. Then we prove convergence of any family of weak solutions of the parabolic-hyperbolic model to a weak solution of the standard Penrose-Fife model as μ ↘ 0. However, the main novelty of this paper consists in proving some regularity results on solutions of the parabolic-hyperbolic system (including also estimates of Moser type) that could be useful for the study of the longterm dynamics.     

A conservative SPH scheme using exact projection with semi-analytical boundary method for free-surface flows

This paper presents a novel SPH scheme for modelling incompressible and divergence-free flow with a free surface (IDFSPH) associated with semi-analytical wall boundary conditions. In line with the projection method, the velocity field is decoupled from the pressure field in the momentum equation. A Poisson equation, serving as the pressure solver, is obtained by which pressure field is decoupled completely from the velocity field. In particular, an exact projection scheme is deployed to fulfil the requirement of the divergence-free velocity field. The condition of incompressibility is satisfied by iteratively updating the density field till the convergence. The two-equation k–ε model is employed to describe the turbulence effects in Newtonian flows. It is shown that the discretised SPH schemes have the feature of both linear and angular momentum conservations. The semi-analytical wall method implements the appropriate integrals to evaluate the boundary contributions to the mass and momentum equations. In comparison to the boundary particle methods, it can greatly enhance the feasibility and efficiency with the complex geometries. The algorithm presented within this paper is applied to several academic test cases for which either analytical results or simulations with other methods are available. The comparisons verify that this scheme is provided with convincing efficiency and extensive applicability.     

A Phase Field Approach for Martensitic Transformations and Crystal Plasticity

This work is motivated by cryogenic turning which allows end shape machining and simultaneously attaining a hardened surface due to deformation induced martensitic transformations. To study the process on the microscale, a multivariant phase field model for martensitic transformations in conjunction with a crystal plastic material model is introduced. The evolution of microstructure is assumed to follow a time-dependent Ginzburg-Landau equation. To solve the field equations the finite element method is used. Time integration is performed with Euler backward schemes, on the global level for the evolution equation of the phase field, and on the element level for the crystal plastic material law. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

轨道结构随机场模型与车辆-轨道耦合随机动力分析

将轨道结构视为一个参数随机系统,提出并建立了轨道结构的随机场模型.利用车辆-轨道耦合动力学的基本方法,将轨道系统有限单元模型与多刚体车辆模型相结合,建立了考虑铁路线路参数空-时随机变化的车辆-轨道动力计算模型.算例表明:所提出的方法较为可靠且高效;线路参数随机性对车辆-轨道系统的动力响应有明显的影响,随线路参数离散程度的增加,可能造成行车不安全、轨道损伤加剧等一些问题.     

SPH method applied to compression of solid materials for a variety of loading conditions

Smoothed Particle Hydrodynamics (SPH) is a numerical method that does not use a mesh or grid when solving a set of partial differential equations. This makes it particularly useful in application to solid mechanics problems where the sample undergoes large deformation. Whereas mesh-based methods have difficulty when the sample becomes severely distorted, SPH naturally deals with this important engineering scenario. We implement the SPH method for compressional deformation of solid samples and focus on uniaxial, biaxial and triaxial loading. We develop numerical procedures that naturally deal with these three different sets of boundary conditions and apply it to both small and larger strains in elastic and more complex materials. The method is shown to be robust up to large strains of 30%. Under uniaxial loading, a cylindrical sample tends to deform by bulging while under triaxial loading the cylindrical sample will remain cylindrical, but the diameter of the sample increases accordingly.     

Investigation of a Phase Field Model for Elasto-plastic Fracture

Phase field modelling of brittle fracture is very well understood today. However, the attempts of investigation of elasto-plastic fracture by the phase field approach are limited. This contribution deals with the investigation of a phase field model for elasto-plastic fracture. Based on a free energy density comprising elastic, fracture and plastic contributions, the model describes an extension of the linear elastic model towards von Mises plasticity. In this work it is analyzed numerically to which extend analytical findings concerning the interpretation of the model parameters in 1D are transferable to 2D scenarios. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Exponential Finite Elements for a Phase Field Fracture Model

Sharp interface material models can be related to phase field models by introducing an order parameter, whose value is assigned to the different phases of a material. The elastic material law is coupled to the evolution equation of the order parameter and cracking is addressed as a phase transition problem instead of a moving boundary value problem. A regularization parameter ϵ controls the width of the diffuse cracks represented by the order parameter and the underlying sharp interface model can be recovered from the phase field model by the limit ϵ → 0. However, in numerical simulations using standard finite elements with linear shape functions, the minimum value of ϵ is restricted by the grid size and therefore the discretization of the crack field requires extensive mesh refinement for small values of ϵ. In this work, we construct special 2d shape functions which take into account the exponential character of the crack field and its dependence on the parameter ϵ. Especially in simulations with small values of ϵ and a rather coarse mesh, the elements with exponential shape functions perform significantly better than standard linear elements. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

非线性耦合标量场方程显式解析解的研究

利用两种不同的变换,获得了一类非线性耦合标量场方程的若干类型的精确解析解,其中包括孤子解、奇性孤波解和三角函数解,从而丰富了方程解的内容。这些结论可以应用于其它的非线性方程。此外还纠正了一些文献的部分结论。     

One Non-Relativistic Particle Coupled to a Photon Field

We investigate the ground state energy of an electron coupled to a photon field. First, we regard the self-energy of a “free” electron, which we describe by the Pauli-Fierz Hamiltonian. We show that, in the case of small values of the coupling constant a \alpha , the leading order term is represented by 2ap^-1(L - ln[1 + L]) 2\alpha\pi^{-1}(\Lambda - \textrm{ln}[1 + \Lambda]) .¶ Next we put the electron in the field of an arbitrary external potential V , such that the corresponding Schrödinger operator p² + V has at least one eigenvalue, and show that by coupling to the radiation field the binding energy increases, at least for small enough values of the coupling constant a \alpha . Moreover, we provide concrete numbers for a \alpha , the ultraviolet cut-off, and the radiative correction for which our procedure works.     

Explicit Traveling Wave Solutions and Their Dynamical Behaviors for the Coupled Higgs Field Equation

In this paper, we focus on the traveling wave solutions of the coupled Higgs field equation from the perspective of dynamical systems. Through the phase portraits, in addition to periodic wave solutions and solitary wave solutions, we also obtain explicit periodic singular wave solutions, singular wave solutions and kink wave solutions, which were not found in the previous works. The dynamical behavior of these solutions and their internal relations are revealed through asymptotic analysis. The results will help supplement the study of field equation.     

Coupled boundary element method and finite difference method for the heat conduction in laser processing

This paper is presented as a way to model transient heat conduction in a 3-D axisymmetric case where large rates of heat fluxes are applied on the surfaces as done in the case of laser processing. This would result in large temperature gradients in a small area irradiated by the laser on the incident surface that could also reach melting and subsequent vaporization. BEM can handle large fluxes very easily and it also can be formulated if needed to incorporate the moving boundary problem in a unique manner while on the other hand FDM is a fast and efficient method. For these reasons a coupled BEM–FDM method is formulated to simulate the heat conduction process. In the BEM method linear elements for the boundary and quadratic elements for the domain were used. The integrals in BEM were integrated in time using the asymptotic expansion for the modified Bessel functions in the Green’s function. To further improve the accuracy, special techniques were employed in the spatial integration. As for the FDM formulation, a flux conservation scheme with a 4th order formula for the fluxes was used. The FDM and BEM were coupled at the interface by the temperature from the FDM formulation being imposed on the BEM and the flux from the BEM being utilized by the FDM elements near to the interface. To advance in time, the Crank–Nicholson scheme was used on the FDM directly and due to coupling indirectly on the BEM. The relative errors for the simulation of constant and variable flux cases demonstrate the successful nature of the numerical model.     

A Phase Field Approach for Martensitic Transformations in Elastoplastic Materials

A multivariant phase field model for martensitic transformations in elastoplastic materials is introduced which is in mathematical terms the regularization of a sharp interface approach. The evolution of microstructure is assumed to follow a time dependent Ginzburg-Landau equation. The coupled problem of the mechanical balance equation and the evolution equations is solved using finite elements and an implicit time integration scheme. In this work, plasticity is considered for the austenitic phase which influences the martensitic evolution. With aid of the model these interactions are studied in detail. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Hydraulic fracturing modeling using the enriched numerical manifold method

Recent attempts to solve rock mechanics problems using the numerical manifold method (NMM) have been regarded as fruitful. In this paper, a coupled hydro-mechanical (HM) model is incorporated into the enriched NMM to simulate fluid driven fracturing in rocks. In this HM model, a “cubic law” is employed to model fluid flow through fractures. Several benchmark problems are investigated to verify the coupled HM model. The simulation results agree well with the analytical and experimental results, indicating that the coupled HM model is able to simulate the hydraulic fracturing process reliably and correctly.     

Heat Conduction in a Phase Field Model for Martensitic Transformation

We study the martensitic transformation with a phase field model, where we consider the Bain transformation path in a small strain setting. For the order parameter, interpolating between an austenitic parent phase and martensitic phases, we use a Ginzburg-Landau evolution equation, assuming a constant mobility. In [1], a temperature dependent separation potential is introduced. We use this potential to extend the model in [2], by considering a transient temperature field, where the temperature is introduced as an additional degree of freedom. This leads to a coupling of both the evolution equation of the order parameter and the mechanical field equations (in terms of thermal expansion) with the heat equation. The model is implemented in FEAP as a 4-node element with bi-linear shape functions. Numerical examples are given to illustrate the influence of the temperature on the evolution of the martensitic phase. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)