Home  About us  Contact
 

Computational Physics (computational + physics)

Distribution by Scientific Domains
Engineering66%

Selected Abstracts

Comparison of three second-order accurate reconstruction schemes for 2D Euler and Navier,Stokes compressible flows on unstructured grids

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 5 2001
N. P. C. Marques
Abstract This paper reports an intercomparison of three second-order accurate reconstruction schemes to predict 2D steady-state compressible Euler and Navier,Stokes flows on unstructured meshes. The schemes comprise one monotone slope limiter (Barth and Jespersen, A1AA Paper 89-0366, 1989) and two approximately monotone methods: the slope limiter due to Venkatakrishnan and a data-dependent weighting least-squares procedure (Gooch, Journal of Computational Physics, 1997; 133:6,17). In addition to the 1D scalar wave problem, comparisons were performed under two inviscid test cases: a supersonic 10° ramp and a supersonic bump; and two viscous laminar compressible flow cases: the Blasius boundary layer and a double-throated nozzle. The data-dependent oscillatory behaviour is found to be dependent on a user-supplied constant. The three schemes are compared in terms of accuracy and computational efficiency. The results show that the data-dependent procedure always returns a numerical steady-state solution, more accurate than the ones returned by the slope limiters. Its use for Navier,Stokes flow calculations is recommended. Copyright © 2001 John Wiley & Sons, Ltd. [source]

Generation of Arbitrary Lagrangian,Eulerian (ALE) velocities, based on monitor functions, for the solution of compressible fluid equations

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10-11 2005
B. V. Wells
Abstract A moving mesh method is outlined based on the use of monitor functions. The method is developed from a weak conservation principle. From this principle a conservation law for the mesh position is derived. Using the Helmholtz decomposition theorem, this conservation law can be converted into an elliptic equation for a mesh velocity potential. The moving mesh method is discretized using standard finite elements. Once the mesh velocities are obtained an arbitrary Lagrangian,Eulerian (ALE) (Journal of Computational Physics 1974; 14:227) fluid solver is used to update the solution on the adaptive mesh. Results are shown for the compressible Euler equations of gas dynamics in one and two spatial dimensions. Two monitor functions are used, the fluid density (which corresponds to a Lagrangian description), and a function which includes the density gradient. A variety of test problems are considered. Copyright © 2005 John Wiley & Sons, Ltd. [source]

Absorbing boundary condition on elliptic boundary for finite element analysis of water wave diffraction by large elongated bodies

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2001
Subrata Kumar Bhattacharyya
Abstract In a domain method of solution of exterior scalar wave equation, the radiation condition needs to be imposed on a truncation boundary of the modelling domain. The Bayliss, Gunzberger and Turkel (BGT) boundary dampers of first- and second-orders, which require a circular cylindrical truncation boundary in the diffraction-radiation problem of water waves, have been particularly successful in this task. However, for an elongated body, an elliptic cylindrical truncation boundary has the potential to reduce the modelling domain and hence the computational effort. Grote and Keller [On non-reflecting boundary conditions. Journal of Computational Physics 1995; 122: 231,243] proposed extension of the first- and second-order BGT dampers for the elliptic radiation boundary and used these conditions to the acoustic scattering by an elliptic scatterer using the finite difference method. In this paper, these conditions are implemented for the problem of diffraction of water waves using the finite element method. Also, it is shown that the proposed extension works well only for head-on wave incidence. To remedy this, two new elliptic dampers are proposed, one for beam-on incidence and the other for general wave incidence. The performance of all the three dampers is studied using a numerical example of diffraction by an elliptic cylinder. Copyright © 2001 John Wiley & Sons, Ltd. [source]

Replica Exchange Light Transport

COMPUTER GRAPHICS FORUM, Issue 8 2009
Shinya Kitaoka
I.3.7 [Computer Graphics]: Three-Dimensional Graphics and Realism; I.3.3 [Computer Graphics]: Picture/Image Generation Abstract We solve the light transport problem by introducing a novel unbiased Monte Carlo algorithm called replica exchange light transport, inspired by the replica exchange Monte Carlo method in the fields of computational physics and statistical information processing. The replica exchange Monte Carlo method is a sampling technique whose operation resembles simulated annealing in optimization algorithms using a set of sampling distributions. We apply it to the solution of light transport integration by extending the probability density function of an integrand of the integration to a set of distributions. That set of distributions is composed of combinations of the path densities of different path generation types: uniform distributions in the integral domain, explicit and implicit paths in light (particle/photon) tracing, indirect paths in bidirectional path tracing, explicit and implicit paths in path tracing, and implicit caustics paths seen through specular surfaces including the delta function in path tracing. The replica-exchange light transport algorithm generates a sequence of path samples from each distribution and samples the simultaneous distribution of those distributions as a stationary distribution by using the Markov chain Monte Carlo method. Then the algorithm combines the obtained path samples from each distribution using multiple importance sampling. We compare the images generated with our algorithm to those generated with bidirectional path tracing and Metropolis light transport based on the primary sample space. Our proposing algorithm has better convergence property than bidirectional path tracing and the Metropolis light transport, and it is easy to implement by extending the Metropolis light transport. [source]

Institute of Physics, CAS: (Adv. Mater.

ADVANCED MATERIALS, Issue 45 2009
45/2009)
Founded in 1928, the Institute of Physics of the Chinese Academy of Sciences has become one of China's leading research institutions. In 2003, the Institute acquired the Beijing National Laboratory for Condensed Matter Physics, one of the first six national laboratories in China. Its current programs focus on condensed matter physics and its theory, optical physics, atomic and molecular physics, soft matter, plasma physics, and computational physics. [source]

On generalized stochastic perturbation-based finite element method

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 1 2006
Marcin Kami
Abstract Generalized nth order stochastic perturbation technique, that can be applied to solve some boundary value or boundary initial problems in computational physics and/or engineering with random parameters is proposed here. This technique is demonstrated in conjunction with the finite element method (FEM) to model 1D linear elastostatics problem with a single random variable. The symbolic computer program is employed to perform computational studies on convergence of the first two probabilistic moments for simple unidirectional tension of a bar. These numerical studies verify the influence of coefficient of variation of the random input and, at the same time, of the perturbation parameter on the first two probabilistic moments of the final solution vector. Copyright © 2005 John Wiley & Sons, Ltd. [source]