474 Matching Results

Search Results

Advanced search parameters have been applied.

NUMERICAL SIMULATION OF VORTEX BREAKDOWN

Description: Vortex breakdown is simulated by a three dimensional Lagrangian method using vortex filaments. The filaments are approximated by vortex elements and their velocity is computed by a Biot-Savart type law of interaction. The numerical calculations show the development of an axisymmetric bubble with a recirculation zone and resemble in many respects the results obtained in the physical experiments on vortex breakdown.
Date: October 1, 1978
Creator: Prete, Vincenza Del
Partner: UNT Libraries Government Documents Department

1-D Van der Waals Foams Heated by Ion Beam Energy Deposition

Description: One dimensional simulations of various initial average density aluminum foams (modeled as slabs of solid metal separated by low density regions) heated by volumetric energy deposition are conducted with a Lagrangian hydrodynamics code using a van der Waals equation of tate (EOS). The resulting behavior is studied to facilitate the design of future warm dense matter (WDM) experiments at LBNL. In the simulations the energy deposition ranges from 10 to 30 kJ/g and from 0.075 to 4.0 ns total pulse length, resulting in temperatures from approximately 1 o 4 eV. We study peak pressures and temperatures in the foams, expansion velocity, and the phase evolution. Five relevant time scales in the problem are identified. Additionally, we present a method for characterizing the level of inhomogeneity in a foam target as it is heated and the time it takes for a foam to homogenize.
Date: December 23, 2009
Creator: Zylstra, A. B.; Barnard, J. J. & More, R. M.
Partner: UNT Libraries Government Documents Department

Species doubling and effective Lagrangians

Description: Coupling gauge fields to the chiral currents from an effective Lagrangian for pseudoscalar mesons naturally gives rise to a species doubling phenomenon similar to that seen with fermionic fields in lattice gauge theory. 17 refs.
Date: September 1, 1996
Creator: Creutz, M. & Tytgat, M.
Partner: UNT Libraries Government Documents Department

Curved mesh generation and mesh refinement using Lagrangian solid mechanics

Description: We propose a method for generating well-shaped curved unstructured meshes using a nonlinear elasticity analogy. The geometry of the domain to be meshed is represented as an elastic solid. The undeformed geometry is the initial mesh of linear triangular or tetrahedral elements. The external loading results from prescribing a boundary displacement to be that of the curved geometry, and the final configuration is determined by solving for the equilibrium configuration. The deformations are represented using piecewise polynomials within each element of the original mesh. When the mesh is sufficiently fine to resolve the solid deformation, this method guarantees non-intersecting elements even for highly distorted or anisotropic initial meshes. We describe the method and the solution procedures, and we show a number of examples of two and three dimensional simplex meshes with curved boundaries. We also demonstrate how to use the technique for local refinement of non-curved meshes in the presence of curved boundaries.
Date: December 31, 2008
Creator: Persson, P.-O. & Peraire, J.
Partner: UNT Libraries Government Documents Department

SPH and Material Failure: Progress Report

Description: Smoothed Particle Hydrodynamics (SPH) is a meshless Lagrangian technique for modeling hydrodynamics, and as such offers some unique advantages when applied to problems of material failure and breakup. The two most important of these advantages are: (1) SPH is Lagrangian and robust--i.e., it is never necessary to advect or remap. Damage models typically involve a number of complex history variables (such as the damage associated with the Lagrangian mass, crack orientations, etc.), and advecting these quantities as is required in a mesh based algorithm is a very challenging problem. (2) SPH allows the Lagrangian points to move about, reconnect, or separate as dictated by the material flow. This naturally allows for the points to move apart as distinct fragments of material form, resulting in gaps or cracks between the fragments. Typically mesh based algorithms represent the ''cracks'' between fragments as zones of failed material, which is quite different than allowing voids devoid of material to form.
Date: April 12, 2005
Creator: Owen, J M
Partner: UNT Libraries Government Documents Department

A METHOD OF TREATING UNSTRUCTURED CONCAVE CELLS IN STAGGERED-GRID LAGRANGIAN HYDRODYNAMICS

Description: A method is proposed for the treatment of concave cells in staggered-grid Lagrangian hydrodynamics. The method is general enough to be applied to two- and three-dimensional unstructured cells. Instead of defining a cell-point as the geometric average of its nodes (a cell-center), the cell-point is that which equalizes the triangular/tetrahedral area/volume in two/three dimensions. Examples are given.
Date: December 1, 2000
Creator: ROUSCULP, C. & BURTON, D.
Partner: UNT Libraries Government Documents Department

Solving the SUSY CP problem with flavor breaking F-terms

Description: Supersymmetric flavor models for the radiative generation of fermion masses offer an alternative way to solve the SUSY-CP problem. We assume that the supersymmetric theory is flavor and CP conserving. CP violating phases are associated to the vacuum expectation values of flavor violating susy-breaking fields. As a consequence, phases appear at tree level only in the soft supersymmetry breaking matrices. Using a U(2) flavor model as an example we show that it is possible to generate radiatively the first and second generation of quark masses and mixings as well as the CKM CP phase. The one-loop supersymmetric contributions to EDMs are automatically zero since all the relevant parameters in the lagrangian are flavor conserving and as a consequence real. The size of the flavor and CP mixing in the susy breaking sector is mostly determined by the fermion mass ratios and CKM elements. We calculate the contributions to {epsilon}, {epsilon}' and to the CP asymmetries in the B decays to {psi}K{sub s}, {phi}K{sub s}, {eta}'K{sub s} and X{sub s}{gamma}. We analyze a case study with maximal predictivity in the fermion sector. For this worst case scenario the measurements of {Delta}m{sub K}, {Delta}m{sub B} and {epsilon} constrain the model requiring extremely heavy squark spectra.
Date: May 11, 2005
Creator: Diaz-Cruz, Lorenzo J. & Ferrandis, Javier
Partner: UNT Libraries Government Documents Department

The Fidelity of Ocean Models With Explicit Eddies (Chapter 17)

Description: Current practices within the oceanographic community have been reviewed with regard to the use of metrics to assess the realism of the upper-ocean circulation, ventilation processes diagnosed by time-evolving mixed layer depth and mode water formation, and eddy heat fluxes in large-scale fine resolution ocean model simulations. We have striven to understand the fidelity of these simulations in the context of their potential use in future fine-resolution coupled climate system studies. A variety of methodologies are used to assess the veracity of the numerical simulations. Sea surface height variability and the location of western boundary current paths from altimetry have been used routinely as basic indicators of fine-resolution model performance. Drifters and floats have also been used to provide pseudo-Eulerian measures of the mean and variability of surface and sub-surface flows, while statistical comparisons of observed and simulated means have been carried out using James tests. Probability density functions have been used to assess the Gaussian nature of the observed and simulated flows. Length and time scales have been calculated in both Eulerian and Lagrangian frameworks from altimetry and drifters, respectively. Concise measures of multiple model performance have been obtained from Taylor diagrams. The time-evolution of the mixed layer depth at monitoring stations has been compared with simulated time series. Finally, eddy heat fluxes are compared to climatological inferences.
Date: August 1, 2007
Creator: McClean, J.; Jayne, S.; Maltrud, M. & Ivanova, D.
Partner: UNT Libraries Government Documents Department

Presto 4.18 user's guide.

Description: Presto is a Lagrangian, three-dimensional explicit, transient dynamics code that is used to analyze solids subjected to large, suddenly applied loads. The code is designed for a parallel computing environment and for problems with large deformations, nonlinear material behavior, and contact. Presto also has a versatile element library that incorporates both continuum elements and structural elements. This user's guide describes the input for Presto that gives users access to all the current functionality in the code. The environment in which Presto is built allows it to be coupled with other engineering analysis codes. Using a concept called scope, the input structure reflects the fact that Presto can be used in a coupled environment. The user's guide describes how scope is implemented from the outermost to the innermost scopes. Within a given scope, the descriptions of input commands are grouped based on functionality of the code. For example, all material input command lines are described in a chapter of the user's guide for all the material models that can be used in Presto.
Date: September 1, 2010
Creator: Spencer, Benjamin Whiting
Partner: UNT Libraries Government Documents Department

Dynamic Data-Driven Event Reconstruction for Atmospheric Releases

Description: Accidental or terrorist releases of hazardous materials into the atmosphere can impact large populations and cause significant loss of life or property damage. Plume predictions have been shown to be extremely valuable in guiding an effective and timely response. The two greatest sources of uncertainty in the prediction of the consequences of hazardous atmospheric releases result from poorly characterized source terms and lack of knowledge about the state of the atmosphere as reflected in the available meteorological data. In this report, we discuss the development of a new event reconstruction methodology that provides probabilistic source term estimates from field measurement data for both accidental and clandestine releases. Accurate plume dispersion prediction requires the following questions to be answered: What was released? When was it released? How much material was released? Where was it released? We have developed a dynamic data-driven event reconstruction capability which couples data and predictive models through Bayesian inference to obtain a solution to this inverse problem. The solution consists of a probability distribution of unknown source term parameters. For consequence assessment, we then use this probability distribution to construct a ''''composite'' forward plume prediction which accounts for the uncertainties in the source term. Since in most cases of practical significance it is impossible to find a closed form solution, Bayesian inference is accomplished by utilizing stochastic sampling methods. This approach takes into consideration both measurement and forward model errors and thus incorporates all the sources of uncertainty in the solution to the inverse problem. Stochastic sampling methods have the additional advantage of being suitable for problems characterized by a non-Gaussian distribution of source term parameters and for cases in which the underlying dynamical system is non-linear. We initially developed a Markov Chain Monte Carlo (MCMC) stochastic methodology and demonstrated its effectiveness by reconstructing a wide range ...
Date: February 22, 2007
Creator: Kosovic, B; Belles, R; Chow, F K; Monache, L D; Dyer, K; Glascoe, L et al.
Partner: UNT Libraries Government Documents Department

Two-way coupling of Presto v2.8 and CTH v8.1.

Description: A loose two-way coupling of SNL's Presto v2.8 and CTH v8.1 analysis code has been developed to support the analysis of explosive loading of structures. Presto is a Lagrangian, three-dimensional explicit, transient dynamics code in the SIERRA mechanics suite for the analysis of structures subjected to impact-like loads. CTH is a hydro code for modeling complex multi-dimensional, multi-material problems that are characterized by large deformations and/or strong shocks. A fundamental assumption in this loose coupling is that the compliance of the structure modeled with Presto is significantly smaller than the compliance of the surrounding medium (e.g. air) modeled with CTH. A current limitation of the coupled code is that the interaction between CTH and thin structures modeled in Presto (e.g. shells) is not supported. Research is in progress to relax this thin-structure limitation.
Date: July 1, 2008
Creator: Bruner, Christopher William Stuteville; Bishop, Joseph E.; Crawford, David A. & Edwards, Harold Carter
Partner: UNT Libraries Government Documents Department

Adagio 4.18 user's guide.

Description: Adagio is a Lagrangian, three-dimensional, implicit code for the analysis of solids and structures. It uses a multi-level iterative solver, which enables it to solve problems with large deformations, nonlinear material behavior, and contact. It also has a versatile library of continuum and structural elements, and an extensive library of material models. Adagio is written for parallel computing environments, and its solvers allow for scalable solutions of very large problems. Adagio uses the SIERRA Framework, which allows for coupling with other SIERRA mechanics codes. This document describes the functionality and input structure for Adagio.
Date: September 1, 2010
Creator: Spencer, Benjamin Whiting
Partner: UNT Libraries Government Documents Department

Numerical Tests and Properties of Waves in Radiating Fluids

Description: We discuss the properties of an analytical solution for waves in radiating fluids, with a view towards its implementation as a quantitative test of radiation hydrodynamics codes. A homogeneous radiating fluid in local thermodynamic equilibrium is periodically driven at the boundary of a one-dimensional domain, and the solution describes the propagation of the waves thus excited. Two modes are excited for a given driving frequency, generally referred to as a radiative acoustic wave and a radiative diffusion wave. While the analytical solution is well known, several features are highlighted here that require care during its numerical implementation. We compare the solution in a wide range of parameter space to a numerical integration with a Lagrangian radiation hydrodynamics code. Our most significant observation is that flux-limited diffusion does not preserve causality for waves on a homogeneous background.
Date: September 3, 2009
Creator: Johnson, B M & Klein, R I
Partner: UNT Libraries Government Documents Department

1-D Van der Waals Foams Heated by Ion Beam Energy Deposition

Description: One dimensional simulations of various initial average density aluminum foams (modeled as slabs of solid metal separated by low density regions) heated by volumetric energy deposition are conducted with a Lagrangian hydrodynamics code using a van der Waals equation of state (EOS). The resulting behavior is studied to facilitate the design of future warm dense matter (WDM) experiments at LBNL. In the simulations the energy deposition ranges from 10 to 30 kJ/g and from 0.075 to 4.0 ns total pulse length, resulting in temperatures from approximately 1 to 4 eV. We study peak pressures and temperatures in the foams, expansion velocity, and the phase evolution. Five relevant time scales in the problem are identified. Additionally, we present a method for characterizing the level of inhomogeneity in a foam target as it is heated and the time it takes for a foam to homogenize.
Date: March 19, 2010
Creator: Zylstra, A; Barnard, J J & More, R M
Partner: UNT Libraries Government Documents Department

Topological Twisted Sigma Model with H-flux Revisited

Description: In this paper we revisit the topological twisted sigma model with H-flux. We explicitly expand and then twist the worldsheet Lagrangian for bi-Hermitian geometry. we show that the resulting action consists of a BRST exact term and pullback terms, which only depend on one of the two generalized complex structures and the B-field. We then discuss the topological feature of the model.
Date: August 18, 2006
Creator: Chuang, Wu-yen
Partner: UNT Libraries Government Documents Department

Modeling pollutant transport using a meshless-lagrangian particle model

Description: A combined meshless-Lagrangian particle transport model is used to predict pollutant transport over irregular terrain. The numerical model for initializing the velocity field is based on a meshless approach utilizing multiquadrics established by Kansa. The Lagrangian particle transport technique uses a random walk procedure to depict the advection and dispersion of pollutants over any type of surface, including street and city canyons
Date: January 1, 2002
Creator: Carrington, D. B. (David B.) & Pepper, D. W. (Darrell W.)
Partner: UNT Libraries Government Documents Department

Sierra/solid mechanics 4.22 user's guide.

Description: Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for the analysis of solids and structures. It provides capabilities for explicit dynamic and implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and an extensive library of material models. The code is written for parallel computing environments, and it allows for scalable solutions of very large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which allows for coupling with other SIERRA mechanics codes. This document describes the functionality and input structure for Sierra/SM.
Date: October 1, 2011
Creator: Thomas, Jesse David
Partner: UNT Libraries Government Documents Department

A pressure relaxation closure model for one-dimensional, two-material Lagrangian hydrodynamics based on the Riemann problem

Description: Despite decades of development, Lagrangian hydrodynamics of strengthfree materials presents numerous open issues, even in one dimension. We focus on the problem of closing a system of equations for a two-material cell under the assumption of a single velocity model. There are several existing models and approaches, each possessing different levels of fidelity to the underlying physics and each exhibiting unique features in the computed solutions. We consider the case in which the change in heat in the constituent materials in the mixed cell is assumed equal. An instantaneous pressure equilibration model for a mixed cell can be cast as four equations in four unknowns, comprised of the updated values of the specific internal energy and the specific volume for each of the two materials in the mixed cell. The unique contribution of our approach is a physics-inspired, geometry-based model in which the updated values of the sub-cell, relaxing-toward-equilibrium constituent pressures are related to a local Riemann problem through an optimization principle. This approach couples the modeling problem of assigning sub-cell pressures to the physics associated with the local, dynamic evolution. We package our approach in the framework of a standard predictor-corrector time integration scheme. We evaluate our model using idealized, two material problems using either ideal-gas or stiffened-gas equations of state and compare these results to those computed with the method of Tipton and with corresponding pure-material calculations.
Date: January 1, 2009
Creator: Kamm, James R & Shashkov, Mikhail J
Partner: UNT Libraries Government Documents Department

Magnetic gauge measurements on the two-stage gun : homogeneous and heterogeneous initiation of high explosives /

Description: One of the reasons for building our gas-driven two-stage gun at Los Alamos was to be able to do shock initiation experiments on high explosives that were too insensitive to initiate with the single-stage gun. In past ARA meetings we have discussed the operation of the gun and the magnetic gauge measurement method. During the past couple of years we have done a number of magnetic gauge experiments on both liquid and solid high explosives. Shock initiation of high explosives depends on the nature of the material - whether it is homogeneous (liquid) or heterogeneous (pressed solid). In the solid explosives, mostly heterogeneous behavior has been measured. In the liquid explosive isopropyl nitrate, classic homogeneous initiation has been measured including the formation of a superdetonation in the shocked liquid. Experiments in both materials are discussed including the particle (mass) velocity profiles at a number of Lagrangian positions in the flow, progress of the shock front as measured by shock tracker gauges, and the position when the reactive wave reaches a detonation condition. The two-stage gun, in conjunction with the multiple magnetic gauging method, has proven very useful for generating new information in initiation experiments. Information from these experiments is of great value to modelers trying to determine the proper reaction rate models to use in simulations of the shock initiation process.
Date: January 1, 2001
Creator: Gustavsen, R. L. (Richard L.); Alcon, R. R. (Robert R.); Davis, L. L. (Lloyd L.) & Sheffield, S. A. (Stephen A.)
Partner: UNT Libraries Government Documents Department

Convergent, compressible Richtmyer-Meshkov experiment-zero order hydrodynamics

Description: Cylindrical experiments were performed on the OMEGA laser at the University of Rochester to study unstable interfaces in single and double shells. For single shells a marker layer of more opaque and higher density material is placed between foam and an outside ablator. The marker is either smooth or with a well defined surface roughness. For double shells an inner cylinder is placed along the outer cylinder axis. The outer cylinder is irradiated directly with 50 laser beams which produces a strong shock (mach number 5-15) that passes through the unstable marker interf'ace creating a Richtmyer-Meshkov (RM) instability. For double shells this shock bounces off the inner cylinder back to the incoming marker layer causing it to decelerate. We present comparisons of the measured smooth marker layer hydrodynamics with computer simulations using both Lagrangian and Eulerian codes.
Date: January 1, 2002
Creator: R.), Magelssen G. R. (Glenn; Scott, J. M. (John M.); Lanier, N. E. (Nicholas E.); Batha, S. H. (Steven H.); Balkey, M. M. (Matthew M.); Barnes, C. W. (Christopher, W.) et al.
Partner: UNT Libraries Government Documents Department

Atmospheric Inverse Estimates of Methane Emissions from Central California

Description: Methane mixing ratios measured at a tall-tower are compared to model predictions to estimate surface emissions of CH{sub 4} in Central California for October-December 2007 using an inverse technique. Predicted CH{sub 4} mixing ratios are calculated based on spatially resolved a priori CH{sub 4} emissions and simulated atmospheric trajectories. The atmospheric trajectories, along with surface footprints, are computed using the Weather Research and Forecast (WRF) coupled to the Stochastic Time-Inverted Lagrangian Transport (STILT) model. An uncertainty analysis is performed to provide quantitative uncertainties in estimated CH{sub 4} emissions. Three inverse model estimates of CH{sub 4} emissions are reported. First, linear regressions of modeled and measured CH{sub 4} mixing ratios obtain slopes of 0.73 {+-} 0.11 and 1.09 {+-} 0.14 using California specific and Edgar 3.2 emission maps respectively, suggesting that actual CH{sub 4} emissions were about 37 {+-} 21% higher than California specific inventory estimates. Second, a Bayesian 'source' analysis suggests that livestock emissions are 63 {+-} 22% higher than the a priori estimates. Third, a Bayesian 'region' analysis is carried out for CH{sub 4} emissions from 13 sub-regions, which shows that inventory CH{sub 4} emissions from the Central Valley are underestimated and uncertainties in CH{sub 4} emissions are reduced for sub-regions near the tower site, yielding best estimates of flux from those regions consistent with 'source' analysis results. The uncertainty reductions for regions near the tower indicate that a regional network of measurements will be necessary to provide accurate estimates of surface CH{sub 4} emissions for multiple regions.
Date: November 21, 2008
Creator: Zhao, Chuanfeng; Andrews, Arlyn E.; Bianco, Laura; Eluszkiewicz, Janusz; Hirsch, Adam; MacDonald, Clinton et al.
Partner: UNT Libraries Government Documents Department