26,776 research outputs found
The diffuse Nitsche method: Dirichlet constraints on phase-field boundaries
We explore diffuse formulations of Nitsche's method for consistently imposing Dirichlet boundary conditions on phase-field approximations of sharp domains. Leveraging the properties of the phase-field gradient, we derive the variational formulation of the diffuse Nitsche method by transferring all integrals associated with the Dirichlet boundary from a geometrically sharp surface format in the standard Nitsche method to a geometrically diffuse volumetric format. We also derive conditions for the stability of the discrete system and formulate a diffuse local eigenvalue problem, from which the stabilization parameter can be estimated automatically in each element. We advertise metastable phase-field solutions of the Allen-Cahn problem for transferring complex imaging data into diffuse geometric models. In particular, we discuss the use of mixed meshes, that is, an adaptively refined mesh for the phase-field in the diffuse boundary region and a uniform mesh for the representation of the physics-based solution fields. We illustrate accuracy and convergence properties of the diffuse Nitsche method and demonstrate its advantages over diffuse penalty-type methods. In the context of imaging based analysis, we show that the diffuse Nitsche method achieves the same accuracy as the standard Nitsche method with sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field, the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method by analyzing stresses in a human vertebral body
Numerical methods for time-fractional evolution equations with nonsmooth data: a concise overview
Over the past few decades, there has been substantial interest in evolution
equations that involving a fractional-order derivative of order
in time, due to their many successful applications in
engineering, physics, biology and finance. Thus, it is of paramount importance
to develop and to analyze efficient and accurate numerical methods for reliably
simulating such models, and the literature on the topic is vast and fast
growing. The present paper gives a concise overview on numerical schemes for
the subdiffusion model with nonsmooth problem data, which are important for the
numerical analysis of many problems arising in optimal control, inverse
problems and stochastic analysis. We focus on the following aspects of the
subdiffusion model: regularity theory, Galerkin finite element discretization
in space, time-stepping schemes (including convolution quadrature and L1 type
schemes), and space-time variational formulations, and compare the results with
that for standard parabolic problems. Further, these aspects are showcased with
illustrative numerical experiments and complemented with perspectives and
pointers to relevant literature.Comment: 24 pages, 3 figure
Phase-field boundary conditions for the voxel finite cell method: surface-free stress analysis of CT-based bone structures
The voxel finite cell method employs unfitted finite element meshes and voxel quadrature rules to seamlessly
transfer CT data into patient-specific bone discretizations. The method, however, still requires the explicit
parametrization of boundary surfaces to impose traction and displacement boundary conditions, which
constitutes a potential roadblock to automation. We explore a phase-field based formulation for imposing
traction and displacement constraints in a diffuse sense. Its essential component is a diffuse geometry model
generated from metastable phase-field solutions of the Allen-Cahn problem that assumes the imaging data as
initial condition. Phase-field approximations of the boundary and its gradient are then employed to transfer
all boundary terms in the variational formulation into volumetric terms. We show that in the context of the
voxel finite cell method, diffuse boundary conditions achieve the same accuracy as boundary conditions
defined over explicit sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field,
the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method
by analyzing stresses in a human femur and a vertebral body
Space-modulated Stability and Averaged Dynamics
In this brief note we give a brief overview of the comprehensive theory,
recently obtained by the author jointly with Johnson, Noble and Zumbrun, that
describes the nonlinear dynamics about spectrally stable periodic waves of
parabolic systems and announce parallel results for the linearized dynamics
near cnoidal waves of the Korteweg-de Vries equation. The latter are expected
to contribute to the development of a dispersive theory, still to come.Comment: Proceedings of the "Journ\'ees \'Equations aux d\'eriv\'ees
partielles", Roscoff 201
A Multiscale Thermo-Fluid Computational Model for a Two-Phase Cooling System
In this paper, we describe a mathematical model and a numerical simulation
method for the condenser component of a novel two-phase thermosyphon cooling
system for power electronics applications. The condenser consists of a set of
roll-bonded vertically mounted fins among which air flows by either natural or
forced convection. In order to deepen the understanding of the mechanisms that
determine the performance of the condenser and to facilitate the further
optimization of its industrial design, a multiscale approach is developed to
reduce as much as possible the complexity of the simulation code while
maintaining reasonable predictive accuracy. To this end, heat diffusion in the
fins and its convective transport in air are modeled as 2D processes while the
flow of the two-phase coolant within the fins is modeled as a 1D network of
pipes. For the numerical solution of the resulting equations, a Dual
Mixed-Finite Volume scheme with Exponential Fitting stabilization is used for
2D heat diffusion and convection while a Primal Mixed Finite Element
discretization method with upwind stabilization is used for the 1D coolant
flow. The mathematical model and the numerical method are validated through
extensive simulations of realistic device structures which prove to be in
excellent agreement with available experimental data
- …