25 Matching Results

Search Results

Advanced search parameters have been applied.


Description: Multimaterial Eulerian and Arbitrary Lagragian-Eulerian (ALE) codes usually use volume fractions of materials to track individual components in mixed cells. Material advection usually is calculated either by interface capturing, where a high-order van Leer-like slope reconstruction technique is applied, or interface tracking, where a normal reconstruction technique is applied. The former approach is more appropriate for gas-like substances, and the latter is ideal for solids and liquids, since it does not smear out material interfaces. A wide range of problems involves both diffuse and sharp interfaces between substances and demands a combination of these techniques. It is possible to treat all substances that can diffuse into each other as a single material and only keep mass fractions of the individual components of the mixture. The material response can be determined based on the assumption of pressure and temperature equilibrium between components of the mixture. Unfortunately, it is extremely difficult to solve the corresponding system of equations. In order to avoid these problems one can introduce an effective gamma and employ the ideal gas approximation to calculate mixture response. This method provides reliable results, is able to compute strong shock waves, and deals with complex equations of state. Results from a number of simulations using this scheme are presented.
Date: September 29, 2005
Creator: Lomov, I & Liu, B
Partner: UNT Libraries Government Documents Department


Description: Geothermal power promises clean, renewable, reliable and potentially widely-available energy, but is limited by high initial capital costs. New drilling technologies are required to make geothermal power financially competitive with other energy sources. One potential solution is offered by Thermal Spallation Drilling (TSD) - a novel drilling technique in which small particles (spalls) are released from the rock surface by rapid heating. While TSD has the potential to improve drilling rates of brittle granitic rocks, the coupled thermomechanical processes involved in TSD are poorly described, making system control and optimization difficult for this drilling technology. In this paper, we discuss results from a new modeling effort investigating thermal spallation drilling. In particular, we describe an explicit model that simulates the grain-scale mechanics of thermal spallation and use this model to examine existing theories concerning spalling mechanisms. We will report how borehole conditions influence spall production, and discuss implications for macro-scale models of drilling systems.
Date: January 19, 2012
Creator: Walsh, S C; Lomov, I & Roberts, J J
Partner: UNT Libraries Government Documents Department

Numerical Modeling of Mixing and Venting from Explosions in Underground Chambers

Description: 2D and 3D numerical simulations were performed to study the dynamic interaction of explosion products in an underground concrete chamber with ambient air, barrels of water, and the surrounding walls and structure. The simulations were carried out with GEODYN, a multi-material, Godunov-based Eulerian code that employs adaptive mesh refinement and runs efficiently on massively parallel computer platforms. Tabular equations of state were used to model materials under shock loading. An appropriate constitutive model was used to describe the concrete. Interfaces between materials were either tracked with a volume-of-fluid method that used high-order reconstruction to specify the interface location and orientation, or a capturing approach was employed with the assumption of local thermal and mechanical equilibrium. A major focus of the study was to estimate the extent of water heating that could be obtained prior to venting of the chamber. Parameters investigated included the chamber layout, energy density in the chamber and the yield-to-water mass ratio. Turbulent mixing was found to be the dominant heat transfer mechanism for heating the water.
Date: June 22, 2005
Creator: Liu, B T; Lomov, I & Glenn, L A
Partner: UNT Libraries Government Documents Department

Simulation of Comet Impact and Survivability of Organic Compounds

Description: Comets have long been proposed as a potential means for the transport of complex organic compounds to early Earth. For this to be a viable mechanism, a significant fraction of organic compounds must survive the high temperatures due to impact. We have undertaken three-dimensional numerical simulations to track the thermodynamic state of a comet during oblique impacts. The comet was modeled as a 1-km water-ice sphere impacting a basalt plane at 11.2 km/s; impact angles of 15{sup o} (from horizontal), 30{sup o}, 45{sup o}, 65{sup o}, and 90{sup o} (normal impact) were examined. The survival of organic cometary material, modeled as water ice for simplicity, was calculated using three criteria: (1) peak temperatures, (2) the thermodynamic phase of H{sub 2}O, and (3) final temperature upon isentropic unloading. For impact angles greater than or equal to 30{sup o}, no organic material is expected to survive the impact. For the 15{sup o} impact, most of the material survives the initial impact and significant fractions (55%, 25%, and 44%, respectively) satisfy each survival criterion at 1 second. Heating due to deceleration, in addition to shock heating, plays a role in the heating of the cometary material for nonnormal impacts. This effect is more noticeable for more oblique impacts, resulting in significant deviations from estimates using scaling of normal impacts. The deceleration heating of the material at late times requires further modeling of breakup and mixing.
Date: July 18, 2007
Creator: Liu, B T; Lomov, I N; Blank, J G & Antoun, T H
Partner: UNT Libraries Government Documents Department

Patch-based Adaptive Mesh Refinement for Multimaterial Hydrodynamics

Description: We present a patch-based direct Eulerian adaptive mesh refinement (AMR) algorithm for modeling real equation-of-state, multimaterial compressible flow with strength. Our approach to AMR uses a hierarchical, structured grid approach first developed by (Berger and Oliger 1984), (Berger and Oliger 1984). The grid structure is dynamic in time and is composed of nested uniform rectangular grids of varying resolution. The integration scheme on the grid hierarchy is a recursive procedure in which the coarse grids are advanced, then the fine grids are advanced multiple steps to reach the same time, and finally the coarse and fine grids are synchronized to remove conservation errors during the separate advances. The methodology presented here is based on a single grid algorithm developed for multimaterial gas dynamics by (Colella et al. 1993), refined by(Greenough et al. 1995), and extended to the solution of solid mechanics problems with significant strength by (Lomov and Rubin 2003). The single grid algorithm uses a second-order Godunov scheme with an approximate single fluid Riemann solver and a volume-of-fluid treatment of material interfaces. The method also uses a non-conservative treatment of the deformation tensor and an acoustic approximation for shear waves in the Riemann solver. This departure from a strict application of the higher-order Godunov methodology to the equation of solid mechanics is justified due to the fact that highly nonlinear behavior of shear stresses is rare. This algorithm is implemented in two codes, Geodyn and Raptor, the latter of which is a coupled rad-hydro code. The present discussion will be solely concerned with hydrodynamics modeling. Results from a number of simulations for flows with and without strength will be presented.
Date: October 18, 2005
Creator: Lomov, I; Pember, R; Greenough, J & Liu, B
Partner: UNT Libraries Government Documents Department


Description: Motion along joints and fractures in the rock has been proposed as one of the sources of near-source shear wave generation, and demonstrating the validity of this hypothesis is a focal scientific objective of the source physics experimental campaign in the Climax Stock granitic outcrop. A modeling effort has been undertaken by LLNL to complement the experimental campaign, and over the long term provide a validated computation capability for the nuclear explosion monitoring community. The approach involves performing the near-field nonlinear modeling with hydrodynamic codes (e.g., GEODYN, GEODYN-L), and the far-field seismic propagation with an elastic wave propagation code (e.g., WPP). the codes will be coupled together to provide a comprehensive source-to-sensor modeling capability. The technical approach involves pre-test predictions of each of the SPE experiments using their state of the art modeling capabilities, followed by code improvements to alleviate deficiencies identified in the pre-test predictions. This spiral development cycle wherein simulations are used to guide experimental design and the data from the experiment used to improve the models is the most effective approach to enable a transition from the descriptive phenomenological models in current use to the predictive, hybrid physics models needed for a science-based modeling capability for nuclear explosion monitoring. The objective of this report is to describe initial results of non-linear motion predictions of the first two SPE shots in the Climax Stock: a 220-lb shot at a depth of 180 ft (SPE No.1), and a 2570-lb shot at a depth of 150 ft (SPE No.2). The simulations were performed using the LLNL ensemble granite model, a model developed to match velocity and displacement attenuation from HARDHAT, PILE DRIVER, and SHOAL, as well as Russian and French nuclear test data in granitic rocks. This model represents the state of the art modeling capabilities as they existed when ...
Date: October 20, 2011
Creator: Antoun, T; Xu, H; Vorobiev, O & Lomov, I
Partner: UNT Libraries Government Documents Department

Simulation of penetration into porous geologic media

Description: We present a computational study on the penetration of steel projectiles into porous geologic materials. The purpose of the study is to extend the range of applicability of a recently developed constitutive model to simulations involving projectile penetration into geologic media. The constitutive model is non-linear, thermodynamically consistent, and properly invariant under superposed rigid body motions. The equations are valid for large deformations and they are hyperelastic in the sense that the stress tensor is related to a derivative of the Helmholtz free energy. The model uses the mathematical structure of plasticity theory to capture the basic features of the mechanical response of geological materials including the effects of bulking, yielding, damage, porous compaction and loading rate on the material response. The new constitutive model has been successfully used to simulate static laboratory tests under a wide range of triaxial loading conditions, and dynamic spherical wave propagation tests in both dry and saturated geologic media.
Date: May 31, 2005
Creator: Vorobiev, O Y; Liu, B T; Lomov, I N & Antoun, T
Partner: UNT Libraries Government Documents Department

Simulation of Hypervelocity Penetration in Limestone

Description: A parameter study was performed to examine the (shock) damage obtained with long-rod and spherical mono-material penetrators impacting two varieties of limestone. In all cases, the impacts were assumed to be normal to the plane of the rock and at zero angle of attack (in the case of the rods). Impact velocities ranged to 15 km/s but most calculations were performed at 4 and 6 km/s and the penetrator mass was fixed at 1000 kg. For unlined underground structures, incipient damage was defined to occur when the peak stress, {sigma}{sub pk}, exceeds 1 kb (100 MPa) and the applied impulse per unit area, I{sub pk}, exceeds 1 ktap (1 kb-{micro}s). Severe damage was assumed to occur when {sigma}{sub pk} exceeds 1 kb and I{sub pk} exceeds 1000 ktaps. Using the latter definition it was found that severe damage in hard, non-porous limestone with spherical impactors extended to a depth of 9 m on-axis for an impact velocity of 4 km/s and 12 m at 6 km/s. Cylinders with length-to-diameter (L/D) ratio of 8.75 achieved depth to severe damage of 23 m and 40 m, respectively under the same conditions. For a limestone medium with 2% initial gas porosity, the latter numbers were reduced to 12 m and 18 m.
Date: May 31, 2005
Creator: Antoun, T; Glenn, L; Walton, O; Goldstein, P; Lomov, I & Liu, B
Partner: UNT Libraries Government Documents Department

Numerical Simulation of Interaction of Hypervelocity Particle Stream with a Target

Description: We present results of direct numerical simulations of impact of hypervelocity particle stream with a target. The stream of interest consists of submillimeter (30-300 micron) brittle ceramic particles. Current supercomputer capabilities make it possible to simulate a realistic size of streams (up to 20 mm in diameter and 500 mm in length) while resolving each particle individually. Such simulations make possible to study the damage of the target from synergistic effects of individual impacts. In our research we fixed the velocity distribution along the axis of the stream (1-4 km/s) and volume fraction of the solid material (1-10%) and study effects of particle size variation, particle and target material properties and surrounding air properties. We ran 3D calibration simulations with up to 10 million individual particles and conducted sensitivity studies with 2D cylindrically symmetric simulations. We used an Eulerian Godunov hydrocode with adaptive mesh refinement. The particles, target material and air are represented with volume-of-fluid approach. Brittle particle and target material has been simulated with pressure-dependent yield strength and Steinberg model has been used for metal targets. Simulations demonstrated penetration depth and a hole diameter similar to experimental observations and can explain the influence of parameters of the stream on the character of the penetration.
Date: July 31, 2007
Creator: Lomov, I; Liu, B; Georgevich, V & Antoun, T
Partner: UNT Libraries Government Documents Department

Numerical Simulation of Gas-Solid Interfaces with Large Deformations

Description: A method of treatment of multimaterial interfaces on Eulerian grids is developed which works well for mixtures of materials with diverse compressibilities and shear moduli. This makes it possible to use this method not only for problems of gas dynamics and solid mechanics but also to model fluid-structure interaction problems.
Date: February 1, 2000
Creator: Vorobiev, O.Y. & Lomov, I.N.
Partner: UNT Libraries Government Documents Department

A strength and damage model for rock under dynamic loading

Description: A thermodynamically consistent strength and failure model for granite under dynamic loading has been developed and evaluated. The model agrees with static strength measurements and describes the effects of pressure hardening, bulking, porous compaction, porous dilation, tensile failure, and failure under compression due to distortional deformations. This paper briefly describes the model and the sensitivity of the simulated response to variations in the model parameters and in the inelastic deformation processes used in different simulations. 1D simulations of an underground explosion in granite are used in the sensitivity study.
Date: June 14, 1999
Creator: Antoun, T H; Glenn, L A; Lomov, I N & Vorobiev, O
Partner: UNT Libraries Government Documents Department

Simulations of an underground explosion in granite

Description: This paper describes the results of a computational study performed to investigate the behavior of granite under shock wave loading conditions. A thermomechanically consistent constitutive model that includes the effects of bulking, yielding, material damage, and porous compaction on the material response was used in the simulations. The model parameters were determined based on experimental data, and the model was then used in a series of one-dimensional simulations of PILE DRIVER, a deeply-buried explosion in a granite formation at the Nevada Test Site. Particle velocity histories, peak velocity and peak displacement as a function of slant range, and the cavity radius obtained from the code simulations compared favorably with PILE DRIVER data.
Date: June 14, 1999
Creator: Antoun, T; Glenn, L A; Lomov, I N & Vorobiev, O Y
Partner: UNT Libraries Government Documents Department

Simulating Perforation Permeability Damage and Cleanup

Description: Completion of cased and cemented wells by shaped charge perforation causes its own damage to the formation, potentially reducing well productivity. In practice it is found that underbalance conditions clean up the damaged zone to some extent, however, the mechanisms of these processes are poorly understood. Most hydrocodes typically used to simulate rock response to shaped charge penetration do not provide permeability estimates. Furthermore, the time scales for formation clean up are potentially much longer than the period of jet penetration. We have developed a simple, yet accurate model for the evolution of porosity and permeability which can easily be incorporated into existing hydrocodes using information from the history of each cell. In addition, we have developed a code that efficiently simulates fines migration during the post-shot surge period using initial conditions taken directly from hydrocode simulations of jet penetration. Results from a one-dimensional model simulation are in excellent agreement with measured fines and permeability distributions. We also present two-dimensional numerical results which qualitatively reproduce experimentally obtained permeability maps for different values of underbalance. Although initial results have been promising, further comparison with experiment is essential to tune the coupling between the hydrocode and fines migration simulator. Currently the permeability model is most appropriate for high permeability sandstones (such as Berea), but with little effort, the model can be extended to other rock types, given sufficient experimental data.
Date: September 1, 2000
Creator: Morris, J.P.; Lomov, I.N. & Glenn, L.A.
Partner: UNT Libraries Government Documents Department

Numerical Integration of Elastoviscoplasticity Model with Stiff Hardening and Softening

Description: The constitutive equations for viscoplasticity typically are stiff differential equations and require special numerical methods to integrate them efficiently. The objective of this paper is to propose a class of rate-dependent viscoplastic constitutive equations which can be integrated by an efficient explicit scheme that includes the first order effect of pressure and plastic strain hardening.
Date: February 1, 2000
Creator: Vorobiev, O.Y.; Lomov, I.N; Glenn, L.A. & Rubin, M.B.
Partner: UNT Libraries Government Documents Department

A Strength and Damage Model for Rock Under Dynamic Loading

Description: A thermodynamically consistent strength and failure model for granite under dynamic loading has been developed and evaluated. The model agrees with static strength measurements and describes the effects of pressure hardening, bulking, shear-enhanced compaction, porous dilation, tensile failure, and failure under compression due to distortional deformations. This paper briefly describes the model and the sensitivity of the simulated response to variations in the model parameters and in the inelastic deformation processes used in different simulations. Numerical simulations of an underground explosion in granite are used in the sensitivity study.
Date: December 1, 1999
Creator: Vorobiev, O.Y.; Antoun, T.H.; Lomov, I.N. & Glenn, L.A.
Partner: UNT Libraries Government Documents Department