Home About us Contact
|
Home About us Contact | ||
|
|
||
Benchmark Tests (benchmark + test)
Selected Abstracts2-D/3-D multiply transmitted, converted and reflected arrivals in complex layered media with the modified shortest path methodGEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2009Chao-Ying Bai SUMMARY Grid-cell based schemes for tracing seismic arrivals, such as the finite difference eikonal equation solver or the shortest path method (SPM), are conventionally confined to locating first arrivals only. However, later arrivals are numerous and sometimes of greater amplitude than the first arrivals, making them valuable information, with the potential to be used for precise earthquake location, high-resolution seismic tomography, real-time automatic onset picking and identification of multiple events on seismic exploration data. The purpose of this study is to introduce a modified SPM (MSPM) for tracking multiple arrivals comprising any kind of combination of transmissions, conversions and reflections in complex 2-D/3-D layered media. A practical approach known as the multistage scheme is incorporated into the MSPM to propagate seismic wave fronts from one interface (or subsurface structure for 3-D application) to the next. By treating each layer that the wave front enters as an independent computational domain, one obtains a transmitted and/or converted branch of later arrivals by reinitializing it in the adjacent layer, and a reflected and/or converted branch of later arrivals by reinitializing it in the incident layer. A simple local grid refinement scheme at the layer interface is used to maintain the same accuracy as in the one-stage MSPM application in tracing first arrivals. Benchmark tests against the multistage fast marching method are undertaken to assess the solution accuracy and the computational efficiency. Several examples are presented that demonstrate the viability of the multistage MSPM in highly complex layered media. Even in the presence of velocity variations, such as the Marmousi model, or interfaces exhibiting a relatively high curvature, later arrivals composed of any combination of the transmitted, converted and reflected events are tracked accurately. This is because the multistage MSPM retains the desirable properties of a single-stage MSPM: high computational efficiency and a high accuracy compared with the multistage FMM scheme. [source] Optimal transportation meshfree approximation schemes for fluid and plastic flowsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 12 2010B. Li Abstract We develop an optimal transportation meshfree (OTM) method for simulating general solid and fluid flows, including fluid,structure interaction. The method combines concepts from optimal transportation theory with material-point sampling and max-ent meshfree interpolation. The proposed OTM method generalizes the Benamou,Brenier differential formulation of optimal mass transportation problems to problems including arbitrary geometries and constitutive behavior. The OTM method enforces mass transport and essential boundary conditions exactly and is free from tension instabilities. The OTM method exactly conserves linear and angular momentum and its convergence characteristics are verified in standard benchmark problems. We illustrate the range and scope of the method by means of two examples of application: the bouncing of a gas-filled balloon off a rigid wall; and the classical Taylor-anvil benchmark test extended to the hypervelocity range. Copyright © 2010 John Wiley & Sons, Ltd. [source] A stabilized SPH method for inviscid shallow water flows,INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 2 2005Riadh Ata Abstract In this paper, the smoothed particle hydrodynamics (SPH) method is applied to the solution of shallow water equations. A brief review of the method in its standard form is first described then a variational formulation using SPH interpolation is discussed. A new technique based on the Riemann solver is introduced to improve the stability of the method. This technique leads to better results. The treatment of solid boundary conditions is discussed but remains an open problem for general geometries. The dam-break problem with a flat bed is used as a benchmark test. Copyright © 2004 John Wiley & Sons, Ltd. [source] Development of a convection,diffusion-reaction magnetohydrodynamic solver on non-staggered gridsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2004Tony W. H. Sheu Abstract This paper presents a convection,diffusion-reaction (CDR) model for solving magnetic induction equations and incompressible Navier,Stokes equations. For purposes of increasing the prediction accuracy, the general solution to the one-dimensional constant-coefficient CDR equation is employed. For purposes of extending this discrete formulation to two-dimensional analysis, the alternating direction implicit solution algorithm is applied. Numerical tests that are amenable to analytic solutions were performed in order to validate the proposed scheme. Results show good agreement with the analytic solutions and high rate of convergence. Like many magnetohydrodynamic studies, the Hartmann,Poiseuille problem is considered as a benchmark test to validate the code. Copyright © 2004 John Wiley & Sons, Ltd. [source] New algorithms and an in silico benchmark for computational enzyme designPROTEIN SCIENCE, Issue 12 2006Alexandre Zanghellini Abstract The creation of novel enzymes capable of catalyzing any desired chemical reaction is a grand challenge for computational protein design. Here we describe two new algorithms for enzyme design that employ hashing techniques to allow searching through large numbers of protein scaffolds for optimal catalytic site placement. We also describe an in silico benchmark, based on the recapitulation of the active sites of native enzymes, that allows rapid evaluation and testing of enzyme design methodologies. In the benchmark test, which consists of designing sites for each of 10 different chemical reactions in backbone scaffolds derived from 10 enzymes catalyzing the reactions, the new methods succeed in identifying the native site in the native scaffold and ranking it within the top five designs for six of the 10 reactions. The new methods can be directly applied to the design of new enzymes, and the benchmark provides a powerful in silico test for guiding improvements in computational enzyme design. [source] A global Eta model on quasi-uniform gridsTHE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 623 2007Hai Zhang Abstract The application of quasi-uniform grids in global models of the atmosphere is an attempt to increase the computational efficiency by a more cost-effective exploitation of the computing infrastructure. This paper describes the development of a global version of NCEP's regional, step-coordinate, Eta model on two quasi-uniform grids: cubic and octagonal. The governing equations are expressed in a general curvilinear form, so that the cubic and the octagonal versions of the model share the same code in spite of different mapping of the computational domain. The dynamical core of the derived global Eta model is successfully tested in the benchmark test of Held and Suarez. The model with the step-wise formulation of the terrain and full physics is integrated in a series of tests with real data, and the results are compared both with the analysis and the results of the regional Eta model. Copyright © 2007 Royal Meteorological Society [source] OpenUH: an optimizing, portable OpenMP compilerCONCURRENCY AND COMPUTATION: PRACTICE & EXPERIENCE, Issue 18 2007Chunhua Liao Abstract OpenMP has gained wide popularity as an API for parallel programming on shared memory and distributed shared memory platforms. Despite its broad availability, there remains a need for a portable, robust, open source, optimizing OpenMP compiler for C/C++/Fortran 90, especially for teaching and research, for example into its use on new target architectures, such as SMPs with chip multi-threading, as well as learning how to translate for clusters of SMPs. In this paper, we present our efforts to design and implement such an OpenMP compiler on top of Open64, an open source compiler framework, by extending its existing analysis and optimization and adopting a source-to-source translator approach where a native back end is not available. The compilation strategy we have adopted and the corresponding runtime support are described. The OpenMP validation suite is used to determine the correctness of the translation. The compiler's behavior is evaluated using benchmark tests from the EPCC microbenchmarks and the NAS parallel benchmark. Copyright © 2007 John Wiley & Sons, Ltd. [source] A non-linear triangular curved shell elementINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 4 2004T. Wenzel Abstract The objective of this paper is to present and test a simple triangular finite shell element that uses five degrees of freedom at each node. The element is characterized by three position vectors and three unit directors. It depicts the plane stress state version of the element presented (Comput. Struct. 1989; 32(2):379). The element is of the ANS-type (assumed natural strain (J. Appl. Mech. 1981; 48:587). All strains inside the element contain dot products of the six actual element nodal vectors. The construction of the element also allows non-linear material behaviour. Since an enhancement of the membrane strains by the EAS (enhanced assumed strain method) is not possible inside a three node triangle element, the membrane strains perform poor. But via the DKT (discrete Kirchhoff theory) the three directors reveal an excellent bending behaviour for thin shells. The main concern of this paper is to test, if superimposing the CST (constant srain) with the classic DKT leads to good results in standard benchmark tests. Copyright © 2004 John Wiley & Sons, Ltd. [source] On singularities in the solution of three-dimensional Stokes flow and incompressible elasticity problems with cornersINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 4 2004A. Dimitrov Abstract In this paper, a numerical procedure is presented for the computation of corner singularities in the solution of three-dimensional Stokes flow and incompressible elasticity problems near corners of various shape. For obtaining the order and mode of singularity, a neighbourhood of the singular point is considered with only local boundary conditions. The weak formulation of this problem is approximated using a mixed u, p Galerkin,Petrov finite element method. Additionally, a separation of variables is used to reduce the dimension of the original problem. As a result, the quadratic eigenvalue problem (P+,Q+,2R)d=0 is obtained, where the saddle-point-type matrices P, Q, R are defined explicitly. For a numerical solution of the algebraic eigenvalue problem an iterative technique based on the Arnoldi method in combination with an Uzawa-like scheme is used. This technique needs only one direct matrix factorization as well as few matrix,vector products for finding all eigenvalues in the interval ,,(,) , (,0.5, 1.0), as well as the corresponding eigenvectors. Some benchmark tests show that this technique is robust and very accurate. Problems from practical importance are also analysed, for instance the surface-breaking crack in an incompressible elastic material and the three-dimensional viscous flow of a Newtonian fluid past a trihedral corner. Copyright © 2004 John Wiley & Sons, Ltd. [source] Hybrid-stress solid elements for shell structures based upon a modified variational functionalINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 12 2002K. Y. Sze Abstract In this paper, we start with a modified generalized laminate stiffness matrix that serves as a remedy to resolve the thickness locking and some abnormalities encountered by solid-shell elements in laminate analyses. A modified Hellinger,Reissner functional having displacement and a set of generalized stresses as independent fields is devised. Based upon the functional, eight-node and 18-node hybrid-stress solid-shell elements are proposed. A number of benchmark tests on homogenous and laminated plates/shells are conducted. The accuracy of the elements is promising. Copyright © 2002 John Wiley & Sons, Ltd. [source] The natural volume method (NVM): Presentation and application to shallow water inviscid flowsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 1 2009R. Ata Abstract In this paper a fully Lagrangian formulation is used to simulate 2D shallow water inviscid flows. The natural element method (NEM), which has been used successfully with several solid and fluid mechanics applications, is used to approximate the fluxes over Voronoi cells. This particle-based method has shown huge potential in terms of handling problems involving large deformations. Its main advantage lies in the interpolant character of its shape function and consequently the ease it allows with respect to the imposition of Dirichlet boundary conditions. In this paper, we use the NEM collocationally, and in a Lagrangian kinematic description, in order to simulate shallow water flows that are boundary moving problems. This formulation is ultimately shown to constitute a finite-volume methodology requiring a flux computation on Voronoi cells rather than the standard elements, in a triangular or quadrilateral mesh. St Venant equations are used as the mathematical model. These equations have discontinuous solutions that physically represent the existence of shock waves, meaning that stabilization issues have thus been considered. An artificial viscosity deduced from an analogy with Riemann solvers is introduced to upwind the scheme and therefore stabilize the method. Some inviscid bidimensional flows were used as preliminary benchmark tests, which produced decent results, leading to well-founded hopes for the future of this method in real applications. Copyright © 2008 John Wiley & Sons, Ltd. [source] A further work on multi-phase two-fluid approach for compressible multi-phase flowsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 8 2008Yang-Yao Niu Abstract This paper is to continue our previous work Niu (Int. J. Numer. Meth. Fluids 2001; 36:351,371) on solving a two-fluid model for compressible liquid,gas flows using the AUSMDV scheme. We first propose a pressure,velocity-based diffusion term originally derived from AUSMDV scheme Wada and Liou (SIAM J. Sci. Comput. 1997; 18(3):633,657) to enhance its robustness. The scheme can be applied to gas and liquid fluids universally. We then employ the stratified flow model Chang and Liou (J. Comput. Physics 2007; 225:240,873) for spatial discretization. By defining the fluids in different regions and introducing inter-phasic force on cell boundary, the stratified flow model allows the conservation laws to be applied on each phase, and therefore, it is able to capture fluid discontinuities, such as the fluid interfaces and shock waves, accurately. Several benchmark tests are studied, including the Ransom's Faucet problem, 1D air,water shock tube problems, 2D shock-water column and 2D shock-bubble interaction problems. The results indicate that the incorporation of the new dissipation into AUSM+ -up scheme and the stratified flow model is simple, accurate and robust enough for the compressible multi-phase flows. Copyright © 2008 John Wiley & Sons, Ltd. [source] The shallow flow equations solved on adaptive quadtree gridsINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2001A. G. L. Borthwick Abstract This paper describes an adaptive quadtree grid-based solver of the depth-averaged shallow water equations. The model is designed to approximate flows in complicated large-scale shallow domains while focusing on important smaller-scale localized flow features. Quadtree grids are created automatically by recursive subdivision of a rectangle about discretized boundary, bathymetric or flow-related seeding points. It can be fitted in a fractal-like sense by local grid refinement to any boundary, however distorted, provided absolute convergence to the boundary is not required and a low level of stepped boundary can be tolerated. Grid information is stored as a tree data structure, with a novel indexing system used to link information on the quadtree to a finite volume discretization of the governing equations. As the flow field develops, the grids may be adapted using a parameter based on vorticity and grid cell size. The numerical model is validated using standard benchmark tests, including seiches, Coriolis-induced set-up, jet-forced flow in a circular reservoir, and wetting and drying. Wind-induced flow in the Nichupté Lagoon, México, provides an illustrative example of an application to flow in extremely complicated multi-connected regions. Copyright © 2001 John Wiley & Sons, Ltd. [source] Modeling UHMWPE wear debris generationJOURNAL OF BIOMEDICAL MATERIALS RESEARCH, Issue 2 2007H. Baudriller Abstract It is widely recognized that polyethylene wear debris is one of the main causes of long-term prosthesis loosening. The noxious bioreactivity associated with this debris is determined by its size, shape, and quantity. The aim of this study was to develop a numerical tool that can be used to investigate the primary polyethylene wear mechanisms involved. This model illustrates the formation of varying flow of polyethylene debris with various shapes and sizes caused by elementary mechanical processes. Instead of using the classical continuum mechanics formulation for this purpose, we used a divided materials approach to simulate debris production and release. This approach involves complex nonlinear bulk behaviors, frictional adhesive contact, and characterizes material damage as a loss of adhesion. All the associated models were validated with various benchmark tests. The examples given show the ability of the numerical model to generate debris of various shapes and sizes such as those observed in implant retrieval studies. Most of wear mechanisms such as abrasion, adhesion, and the shearing off of micro-asperities can be described using this approach. Furthermore, it could be applied to study the effects of friction couples, macroscopic geometries, and material processing (e.g. irradiation) on wear. © 2006 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater, 2007 [source] Validation of numerical codes for impact and explosion cratering: Impacts on strengthless and metal targetsMETEORITICS & PLANETARY SCIENCE, Issue 12 2008E. PIERAZZO When properly benchmarked and validated against observation, computer models offer a powerful tool for understanding the mechanics of impact crater formation. This work presents results from the first phase of a project to benchmark and validate shock codes. A variety of 2D and 3D codes were used in this study, from commercial products like AUTODYN, to codes developed within the scientific community like SOVA, SPH, ZEUS-MP, iSALE, and codes developed at U.S. National Laboratories like CTH, SAGE/RAGE, and ALE3D. Benchmark calculations of shock wave propagation in aluminum-on-aluminum impacts were performed to examine the agreement between codes for simple idealized problems. The benchmark simulations show that variability in code results is to be expected due to differences in the underlying solution algorithm of each code, artificial stability parameters, spatial and temporal resolution, and material models. Overall, the inter-code variability in peak shock pressure as a function of distance is around 10 to 20%. In general, if the impactor is resolved by at least 20 cells across its radius, the underestimation of peak shock pressure due to spatial resolution is less than 10%. In addition to the benchmark tests, three validation tests were performed to examine the ability of the codes to reproduce the time evolution of crater radius and depth observed in vertical laboratory impacts in water and two well-characterized aluminum alloys. Results from these calculations are in good agreement with experiments. There appears to be a general tendency of shock physics codes to underestimate the radius of the forming crater. Overall, the discrepancy between the model and experiment results is between 10 and 20%, similar to the inter-code variability. [source] | |||||||||||||