We develop a special phase field/diffusive interface method to model the nuclear architecture reorganization process. In particular, we use a Lagrange multiplier approach in the phase field model to preserve the specific physical and geometrical constraints for the biological events. We develop several efficient and robust linear and weakly nonlinear schemes for this new model. To validate the model and numerical methods, we present ample numerical simulations which in particular reproduce several processes of nuclear architecture reorganization from the experiment literature.
We introduce a phase-field method for continuous modeling of cracks with frictional contacts. Compared with standard discrete methods for frictional contacts, the phase-field method has two attractive features: (1) it can represent arbitrary crack geometry without an explicit function or basis enrichment, and (2) it does not require an algorithm for imposing contact constraints. The first feature, which is common in phase-field models of fracture, is attained by regularizing a sharp interface geometry using a surface density functional. The second feature, which is a unique advantage for contact problems, is achieved by a new approach that calculates the stress tensor in the regularized interface region depending on the contact condition of the interface. Particularly, under a slip condition, this approach updates stress components in the slip direction using a standard contact constitutive law, while making other stress components compatible with stress in the bulk region to ensure non-penetrating deformation in other directions. We verify the proposed phase-field method using stationary interface problems simulated by discrete methods in the literature. Subsequently, by allowing the phase field to evolve according to brittle fracture theory, we demonstrate the proposed methods capability for modeling crack growth with frictional contact.
We compare the accuracy, convergence rate and computational cost of eigenerosion (EE) and phase-field (PF) methods. For purposes of comparison, we specifically consider the standard test case of a center-crack panel loaded in biaxial tension and assess the convergence of the energy error as the length scale parameter and mesh size tend to zero simultaneously. The panel is discretized by means of a regular mesh consisting of standard bilinear or Q1 elements. The exact stresses from the known analytical linear elastic solution are applied to the boundary. All element integrals over the interior and the boundary of the domain are evaluated exactly using the symbolic computation program Mathematica. When the EE inelastic energy is enhanced by means of Richardson extrapolation, EE is found to converge at twice the rate of PF and to exhibit much better accuracy. In addition, EE affords a one-order-of-magnitude computational speed-up over PF.
Techniques from numerical bifurcation theory are very useful to study transitions between steady fluid flow patterns and the instabilities involved. Here, we provide computational methodology to use parameter continuation in determining probability density functions of systems of stochastic partial differential equations near fixed points, under a small noise approximation. Key innovation is the efficient solution of a generalized Lyapunov equation using an iterative method involving low-rank approximations. We apply and illustrate the capabilities of the method using a problem in physical oceanography, i.e. the occurrence of multiple steady states of the Atlantic Ocean circulation.
In this paper, we present an efficient numerical algorithm for solving the time-dependent Cahn--Hilliard--Navier--Stokes equations that model the flow of two phases with different densities. The pressure-correction step in the projection method consists of a Poisson problem with a modified right-hand side. Spatial discretization is based on discontinuous Galerkin methods with piecewise linear or piecewise quadratic polynomials. Flux and slope limiting techniques successfully eliminate the bulk shift, overshoot and undershoot in the order parameter, which is shown to be bound-preserving. Several numerical results demonstrate that the proposed numerical algorithm is effective and robust for modeling two-component immiscible flows in porous structures and digital rocks.
Several aspects influence corrosive processes in RC structures, such as environmental conditions, structural geometry, and mechanical properties. Since these aspects present large randomnesses, probabilistic models allow a more accurate description of the corrosive phenomena. On the other hand, the definition of limit states, applied in the reliability assessment, requires a proper mechanical model. In this context, this study proposes an accurate methodology for the mechanical-probabilistic modelling of RC structures subjected to reinforcements corrosion. To this purpose, an improved damage approach is proposed to define the limit states for the probabilistic modelling, considering three main degradation phenomena: concrete cracking, rebar yielding, and rebar corrosion caused either by chlorides or carbonation process. The stochastic analysis is evaluated by the Monte Carlo simulation method due to the computational efficiency of the LDMC. The proposed mechanical-probabilistic methodology is implemented in a computational framework and applied to the analysis of a simply supported RC beam, and a 2D RC frame. Curves illustrate the probability of failure over a service life of 50 years. Moreover, the proposed model allows drawing the probability of failure map and then identify the critical failure path for progressive collapse analysis. Collapse path changes caused by the corrosion phenomena are observed.
Qing Cheng
,Pourya Delafrouz
,Jie Liang
.
(2021)
.
"Modeling and simulation of nuclear architecture reorganization process using a phase field approach"
.
Qing Cheng
هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا