Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Engineering Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Computational Flame Diagnostics for Direct Numerical Simulations with Detailed Chemistry of Transportation Fuels

The goal of the proposed research is to create computational flame diagnostics (CFLD) that are rigorous numerical algorithms for systematic detection of critical flame features, such as ignition, extinction, and premixed and non-premixed flamelets, and to understand the underlying physicochemical processes controlling limit flame phenomena, flame stabilization, turbulence-chemistry interactions and pollutant emissions etc. The goal has been accomplished through an integrated effort on mechanism reduction, direct numerical simulations (DNS) of flames at engine conditions and a variety of turbulent flames with transport fuels, computational diagnostics, turbulence modeling, and DNS data mining and data reduction. The computational diagnostics are primarily based on the chemical explosive mode analysis (CEMA) and a recently developed bifurcation analysis using datasets from first-principle simulations of 0-D reactors, 1-D laminar flames, and 2-D and 3-D DNS (collaboration with J.H. Chen and S. Som at Argonne, and C.S. Yoo at UNIST). Non-stiff reduced mechanisms for transportation fuels amenable for 3-D DNS are developed through graph-based methods and timescale analysis. The flame structures, stabilization mechanisms, local ignition and extinction etc., and the rate controlling chemical processes are unambiguously identified through CFLD. CEMA is further employed to segment complex turbulent flames based on the critical flame features, such as premixed reaction fronts, and to enable zone-adaptive turbulent combustion modeling.

42 ENGINEERING↗
DOE OSTI2012

EPISODIC JETS AS THE CENTRAL ENGINE OF GAMMA-RAY BURSTS

Most gamma-ray bursts (GRBs) have erratic light curves, which demand that the GRB central engine launches an episodic outflow. Recent Fermi observations of some GRBs indicate a lack of the thermal photosphere component as predicted by the baryonic fireball model, which suggests a magnetic origin of GRBs. Given that powerful episodic jets have been observed along with continuous jets in other astrophysical black hole systems, here we propose an intrinsically episodic, magnetically dominated jet model for the GRB central engine. Accumulation and eruption of free magnetic energy in the corona of a differentially rotating, turbulent accretion flow around a hyperaccreting black hole lead to ejections of episodic, magnetically dominated plasma blobs. These blobs are accelerated magnetically, collide with each other at large radii, trigger rapid magnetic reconnection and turbulence, efficient particle acceleration, and radiation, and power the observed episodic prompt gamma-ray emission from GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

CFD Predictions of Severe Accident Natural Circulation Flows in a Combustion Engineering PWR

Computational fluid dynamics (CFD) is used to predict natural circulation flows between a reactor vessel upper plenum and a steam generator during specific severe accidents in a pressurized-water reactor (PWR). The predictions are based on a model for a simplified Combustion Engineering (CE) design. This work is part of a larger study to investigate the potential for severe accident induced steam generator (SG) tube failures in PWRs being conducted by the U.S. Nuclear Regulatory Commission Office of Nuclear Regulatory Research. The CFD model builds upon the lessons learned from previous studies of a generic Westinghouse PWR under similar conditions. The CE SG is significantly larger than the Westinghouse design and geometric differences are expected to reduce the mixing and increase the potential thermal challenge to the SG tubes. The CFD predictions are utilized to improve the basis for simplified flow and mixing models applied in one-dimensional severe accident codes. Parameters are determined to support system code modeling of hot leg and SG inlet plenum mixing, hot leg and SG tube bundle flow rates, and tube inlet temperatures. This updated modeling uses boundary conditions obtained from system analysis code predictions of severe accident behavior in a CE PWR. Model refinements include an improved tube bundle representation and improvements to the thermal boundary conditions. Results indicate less mixing and higher tube inlet temperatures for the CE design compared to previous Westinghouse studies. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Bayesian analysis for estimating statistical parameter distributions of elasto-viscoplastic material models

High temperature design methods rely on constitutive models for inelastic deformation and failure typically calibrated against the mean of experimental data without considering the associated scatter. Variability may arise from the experimental data acquisition process, from heat-to-heat material property variations, or both and need to be accurately captured to predict parameter bounds leading to efficient component design. Applying the Bayesian Markov Chain Monte Carlo (MCMC) method to produce statistical models capturing the underlying uncertainty in the experimental data is an area of ongoing research interest. This work varies aspects of the Bayesian MCMC method and explores their effect on the posterior parameter distributions for a uniaxial elasto-viscoplastic damage model using synthetically generated reference data. From our analysis with the uniaxial inelastic model we determine that an informed prior distribution including different types of test conditions results in more accurate posterior parameter distributions. The parameter posterior distributions, however, do not improve when increasing the number of similar experimental data. Additionally, changing the amount of scatter in the data affects the quality of the posterior distributions, especially for the less sensitive model parameters. Moreover, we perform a sensitivity study of the model parameters against the likelihood function prior to the Bayesian analysis. The results of the sensitivity analysis help to determine the reliability of the posterior distributions and reduce the dimensionality of the problem by fixing the insensitive parameters. The comprehensive study described in this work demonstrates how to efficiently apply the Bayesian MCMC methodology to capture parameter uncertainties in high temperature inelastic material models. Quantifying these uncertainties in inelastic models will improve high temperature engineering design practices and lead to safer, more effective component designs.

42 ENGINEERING↗
DOE OSTI2024

Computational analysis of flame initiation, quenching, and re-ignition in a prechamber natural gas engine under varying EGR-dilution levels

The on-road natural-gas (NG) fueled transportation relies on stoichiometric spark-ignition engines for the advantages of simple after-treatment system despite the efficiency penalty relative to lean-burn combustion strategies. Exhaust gas recirculation (EGR) has the potential to reduce this efficiency gap at low to moderate loads without the need for complex lean-exhaust aftertreatment systems. However, EGR dilution leads to reduced combustion stability and increased cycle-to-cycle variability. A promising technology that has the potential to achieve reliable operation under diluted conditions is the prechamber ignition (or turbulent jet ignition) which uses chemically active turbulent jets generated from combustion inside a prechamber to initiate, stabilize and accelerate combustion of the mixture inside the main chamber. The present work focusses on developing a RANS-based CFD approach to accurately reproduce in-cylinder phenomena in a stoichiometric NG prechamber-assisted heavy-duty engine without relying on complex combustion models that account for turbulence-chemistry interactions. This is necessary because reactive prechamber jets at high EGR dilution tend to extinguish while emerging into the main chamber, which is followed by a phase of re-ignition — a phenomenon that conventional G-equation or well-stirred reactor combustion models cannot reproduce. With addition of a damping multiplier to the well-stirred reactor model, the predictions are seen to show good agreement with experimental pressure evolution and combustion images acquired from a single cylinder Cummins N-14 optical diesel engine retrofitted with a prechamber ignition system. Model predictions of local heat release in the flame and temperature evolution inside the flame are used to investigate combustion dynamics in the prechamber and the main chamber. It is seen that the well-stirred reactor model with the inclusion of damping is able to reproduce the temporary reduction in heat release within the flame, which can be considered equivalent to quenching of jets, and the subsequent re-ignition of the flame inside the main chamber. The delay between quenching and re-ignition depends on the amount of dilution, as explained by an illustration of flame evolution in a Borghi diagram.

Prechamber ignition↗
DOE OSTI2022

Overview of Base Model in Parametric Studies Specific to Performance of U-Mo Plates

This paper provides an overview of the base model specifically developed to perform parametric sensitivity studies on the U-10Mo monolithic fuel system. U-Mo monolithic fuels are being considered for the conversion of test reactors into high-performance research reactors that operate using proliferation-resistant, low-enriched uranium (LEU) fuels. These plate-type fuels contain a high-density, low-enrichment fuel sandwiched between zirconium diffusion barriers and encapsulated in aluminum claddings. All U.S. high-performance research reactors have released the designs of their LEU monolithic fuel reactor cores. These designs include nearly 50 distinct fuel plate geometries with different operational parameters. Consequently, a single generic plate geometry representing all the extreme points in this design matrix is unrealistic. To evaluate the performance for various parameters, a set of sensitivity studies was performed. These studies considered various input parameters (i.e., geometric, operational, and material property-related). The results revealed valuable information about plate performance and the sensitivity of this performance to various modeling inputs. To establish a reference state for comparing these result, base model featuring representative irradiation conditions was developed. To capture in-reactor behavior accurately, incorporation of representative constitutive models capable of evolving properties with respect to temperature, irradiation time, and burnup was needed. The behavioral models considered burnup-dependent properties, swelling, creep, and degradation. This paper introduces the base model created for the parametric sensitivity studies. The detailed description of the procedure includes the model geometry, model discretization, thermo-mechanical coupling, material properties and behavioral models. This paper also provides selected results and assesses the performance of the base model.

42 ENGINEERING↗
DOE OSTI2018

Efficiency and Emissions Characteristics of an HCCI Engine Fueled by Primary Reference Fuels

Here, this article investigates the effects of various primary reference fuel (PRF) blends, compression ratios, and intake temperatures on the thermodynamics and performance of homogeneous charge compression ignition (HCCI) combustion in a Cooperative Fuels Research (CFR) engine. Combustion phasing was kept constant at a CA50 phasing of 5° after top dead center (aTDC) and the equivalence ratio was kept constant at 0.3. Meanwhile, the compression ratio varied from 8:1 to 15:1 as the PRF blends ranged from pure n-heptane to nearly pure isooctane. The intake temperature was used to match CA50 phasing. In addition to the experimental results, a GT-Power model was constructed to simulate the experimental engine and the model was validated against the experimental data. The GT-Power model and simulation results were used to help analyze the energy flows and thermodynamic conditions tested in the experiment. The results indicate that an increase of compression ratio causes higher thermal efficiency and fuel conversion efficiency; however, at the same compression ratio, an increase in PRF number results in lower efficiency due to the required increase in intake temperature and the associated decrease in charge density. While the efficiency does increase with compression ratio, the results show that the effect of increased expansion work is partially offset by higher heat transfer losses and lower ratios of specific heats at higher compression ratios. The results indicate that the maximum pressure rise rate (MPRR) in HCCI significantly increases with compression ratio. Combustion efficiency shows a strong trend with peak temperature regardless of the PRF number or compression ratio, indicating that the CO-to-CO2 conversion is independent of the parent fuel chemistry in the case of the PRFs, whereas the unburned hydrocarbon emissions showed the opposite trend, depending mostly on the parent fuel’s autoignition tendency.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2012

Benchmarking a new closed-form thermal analysis technique against a traditional lumped parameter, finite-difference method

A benchmarking effort was conducted to determine the accuracy of a new analytic generic geology thermal repository model developed at LLNL relative to a more traditional, numerical, lumped parameter technique. The fast-running analytical thermal transport model assumes uniform thermal properties throughout a homogenous storage medium. Arrays of time-dependent heat sources are included geometrically as arrays of line segments and points. The solver uses a source-based linear superposition of closed form analytical functions from each contributing point or line to arrive at an estimate of the thermal evolution of a generic geologic repository. Temperature rise throughout the storage medium is computed as a linear superposition of temperature rises. It is modeled using the MathCAD mathematical engine and is parameterized to allow myriad gridded repository geometries and geologic characteristics [4]. It was anticipated that the accuracy and utility of the temperature field calculated with the LLNL analytical model would provide an accurate 'birds-eye' view in regions that are many tunnel radii away from actual storage units; i.e., at distances where tunnels and individual storage units could realistically be approximated as physical lines or points. However, geometrically explicit storage units, waste packages, tunnel walls and close-in rock are not included in the MathCAD model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Reduced-order modeling through machine learning and graph-theoretic approaches for brittle fracture applications

Typically, thousands of computationally expensive micro-scale simulations of brittle crack propagation are needed to upscale lower length scale phenomena to the macro-continuum scale. Running such a large number of crack propagation simulations presents a significant computational challenge, making reduced-order models (ROMs) attractive for this task. Here, the ultimate goal of this research is to develop ROMs that have sufficient accuracy and low computational cost so that these upscaling simulations can be readily performed. However, constructing ROMs for these complex simulations presents its own challenge. Here, we present and compare four different approaches for reduced-order modeling of brittle crack propagation in geomaterials. These methods rely on machine learning (ML) and graph-theoretic algorithms to approximate key aspects of the brittle crack problem. These methods also incorporate different physics-based assumptions in order to reduce the training requirements while maintaining accurate physics as much as possible. Results from the ROMs are directly compared against a high-fidelity model of brittle crack propagation. Further, the strengths and weaknesses of the ROMs are discussed, and we conclude that combining smart physics-informed feature engineering with highly trainable ML models provides the best performance. The ROMs considered here have computational costs that are orders-of-magnitude less than the cost associated with high-fidelity physical models while maintaining good accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Multi-level computational exploration of advanced combustion engine operating strategies

Advanced combustion engine (ACE) research is typically carried out on single-cylinder research engines. These engines are designed to tightly control fueling and conditions at intake valve closure (IVC) and to precisely measure in-cylinder conditions and emissions. However, to be able to measure and control engine operation so precisely, these research engines typically do not feature intake and exhaust tracts that resemble those in production engines, specifically in regards to turbomachinery, heat exchangers, and exhaust gas recirculation (EGR) systems. For this reason, these research engines are effective for understanding in-cylinder combustion parameters such as heat release rate, burn duration, combustion efficiency, pollutant formation, and exhaust valve opening (EVO) conditions. Here, this paper applies high fidelity simulations to determine the feasibility of achieving a chosen single cylinder engine operating point on a production type homogeneous charge compression ignition (HCCI) engine, using a partial fuel stratification (PFS) strategy. To accomplish this, a Converge 3 dimensional (3D) – computational fluid dynamics (CFD) model of the experimental combustion chamber and intake and exhaust runners was created to simulate the experimental engine. This model was used to simulate an operating point achieved experimentally, as well as to determine the sensitivity of the operating point to variations in intake pressure, intake temperature, injection timing, injected mass, and EGR fraction. The results from these simulations were fed into a 1-dimensional engine simulation created in AVL Boost, featuring production-type intake and exhaust systems, including turbomachinery and heat exchangers necessary to create the required IVC conditions. This full engine simulation was used to assess the cycle efficiency of the engine at the experimental operating condition, and to assess whether changes to this operating point in intake temperature, intake pressure, direct injection timing, or fueling are beneficial to the cycle efficiency and engine-out emissions. In addition, the sensitivity of promising engine operating points to injection timing and injection mass are determined to evaluate the potential stability of these operating points.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Scalable modular dynamic molten salt reactor system model with decay heat

Nodal dynamic models allow rapid reduced-order simulations of complex systems. This paper presents a publicly available system dynamic model of a molten salt reactor system with a topology of the Molten Salt Reactor Experiment (MSRE), in MATLAB/Simulink, building upon a model previously published and validated with MSRE data. It adds scaling of nominal power, dynamic representation of decay heat generation, a generic decay heat removal system, improved documentation, modularity, and test cases. One hour of system transient runs in less than 10 min on a laptop. This model is intended for first order engineering calculations, developing insights into system transient behavior, equipment sizing, parameter sensitivity studies. The model contains an easy to use toolkit which can be customized and expanded to describe any molten salt reactor system. Finally, scenarios involving off-normal transients are presented. Relation between shutdown reactivity insertions, decay heat removal, core recriticality, and safe shutdown conditions are discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

SolarPILOT: A power tower solar field layout and characterization tool

This paper develops and demonstrates a new SolarPower tower Integrated Layout and Optimization Tool (SolarPILOT). The tool uses the analytical flux image Hermite series approximation originally implemented in the DELSOL3 software developed by Sandia National Laboratory in the early 1980s. By applying the analytical model to individual heliostat images rather than large groups or zones of heliostats, SolarPILOT can characterize a wide variety of heliostat field layouts. The individual heliostat modeling approach increases computational expense in comparison with DELSOL3, so SolarPILOT implements a number of improvements to the analytical approximation method to improve model accuracy and computational efficiency. This paper discusses several of these methods, including dynamic heliostat grouping to reduce the expense of intercept factor evaluation, approximation of annual productivity with a subset of time steps throughout the year, polygon clipping to accurately calculate inter-heliostat shadowing and blocking, receiver and tower geometry optimization, and a trigonometric image transformation technique that ensures model accuracy for small heliostats. SolarPILOT also integrates a Monte-Carlo ray tracing engine (SolTrace), providing improved receiver optical modeling capability, a user-friendly front end for geometry definition, and side-by-side validation of the analytical algorithms.

14 SOLAR ENERGY↗
DOE OSTI2016

Dynamic order reduction of thin-film deposition kinetics models: A reaction factorization approach

A set of numerical tools for the analysis and dynamic dimension reduction of chemical vapor and atomic layer deposition (ALD) surface reaction models is developed in this work. The approach is based on a two-step process where in the first, the chemical species surface balance dynamic equations are factored to effectively decouple the (nonlinear) reaction rates, a process that eliminates redundant dynamic modes and that identifies conserved quantities. If successful, the second phase is implemented to factor out redundant dynamic modes when species relatively minor in concentration are omitted; if unsuccessful, the technique points to potential model structural problems. An alumina ALD process is used for an example consisting of 19 reactions and 23 surface and gas-phase species. Using the approach developed, the model is reduced by nineteen modes to a four-dimensional dynamic system without any knowledge of the reaction rate values. Results are interpreted in the context of potential model validation studies.

42 ENGINEERING↗
DOE OSTI2018

High-Throughput Material Characterization using the Virtual Fields Method

Modeling material and component behavior using finite element analysis (FEA) is critical for modern engineering. One key to a credible model is having an accurate material model, with calibrated model parameters, which describes the constitutive relationship between the deformation and the resulting stress in the material. As such, identifying material model parameters is critical to accurate and predictive FEA. Traditional calibration approaches use only global data (e.g. extensometers and resultant force) and simplified geometries to find the parameters. However, the utilization of rapidly maturing full-field characterization techniques (e.g. Digital Image Correlation (DIC)) with inverse techniques (e.g. the Virtual Feilds Method (VFM)) provide a new, novel and improved method for parameter identification. This LDRD tested that idea: in particular, whether more parameters could be identified per test when using full-field data. The research described in this report successfully proves this hypothesis by comparing the VFM results with traditional calibration methods. Important products of the research include: verified VFM codes for identifying model parameters, a new look at parameter covariance in material model parameter estimation, new validation techniques to better utilize full-field measurements, and an exploration of optimized specimen design for improved data richness.

36 MATERIALS SCIENCE↗
DOE OSTI2017

The Institute for Sustained Performance, Energy, and Resilience Final Project Report

Over the course of the five years of the Scientific Discovery through Advanced Computation (SciDAC-3) program, computational scientists working on behalf of the Department of Energy's Office of Science (DOE SC) exploited a new generation of petascale computing resources to make previously inaccessible discoveries in a broad range of disciplines including physics, chemistry, and materials science. The computing systems underpinning this work increased in performance potential from tens to hundreds of PFlop/s, but in the process evolving significantly from architectures in 2011, when SciDAC-3 began. Although Moore's law continues unabated, the end of Dennard scaling has necessitated a fundamental shift in computer architecture focused on power efficiency. To that end, processors are increasingly varied as they strive to satisfy performance, productivity, reliability, and energy efficiency in the face of divergent computational requirements. The diversity among these machines presents a number of challenges to merely porting today's scientific applications, much less achieving good performance. Therefore, at the beginning of SciDAC-3, USC and its collaborators anticipated that both vastly increased scale (e.g., more chips, more cores per chip, wider SIMD) and heterogeneity would exacerbate performance optimization challenges while simultaneously promoting the issues of energy and heterogeneity will exacerbate performance optimization challenges while simultaneously promoting the issues of energy consumption and resilience to the forefront. To ensure that DOE's computational scientists could successfully exploit emerging high performance computing (HPC) systems, DOE SC’s Advanced Scientific Computing Research Program (ASCR) funded the SciDAC-3 Institute for Sustained Performance, Energy, and Resilience (SUPER). Led by USC, this was a broadly based project with expertise in compilers and system tools, performance engineering, energy management, and resilience. SUPER followed the successful model developed in the previous SciDAC-2 Performance Engineering Research Institute (PERI) of leveraging the research investments DOE and others have made and integrating the results to create new capabilities beyond the reach of any one group. To accomplish this, SUPER was organized into multi-institution teams conducting the research activities in performance portability, energy efficiency, exascale resilience, and overall optimization. SUPER researchers also collaborated with other SciDAC-3 institutes, 13 of the 19 Scientific Computation Application Partnerships (SAP), as well as the broader DOE computational sciences community. These collaborations both focused our research on the real challenges facing petascale scientific computing as well as ensured the broad and immediate impact of the results of our research.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Modeling the Physical and Biochemical Influence of Ocean Thermal Energy Conversion Plant Discharges into their Adjacent Waters

This paper describes the modeling work by Makai Ocean Engineering, Inc. to simulate the biochemical effects of of the nutrient-enhanced seawater plumes that are discharged by one or several 100 megawatt OTEC plants. The modeling is needed to properly design OTEC plants that can operate sustainably with acceptably low biological impact.

16 TIDAL AND WAVE POWER↗
DOE OSTI2017

Air Flow Modeling in the Wind Tunnel of the FHWA Aerodynamics Laboratory at Turner-Fairbank Highway Research Center

Computational fluid dynamics (CFD) modeling is widely used in industry for design and in the research community to support, compliment, and extend the scope of experimental studies. Analysis of transportation infrastructure using high performance cluster computing with CFD and structural mechanics software is done at the Transportation Research and Analysis Computing Center (TRACC) at Argonne National Laboratory. These resources, available at TRACC, were used to perform advanced three-dimensional computational simulations of the wind tunnel laboratory at the Turner-Fairbank Highway Research Center (TFHRC). The goals were to verify the CFD model of the laboratory wind tunnel and then to use versions of the model to provide the capability to (1) perform larger parametric series of tests than can be easily done in the laboratory with available budget and time, (2) to extend testing to wind speeds that cannot be achieved in the laboratory, and (3) to run types of tests that are very difficult or impossible to run in the laboratory. Modern CFD software has many physics models and domain meshing options. Models, including the choice of turbulence and other physics models and settings, the computational mesh, and the solver settings, need to be validated against measurements to verify that the results are sufficiently accurate for use in engineering applications. The wind tunnel model was built and tested, by comparing to experimental measurements, to provide a valuable tool to perform these types of studies in the future as a complement and extension to TFHRC’s experimental capabilities. Wind tunnel testing at TFHRC is conducted in a subsonic open-jet wind tunnel with a 1.83 m (6 foot) by 1.83 m (6 foot) cross section. A three component dual force-balance system is used to measure forces acting on tested models, and a three degree of freedom suspension system is used for dynamic response tests. Pictures of the room are shown in Figure 1-1 to Figure 1-4. A detailed CAD geometry and CFD model of the wind tunnel laboratory at TFHRC was built and tested. Results were compared against experimental wind velocity measurements at a large number of locations around the room. This testing included an assessment of the air flow uniformity provided by the tunnel to the test zone and assessment of room geometry effects, such as influence of the proximity the room walls, the non-symmetrical position of the tunnel in the room, and the influence of the room setup on the air flow in the room. This information is useful both for simplifying the computational model and in deciding whether or not moving, or removing, some of the furniture or other movable objects in the room will change the flow in the test zone.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Failure Formulations in Modular Plasticity Models

Computational prediction of ductile failure remains a challenging and important problem as demonstrated by the recent Sandia Fracture Challenges. In addition to emphasizing the complexity of such problems, the variety of solution strategies also highlighted the number of possible approaches to this problem. A common engineering approach for such efforts is to use a failure model in conjunction with element deletion. In the second Sandia Fracture Challenge, for instance, nine of the fourteen teams used some form of element deletion. For such schemes, a critical decision pertains to the selection of the appropriate failure model; of which many may be found in the literature (see the review of Corona and Reedlunn). The variety may also be observed in the aforementioned second Sandia Fracture Challenge in which at least eight different failure criteria are listed for the nine element deletion based approaches. The selection of the appropriate failure model is a difficult challenge depending on the material being considered and such criteria can variously depend on stress state (i.e. triaxiality, Lode angle) and loading conditions (i.e. strain rate, temperature). Separate implementations of each criteria with different plasticity models can be a repetitive and cumbersome process which may limit available models for an engineering analyst. To mitigate this issue, an effort was pursued to flexibly implement failure models in which different failure models could be specified and utilized within the same elastic-plastic constitutive routine by simply changing the input syntax. Similarly, the same models are implemented across a suite of elastic-plastic formulations enabling consistent definitions. As will be discussed later, a specific "modular failure" model is also implemented which allows for the selection or specification of different dependencies depending on the current need. At this stage, this effort is limited to defining failure models; progression/damage evolution in the constitutive model is not treated and left to future efforts.

42 ENGINEERING↗