10,282 research outputs found
Non-Linear Model Predictive Control with Adaptive Time-Mesh Refinement
In this paper, we present a novel solution for real-time, Non-Linear Model
Predictive Control (NMPC) exploiting a time-mesh refinement strategy. The
proposed controller formulates the Optimal Control Problem (OCP) in terms of
flat outputs over an adaptive lattice. In common approximated OCP solutions,
the number of discretization points composing the lattice represents a critical
upper bound for real-time applications. The proposed NMPC-based technique
refines the initially uniform time horizon by adding time steps with a sampling
criterion that aims to reduce the discretization error. This enables a higher
accuracy in the initial part of the receding horizon, which is more relevant to
NMPC, while keeping bounded the number of discretization points. By combining
this feature with an efficient Least Square formulation, our solver is also
extremely time-efficient, generating trajectories of multiple seconds within
only a few milliseconds. The performance of the proposed approach has been
validated in a high fidelity simulation environment, by using an UAV platform.
We also released our implementation as open source C++ code.Comment: In: 2018 IEEE International Conference on Simulation, Modeling, and
Programming for Autonomous Robots (SIMPAR 2018
Real-time Error Control for Surgical Simulation
Objective: To present the first real-time a posteriori error-driven adaptive
finite element approach for real-time simulation and to demonstrate the method
on a needle insertion problem. Methods: We use corotational elasticity and a
frictional needle/tissue interaction model. The problem is solved using finite
elements within SOFA. The refinement strategy relies upon a hexahedron-based
finite element method, combined with a posteriori error estimation driven local
-refinement, for simulating soft tissue deformation. Results: We control the
local and global error level in the mechanical fields (e.g. displacement or
stresses) during the simulation. We show the convergence of the algorithm on
academic examples, and demonstrate its practical usability on a percutaneous
procedure involving needle insertion in a liver. For the latter case, we
compare the force displacement curves obtained from the proposed adaptive
algorithm with that obtained from a uniform refinement approach. Conclusions:
Error control guarantees that a tolerable error level is not exceeded during
the simulations. Local mesh refinement accelerates simulations. Significance:
Our work provides a first step to discriminate between discretization error and
modeling error by providing a robust quantification of discretization error
during simulations.Comment: 12 pages, 16 figures, change of the title, submitted to IEEE TBM
Simulating radiative shocks in nozzle shock tubes
We use the recently developed Center for Radiative Shock Hydrodynamics
(CRASH) code to numerically simulate laser-driven radiative shock experiments.
These shocks are launched by an ablated beryllium disk and are driven down
xenon-filled plastic tubes. The simulations are initialized by the
two-dimensional version of the Lagrangian Hyades code which is used to evaluate
the laser energy deposition during the first 1.1ns. The later times are
calculated with the CRASH code. This code solves for the multi-material
hydrodynamics with separate electron and ion temperatures on an Eulerian
block-adaptive-mesh and includes a multi-group flux-limited radiation diffusion
and electron thermal heat conduction. The goal of the present paper is to
demonstrate the capability to simulate radiative shocks of essentially
three-dimensional experimental configurations, such as circular and elliptical
nozzles. We show that the compound shock structure of the primary and wall
shock is captured and verify that the shock properties are consistent with
order-of-magnitude estimates. The produced synthetic radiographs can be used
for comparison with future nozzle experiments at high-energy-density laser
facilities.Comment: submitted to High Energy Density Physic
Monolithic simulation of convection-coupled phase-change - verification and reproducibility
Phase interfaces in melting and solidification processes are strongly
affected by the presence of convection in the liquid. One way of modeling their
transient evolution is to couple an incompressible flow model to an energy
balance in enthalpy formulation. Two strong nonlinearities arise, which account
for the viscosity variation between phases and the latent heat of fusion at the
phase interface.
The resulting coupled system of PDE's can be solved by a single-domain
semi-phase-field, variable viscosity, finite element method with monolithic
system coupling and global Newton linearization. A robust computational model
for realistic phase-change regimes furthermore requires a flexible
implementation based on sophisticated mesh adaptivity. In this article, we
present first steps towards implementing such a computational model into a
simulation tool which we call Phaseflow.
Phaseflow utilizes the finite element software FEniCS, which includes a
dual-weighted residual method for goal-oriented adaptive mesh refinement.
Phaseflow is an open-source, dimension-independent implementation that, upon an
appropriate parameter choice, reduces to classical benchmark situations
including the lid-driven cavity and the Stefan problem. We present and discuss
numerical results for these, an octadecane PCM convection-coupled melting
benchmark, and a preliminary 3D convection-coupled melting example,
demonstrating the flexible implementation. Though being preliminary, the latter
is, to our knowledge, the first published 3D result for this method. In our
work, we especially emphasize reproducibility and provide an easy-to-use
portable software container using Docker.Comment: 20 pages, 8 figure
On Validating an Astrophysical Simulation Code
We present a case study of validating an astrophysical simulation code. Our
study focuses on validating FLASH, a parallel, adaptive-mesh hydrodynamics code
for studying the compressible, reactive flows found in many astrophysical
environments. We describe the astrophysics problems of interest and the
challenges associated with simulating these problems. We describe methodology
and discuss solutions to difficulties encountered in verification and
validation. We describe verification tests regularly administered to the code,
present the results of new verification tests, and outline a method for testing
general equations of state. We present the results of two validation tests in
which we compared simulations to experimental data. The first is of a
laser-driven shock propagating through a multi-layer target, a configuration
subject to both Rayleigh-Taylor and Richtmyer-Meshkov instabilities. The second
test is a classic Rayleigh-Taylor instability, where a heavy fluid is supported
against the force of gravity by a light fluid. Our simulations of the
multi-layer target experiments showed good agreement with the experimental
results, but our simulations of the Rayleigh-Taylor instability did not agree
well with the experimental results. We discuss our findings and present results
of additional simulations undertaken to further investigate the Rayleigh-Taylor
instability.Comment: 76 pages, 26 figures (3 color), Accepted for publication in the ApJ
- …