97,635 research outputs found
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
Wavelet-based Adaptive Techniques Applied to Turbulent Hypersonic Scramjet Intake Flows
The simulation of hypersonic flows is computationally demanding due to large
gradients of the flow variables caused by strong shock waves and thick boundary
or shear layers. The resolution of those gradients imposes the use of extremely
small cells in the respective regions. Taking turbulence into account
intensives the variation in scales even more. Furthermore, hypersonic flows
have been shown to be extremely grid sensitive. For the simulation of
three-dimensional configurations of engineering applications, this results in a
huge amount of cells and prohibitive computational time. Therefore, modern
adaptive techniques can provide a gain with respect to computational costs and
accuracy, allowing the generation of locally highly resolved flow regions where
they are needed and retaining an otherwise smooth distribution. An h-adaptive
technique based on wavelets is employed for the solution of hypersonic flows.
The compressible Reynolds averaged Navier-Stokes equations are solved using a
differential Reynolds stress turbulence model, well suited to predict
shock-wave-boundary-layer interactions in high enthalpy flows. Two test cases
are considered: a compression corner and a scramjet intake. The compression
corner is a classical test case in hypersonic flow investigations because it
poses a shock-wave-turbulent-boundary-layer interaction problem. The adaptive
procedure is applied to a two-dimensional confguration as validation. The
scramjet intake is firstly computed in two dimensions. Subsequently a
three-dimensional geometry is considered. Both test cases are validated with
experimental data and compared to non-adaptive computations. The results show
that the use of an adaptive technique for hypersonic turbulent flows at high
enthalpy conditions can strongly improve the performance in terms of memory and
CPU time while at the same time maintaining the required accuracy of the
results.Comment: 26 pages, 29 Figures, submitted to AIAA Journa
POD model order reduction with space-adapted snapshots for incompressible flows
We consider model order reduction based on proper orthogonal decomposition
(POD) for unsteady incompressible Navier-Stokes problems, assuming that the
snapshots are given by spatially adapted finite element solutions. We propose
two approaches of deriving stable POD-Galerkin reduced-order models for this
context. In the first approach, the pressure term and the continuity equation
are eliminated by imposing a weak incompressibility constraint with respect to
a pressure reference space. In the second approach, we derive an inf-sup stable
velocity-pressure reduced-order model by enriching the velocity reduced space
with supremizers computed on a velocity reference space. For problems with
inhomogeneous Dirichlet conditions, we show how suitable lifting functions can
be obtained from standard adaptive finite element computations. We provide a
numerical comparison of the considered methods for a regularized lid-driven
cavity problem
Reliability-based design optimization using kriging surrogates and subset simulation
The aim of the present paper is to develop a strategy for solving
reliability-based design optimization (RBDO) problems that remains applicable
when the performance models are expensive to evaluate. Starting with the
premise that simulation-based approaches are not affordable for such problems,
and that the most-probable-failure-point-based approaches do not permit to
quantify the error on the estimation of the failure probability, an approach
based on both metamodels and advanced simulation techniques is explored. The
kriging metamodeling technique is chosen in order to surrogate the performance
functions because it allows one to genuinely quantify the surrogate error. The
surrogate error onto the limit-state surfaces is propagated to the failure
probabilities estimates in order to provide an empirical error measure. This
error is then sequentially reduced by means of a population-based adaptive
refinement technique until the kriging surrogates are accurate enough for
reliability analysis. This original refinement strategy makes it possible to
add several observations in the design of experiments at the same time.
Reliability and reliability sensitivity analyses are performed by means of the
subset simulation technique for the sake of numerical efficiency. The adaptive
surrogate-based strategy for reliability estimation is finally involved into a
classical gradient-based optimization algorithm in order to solve the RBDO
problem. The kriging surrogates are built in a so-called augmented reliability
space thus making them reusable from one nested RBDO iteration to the other.
The strategy is compared to other approaches available in the literature on
three academic examples in the field of structural mechanics.Comment: 20 pages, 6 figures, 5 tables. Preprint submitted to Springer-Verla
- …