Search results for "Computation"
showing 10 items of 7362 documents
High Order Extrapolation Techniques for WENO Finite-Difference Schemes Applied to NACA Airfoil Profiles
2017
Finite-difference WENO schemes are capable of approximating accurately and efficiently weak solutions of hyperbolic conservation laws. In this context high order numerical boundary conditions have been proven to increase significantly the resolution of the numerical solutions. In this paper a finite-difference WENO scheme is combined with a high order boundary extrapolation technique at ghost cells to solve problems involving NACA airfoil profiles. The results obtained are comparable with those obtained through other techniques involving unstructured meshes.
Computation of travelling wave solutions of scalar conservation laws with a stiff source term
2003
Abstract In this paper we propose a nonoscillatory numerical technique to compute the travelling wave solution of scalar conservation laws with a stiff source term. This procedure is based on the dynamical behavior described by the associated stationary ODE and it reduces/avoids numerical errors usually encountered with these problems, i.e., spurious oscillations and incorrect wave propagation speed. We combine this treatment with either the first order Lax–Friedrichs scheme or the second order Nessyahu–Tadmor scheme. We have tested several model problems by LeVeque and Yee for which the stiffness coefficient can be increased. We have also tested a problem with a nonlinear flux and a discon…
Approximate Lax–Wendroff discontinuous Galerkin methods for hyperbolic conservation laws
2017
Abstract The Lax–Wendroff time discretization is an alternative method to the popular total variation diminishing Runge–Kutta time discretization of discontinuous Galerkin schemes for the numerical solution of hyperbolic conservation laws. The resulting fully discrete schemes are known as LWDG and RKDG methods, respectively. Although LWDG methods are in general more compact and efficient than RKDG methods of comparable order of accuracy, the formulation of LWDG methods involves the successive computation of exact flux derivatives. This procedure allows one to construct schemes of arbitrary formal order of accuracy in space and time. A new approximation procedure avoids the computation of ex…
Kinematic splitting algorithm for fluid–structure interaction in hemodynamics
2013
Abstract In this paper we study a kinematic splitting algorithm for fluid–structure interaction problems. This algorithm belongs to the class of loosely-coupled fluid–structure interaction schemes. We will present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Fluid flow is described by the conservation laws with nonlinearities in convective and diffusive terms. For simplicity of presentation the structure is modelled by the generalized string equation, but the results presented in the paper may be generalized to more complex structure models. The arbitrary Lagrangian–Eulerian approach is used in order to take…
Riemann solvers in relativistic astrophysics
1999
AbstractOur contribution reviews High Resolution Shock Capturing methods (HRSC) in the field of relativistic hydrodynamics with special emphasis on Riemann solvers. HRSC techniques achieve highly accurate numerical approximations (formally second order or better) in smooth regions of the flow, and capture the motion of unresolved steep gradients without creating spurious oscillations. One objective of our contribution is to show how these techniques have been extended to relativistic hydrodynamics, making it possible to explore some challenging astrophysical scenarios. We will review recent literature concerning the main properties of different special relativistic Riemann solvers, and disc…
Adaptive mesh refinement techniques for high-order shock capturing schemes for multi-dimensional hydrodynamic simulations
2006
The numerical simulation of physical phenomena represented by non-linear hyperbolic systems of conservation laws presents specific difficulties mainly due to the presence of discontinuities in the solution. State of the art methods for the solution of such equations involve high resolution shock capturing schemes, which are able to produce sharp profiles at the discontinuities and high accuracy in smooth regions, together with some kind of grid adaption, which reduces the computational cost by using finer grids near the discontinuities and coarser grids in smooth regions. The combination of both techniques presents intrinsic numerical and programming difficulties. In this work we present a …
A Polynomial Approach to the Piecewise Hyperbolic Method
2003
In this paper, a local (third-order accurate) shock capturing method for hyperbolic conservation laws is presented. The method has been made with the same idea as the PHM method, but with a simpler reconstruction. A comparison with the classic high order methods is discussed.
Constraint qualifications and Lagrange multipliers in nondifferentiable programming problems
1994
In this paper, we present several constraint qualifications, and we show that these conditions guarantee the nonvacuity and the boundedness of the Lagrange multiplier sets for general nondifferentiable programming problems. The relationships with various constraint qualifications are investigated.
Density as a constraint and the separation of internal excitation energy in TDHF
1985
We present a fast and efficient constrained Hartree-Fock iteration scheme which constraints the complete density distribution to remain constant. The scheme is particularly suited to a coordinate- or momentum-space representation. The technique is applied to separate the collective and the internal energy in a propagating TDHF state. We study the behavior of these two energies in an16O+16O collision.
How to apply the Kirkwood–Buff theory to individual species in salt solutions
2013
It is generally assumed that the Kirkwood–Buff (KB) theory cannot be applied to anions and cations individually in a solution, as one cannot simulate this system in an open ensemble due to the electroneutrality constraint. By applying our recently derived KB theory for closed systems, we show that one does have access to single-ion properties in Molecular Dynamics. Our findings are supported by simulations for a model of a salt solution in which particles interact with WCA potentials, as well as for the NaCl/water system using the Particle Mesh Ewald technique for electrostatics.