Search NASA⌕ Search

SEARCH · Search NASA

Results for “Contact prediction”

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 361 records · Page 20

A Bayesian Multi-fidelity Neural Network to Predict Nonlinear Frequency Backbone Curves

The use of structural mechanics models during the design process often leads to the development of models of varying fidelity. Often low-fidelity models are efficient to simulate but lack accuracy, while the high-fidelity counterparts are accurate with less efficiency. Here, this paper presents a multi-fidelity surrogate modeling approach that combines the accuracy of a high-fidelity finite element model with the efficiency of a low-fidelity model to train an even faster surrogate model that parameterizes the design space of interest. The objective of these models is to predict the nonlinear frequency backbone curves of the Tribomechadynamics Research Challenge benchmark structure which exhibits simultaneous nonlinearities from frictional contact and geometric nonlinearity. The surrogate model consists of an ensemble of neural networks that learn the mapping between low and high-fidelity data through nonlinear transformations. Bayesian neural networks are used to assess the surrogate model's uncertainty. Once trained, the multi-fidelity neural network is used to perform sensitivity analysis to assess the influence of the design parameters on the predicted backbone curves. Additionally, Bayesian calibration is performed to update the input parameter distributions to correlate the model parameters to the collection of experimentally measured backbone curves.

42 ENGINEERING↗

Planetary Protection Bioburden Analysis Program

This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware. This work was done by David Shores and Sharon P. Goza of Johnson Space Center; Cheyenne McKeegan, Rick Easley, Janet Way, and Shonn Everett of MEI Technologies; Mark Manning of PTI; and Mark Guerra, Ray Kraesig, and William Leu of Tietronix Software, Inc. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24211-1 Spitzer Telemetry Processing System NASA's Jet Propulsion Laboratory, Pasadena, California The Spitzer Telemetry Processing System (SirtfTlmProc) was designed to address objectives of JPL's Multi-mission Image Processing Lab (MIPL) in processing spacecraft telemetry and distributing the resulting data to the science community. To minimize costs and maximize operability, the software design focused on automated error recovery, performance, and information management. The system processes telemetry from the Spitzer spacecraft and delivers Level 0 products to the Spitzer Science Center. SirtfTlmProc is a unique system with automated error notification and recovery, with a real-time continuous service that can go quiescent after periods of inactivity. The software can process 2 GB of telemetry and deliver Level 0 science products to the end user in four hours. It provides analysis tools so the operator can manage the system and troubleshoot problems. It automates telemetry processing in order to reduce staffing costs. This work was done by Alice Stanboli, Elmain M. Martinez, and James M. McAuley of Caltech for NASA's Jet Propulsion Laboratory. For more information, contact iaoffice @jpl.nasa.gov. This software is available for commercial licensing. Please contact Dan Broderick at Daniel.F. Broderick@jpl.nasa.gov. Refer to NPO-47803. NASA Tech Briefs, September 2013 29 This rapid response computer program predicts Orbiter Wing Leading Edge (WLE) damage caused by ice or foam impact during a Space Shuttle launch (Program "IMPACT2"). The program was developed after the Columbia accident in order to assess quickly WLE damage due to ice, foam, or metal impact (if any) during a Shuttle launch. IMPACT2 simulates an impact event in a few minutes for foam impactors, and in seconds for ice and metal impactors. The damage criterion is derived from results obtained from one sophisticated commercial program, which requires hours to carry out simulations of the same impact events. The program was designed to run much faster than the commercial program with prediction of projectile threshold velocities within 10 to 15% of commercial-program values. The mathematical model involves coupling of Orbiter wing normal modes of vibration to nonlinear or linear springmass models. IMPACT2 solves nonlinear or linear impact problems using classical normal modes of vibration of a target, and nonlinear/ linear time-domain equations for the projectile. Impact loads and stresses developed in the target are computed as functions of time. This model is novel because of its speed of execution. A typical model of foam, or other projectile characterized by material nonlinearities, impacting an RCC panel is executed in minutes instead of hours needed by the commercial programs. Target damage due to impact can be assessed quickly, provided that target vibration modes and allowable stress are known. This work was done by Robert Clark, Jr., Paul Cotter, and Constantine Michalopoulos of The Boeing Company for Johnson Space Center. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24988-1 Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2) Lyndon B. Johnson Space Center, Houston, Texas (3) the MSL statistics for only the accountable assays. Other options on the main menu include a data editing form and utility programs that produce various reports requested by the microbiologists and the project, and tools to generate the groupings for the final analyses. The analyses can be carried out in three ways: Each assay can be treated separately, the assays can be collectively treated for the whole zone as a group, or the assays can be collected in groups designated by the JPL Planetary Protection Manager. The latter approach was used to generate the final report because assays on the same equipment or similar equipment can be assumed to have been exposed to the same environment and cleaning. Thus, the statistics are improved by having a larger population, thereby reducing the standard deviation by the square root of N. For each method mentioned above, three reports are available. The first is a detailed report including all the data. This version was very useful in verifying the calculations. The second is a brief report that is similar to the full detailed report, but does not print out the data. The third is a grand total and summary report in which each assay requires only one line. For the first and second reports, most of the calculations are performed in the report section itself. For the third, all the calculations are performed directly in the query bound to the report. All the numerical results were verified by comparing them with Excel templates, then exporting the data from the Planetary Protection Analysis program to Excel.

Beaudet, Robert A.↗

Attention Network Forecasts Time-to-Failure in Laboratory Shear Experiments

Rocks under stress deform by creep mechanisms that include formation and slip on small-scale internal cracks. Intragranular cracks and slip along grain contacts release energy as elastic waves termed acoustic emissions (AE). AEs are thought to contain predictive information that can be used for fault failure forecasting. Here, we present a method using unsupervised classification and an attention network to forecast labquakes using AE waveform features. Our data were generated in a laboratory setting using a biaxial shearing device with granular fault gouge intended to mimic the conditions of tectonic faults. Here, we analyzed the temporal evolution of AEs generated throughout several hundred laboratory earthquake cycles. We used a Conscience Self-Organizing Map (CSOM) to perform topologically ordered vector quantization based on waveform properties. The resulting map was used to interactively cluster AEs. We examined the clusters over time to identify those with predictive ability. Finally, we used a variety of LSTM and attention-based networks to test the predictive power of the AE clusters. By tracking cumulative waveform features over the seismic cycle, the network is able to forecast the time-to-failure (TTF) of lab earthquakes. Our results show that analyzing the data to isolate predictive signals and using a more sophisticated network architecture are key to robustly forecasting labquakes. In the future, this method could be applied on tectonic faults to monitor earthquakes and augment early warning systems.

58 GEOSCIENCES↗

A particle finite element-based model for droplet spreading analysis

Here, a particle finite element method-based model is proposed to analyze droplet dynamics problems, particularly droplet spreading on solid substrates (wetting). The model uses an updated Lagrangian framework to formulate the governing equations of the liquid. The curvature of the liquid surface is tracked accurately using a deforming boundary mesh. In order to predict the spreading rate of the droplet on the solid substrate and track the corresponding contact angle evolution, dissipative forces at the contact line are included in the formulation in addition to the Navier-slip boundary conditions at the solid–liquid interface. The inclusion of these boundary conditions makes it possible to account for the induced Young’s stress at the contact line and for the viscous dissipation along the solid–liquid interfacial region. These are found to be essential to obtain a mesh-independent physical solution. The temporal evolution of the contact angle and the contact line velocity of the proposed model are compared with spreading droplets and micro-sessile droplet injection experiments and are shown to be in good agreement.

42 ENGINEERING↗

Evaluating MetGlasTM Solderability with Tin-Silver-Copper and Tin-Silver-Bismuth Solder Alloys

The solderability of MetGlas TM (subsidiary of Hitachi Metals, America Ltd) 2826 MB, a rapidly solidified metallic foil, was evaluated by the meniscus height/wetting force method for tin-silver-copper (SnAgCu) and tin-silver-bismuth (SnAgBi) solders to understand the effects of the extreme non-equilibrium condition of the MetglasTM surface on solderability performance. Of the variables studied here (solder temperature, heat treatment, and solder composition), solder composition had the largest impact on contact angle. Flux and foil composition remained constant throughout; but, these factors would also be predicted to significantly affect solderability. A greater understanding of the manner whereby non-equilibrium cooling affects solderability of these foils will broaden the application of soldering technology of structures fabricated by rapid cooling process (i.e. additively manufactured coatings and parts). Developing a robust database for Pb-free solderability behavior is also necessary, as industry transitions from tin-lead (SnPb) to lead free (Pb-free) solders.

36 MATERIALS SCIENCE↗

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

FY23 Status Report on MOOSE-Based Approaches to Modeling Core Bowing in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY23, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21 and FY22. Several benchmark verification examples were selected from the IAEA verification and validation report with increased number of ducts and more complex contact interaction behavior. These examples involve a full symmetric sector with restraint rings at multiple load pad locations to simulate a limited free-bow restraint system concept, as well as irradiation induced swelling and creep effects in a sector. In addition, improvements to the contact module sideset assignment were assessed and compared with previous MOOSE results to verify the contact behavior. The MOOSE-based results were compared to IAEA benchmark participants’ results. In every case, the MOOSE results agreed with other simulations results for estimating bowing behavior, providing additional verification basis of the tools for this particular physics application. Estimation of contact forces were mostly in agreement, with a few outlier results. A plan was suggested for dealing with the discrepancies with estimating contact force values. In addition, a 1-way coupling demonstration was performed using subchannel analysis code Pronghorn-SC and MOOSE on an ABR-1000-design sodium-cooled fast reactor assembly to evaluate coolant and duct wall temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Accelerating computational fluid dynamics simulation of post-combustion carbon capture modeling with MeshGraphNets

Packed columns are commonly used in post-combustion processes to capture CO 2 emissions by providing enhanced contact area between a CO 2 -laden gas and CO 2 -absorbing solvent. To study and optimize solvent-based post-combustion carbon capture systems (CCSs), computational fluid dynamics (CFD) can be used to model the liquid–gas countercurrent flow hydrodynamics in these columns and derive key determinants of CO 2 -capture efficiency. However, the large design space of these systems hinders the application of CFD for design optimization due to its high computational cost. In contrast, data-driven modeling approaches can produce fast surrogates to study large-scale physics problems. We build our surrogates using MeshGraphNets (MGN), a graph neural network framework that efficiently learns and produces mesh-based simulations. We apply MGN to a random packed column modeled with over 160K graph nodes and a design space consisting of three key input parameters: solvent surface tension, inlet velocity, and contact angle. Our models can adapt to a wide range of these parameters and accurately predict the complex interactions within the system at rates over 1700 times faster than CFD, affirming its practicality in downstream design optimization tasks. This underscores the robustness and versatility of MGN in modeling complex fluid dynamics for large-scale CCS analyses.

97 MATHEMATICS AND COMPUTING↗

Thermal Analysis of Axially Grooved Heat Pipes

The liquid viscous drag and the evaporator film coefficient were determined analytically for a heat pipe with axial rectangular grooves under zero-g conditions. The liquid pressure drop and the evaporator film coefficient for an axially grooved heat pipe were found to be expressed by the functions of the meniscus angle, which makes it possible to predict its performance analytically. The performance near its capillary pumping limit is influenced by the contact angle.

Kamotani, Y.↗

ERBS orbit ascent utilizing continuous low thrust maneuvers

The Earth Radiation Budget Satellite (ERBS) is a National Aeronautics and Space Administration (NASA) spacecraft whose purpose is to investigate the components of the earth's radiation budget. The ERBS was deployed in a 57 degree inclined, 352.2 kilometer altitude circular orbit by the NASA Space Transportation System (STS) on October 5, 1984. The spacecraft then ascended to its 603 kilometer, near-circular mission orbit by a series of continuous low-thrust maneuvers. The ERBS was the first free-flyer mission to rely on continuous low thrust to spiral from one circular orbit to another. Careful planning and monitoring of these maneuvers were essential to their successful execution. Errors in the prediction of the propulsion system performance or in burn duration would result in loss of contact with the spacecraft during a maneuver and could result in a premature end to the maneuver and difficulty in computing a definitive orbit.

Hoge, S. L.↗

Characterizing the wake vortex signature for an active line of sight remote sensor

A recurring phenomenon, described as a wake vortex, develops as an aircraft approaches the runway to land. As the aircraft moves along the runway, each of the wing tips generates a spiraling and expanding cone of air. During the lifetime of this turbulent event, conditions exist over the runway which can be hazardous to following aircraft, particularly when a small aircraft is following a large aircraft. Left to themselves, these twin vortex patterns will converge toward each other near the center of the runway, harmlessly dissipating through interaction with each other or by contact with the ground. Unfortunately, the time necessary to disperse the vortex is often not predictable, and at busy airports can severely impact terminal area productivity. Rudimentary methods of avoidance are in place. Generally, time delays between landing aircraft are based on what is required to protect a small aircraft. Existing ambient wind conditions can complicate the situation. Reliable detection and tracking of a wake vortex hazard is a major technical problem which can significantly impact runway productivity. Landing minimums could be determined on the basis of the actual hazard rather than imposed on the basis of a worst case scenario. This work focuses on using a windfield description of a wake vortex to generate line-of-sight Doppler velocity truth data appropriate to an arbitrarily located active sensor such as a high resolution radar or lidar. The goal is to isolate a range Doppler signature of the vortex phenomenon that can be used to improve detection. Results are presented based on use of a simplified model of a wake vortex pattern. However, it is important to note that the method of analysis can easily be applied to any vortex model used to generate a windfield snapshot. Results involving several scan strategies are shown for a point sensor with a range resolution of 1 to 4 meters. Vortex signatures presented appear to offer potential for detection and tracking.

Heil, Robert Milton↗

Computational Virtual Reality (VR) as a human-computer interface in the operation of telerobotic systems

This presentation focuses on the application of computer graphics or 'virtual reality' (VR) techniques as a human-computer interface tool in the operation of telerobotic systems. VR techniques offer very valuable task realization aids for planning, previewing and predicting robotic actions, operator training, and for visual perception of non-visible events like contact forces in robotic tasks. The utility of computer graphics in telerobotic operation can be significantly enhanced by high-fidelity calibration of virtual reality images to actual TV camera images. This calibration will even permit the creation of artificial (synthetic) views of task scenes for which no TV camera views are available.

Bejczy, Antal K.↗

Bubble Formation and Detachment in Reduced Gravity Under the Influence of Electric Fields

The objective of the study is to investigate the behavior of individual air bubbles injected through an orifice into an electrically insulating liquid under the influence of a static electric field. Both uniform and nonuniform electric field configurations were considered. Bubble formation and detachment were recorded and visualized in reduced gravity (corresponding to gravity levels on Mars, on the Moon as well as microgravity) using a high-speed video camera. Bubble volume, dimensions and contact angle at detachment were measured. In addition to the experimental studies, a simple model, predicting bubble characteristics at detachment was developed. The model, based on thermodynamic considerations, accounts for the level of gravity as well as the magnitude of the uniform electric field. Measured data and model predictions show good agreement and indicate that the level of gravity and the electric field magnitude significantly affect bubble shape, volume and dimensions.

Herman, Cila↗

Residual Gas Effects on Detached Solidification in Microgravity

Our long term goal has been to make detached solidification reproducible, which requires a full understanding of the mechanisms underlying it. Our Moving Meniscus Model of steady-state detachment predicts that it depends strongly on the surface tension of the melt and the advancing contact angle with the ampoule wall. Thus, the objective of the current project was to determine the influence of residual gases on the surface tension and contact angle of molten semiconductors on typical ampoule materials. Our focus was on the influence of oxygen on indium antimonide on clean silica ("quartz"). The research was performed by three chemical engineering graduate students, the third of whom will complete his research in the summer of 2005. Originally, we had planned to use a sealed silica cell containing a zirconia electrochemical element to control the oxygen partial pressure. However, zirconia requires an operating temperature above the 530 C melting point of InSb and is difficult to form a gas-tight seal with silica. Thus, we decided instead to flow an oxygen-containing gas through the cell. A special apparatus was designed, built and perfected. A piece of InSb was placed on a horizontal silica plate in a quartz cell. High purity argon, helium or hydrogen-containing gas is passed continuously through the cell while the oxygen concentration in the effluent gas is measured. The shape of the resulting drop was used to determine contact angle and surface tension of Ga-doped and high purity InSb. Oxygen appeared to decrease the contact angle, and definitely did not increase it. The following section gives the background for the research. Section 2 summarizes the results obtained on Ga-doped InSb with relatively high oxygen concentrations. Section 3 describes recent improvements made to the apparatus and methods of analysis. Section 4 gives recent results for high-purity InSb at low oxygen concentrations. Final results will be obtained only this summer (2005). Each section has its own references.

Wilcox, William R.↗

Bubble Detachment in Variable Gravity Under the Influence of Electric Fields

The objective of the research is to investigate the behavior of individual air bubbles injected through an orifice into an electrically insulating liquid under the influence of a static electric field. Situations were considered with both uniform and nonuniform electric fields. Bubble formation and detachment were visualized in terrestrial gravity as well as for several levels of reduced gravity (lunar, martian and microgravity) using a high-speed video camera. Bubble volume, dimensions and contact angles at detachment were measured. In addition to the experimental studies, a simple model, predicting bubble characteristics at detachment in an initially uniform electric field was developed. The model, based on thermodynamic considerations, accounts for the level of gravity as well as the magnitude of the uniform electric field. The results of the study indicate that the level of gravity and the electric field magnitude significantly affect bubble behavior as well as shape, volume and dimensions.

Herman, Cila↗

Electric Field Effects on an Injected Air Bubble at Detachment in a Low Gravity Environment

The objective of the study is to investigate the behavior of individual air bubbles injected through an orifice into an electrically insulating liquid under the influence of a static and uniform electric field. Bubble formation and detachment were visualized and recorded in microgravity using a high-speed video camera. Bubble volume, dimensions and contact angle at detachment were measured. In addition to the experimental studies, a simple model, predicting bubble characteristics at detachment was developed. The model, based on thermodynamic considerations, accounts for the level of gravity as well as the magnitude of the uniform electric field. Measured data and model predictions show good agreement, and indicate that the level of gravity and the electric field magnitude significantly affect bubble shape, volume and dimensions.

Iacona, Estelle↗

Non-contact velocity measurement instruments and systems, and related methods

Methods for calibration of non-contact velocity measurements and systems for implementing the same are described. Generally, the method comprises inducing a shock wave into a sample at a stress intensity that varies across the sample's elastic limit, which corresponds to the elastic-plastic state transition of the sample. That transition state may be at the sample's Hugoniot elastic limit. The velocity of the sample is measured using a non-contact velocity measurement instrument such as a velocimeter. The measurement may be compared to a predicted velocity or a velocity measurement made by another system to determine calibration parameters.

47 OTHER INSTRUMENTATION↗

Big Hill Geomechanical Analysis Using M-D Viscoplastic Material Model

The integrity of wellbores at the interbed between the caprock and salt is a serious concern in the Big Hill site. For the remediation and life extension of wellbores, more accurate predictions from the global model are needed. The Big Hill global model is improved using the M-D viscoplastic contact surface model and the mesh containing the interbed layer with contact surfaces between the salt and caprock layers, and fault blocks in overburden and caprock layers. The model calibration has been performed based on the cavern volumetric closures obtained from the Caveman calculations. The results agree well from 1991 to the early 2000s. The difference starts to widen after that, it might be because of frequent fluid movement and raw water injection. Therefore, the predictions from this improved model could be used to examine the structural integrity of caverns in Big Hill salt dome.

42 ENGINEERING↗