14,048 research outputs found
A least-squares implicit RBF-FD closest point method and applications to PDEs on moving surfaces
The closest point method (Ruuth and Merriman, J. Comput. Phys.
227(3):1943-1961, [2008]) is an embedding method developed to solve a variety
of partial differential equations (PDEs) on smooth surfaces, using a closest
point representation of the surface and standard Cartesian grid methods in the
embedding space. Recently, a closest point method with explicit time-stepping
was proposed that uses finite differences derived from radial basis functions
(RBF-FD). Here, we propose a least-squares implicit formulation of the closest
point method to impose the constant-along-normal extension of the solution on
the surface into the embedding space. Our proposed method is particularly
flexible with respect to the choice of the computational grid in the embedding
space. In particular, we may compute over a computational tube that contains
problematic nodes. This fact enables us to combine the proposed method with the
grid based particle method (Leung and Zhao, J. Comput. Phys. 228(8):2993-3024,
[2009]) to obtain a numerical method for approximating PDEs on moving surfaces.
We present a number of examples to illustrate the numerical convergence
properties of our proposed method. Experiments for advection-diffusion
equations and Cahn-Hilliard equations that are strongly coupled to the velocity
of the surface are also presented
Towards a unified linear kinetic transport model with the trace ion module for EIRENE
Linear kinetic Monte Carlo particle transport models are frequently employed
in fusion plasma simulations to quantify atomic and surface effects on the main
plasma flow dynamics. Separate codes are used for transport of neutral
particles (incl. radiation) and charged particles (trace impurity ions).
Integration of both modules into main plasma fluid solvers provides then self
consistent solutions, in principle. The required interfaces are far from
trivial, because rapid atomic processes in particular in the edge region of
fusion plasmas require either smoothing and resampling, or frequent transfer of
particles from one into the other Monte Carlo code. We propose a different
scheme here, in which despite the inherently different mathematical form of
kinetic equations for ions and neutrals (e.g. Fokker-Planck vs. Boltzmann
collision integrals) both types of particle orbits can be integrated into one
single code. We show that the approximations and shortcomings of this "single
sourcing" concept (e.g., restriction to explicit ion drift orbit integration)
can be fully tolerable in a wide range of typical fusion edge plasma
conditions, and be overcompensated by the code-system simplicity, as well as by
inherently ensured consistency in geometry (one single numerical grid only) and
(the common) atomic and surface process modulesComment: 15 pages, 7 figure
A level-set method for the evolution of cells and tissue during curvature-controlled growth
Most biological tissues grow by the synthesis of new material close to the
tissue's interface, where spatial interactions can exert strong geometric
influences on the local rate of growth. These geometric influences may be
mechanistic, or cell behavioural in nature. The control of geometry on tissue
growth has been evidenced in many in-vivo and in-vitro experiments, including
bone remodelling, wound healing, and tissue engineering scaffolds. In this
paper, we propose a generalisation of a mathematical model that captures the
mechanistic influence of curvature on the joint evolution of cell density and
tissue shape during tissue growth. This generalisation allows us to simulate
abrupt topological changes such as tissue fragmentation and tissue fusion, as
well as three dimensional cases, through a level-set-based method. The
level-set method developed introduces another Eulerian field than the level-set
function. This additional field represents the surface density of tissue
synthesising cells, anticipated at future locations of the interface. Numerical
tests performed with this level-set-based method show that numerical
conservation of cells is a good indicator of simulation accuracy, particularly
when cusps develop in the tissue's interface. We apply this new model to
several situations of curvature-controlled tissue evolutions that include
fragmentation and fusion.Comment: 15 pages, 10 figures, 3 supplementary figure
A Model for Optimal Human Navigation with Stochastic Effects
We present a method for optimal path planning of human walking paths in
mountainous terrain, using a control theoretic formulation and a
Hamilton-Jacobi-Bellman equation. Previous models for human navigation were
entirely deterministic, assuming perfect knowledge of the ambient elevation
data and human walking velocity as a function of local slope of the terrain.
Our model includes a stochastic component which can account for uncertainty in
the problem, and thus includes a Hamilton-Jacobi-Bellman equation with
viscosity. We discuss the model in the presence and absence of stochastic
effects, and suggest numerical methods for simulating the model. We discuss two
different notions of an optimal path when there is uncertainty in the problem.
Finally, we compare the optimal paths suggested by the model at different
levels of uncertainty, and observe that as the size of the uncertainty tends to
zero (and thus the viscosity in the equation tends to zero), the optimal path
tends toward the deterministic optimal path
Simulation of 2-dimensional viscous flow through cascades using a semi-elliptic analysis and hybrid C-H grids
A semi-elliptic formulation, termed the interacting parabolized Navier-Stokes (IPNS) formulation, is developed for the analysis of a class of subsonic viscous flows for which streamwise diffusion is neglible but which are significantly influenced by upstream interactions. The IPNS equations are obtained from the Navier-Stokes equations by dropping the streamwise viscous-diffusion terms but retaining upstream influence via the streamwise pressure-gradient. A two-step alternating-direction-explicit numerical scheme is developed to solve these equations. The quasi-linearization and discretization of the equations are carefully examined so that no artificial viscosity is added externally to the scheme. Also, solutions to compressible as well as nearly compressible flows are obtained without any modification either in the analysis or in the solution process. The procedure is applied to constricted channels and cascade passages formed by airfoils of various shapes. These geometries are represented using numerically generated curilinear boundary-oriented coordinates forming an H-grid. A hybrid C-H grid, more appropriate for cascade of airfoils with rounded leading edges, was also developed. Satisfactory results are obtained for flows through cascades of Joukowski airfoils
- …