Results for “geometric engineering” — Page 56

NASA NTRSDate not supplied

Integral Channel Nozzles and Heat Exchangers using Additive Manufacturing Directed Energy Deposition NASA HR-1 Alloy

Heat exchangers for use in propulsion applications are very critical components because they must be efficient, compact and light and often operate with working fluids at extreme temperatures or pressures or both. Various components and systems use heat exchangers such as combustion chambers of gas turbines and internal combustion engines, fuel cells (air supply and thermal management), electric batteries (thermal management), evaporators and recuperators of waste-heat-to-power systems, and rocket engines. Even if the results are more generally applicable, the heat exchangers applications to which this study is more closely related are regeneratively cooled rocket nozzles and chambers, and repressurization systems for the launch vehicles. These components are often thin-walled and contain pressurized fluids, like propellants at cryogenic or elevated temperatures. Given that the environments that these propulsion components must endure are challenging, the manufacturing to meet these specifications often require long lead times due to specialty processes and unique tooling associated with the combined thin-wall integral channel and large-scale structures. Additive manufacturing (AM) offers programmatic advantages for reduction in processing time and cost in addition to various technical advantages, including the possibility to achieve enhanced hardware complexity targeted to superior performance, part consolidation, and the capability of processing of novel alloys. While AM is already being utilized for heat exchanger components in propulsion applications, almost all these AM components are made by means of Laser Powder Bed Fusion (L-PBF). L-PBF allows for fine features but is rather limited with respect to the overall size of the components that can be manufactured. Recent developments are maturing the Laser Powder Directed Energy Deposition (LP-DED) process which may be used, for example, to make integral channel thin-wall regeneratively-cooled rocket nozzles with diameters greater than 1 m. This paper highlights some integral channel heat exchanger demonstrator hardware applications of LP-DED, as well as the characterization of this process in combination with the use of the NASA HR-1 alloy. To properly utilize LP-DED for heat exchanger manufacturing, various aspects are being characterized such as geometry limitations, measurement of surface texture and geometric angled surfaces, surface enhancements for internal channels, and material evaluation. NASA HR-1 (FeNi-Cr) is a high strength hydrogen resistant superalloy developed for use in aerospace applications, such as heat exchangers. Some aspects and considerations about the design of heat exchangers are summarized together with data relevant to LP-DED manufacturing in combination with the NASA HR-1 alloy. Microchannels were successful deposited down to 2.54 mm and 1 mm wall thickness, wall angles of 30°, both with high reproducibility. It was also found that the areal surface roughness is highly dependent on the size of the powder feedstock used for deposition. The characterization of these LP-DED features is critical for fluid flow and heat transfer predictions as it can be exploited to enhance heat transfer at the cost of increased pressure drop.

Additive Manufacturing↗

NASA NTRSDate not supplied

Integral Channel Nozzles and Heat Exchangers using Additive Manufacturing Directed Energy Deposition NASA HR-1 Alloy

Heat exchangers for use in propulsion applications are very critical components because they must be efficient, compact and light and often operate with working fluids at extreme temperatures or pressures or both. Various components and systems use heat exchangers such as combustion chambers of gas turbines and internal combustion engines, fuel cells (air supply and thermal management), electric batteries (thermal management), evaporators and recuperators of waste-heat-to-power systems, and rocket engines. Even if the results are more generally applicable, the heat exchangers applications to which this study is more closely related are regeneratively cooled rocket nozzles and chambers, and repressurization systems for the launch vehicles. These components are often thin-walled and contain pressurized fluids, like propellants at cryogenic or elevated temperatures. Given that the environments that these propulsion components must endure are challenging, the manufacturing to meet these specifications often require long lead times due to specialty processes and unique tooling associated with the combined thin-wall integral channel and large-scale structures. Additive manufacturing (AM) offers programmatic advantages for reduction in processing time and cost in addition to various technical advantages, including the possibility to achieve enhanced hardware complexity targeted to superior performance, part consolidation, and the capability of processing of novel alloys. While AM is already being utilized for heat exchanger components in propulsion applications, almost all these AM components are made by means of Laser Powder Bed Fusion (L-PBF). L-PBF allows for fine features but is rather limited with respect to the overall size of the components that can be manufactured. Recent developments are maturing the Laser Powder Directed Energy Deposition (LP-DED) process which may be used, for example, to make integral channel thin-wall regeneratively-cooled rocket nozzles with diameters greater than 1 m. This paper highlights some integral channel heat exchanger demonstrator hardware applications of LP-DED, as well as the characterization of this process in combination with the use of the NASA HR-1 alloy. To properly utilize LP-DED for heat exchanger manufacturing, various aspects are being characterized such as geometry limitations, measurement of surface texture and geometric angled surfaces, surface enhancements for internal channels, and material evaluation. NASA HR-1 (FeNi-Cr) is a high strength hydrogen resistant superalloy developed for use in aerospace applications, such as heat exchangers. Some aspects and considerations about the design of heat exchangers are summarized together with data relevant to LP-DED manufacturing in combination with the NASA HR-1 alloy. Microchannels were successful deposited down to 2.54 mm and 1 mm wall thickness, wall angles of 30°, both with high reproducibility. It was also found that the areal surface roughness is highly dependent on the size of the powder feedstock used for deposition. The characterization of these LP-DED features is critical for fluid flow and heat transfer predictions as it can be exploited to enhance heat transfer at the cost of increased pressure drop.

additive manufacturing↗

DOE OSTI2020

An experimental investigation of the effectiveness of Ar-CO 2 shielding gas mixture for the wire arc additive process

Wire arc additive manufacturing (AM) is the process by which a large, metallic structure is built layer-by-layer using a welding arc to melt a wire feedstock. A novel opportunity exists to alter the shielding gas composition to fabricate distinct geometrical features without altering the other AM parameters. In this paper, shielding gases with varied concentrations of Argon (Ar) and CO 2 was used to deposit three distinct geometric shapes (walls, infill, and overhang) using a wire-based additive manufacturing system utilizing the gas metal arc welding (GMAW) surface tension transfer (STT) process. Computer-aided design (CAD) models were sliced with a custom-built slicer, and the sliced algorithm was converted into optimal robotic toolpaths. A custom virtual instrument (VI) was built in LabVIEW to compare the temperature profiles on the surface during each deposition process. After each deposition, the geometric features were scanned, and the surface waviness value was evaluated. Tensile and Charpy impact coupons were extracted from the wall geometries in the longitudinal and transverse directions and tested. The results indicated that a higher CO 2 content produced higher melt pool temperatures to an extent, while lower contents of CO 2 resulted in a dimensionally accurate geometry. The data also indicated that the 2%/98% CO 2 /Ar blend produced scatter in tensile strength and the analysis of variance (ANOVA) shows significant difference. However, the intermediate range of CO 2 (5–10%) resulted in uniform tensile properties. Altogether, these results indicate that a 5%/95% CO 2 /Ar blend is the ideal shielding gas for lowering process temperatures and improving mechanical properties in wire arc additive manufacturing using the gas metal arc welding surface tension transfer process. Additionally, varying concentrations of Ar/CO 2 can be used within the same part in order to modify the local properties or process parameters such as strength, toughness, temperature, or dimensional features. This may improve overall manufacturing quality without sacrificing specific properties.

42 ENGINEERING↗

DOE OSTI2020

General Correlations for Laminar Flow Friction Loss and Heat Transfer in Plain Rectangular Plate-Fin Cores

In this paper, new generalized correlations for predicting the average fanning friction factor f and average Nusselt number Nu for laminar flow in plain plate-fin compact cores of rectangular cross section are presented. These are based on extended experimental data, as well as three-dimensional computational simulations, obtained for a broad range of fin density and geometrical attributes. The results indicate that while the fully developed forced convection scales only with the interfin channel cross-sectional ratio α (fin spacing by fin height), the entrance region hydrodynamic and thermal performance is additionally a function of the fin-core length L, flow Reynolds number Re, and fluid Prandtl number Pr. The developing flow and convection is further shown to scale as: (fRe)~(L/d h Re) 1/2 , and Nu ~(L/d h Re) 1/2 Pr 1/3 Φ(α), where f, Re, and Nu are all based on the hydraulic diameter dh of the interfin flow channel. Generalized correlations for both (fRe) and Nu are developed by the corresponding scaling of the forced convection behavior and asymptotic matching of the entrance or developing flow (short fin-core flow length) and the fully developed flow (large fin-core flow length) region performance. Finally, the predictions from these correlations are found to be within ±15% of all available experimental data for air, water, and glycol (0.71 ≤ Pr ≤ 10), and fin cores with 0 < α ≤ 1.

36 MATERIALS SCIENCE↗

NASA NTRS1946

Flight Investigation of Effect of Various Vertical-Tail Modifications on the Directional Stability and Control Characteristics of the P-63A-1 Airplane (AAF No. 42-68889)

Because the results of preliminary flight tests had indicated. the P-63A-1 airplane possessed insufficient directional stability, the NACA and the manufacturer (Bell Aircraft Corporation) suggested three vertical-tail modifications to remedy the deficiencies in the directional characteristics. These modifications included an enlarged vertical tail formed by adding a tip extension to the original vertical tail, a large sharp-edge ventral fin, and a small dorsal fin. The enlarged vertical tail involved only a slight increase in total vertical-tail area from 23.73 to 26.58 square feet but a relatively much larger increase in geometric aspect ratio from 1.24 to 1.73 based on height and area above the horizontal tail. At the request of the Air Material Command, Army Air Forces, flight tests were made to determine the effect of these modifications and of some combinations of these modifications on the directional stability and control characteristics of the airplane, In all, six different vertical-tail. configurations were investigated to determine the lateral and directional oscillation characteristics of the airplane, the sideslip characteristics, the yaw due to ailerons in rudder-fixed rolls from turns and pull-outs, the trim changes due to speed changes; and the trim changes due to power changes. Results of the tests showed that the enlarged vertical tail approximately doubled the directional stability of the airplane and that the pilots considered the directional stability provided by the enlarged vertical tail to be satisfactory. Calculations based on sideslip data obtained at an indicated airspeed of 300 miles per hour showed that the directional stability of the airplane with the original vertical tail corresponded to a value of 0(sub n beta) of -0.00056 whereas for the enlarged vertical tail the estimated va1ue of C(sub n beta) was -0.00130, The ventral fin was found to increase by a moderate amount the directional stability of the airplane with the original vertical tail for smal1 sides1ip angles at low speeds but little consistent change in directional stability was effected by the ventral fin at higher speeds, The effectiveness of the ventral fin was generally much less when used with the enlarged vertical tail than when used with the original vertical tail. The ventral and dorsal fins were found to be very effective in eliminating rudder-force reversals which occurred in low-speed, high-engine-power, sideslipped conditions of flight . Sideslip tests at two altitudes for approximately the sane engine power and indicated airspeed showed that a small decrease in static directional stability occurred with increasing altitude and this decrease in stability was attributed to the increased propeller blade angles required at high altitudes. The variations of rudder pedal force with indicated airspeed using normal rated power and a constant rudder tab setting through the speed range were desirably small for all the configurations tested. The rudder pedal force changed by about 50 pounds for a power change from engine idling power, to normal rated power and this pedal force change was largely independent of airspeed or of vertical-tail configuration for the various configurations tested.

Johnson, Harold I.↗

DOE OSTI2020

Characterization and Scaling of Forced Convective Swirl in Sinusoidal Wavy-Plate-Fin Cores of Compact Heat Exchangers

The characterization and scaling of the thermal-hydraulic performance in wavy plate-fin compact heat exchanger cores, based on the understanding of the physical phenomena and heat transfer enhancement mechanism is delineated. Experimental data are presented for forced convection in air (Pr = 0.71) with flow rates in the range 50 ≤ Re ≤ 4000. A variety of wavy-fin cores that span viable applications, with geometrical attributes described by the cross section aspect ratio α (ratio of fin spacing to height), fin corrugation aspect ratio γ (ratio of 2× corrugation amplitude to wave pitch), and fin spacing ratio ζ (ratio of fin spacing to wave pitch), are considered. To characterize and correlate the vortex-flow mixing in interfin spaces, a Swirl number Sw is introduced from the balance of viscous, inertial and centrifugal forces. It is shown that the laminar, transitional and turbulent flow regimes can be explicitly identified by this Swirl number. Based on the experimental results and extended analysis, new correlations for Fanning friction factor f and Colburn factor j are developed with Sw, α, γ, and ζ as scaling parameters. Requisite expressions are devised by a superposition of both enhancement components due to the corrugated surface area enlargement and induced swirl flow field, and they are combined to cover the laminar, transitional and turbulent regimes by the method of asymptotic matching. The resulting generalized correlations are further shown to predict all available experimental data for f and j factors to within ±20% and ±15%, respectively.

42 ENGINEERING↗

DOE OSTI2025

Modular Hydronic Room Conditioning System (CRADA NFE-24-10120 Final Report)

This project presents the development and validation of high-fidelity numerical models for fin-tube heat exchangers to enable accurate performance prediction and informed design optimization. The modeling framework integrates detailed geometric specifications, thermophysical property data, and system-level constraints to simulate the heat exchanger’s behavior under a range of operating conditions. The model is calibrated using real-world product specifications and validated against experimental data collected from controlled cooling and heating tests. In cooling mode, the model captures the overall trends in capacity and outlet air temperature but tends to underpredict latent effects, especially at lower air flow rates. In heating mode, the simulation consistently overestimates both the thermal capacity and outlet air temperature, indicating the need for refinement in air-side heat transfer assumptions. Despite these deviations, the model provides a solid foundation for optimization, allowing key design variables to be tuned within physical and performance-based constraints. This research advances the ability to simulate, validate, and optimize fin-tube heat exchanger designs with greater confidence and efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Development and Evaluation of a General Drag Model for Gas-Solid Flows via Deep Learning

This project presents the development and evaluation of a general drag model for gas–solid multiphase flows using deep learning techniques. A comprehensive database of more than 4,000 experimental and numerical data points for spherical and non spherical particles was compiled, incorporating geometric features such as sphericity, aspect ratio, and orientation. Several predictive approaches—including traditional em pirical correlations, machine learning, and deep neural networks—were benchmarked, with the proposed Drag Coefficient Correlation-aided Deep Neural Network (DCC DNN) demonstrating superior accuracy. To account for particle–particle interactions, additional drag data were generated using CFD-based simulations of packed and flu idized beds, leading to the development of a retrained model capable of incorporat ing volume fraction effects. Integration of the trained model with the MFiX CFD solver was achieved using FTorch, enabling drag predictions during discrete element method (DEM) simulations. Validation against experimental data for single particles and fluidized beds confirmed the model’s improved predictive ability, particularly for non-spherical geometries. While the model performed strongly under fluidized con ditions, limitations remained in unfluidized regimes, suggesting a need for expanded datasets. Overall, this study demonstrates the feasibility of combining deep learning with physics-informed CFD to improve drag modeling for gas–solid flows, with promis ing implications for scaling multiphase simulations in industrial applications.

42 ENGINEERING↗

DOE OSTI2021

A new procedure for implementing the modified inherent strain method with improved accuracy in predicting both residual stress and deformation for laser powder bed fusion

As a metal additive manufacturing (AM) process, laser powder bed fusion (L-PBF) has been widely used to produce parts with complex geometries. The large thermal gradient caused by the fast, intense, and repeated laser scanning induces significant residual deformation and stress to the as-built parts, which increase manufacturing difficulty and geometrical inaccuracy as a result. The modified inherent strain (MIS) method exploiting multiscale process simulations was developed to simulate residual deformation accurately and efficiently. However, the existing procedure of implementing the MIS method is found to give inaccurate residual stress prediction. Here in this work, a new implementation procedure for the MIS method is proposed to improve the simulation accuracy of residual stress without degrading the residual deformation prediction. The new procedure concerns the application of inherent strains to the part-scale layer-by-layer finite element model to obtain residual stress and deformation field. While the existing implementation of the part-scale MIS model involves only mechanical properties at ambient temperature, the new procedure adds one more solution step employing mechanical properties at an elevated temperature determined from the inherent strain extraction step. Both numerical and experimental studies are conducted to validate the proposed new implementation procedure. It shows that by using the new procedure, the MIS-based simulation can predict both residual stress and deformation of as-built L-PBF metal parts with good accuracy.

36 MATERIALS SCIENCE↗

NASA NTRS2002

SMEX-Lite Modular Solar Array Architecture

For the most part, Goddard solar arrays have been custom designs that are unique to each mission. The solar panel design has been frozen prior to issuing an RFP for their procurement. There has typically been 6-9 months between RFP release and contract award, followed by an additional 24 months for performance of the contract. For Small Explorer (SMEX) missions, with three years between mission definition and launch, this has been a significant problem. The SMEX solar panels have been sufficiently small that the contract performance period has been reduced to 12-15 months. The bulk of this time is used up in the final design definition and fabrication of flight solar cell assemblies. Even so, it has been virtually impossible to have the spacecraft design at a level of maturity sufficient to freeze the solar panel geometry and release the RFP in time to avoid schedule problems with integrating the solar panels to the spacecraft. With that in mind, the SMEX-Lite project team developed a modular architecture for the assembly of solar arrays to greatly reduce the cost and schedule associated with the development of a mission- specific solar array. In the modular architecture, solar cells are fabricated onto small substrate panels. This modular panel (approximately 8.5" x 17" in this case) becomes the building block for constructing solar arrays for multiple missions with varying power requirements and geometrical arrangements. The mechanical framework that holds these modules together as a solar array is the only mission-unique design, changing in size and shape as required for each mission. There are several advantages to this approach. First, the typical solar array development cycle requires a mission unique design, procurement, and qualification including a custom qualification panel. With the modular architecture, a single qualification of the SMEX-Lite modules and the associated mechanical framework in a typical configuration provided a qualification by similarity to multiple missions. It then becomes possible to procure solar array modules in advance of mission definition and respond quickly and inexpensively to a selected mission's unique requirements. The solar array modular architecture allows the procurement of solar array modules before the array geometry has been frozen. This reduces the effect of procurement lead-time on the mission integration and test flow by as much as 50%. Second, by spreading the non-recurring costs over multiple missions, the cost per unit area is also reduced. In the case of the SMEX-Lite procurement, this reduction was by about one third of the cost per unit area compared to previous SMEX mission-unique procurements. Third, the modular architecture greatly facilitates the infusion of new solar cell technologies into flight programs as these technologies become available. New solar cell technologies need only be fabricated onto a standard-sized module to be incorporated into the next available mission. The modular solar array can be flown in a mixed configuration with some new and some standard cell technologies. Since each module has its own wiring terminals, the array can be arranged as desired electrically with little impact to cost and schedule. The solar array modular architecture does impose some additional constraints on systems and subsystem engineers. First, they must work with discrete solar array modules rather than size the array to fit exactly within an available envelope. The array area is constrained to an integer multiple of the module area. Second, the modular design is optimized for space radiation and thermal environments not greatly different from a typical SMEX LEO environment. For example, a mission with a highly elliptical orbit (e.g., Polar, SMEX/FAST) would require thicker coverglasses to protect the solar cells from the more intense radiation environment.

Lyons, John↗

DOE OSTI2024

Universal correlation for falling film evaporation heat transfer coefficients of water and seawater

Horizontal falling film evaporators are widely utilized in desalination industries to increase fresh water supply. However, universal correlations for seawater falling film evaporation under varied operational and geometrical conditions are simply unavailable in open literature. Thus, this study aims to develop such a universal correlation for both plain and enhanced tubes. The detailed heat transfer mechanisms are reviewed, and rational parameters are incorporated to develop the universal correlation. A dataset of 994 data points from 9 sources covering a wide range of conditions was compiled. These conditions include Reynolds numbers from 10 to 7235, heat fluxes from 7.7 to 208 kW/m -2 , saturation temperatures from 278 to 401 K, salinities from 0 to 60 gsalt kg -1 water, tube diameters from 15.8 to 50.8 mm, and liquid feeder height to diameter ratios from 1 to 2.25. Upon analysis, it was found that most of the recommended existing correlations exhibited poor predictive accuracy, as evidenced by larger MADs. The developed correlation in this study demonstrated the best predictive accuracy for the entire dataset, yielding a MAD of 16.8 % and an R 2 of 0.82. Furthermore, the performance of the new correlation was individually assessed across a broader spectrum of operational and design conditions, reflecting the individual conditions’ influences with an overall MAD of 20 %.

42 ENGINEERING↗

DOE OSTI2020

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

42 ENGINEERING↗

DOE OSTI2020

Understanding Microstructural Effects on Dynamic Performance Towards the Development of Shock Metamaterials

With the recent advances in additive manufacturing (AM), long-range periodic lattice assemblies are being developed for vibration and shock mitigation components in aerospace and military applications with unique geometric and topological structures. There has been extensive work in understanding the static properties associated with varying topology of these lattice architectures, but there is almost no understanding of microstructural affects in such structures under high-strain rate dynamic loading conditions. Here we report the shock behavior of lattices with varying intrinsic grain structures achieved by post process annealing. High resolution 316L stainless steel lattices were 3D printed by a laser-powder bed fusion machine and characterized by computed tomography. Subsequent annealing resulted in stress-relieved and recrystallized lattices. Overall the lattices had strong cubic texture aligning with the x-, y- and z-directions of the build with a preference outside the build direction (z). The recrystallized sample had more equiaxed polygonal grains and a layer of BCC ferrite at the surface of the structure approximately 1 grain thick. Upon dynamic compression the as-deposited lattice showed steady compaction behavior while the heat-treated lattices exhibit negative velocity behavior indicative of failure. We attribute this to the stiffer BCC ferrite in the annealed lattices becoming damaged and fragmenting during compression.

36 MATERIALS SCIENCE↗

NASA NTRS1997

A Viscoplastic Constitutive Theory for Monolithic Ceramic Materials

With increasing use of ceramic materials in high temperature structural applications such as advanced heat engine components, the need arises to accurately predict thermomechanical behavior. This paper, which is the first of two in a series, will focus on inelastic deformation behavior associated with these service conditions by providing an overview of a viscoplastic constitutive model that accounts for time-dependent hereditary material deformation (e.g., creep, stress relaxation, etc.) in monolithic structural ceramics. Early work in the field of metal plasticity indicated that inelastic deformations are essentially unaffected by hydrostatic stress. This is not the case, however, for ceramic-based material systems, unless the ceramic is fully dense. The theory presented here allows for fully dense material behavior as a limiting case. In addition, ceramic materials exhibit different time-dependent behavior in tension and compression. Thus, inelastic deformation models for ceramics must be constructed in a fashion that admits both sensitivity to hydrostatic stress and differing behavior in tension and compression. A number of constitutive theories for materials that exhibit sensitivity to the hydrostatic component of stress have been proposed that characterize deformation using time-independent classical plasticity as a foundation. However, none of these theories allow different behavior in tension and compression. In addition, these theories are somewhat lacking in that they are unable to capture creep, relaxation, and rate-sensitive phenomena exhibited by ceramic materials at high temperature. When subjected to elevated service temperatures, ceramic materials exhibit complex thermomechanical behavior that is inherently time-dependent, and hereditary in the sense that current behavior depends not only on current conditions, but also on thermo-mechanical history. The objective of this work is to present the formulation of a macroscopic continuum theory that captures these time-dependent phenomena. Specifically, the overview contained in this paper focuses on the multiaxial derivation of the constitutive model, and examines the scalar threshold function and its attending geometrical implications.

Janosik, Lesley A.↗

DOE OSTI2020

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

Reactor Containment Passive Safety Analysis: Steam Condensation in Presence of Non-condensable Gas Scaled Experiment and Modeling

This study presents steam condensation scaled experiments and semi-empirical models in presence of nitrogen (N)—a noncondensable gas (NCG), simulating air in the reactor containment—to support water-cooled small modular reactors (SMRs) passive containment cooling system (PCCS) design and analysis. Previous experimental studies on PCCS are focused on fixed and smaller tube (mostly 2-in.) geometries and specific test condition variations, bringing challenges with geometric scaling and mismatching with SMR prototypic design. To address these challenges, this study presents steam condensation test dataset obtained from three scaled test sections of 1-, 2-, and 4-in.-diameter steam condensers with an annular/jacket cooling of 2-, 3-, and 6 in.-diameter tubes, respectively. Test data were collected for steam ranges from 58 to 63 kg/hr., and NCG flow of 4.4 to 13.3 kg/hr. Annular cooling water flow was varied to obtain required testing conditions of saturated steam inlet and fully condensed outlet. Axial temperature test data of bulk cooling water, steam and condensate were collected by thermocouples for three test sections and various steam-NCG mixing/testing conditions. A standard data reduction method was adopted—utilizing iterative and nodalized mass and heat transfer calculation—to estimate axial local heat fluxes, heat transfer coefficients (HTCs), condensation rates, film thickness, and Nusselt number. Based on the obtained dataset semi-empirical model results—a ratio of experimental and Nusselt’s theoretical HTC are presented. Such results and findings are supportive of developing scaled-up testing facility, to enable model validations and accelerate next generation of reactors development and deployment

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Characterization of the mud displacement in an enlarged wellbore: An integrated rock-fluid model

Cement-mud displacement plays a crucial role in the sealability of cement sheaths. Irregular geometric features of a wellbore due to washout can have a negative impact on mud and cement mobilization. An unstable interface between two fluids always leads to mud channeling, interfluid mixing, and cement contamination, degrading the cement quality. Many factors, such as mechanical and rheological properties of fluids, annulus geometry, flow pattern, and flow rate, significantly influence the displacement efficiency. This study investigates the characterization of the mud displacement in an irregular horizontal well using a 3D computational fluid dynamics (CFD) model. Mud is displaced in an enlarged wellbore by geopolymer and neat class G cement. The wellbore geometry is developed based on the caliper log data from an unconventional shale well in the Tuscaloosa Marine Shale (TMS) lithology. The effects of pump rate, density difference, and mud contamination are evaluated by numerical simulations. The results present those residual muds mainly exist in the upper annulus of the enlarged section. Geopolymer has a better sealing performance and can resist more water-based mud (WBM) contaminations than neat class G cement. The scenario with a low mud-cement density difference and high cement injection rate results in a high cement volume fraction, mitigating the gas migration.

58 GEOSCIENCES↗