Search NASA⌕ Search

SEARCH · Search NASA

Results for “Casting Modeling Simulation”

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 181 records · Page 10

IMAGING AN EVENT HORIZON: MITIGATION OF SOURCE VARIABILITY OF SAGITTARIUS A*

The black hole in the center of the Galaxy, associated with the compact source Sagittarius A* (Sgr A*), is predicted to cast a shadow upon the emission of the surrounding plasma flow, which encodes the influence of general relativity (GR) in the strong-field regime. The Event Horizon Telescope (EHT) is a Very Long Baseline Interferometry (VLBI) network with a goal of imaging nearby supermassive black holes (in particular Sgr A* and M87) with angular resolution sufficient to observe strong gravity effects near the event horizon. General relativistic magnetohydrodynamic (GRMHD) simulations show that radio emission from Sgr A* exhibits variability on timescales of minutes, much shorter than the duration of a typical VLBI imaging experiment, which usually takes several hours. A changing source structure during the observations, however, violates one of the basic assumptions needed for aperture synthesis in radio interferometry imaging to work. By simulating realistic EHT observations of a model movie of Sgr A*, we demonstrate that an image of the average quiescent emission, featuring the characteristic black hole shadow and photon ring predicted by GR, can nonetheless be obtained by observing over multiple days and subsequent processing of the visibilities (scaling, averaging, and smoothing) before imaging. Moreover, it is shown that this procedure can be combined with an existing method to mitigate the effects of interstellar scattering. Taken together, these techniques allow the black hole shadow in the Galactic center to be recovered on the reconstructed image.

79 ASTRONOMY AND ASTROPHYSICS↗

U-10Mo/Zr Interface Modeling using a Microstructure-Based FEM Approach

The U-10Mo in low enrichments (LEU) has been identified as the most promising alternative to the current highly enriched uranium (HEU) used in the United States’ fleet of high performance research reactors (USHPRRs). The nominal configuration of the new LEU U-10Mo plate fuel comprises a U-10Mo fuel foil enriched to slightly less than 20% U-235 (0.08” to 0.02” thick), a thin Zr interlayer/diffusion barrier (25 m thick) and a relatively thick outer can of 6061 aluminum. Currently, the Zr interlayer is clad by hot roll bonding. Previous studies and observations revealed a thinning of the zirconium (Zr) layer during this fuel fabrication process, which is not desirable from the fuel performance perspective. Coarse UMo grains, dendritic structures, Mo concentration segregation, carbides, and porosity are present in the as-cast material and can lead to a nonuniform UMo/Zr interface. The purpose of the current work is to investigate the effects of these microstructural parameters on the Zr coating variation. A microstructure-based finite-element method model was used in this work, and a study on the effect of homogenization on the interface between U-10Mo and Zr was conducted. The model uses actual backscattered electron–scanning electron microscopy microstructures, Mo concentrations, and mechanical properties to predict the behavior of a representative volume element under compressive loading during the rolling process. The model successfully predicted the experimentally observed thinning of the Zr layer in the as-cast material. The model also uses results from a homogenization model as an input, and a study on the effect of different levels of homogenization on the interface indicated that homogenization helps decrease this thinning. This model can be considered a predictive tool representing a first step for model integration and an input into a larger fuel fabrication performance model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

3D Morphological Analysis and Synthesis of Industrial Materials Surfaces

We present the original results pertaining to understanding topographical features on real industrial alloy surfaces. Quite often, such materials are subjected to different thermomechanical treatments with the goal of improving their yield strength, ductility, corrosion resistance, and other properties. Operations of casting, rolling, extruding, and stretching metal invariably leave their distinct “signatures” on the processed materials’ surfaces. It is highly desirable to be able to perform their characterization not in terms of the numerous surface roughness and/or waviness parameters, but rather in terms of the distinctive contributions of such operations. Indeed, if such a quantitative characterization were possible, it would be much easier to introduce corrective actions into the appropriate chain of production operations. We provide such possibility with the “morphological analysis and synthesis” techniques described in the paper. The results can be quantified either using surface topography measurements or studying light scattering from such surfaces. This last capability also gives a quick check of the “morphological similarity” of the “real” and modeled surfaces (images), which is important, e.g., in phase-field simulations. All calculations were performed using the MorphoHawk© software developed by the author and colleagues.

36 MATERIALS SCIENCE↗

A modified approach for simulating electronically nonadiabatic dynamics via the generalized quantum master equation

We present a modified approach for simulating electronically nonadiabatic dynamics based on the Nakajima-Zwanzig generalized quantum master equation (GQME). The modified approach utilizes the fact that the Nakajima-Zwanzig formalism does not require casting the overall Hamiltonian in system-bath form, which is arguably neither natural nor convenient in the case of the Hamiltonian that governs nonadiabatic dynamics. Within the modified approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density operator is fully captured by a memory kernel super-operator. A methodology for calculating the memory kernel from projection-free inputs is developed. Simulating the electronic dynamics via the modified approach, with a memory kernel obtained using exact or approximate methods, can be more cost effective and/or lead to more accurate results than direct application of those methods. The modified approach is compared to previously proposed GQME-based approaches, and its robustness and accuracy are demonstrated on a benchmark spin-boson model with a memory kernel which is calculated within the Ehrenfest method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deep Reinforcement Learning Based Volt-VAR Optimization in Smart Distribution Systems

This paper develops a model-free volt-VAR optimization (VVO) algorithm via multi-agent deep reinforcement learning (DRL) in unbalanced distribution systems. This method is novel since we cast the VVO problem in distribution networks to an intelligent deep Q-network (DQN) framework, which avoids solving a specific optimization model directly when facing time-varying operating conditions in the systems. We consider statuses/ratios of switchable capacitors, voltage regulators, and smart inverters installed at distributed generators as the action variables of the agents. A delicately designed reward function guides these agents to interact with the distribution system, in the direction of reinforcing voltage regulation and power loss reduction simultaneously. The forward-backward sweep method for radial three-phase distribution systems provides accurate power flow results within a few iterations to the DRL environment. The proposed method realizes the dual goals for VVO. We test this algorithm on the unbalanced IEEE 13-bus and 123-bus systems. Numerical simulations validate the excellent performance of this method in voltage regulation and power loss reduction.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Technetium and Iodine Getters to Improve Cast Stone Performance

To determine the effectiveness of the various getter materials prior to their solidification in Cast Stone, a series of batch sorption experiments was performed at Pacific Northwest National Laboratory. To quantify the effectiveness of the removal of Tc(VII) and I(I) from solution by getters, the distribution coefficient, K d (mL/g), was calculated. Testing involved placing getter material in contact with spiked waste solutions at a 1:100 solid-to-solution ratio for periods up to 45 days with periodic solution sampling. One Tc getter was also tested at a 1:10 solid-to-solution ratio. Two different solution media, 18.2 MΩ deionized water (DI H 2 O) and a 7.8 M Na LAW simulant, were used in the batch sorption tests. Each test was conducted at room temperature in an anoxic chamber containing N2 with a small amount of H 2 (0.7%) to maintain anoxic conditions. Each getter-solution combination was run in duplicate. Three Tc- and I-doping concentrations were used separately in aliquots of both the 18.2 MΩ DI H2O and a 7.8 M Na LAW waste simulant. The 1× concentration was developed based on Hanford Tank Waste Operations Simulator (HTWOS) model runs to support the River Protection Project System Plan Revision 6. The other two concentrations were 5× and 10× of the HTWOS values. The Tc and I tests were run separately (i.e., the solutions did not contain both solutes). Sampling of the solid-solution mixtures occurred nominally after 0.2, 1, 3, 6, 9, 12, 15 days and ~35 to 45 days. Seven getter materials were tested for Tc and five materials were tested for I. The seven Tc getters were blast furnace slag 1 (BFS1) (northwest source), BFS2 (southeast source), Sn(II)-treated apatite, Sn(II) chloride, nano tin phosphate, KMS (a potassium-metal-sulfide), and tin hydroxapatite. The five iodine getters were layered bismuth hydroxide (LBH), argentite mineral, synthetic argentite, silver-treated carbon, and silver-treated zeolite. The Tc Kd values measured from experiments conducted using the 7.8 M Na LAW simulant (the simulant selected to represent LAW) for the first 15 days for four Tc getters (BFS1, BFS2, Sn(II)-treated apatite, and Sn(II) chloride) show no, to a very small, capacity to remove Tc from the LAW simulant. For the Tc-getter experiments in the 7.8 M LAW simulant, the majority of the effluent samples show very small drops in Tc concentrations for the 35-day compared to the 15-day samplings. However, the Tc concentration in the simulant blanks also dropped slightly during this period, so the effect of the getter contacting LAW simulant at 35 days compared to 15 days is minimal; except that the BFS1 1:10 test shows a slow but steady decrease in Tc concentration in the LAW simulant supernatant from the beginning to the 35 day contact at which point about 20% of the original Tc has been removed from solution. Lastly, the KMS getter gives the highest K d value for Tc at 35 days where K d values have increased to 104 mL/g. When considering the different I getters reacting with the 7.8 M LAW simulant, two getters are much more effective than the others: Ag zeolite and Syn Arg. The other getters have calculated iodide distribution coefficients that show very limited effectiveness in the caustic conditions created by the LAW simulant. These are preliminary results that will need more detailed analyses including both pre- and post-batch sorption getter solid-phase characterization using state-of-the-art instrumentation such as synchrotron X-ray absorption spectroscopy, which can delineate the oxidation state of the Tc and likely iodine species as well as some of the getters key major components, sulfur and iron in the BFS, and tin and sulfur in the tin-bearing and sulfur-bearing getters. This report also describes future experimental studies to be performed to better elucidate the mechanisms controlling the Tc and I sequestration processes in the various getters and leach tests of getter-bearing Cast Stone monoliths.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Dual-Purpose Canister Filling Demonstration Project Progress Report at ORNL, 2023

The US DOE Office of Nuclear Energy is investigating the feasibility of direct disposal of dual-purpose canisters (DPCs) in a hypothetical geological repository to offset the potential requirement to repackage spent nuclear fuel (SNF) from existing DPCs into smaller, disposal-ready canisters. Oak Ridge National Laboratory (ORNL) is currently evaluating the feasibility of filling void space in loaded DPCs with an engineered material to prevent a criticality event caused by groundwater/moderator intrusion. Metal alloys are being investigated as a filler material because of their relatively low viscosities when molten, which may facilitate their injection via an existing drainpipe that runs almost the full length of the DPC. ORNL’s strategy for evaluating filler viability includes simulations and physical demonstrations of filling and casting behavior, as well as evaluations of materials for compatibility in the repository environment. Filling of DPCs in this manner is expected to mitigate the risk associated with a post-closure criticality event during the repository performance assessment time frame (10,000 years or greater). Efforts in this fiscal year focused on (1) destructive analysis of experimental filler castings made in FY 2022, (2) a report outlining a conceptual design of a DPC filling facility (Fortner et al., 2023, M3SF 23OR010305044/ ORNL/SPR-2023/2921, May 31, 2023), (3) a report on affected features, events, and processes (FEPs) due to DPC filling (Price et al., 2023, M3SF-23SN010305093, issuance pending), (4) developing and testing a more practical alloy filler based upon a Sn-Al eutectic, and (5) modeling the heating and cooling dynamics of a DPC subjected to molten metal filling. The preliminary results from each of these tasks support the feasibility of filling DPCs with metal as a strategy against the possibility of criticality in the repository. The FY 2022 casting was found to penetrate even very small orifices in the mold and internal structures. Preliminary testing of the Sn-Al eutectic indicate little interaction of the melt with Zircaloy cladding. Thermal modelling shows that a filled DPC will cool to manageable temperatures within 2-3 days, which is likely manageable in a facility.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Variational preparation of the thermofield double state of the Sachdev-Ye-Kitaev model

Here, we provide an algorithm for preparing the thermofield double (TFD) state of the Sachdev-Ye-Kitaev (SYK) model without the need for an auxiliary bath. Following previous work, the TFD can be cast as the approximate ground state of a Hamiltonian, H TFD . Using variational quantum circuits, we propose and implement a gradient-based algorithm for learning parameters that find this ground state, an application of the variational quantum eigensolver. Concretely, we find shallow quantum circuits that prepare the ground state of H TFD for the q = 4 SYK model for N = 8 Majoranas per side. For N = 12, we achieve a variational energy within 1% of the true ground-state energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The characterization of microstructure and chemistry of transmutation fuels

The goal of the fast reactor development programs is fabrication and certification of a nuclear fuel system that performs all functions of a conventional fuel while destroying recycled actinides. Metal fuels, such as U-Pu-Zr fuels, are considered to be candidates for this application. The physical and mechanical properties, as well as behavior laws and correlations of metal alloy fuel, all of which define the behavior and performance of metal alloy fuels has been investigated in detail. The results indicated that the prospect for fueling a future fast spectrum reactor using metal fuels is excellent. There are several challenges associated with the fuel implementation in reactor environment. One of the limiting factors in lifetime and safety of thermal and fast reactor systems is fuel swelling and resulting fission product transport to cladding. Understanding diffusion kinetics and phases formed between metal fuel and cladding upon exposure to high temperatures typical in reactor environment is critical for ensuring integrity, safety, and performance of the material, advancement of this type of fuel in the future, and securing long-term success of the nuclear fleet. The first step to understanding fuel-cladding interaction is investigation of an as-cast fuel microstructure. Despite the importance of these fuels, detailed description of microstructures in these materials is extremely limited. The majority of published results are based on characterization of as-cast microstructures of these fuels in scanning electron microscopes (SEM). In this work, we present the results from microstructural characterization of metallic Pu-based fuels, which is essential in the validation of current and future atomistic and meso-scale modeling efforts. A firm understanding of microstructural evolution in reactor materials is imperative and will lead to deeper understanding of in-reactor degradation behavior, creating a new generation of radiation-resistant materials and fuels with the potential for significant increases in burn-up.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Improving Cost and Efficiency of the Scalable Solid Oxide Fuel Cells Power System

The objective of this project was to design and develop a 20kW range small-scale solid oxide fuel cells (SOFC) power system for applications such as data centers and commercial buildings. The original plan included a 5,000 hours demonstration and a Techno-Economic Analysis (TEA) which were dropped as part of project termination. The original project plan was to use a stack with a cross-flow cell design which had previously been tested for 500 hours at a community college in Malta, NY. However, it was decided to move to the advanced R-SOFC co-flow cell developed under Department of Energy Award DE-FE0031971. The advanced cell design has the advantage of a larger active area for the same manufacturing footprint which results in fewer required cells for the same stack power, hence a higher volumetric power density (kW/L) and lower cost per kW than the original cross-flow cell design. A full SOFC system Simulink model was developed and calibrated with testing data from a fuel cell stack and BOP (balance of plant) components. The simulation results from the calibrated model showed an acceptable match with the experimental data. A structural analysis conducted for various load scenarios indicated no high stress areas for all spatial directions. Major electrical system components were acquired, built and successfully tested. System sensors were verified and validated against controls. Safety checks, a diagnostic check, PID tuning, and control software commissioning tasks were also conducted. The power electronics prototype was delivered and trial testing completed. Balance of Plant component testing and simulation work was conducted to characterize Reformer-Heat Exchanger heat transfer and backpressure and reformer catalyst methane conversion and product selectivity. Simulations were conducted to design the Anode and Cathode fluid passages and size the air-air and fuel-fuel heat exchangers. A Burner operation map was created from test data and the Anode Gas Recirculation blower was tested to evaluate its durability. The SOFC system used a horizontal style design where components sit directly on a casting with a direct connection to the skid. This design has efficient packaging and a small footprint with approximate dimensions of 750 mm x 700 mm x 1700 mm. An SOFC system was built and successfully tested at the Malta, NY facility The system for over 500 hours under load of which over 300 hours was at full load of 20 kW.

30 DIRECT ENERGY CONVERSION↗

On nonlocal problems with Neumann boundary conditions: scaling and convergence for nonlocal operators and solutions

Formulations of Neumann-type boundary conditions for boundary value problems in the nonlocal framework are beset with difficulties, some related to the choice of a proper scaling. Here we identify a space-dependent scaling for a nonlocal Neumann operator, for which we prove linear in δ (δ being the radius for the support for the kernel) convergence of the Neumann operator and $\mathcal{O}$(δ 2 ) convergence of solutions to their classical counterparts. The pointwise-like convergence of the nonlocal normal operator is cast as a new type of two-scale operator-point convergence, which we call condensated convergence . The results hold for general integrable kernels, a setting which is favored in numerical simulations. We support this analysis with numerical convergence studies using a piecewise linear discontinuous Galerkin discretization and show an $\mathcal{O}$(δ 2 ) rate of convergence of solutions, also exhibiting an $\mathcal{O}$(h 2 ) convergence, where h is the mesh size.

97 MATHEMATICS AND COMPUTING↗

Mechanical property degradation and microstructural evolution of cast austenitic stainless steels under short-term thermal aging

Mechanical testing and microstructural characterization were performed on short-term thermally aged cast austenitic stainless steels (CASS) to understand the severity and mechanisms of thermal-aging degradation experienced during extended operation of light water reactor (LWR) coolant systems. Four CASS materials – CF3, CF3M, CF8, and CF8M – were thermally aged for 1500 hours at 290 °C, 330 °C, 360 °C, and 400 °C. All four alloys experienced insignificant change in strength and ductility properties but a significant reduction in absorbed impact energy. The primary microstructural and compositional changes during thermal aging were spinodal decomposition of the δ-ferrite into α/ α`, precipitation of G-phase in the δ-ferrite, segregation of solute to the austenite/ ferrite interphase boundary, and growth of M23C6 carbides on the austenite/ferrite interphase boundary. These changes were shown to be highly dependent on chemical composition, particularly the concentration of C and Mo, and aging temperature. A comprehensive model is being developed to correlate the microstructural evolution with mechanical behavior and simulation for predictive evaluations of LWR coolant system components.

36 MATERIALS SCIENCE↗

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

Modeling of inertial deposition in scaled models of rat and human nasal airways: Towards in vitro regional dosimetry in small animals

Rodents are routinely used in inhalation toxicology tests as human surrogates. However, in vitro dosimetry tests in rodent casts are still scarce due to small rodent airways and in vitro tests to quantify sub-regional dosimetry are still impractical. We hypothesized that for inertial particles whose deposition is dominated by airflow convection (Reynolds number) and particle inertia (Stokes number), the deposition should be similar among airway replicas of different scales if their Reynolds and Stokes numbers are kept the same. In this study, we aimed to (1) numerically test the hypothesis in three airway geometries: a USP induction port, a human nose model, and a Sprague-Dawley rat nose model, and (2) find the range of applicability of this hypothesis. Five variants of the USP and human nose models and three variants of the rat nose model were tested. Inhalation rates and particle sizes were scaled to match the Reynolds number and Stokes numbers. A low-Reynolds-number k–ω model was used to resolve the airflow and a Lagrangian tracking algorithm was used to simulate the particle transport and deposition. Statistical analysis of predicted doses was conducted using ANOVA. For normal inhalation rates and particle dia- meters ranging from 0.5 to 24 mm, the deposition differences between the life-size and scaled models are insignificant for all airway geometries considered (i.e., human nose, USP, and rat nose). Furthermore, the deposition patterns and exit particle profiles also look similar among scaled models. However, deposition rates and patterns start to deviate if inhalation rates are too low, or particle sizes are too large. For the rat nose, the threshold velocity was found to be 0.71 m/s and the threshold Froude number to be 50. Results of this study provide a theoretical foundation for sub-regional in vitro dosimetry tests in small animals and for interpretation of data from inter-species or intra-species with varying body sizes.

54 ENVIRONMENTAL SCIENCES↗

Constitutive Modeling of the Flow Stress of GCr15 Continuous Casting Bloom in the Heavy Reduction Process

According to the calculation results of a 3D thermomechanical-coupled finite-element (FE) model of GCr15 bearing steel bloom during a heavy reduction (HR) process, the variation ranges in the strain rate and strain under HR were described. In addition, the hot deformation behavior of the GCr15 bearing steel was studied over the temperature range from 1023 K to 1573 K (750 °C to 1300 °C) with strain rates of 0.001, 0.01, and 0.1 s{sup −1} in single-pass thermosimulation compression experiments. To ensure the accuracy of the constitutive model, the temperature range was divided into two temperature intervals according to the fully austenitic temperature of GCr15 steel [1173 K (900 °C)]. Two sets of material parameters for the constitutive model were derived based on the true stress–strain curves of the two temperature intervals. A flow stress constitutive model was established using a revised Arrhenius-type constitutive equation, which considers the relationships among the material parameters and true strain. This equation describes dynamic softening during hot compression processes. Considering the effect of glide and climb on the deformation mechanism, the Arrhenius-type constitutive equation was modified by a physically based approach. This model is the most accurate over the temperatures ranging from 1173 K to 1573 K (900 °C to 1300 °C) under HR deformation conditions (ignoring the range from 1273 K to 1573 K (1000 °C to 1300 °C) with a strain rate of 0.1 s{sup −1}). To ensure the convergence of the FE calculation, an approximated method was used to estimate the flow stress at temperatures greater than 1573 K (1300 °C).

36 MATERIALS SCIENCE↗

Optimized Gating System for Steel Castings

Computer modeling using a commercially available software package was used to explore the optimization of filling systems based on the relatively new concepts of avoiding entrainment of air bubbles and oxides by avoiding surface turbulence. The test casting was based on a pattern for a traditional top poured test bar, whose cross section was a tri-lobed clover-like shape. The study illustrates clearly that the detailed design of the filling system has a major influence on the conditions for defect generation during filling. Traditional steel casting systems using the widely popular assembly of preformed refractory tubes were found to behave poorly. Systems were demonstrated which were capable of delivering highly controlled filling behavior. The latest systems to be developed employed (1) a naturally pressurized filling system and (2) the use of filters placed flush on the top of the runner to act as bubble diverters, together with (3) terminal spin traps. These novel filling systems demonstrated excellent performance in simulation, in agreement with practical experience of the capability (of the trident gate in particular) to produce, for the first time in the history of casting, defect-free castings on a routine basis.

36 MATERIALS SCIENCE↗

Wellbore Microannulus Characterization and Modeling.

Subsurface geologic formations used for extracting resources such as oil and gas can subsequently be used as a storage reservoir for the common greenhouse gas CO 2 , a concept known as Carbon Capture and Storage (CCS). Pre-existing wellbores penetrate the reservoirs where supercritical CO 2 is to be injected. These wellbores can potentially be a pathway for contamination if CO 2 leaks through wellbore flaws to an overlying aquifer or the atmosphere. Characterizing wellbore integrity and providing zonal isolation by repairing these wellbore flaws is of critical importance to the long-term isolation of CO 2 and success of CCS. This research aims to characterize the microannulus region of the cement sheath-steel casing interface in terms of its compressibility and permeability. A mock-up of a wellbore system was used for lab-scale testing. Specimens, consisting of a cement sheath cast on a steel casing with microannuli, were subjected to confining pressures and casing pressures in a pressure vessel that allows simultaneous measurement of gas flow along the axis of the specimen. The flow was interpreted as the hydraulic aperture of the microannuli. Numerical models are used to analyze stress and displacement conditions along the casing-cement interface. These numerical results provide good agreement with closed-form elastic solutions. Numerical models incorporating flaws of varying dimensions along the casing-cement interface were then developed to describe the microannulus region. A joint model is used to describe the hydraulic aperture of the microannulus region, whose mechanical stiffness is altered in response to the imposed stress state across the joint interface. The aperture-stress behavior is based upon laboratory measurements of hydraulic aperture as a function of imposed stress conditions. This investigation found that microannulus permeability can satisfactorily be described by a joint model and that the constitutive model imposed in a numerical simulation can play a significant role in the solution behavior and agreement to experimental data. Recommendations for future work include an application of the joint model with a thermally active large-scale reservoir coupled with pore pressure caused by dynamic CO 2 injection and subsequent microannulus region affects.

03 NATURAL GAS↗

Computational Design of a Novel Medium-Carbon, Low-Alloy Steel Microalloyed with Niobium

The design of a new steel with specific properties is always challenging owing to the complex interactions of many variables. In this work, this challenge is dealt with by combining metallurgical principles with computational thermodynamics and kinetics to design a novel steel composition suitable for thermomechanical processing and induction heat treatment to achieve a hardness level in excess of 600 HV with the potential for good fracture toughness. CALPHAD-based packages for the thermodynamics and kinetics of phase transformations and diffusion, namely Thermo-Calc{sup ®} and JMatPro{sup ®}, have been combined with an interdendritic segregation tool (IDS) to optimize the contents of chromium, molybdenum and niobium in a proposed medium-carbon low-manganese steel composition. Important factors taken into account in the modeling and optimization were hardenability and as-quenched hardness, grain refinement and alloying cost. For further investigations and verification, the designed composition, i.e., in wt.% 0.40C, 0.20Si, 0.25Mn, 0.90Cr, 0.50Mo, was cast with two nominal levels of Nb: 0 and 0.012 wt.%. The results showed that an addition of Nb decreases the austenite grain size during casting and after slab reheating prior to hot rolling. Validation experiments showed that the predicted properties, i.e., hardness, hardenability and level of segregation, for the designed composition were realistic. It is also demonstrated that the applied procedure could be useful in reducing the number of experiments required for developing compositions for other new steels.

36 MATERIALS SCIENCE↗