Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SCALE”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Exceptional Mineral Scaling Resistance from the Surface Gas Layer: Impacts of Surface Wetting Properties and the Gas Layer Charging Mechanism

Mineral scaling is a phenomenon that occurs on submerged surfaces in contact with saline solutions. In membrane desalination, heat exchangers, and marine structures, mineral scaling reduces process efficiency and eventually leads to process failure. Therefore, achieving long-term scaling resistance is beneficial to enhancing process performance and reducing operating and maintenance costs. While evidence shows that superhydrophobic surfaces may reduce mineral scaling kinetics, prolonged scaling resistance is limited due to the finite stability of the entrained gas layer present in a Cassie–Baxter wetting state. Additionally, superhydrophobic surfaces are not always feasible for all applications, but strategies for long-term scaling resistance with smooth or even hydrophilic surfaces are often overlooked. In this study, we elucidate the role of interfacial nanobubbles on the scaling kinetics of submerged surfaces of varied wetting properties, including those that do not entrain a gas layer. We show that both solution conditions and surface wetting properties that promote interfacial bubble formation enhances scaling resistance. In the absence of interfacial bubbles, scaling kinetics decrease as surface energy decreases, while the presence of bulk nanobubbles enhances the scaling resistance of the surface with any wetting property. The findings in this study allude to scaling mitigation strategies that are enabled by solution and surface properties that promote the formation and stability of interfacial gas layers and provide insights to surface and process design for greater scaling resistance.

calcite↗
DOE OSTI2020

On the regeneration of fish scales: structure and mechanical behavior

Fish scales serve as a dermal armor that provides protection from physical injury. Owing to a number of outstanding properties, fish scales are inspiring new concepts for layered engineered materials and next-generation flexible armors. Although past efforts have primarily focused on the structure and mechanical behavior of ontogenetic scales, the structure–property relationships of regenerated scales have received limited attention. In the present study, common carp ( Cyprinus carpio ) acquired from the wild were held live in an aquatic laboratory at 10°C and 20°C. Ontogenetic scales were extracted from the fish for analysis, as well as regenerated scales after approximately 1 year of development and growth. Their microstructure was characterized using microscopy and Raman spectroscopy, and the mechanical properties were evaluated in uniaxial tension to failure under hydrated conditions. We note the strength, strain to fracture and toughness of the regenerated scales were significantly lower than those of ontogenetic scales from the same fish, regardless of the water temperature. Scales that regenerated at 20°C exhibited significantly higher strength, strain to fracture and toughness than those regenerated at 10°C. The regenerated scales exhibited a highly mineralized outer layer, but no distinct limiting layer or external elasmodine; they also possessed a significantly lower number of plies in the basal layer than the ontogenetic scales. The results suggest that a mineralized layer develops preferentially during scale regeneration with the topology needed for protection, prior to the development of other qualities.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRSDate not supplied

Adaptive Control and Scaling Approach for the Emulation of Dynamic Subscale Torque Loads

Previous research by the authors proposed a control and scaling approach for emulating dynamic sub-scale torque loads. This approach produced an emulation controller that successfully regulated the dynamic behavior of a sub-scale electro-mechanical system intended to be a dynamical representation of the sub-scale turboelectric powertrain of a single-aisle commercial aircraft. The sub-scale system provides an environment without turbomachinery or rotors for the initial testing of electrified aircraft propulsion (EAP) control algorithms as they would be applied to a full-scale EAP system. The sub-scale turbomachinery/rotor torque loads were produced by electric machines (EMs) driven by the control and scaling approach, a full-scale turboelectric powertrain model, and an advanced EAP control algorithm. Although successfully tested, this approach produces an emulation controller that does not guarantee asymptotic stability. Modifying the original control law and integrating adaptive control techniques into the emulation controller allows the designer to guarantee asymptotic stability. This paper introduces the idea behind the emulation controller modifications, derives the controller, proves asymptotic stability, describes the implementation of the controller on a sub-scale electro-mechanical system intended to represent a parallel hybrid-electric turbofan engine, and describes the testing of a full-scale advanced EAP control algorithm. The turbofan engine and emulation controller performance are compared to results obtained using the previous, non-adaptive control and scaling approach.

adaptive↗
DOE OSTI2016

Fragmentation Prediction for a Brittle Ceramic: a Two-Scale Model Approach

Glass and oxide-ceramic radioactive material forms (e.g., PuO{sub 2}) are susceptible to brittle fracture during handling and transportation. Fracture of these materials may result in the generation of fragments small enough to become airborne, and thus easily dispersed and inhaled. Entrainment and dispersal of particles with dimensions less than 10 μm is generally anticipated under normal aerodynamic conditions. The assessment of the risk to the public and the environment from these materials due to events resulting in their fracture is based on the determination of the source term for environmental release. This source term is typically quantified by an airborne release fraction (ARF) and a respirable fraction (RF) of material generated by a given event. While the recommended empirical bounding values prove satisfactory to comply with the regulatory process for handling and transportation of these materials, they are based on test data for only a few materials subjected to a limited range of loading conditions. Because the final state of fragmentation (in terms of the size and distribution of fragments) is strongly dependent not only on the properties of the material and the energetic loading conditions that lead to breakup, but also on the many intrinsic length scales associated with fragmentation, it is difficult to precisely predict fragmentation characteristics for other brittle solid wastes forms subjected to a wider range of loading conditions based on the test data alone. The purpose of this manuscript is to present a two-scale modeling approach that may be used to simulate the dynamic fragmentation of brittle solid radioactive material forms subjected to general loading and boundary conditions. The method provides the resulting fragment characteristics (average size and size distribution) across the entire fragment range of interest. At the macro-scale, the approach uses a gradient damage elastic continuum mechanics material model within a finite element framework to simulate dynamic fracture under general loading and boundary conditions. Results from the macroscopic calculations are passed to a one-dimensional (1-D) lower length scale model that incorporates the relevant fine scale effects to determine the micro-scale fragmentation. This 1-D micro-scale fragmentation model simulates the internal crack nucleation and opening processes to predict the average fragment size and the fragment size distribution as a function of the material properties, loading conditions, and microscopic characteristic length scale of the material. Results from the lower length scale model are then combined with the macro-scale model fragmentation information to determine the resulting fragment size distribution spanning both length scales. This two-scale model has been exercised on a representative problem drawn from experimental studies, consisting of an impact of a falling weight onto a brittle ceramic. (authors)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2013

A multi-scale metrics approach to forest fragmentation for Strategic Environmental Impact Assessment

Forests are becoming severely fragmented as a result of land development. South Korea has responded to changing community concerns about environmental issues. The nation has developed and is extending a broad range of tools for use in environmental management. Although legally mandated environmental compliance requirements in South Korea have been implemented to predict and evaluate the impacts of land-development projects, these legal instruments are often insufficient to assess the subsequent impact of development on the surrounding forests. It is especially difficult to examine impacts on multiple (e.g., regional and local) scales in detail. Forest configuration and size, including forest fragmentation by land development, are considered on a regional scale. Moreover, forest structure and composition, including biodiversity, are considered on a local scale in the Environmental Impact Assessment process. Recently, the government amended the Environmental Impact Assessment Act, including the SEA, EIA, and small-scale EIA, to require an integrated approach. Therefore, the purpose of this study was to establish an impact assessment system that minimizes the impacts of land development using an approach that is integrated across multiple scales. This study focused on forest fragmentation due to residential development and road construction sites in selected Congestion Restraint Zones (CRZs) in the Greater Seoul Area of South Korea. Based on a review of multiple-scale impacts, this paper integrates models that assess the impacts of land development on forest ecosystems. The applicability of the integrated model for assessing impacts on forest ecosystems through the SEIA process is considered. On a regional scale, it is possible to evaluate the location and size of a land-development project by considering aspects of forest fragmentation, such as the stability of the forest structure and the degree of fragmentation. On a local scale, land-development projects should consider the distances at which impacts occur in the vicinity of the forest ecosystem, and these considerations should include the impacts on forest vegetation and bird species. Impacts can be mitigated by considering the distances at which these influences occur. In particular, this paper presents an integrated environmental impact assessment system to be applied in the SEIA process. The integrated assessment system permits the assessment of the cumulative impacts of land development on multiple scales. -- Highlights: • The model is to assess the impact of forest fragmentation across multiple scales. • The paper suggests the type of forest fragmentation on a regional scale. • The type can be used to evaluate the location and size of a land development. • The paper shows the influence distance of land development on a local scale. • The distance can be used to mitigate the impact at an EIA process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗
DOE OSTI2021

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Scaling patch analysis of planar turbulent wakes

In this work, a scaling patch approach is used to investigate the proper scales in planar turbulent wakes. A proper scale for the mean axial flow is the well-known maximum velocity deficit U ref = U ∞ –U ctr , where U ∞ is the free stream velocity and U ctr is the mean axial velocity at the wake centerline. From an admissible scaling of the mean continuity equation, a proper scale for the mean transverse flow is found as V ref = (dδ/dx)U ref , where dδ/dx is the growth rate of the wake width. From an admissible scaling of the mean momentum equation, a proper scale for the kinematic Reynolds shear stress is found as R uv,ref = U ∞ V ref , which is a mixed scale of the free stream velocity and the mean transverse flow scale. Expressions are derived for the scaled mean transverse velocity and Reynolds shear stress in the far field of planar turbulent wakes. Using a Gaussian function for the mean axial velocity deficit, approximate functions for the scaled mean transverse velocity and Reynolds shear stress are developed and found to agree well with experimental and simulation data. This work reveals that the mean transverse flow, despite its small magnitude, plays an important role in the scaling and understanding of the planar turbulent wake.

42 ENGINEERING↗
DOE OSTI2024

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2003

Comparisons of the Impact Responses of a 1/5-Scale Model and a Full-Scale Crashworthy Composite Fuselage Section

A 25-fps vertical drop test of a 1/5-scale model composite fuselage section was conducted to replicate a previous test of a full-scale fuselage section. The purpose of the test was to obtain experimental data characterizing the impact response of the 1/5-scale model fuselage section for comparison with the corresponding full-scale data. This comparison is performed to assess the scaling procedures and to determine if scaling effects are present. For the drop test, the 1/5-scale model fuselage section was configured in a similar manner as the full-scale section, with lead masses attached to the floor through simulated seat rails. Scaled acceleration and velocity responses are compared and a general assessment of structural damage is made. To further quantify the data correlation, comparisons of the average acceleration data are made as a function of floor location and longitudinal position. Also, the percentage differences in the velocity change (area under the acceleration curve) and the velocity change squared (proportional to kinetic energy) are compared as a function of floor location. Finally, correlation coefficients are calculated for corresponding 1/5- and full-scale data channels and these values are plotted versus floor location. From a scaling perspective, the differences between the 1/5- and full-scale tests are relatively small, indicating that appropriate scaling procedures were used in fabricating the test specimens and in conducting the experiments. The small differences in the scaled test data are attributed to minor scaling anomalies in mass, potential energy, and impact attitude.

Jackson, Karen E.↗
NASA NTRSDate not supplied

TOWARD A METHOD FOR SCALING HUMAN BODY MODELS IN AN IMU-BASED WORKFLOW

BACKGROUND Scaled biomechanical models can more accurately inform crew health decisions when tailored to the wide range of astronaut sizes. One component to improve scaling of existing models to better represent each unique astronaut’s size is the individual length scaling of limbs. Traditionally, these lengths are determined by motion capture or manual measurement. A new method is herein proposed for length scaling which can be done by measuring linear and angular accelerations at a desired point during isolated motion around a point of rotation, then calculating the distance between the desired point and point of rotation. When an Inertial Measurement Unit (IMU) device is placed at the distal point of a limb, the isolated motion is about that limb’s proximal joint. For example, to measure forearm length, an IMU is placed at the wrist, the point of rotation is at the elbow, and the isolated motion is forearm flexion and extension. These calculated lengths are then used to scale models to each unique astronaut’s size, thereby improving the applicability of the model. This method of scaling limb segments can be used for any limb that has an easily defined proximal joint for the limb to rotate around including hands, arms, legs, feet. Utilizing IMUs for data collection also provides the synergistic ability to record data without a dedicated space in a room with many cameras, therefore reducing the data collection footprint, or record data where optical motion capture is not possible, such as inside a spacesuit. METHODS AND RESULTS To test this method, upper body data collection was performed with 5 Xsens DOT IMUs on a single subject. IMUs consist of an accelerometer, a gyroscope, and a magnetometer which collect linear acceleration, angular velocity, and magnetic fluctuations, respectively. Before any ground-based laboratory collection, the magnetic fluctuations are used to correct the heading of the IMU in space relative to the Earth’s magnetic field. The direct measurement of angular velocity is integrated to calculate angular acceleration. Then the linear acceleration ( a ) and angular acceleration (α) are solved using r = at/α to calculate the radius, which in this case is the distance between the IMU and the point of rotation (i.e., segment length). The IMU must be placed at the most distal point of the limb being measured (i.e., ankle if measuring lower leg length) and the test plan must consist of an isolated motion about that limb’s proximal joint (i.e., knee flexion and extension if measuring lower leg length). The distances (radii) calculated at every time interval were filtered (bandpass filter keeping 5th-90th percentile data) to eliminate outliers and spurious data that occur when the isolated motion was stopped or nearly stopped. The remaining distances were averaged, resulting in the calculated limb length. Scaling factors were then computed by dividing the calculated limb length by the unscaled model’s length. These scale factors are plugged into the Scale Tool in OpenSim [1,2] to apply the scaling to the OpenSim Full Body Rajagopal Model [3,4]. Manual measurements of limb lengths were taken before data collection started and used for comparing against the calculated lengths. The Anthropometric Survey of US Army Personnel (ANSUR II) [5] was also used as a third source of reference for limb length measurements. The following measurements were retrieved from the subject before data collection: 34.5 cm from L1 to C7 (thorax), 25.7 cm from C7 Joint Center (JC) to head vertex (neck and head), 36.3 cm from shoulder JC to elbow JC (humerus), 29.5 cm from elbow JC to wrist JC (forearm), and 16.2 cm from clavicle to acromion (clavicle). Of those five, forearm and humerus lengths were calculated using this proposed method to obtain preliminary results. The forearm length after filtering and averaging was calculated to be 37.6 cm. This is a 28% overestimation from the measured forearm length (29.5 cm). The humerus length after filtering and averaging was calculated to be 47.8 cm. This is a 31% difference from the measured subject length (36.3 cm). Sources of error include imperfect isolated motion (method currently expects that motion should be perfectly circular in a 2D plane, include no rotation of the IMU, and be relatively smooth; a more secure IMU attachment method will help), unrefined filter techniques (removed highest and lowest values with 20% high and low pass filters and no smoothing filters), arbitrary removal of stopped or near stopped data (kept data for only a short range before and after the angular velocity peaking), and a more representative method for removing gravitational acceleration is needed (current method is to zero all accelerations against a baseline taken just before the isolated motion which does not account for the gravitational acceleration changed due to IMU rotation during movement). Addressing these error sources will improve the accuracy of the limb length calculation. Next steps include creating a method for whole-body scaling estimation using individual limb scale factors. Continued pursuit of these techniques is expected to enable acquiring anthropometric information using only IMUs in real-time.

E. K. Marecki↗
DOE OSTI2016

Uranium and strontium fate in waste-weathered sediments: Scaling of molecular processes to predict reactive transport

Objectives of the project: 1. Determine the process coupling that occurs between mineral transformation and contaminant (U and Sr) speciation in acid-uranium waste weathered Hanford sediments. 2. Establish linkages between molecular-scale contaminant speciation and meso-scale contaminant lability, release and reactive transport. 3. Make conjunctive use of molecular- to bench-scale data to constrain the development of a mechanistic, reactive transport model that includes coupling of contaminant sorption-desorption and mineral transformation reactions. Hypotheses tested: - Uranium and strontium speciation in legacy sediments from the U-8 and U-12 Crib sites can be reproduced in bench-scale weathering experiments conducted on unimpacted Hanford sediments from the same formations. - Reactive transport modeling of future uranium and strontium releases from the vadose zone of acid-waste weathered sediments can be effectively constrained by combining molecular-scale information on contaminant bonding environment with grain-scale information on contaminant phase partitioning, and meso-scale kinetic data on contaminant release from the waste-weathered porous media. - Although field contamination and laboratory experiments differ in their diagenetic time scales (decades for field vs. months to years for lab), sediment dissolution, neophase nucleation, and crystal growth reactions that occur during the initial disequilibrium induced by waste-sediment interaction leave a strong imprint that persists over subsequent longer-term equilibration time scales and, therefore, give rise to long-term memory effects. Enabling capabilities developed: Our team developed an iterative measure-model approach that is broadly applicable to elucidate the mechanistic underpinnings of reactive contaminant transport in geomedia subject to active weathering. Experimental design: Hypotheses were tested by comparing (with a similar set of techniques) the geochemical transformations and transport behaviors that occured in bench-scale studies of waste-sediment interaction with parallel model systems studies of homogeneous nucleation and neo-phase dissolution. Initial plans were to compare results with core sample extractions from the acid uranium waste impacted U-8 and U-12 Cribs at Hanford (see original proposal and letter of collaboration from J. Zachara). However, this part of the project was impossible because funding for core extractions were eliminated from the DoE budget. Three distinct crib waste aqueous simulants (whose composition is based on the most up-to-date information from field site investigations) were reacted with Hanford sediments in batch and column systems. Coupling of contaminant uptake to mineral weathering was monitored using a suite of methods both during waste-sediment interaction, and after, when waste-weathered sediments were subjected to infusion with circumneutral background pore water solutions. Our research was designed to adapt as needed to maintain a strong dialogue between laboratory and modeling investigations so that model development was increasingly constrained by emergent data and understanding. Potential impact of the project to DOE: Better prediction of contaminant uranium transport was achieved by employing multi-faceted lines of inquiry to build a strong bridge between molecular- and field-scale information. By focusing multiple lines and scales of observation on a common experimental design, our collaborative team revealed non-linear and emergent behavior in contaminated weathering systems. A goal of the current project was to expand our modeling capabilities, originally focused on hyperalkaline legacy waste streams, to include acidic weathering reactions that, as described above, were expected to result in profoundly different products. We were able to achieve this goal, and showed that these products nonetheless undergo analogous silicate and non-silicate transformation, ripening and aging processes. Our prediction that these weathering reactions would vary with waste stimulant chemistry resulted in data that was incorporated directly into a reactive transport model structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Pilot-Scale Testing of an Integrated Circuit for the Extraction of Rare Earth Minerals and Elements from Coal and Coal Byproducts Using Advanced Separation Technologies

The primary objective of this project was to develop and demonstrate an integrated pilot-scale circuitry for recovering high-value rare earth elements (REEs) from coal and coal byproducts. The target performance was to produce a mixed REE product with content of at least two percent by weight on a dry mass basis in a cost-effective and environmentally benign manner. During the first nine months of the project period (Phase 2 Budget Period 2), pilot plant construction was completed including all field site startup activities such as permitting, engineering design, procurement/bidding, unit fabrication, site construction, equipment installation, module assembly, safety training, and circuit shakedown. During the remaining 21 months of project period (Phase 2 Budget Period 3), detailed field-testing activities were performed including feedstock sample collection and preparation, exploratory testing, circuit modification, detailed parametric study, and performance optimization. A detailed techno-economic analysis was performed based on the pilot plant testing findings which provided various scenarios for REE production. The project successfully accomplished the proposed target performance by producing mixed rare earth oxide (REO) with greater than 90% purity by weight in a continuous pilot scale operation from two distinctly different coarse refuse materials (i.e., West Kentucky No. 13 and Fire Clay coal seams), and at least three secondary sources (i.e., heap leach process and naturally formed acid mine drainage system). Project partners included the University of Kentucky, Virginia Tech, West Virginia University, Alliance Coal, Blackhawk Mining, Mineral Refining Company, and Mineral Separation Technologies. The pilot scale test facility was constructed at a former mining complex owned by Alliance Natural Resource Partners (Alliance Coal). The site was rehabilitated to accommodate the equipment installation, construction and fabrication, electrical power requirement, water line management and containment. The process units constructed and installed included X-ray sorting unit, crushing and grinding unit, physical separation unit, acid leaching unit, solvent extraction unit, and wastewater management unit. A rare earth mineral concentration unit was constructed as a standalone unit for flexible operation. A detailed environmental assessment and control plan was carried out to identify and quantify any potential impacts of the pilot-scale processing circuitry on the human and eco-system health and well-being. Corresponding mitigation strategies and control measures were provided. A conceptual flowsheet was developed to effectively remove thorium and uranium from high purity rare earth oxide mix or any potential radionuclide enriched stream. The two distinct feedstock materials were secured from the Blackhawk Mining Complex in eastern Kentucky where the Fire Clay (Hazard No. 4) seam is processed. The West Kentucky No. 13 (Baker) coarse refuse material was collected from an active process stream at an Alliance coal preparation plant located in western Kentucky. Characterization analysis indicated that both of feed materials generated from the two sources contained >300 ppm of TREEs on a dry whole mass basis which met the requirements for a qualified feed stock. The two feedstocks were further upgraded using a dual x-ray sorter to prepare the feed material for hydrometallurgical circuit. Thermal treatment on feed material prior to leaching was found to: 1) improve the leaching recovery of REEs, 2) increase the leaching kinetics, and 3) allow the leaching reaction to occur at lower acidity. Roasting at 600°C was selected as the pre-treatment condition for both West Kentucky No. 13 and Fire Clay coarse refuse material. Over 40% of leaching recovery was achieved by roasting West Kentucky No. 13 material having a top particle size of 3 mm in the pilot scale operation using 1.2M sulfuric acid leaching at 75OC. Initial pilot scale testing involved continuous operation of the pilot plant for 94 hours. The leaching unit was operated at solid-to-liquid ratio of 1 to 10 (w/v) using 0.5M sulfuric acid solution at a temperature of 75°C. The continuous solvent extraction circuit utilized rougher and cleaner units with DEHPA and TBP as the extractants. An innovative stripping circuit was developed to accumulate the REE concentration in the stripping solution to a level above 600 ppm. A bleed stream from the recycled strip solution was treated using oxalic acid precipitation which produced a high grade rare earth oxalate. The oxalate product was roasted to remove the oxalate which produced a rare earth oxide product having a purity greater 90%. Due to high concentrations of contaminant ions in the pregnant leach solution (PLS), a modified flowsheet was developed that involved pre-concentration of the REEs using multiple stages of precipitation and redissolution. The advantage of this process was improved removal of contamination before the downstream purification process and a significant cost reduction relative to the circuit that utilized the solvent extraction process. The modified circuitry included processes involving leaching, multistage precipitation, redissolution, and oxalate precipitation followed by roasting of the oxalate product. The circuit produced a mixed REO that was 92.96% pure from the initial test. A detailed parametric test plan was carried out which involved varying key parameters including solids feed rate, acid flowrate, acid concentration, multistage precipitation pH, redissolution pH, oxalate precipitation dosage and pH. The response variables included REE recovery, contaminant recovery, REO product grade and overall chemical consumption. Test results indicated that the acid-to-solid ratio is the key parameter to leaching efficiency as performance deteriorated with an increase in solids concentration. The optimal pH determined for REE precipitation and redissolution was 6.5 and 2.5, respectively. Additional tests were conducted to further improve the flowsheet. Recirculating a portion of the PLS to the feed of the leach tanks improved the leaching performance by lowering the pH of the leaching system and reducing the contamination recovery by shortening the residence time. Moreover, the removal of Al prior to REE precipitation significantly reduced the oxalic acid consumption in the oxalate precipitation circuit. The modified circuit produced over 90% grade REO by weight from both West Kentucky No. 13 and Fire Clay coarse refuse material in pilot scale continuous test programs. A case analysis model was developed to project the REE and major contaminants concentration in each PLS stream based on the leaching condition, pH cut point, and oxalic acid dosage. A correlation was established using empirical and semi-empirical models. Using the models, chemical consumption required for each stage was predicted based on the projected performance of the hydrometallurgy circuit. After identifying the optimum conditions, validation tests were carried out for the treatment of both West Kentucky No. 13 and Fire Clay coarse refuse materials in the pilot plant. The actual circuit performance and chemical consumptions were very close to the model predictions. Other than the two coarse refuse sources, several secondary feedstocks were also tested in the pilot plant facility. A “heap leach” system was constructed using the coarse refuse material generated from cleaning the West Kentucky No. 13 seam coal. Using the two stage SX rougher and cleaner circuit, a concentrate with a grade >90% REO was produced while recovering >97% of the REEs from the heap leach PLS. Naturally generated acid mine drainage (AMD) from West Kentucky No.13 mine was processed using the multistage precipitation circuit in the pilot plant in a test conducted for a period of 32 hours. The final grade of the mix RE oxide produced from the AMD was 90.84% with an overall circuit recovery of 64%. The primary source of REE was the selective precipitation steps involving iron and aluminum rejection. The hydrophobic-hydrophilic separation (HHS) process was proven to effectively recover coal from fine waste materials. For REM recovery, the HHS process was able to produce concentrates at grades of approximately 1.8% REE on an ash basis; however, recovery values were typically low, <10%, under the optimal conditions determined in the laboratory-scale testing. Staged testing of the pilot-scale HHS process for coal recovery and semi-continuous laboratory testing for REM testing showed that a total concentration ratio of more than 15x was observed for the REM recovery process. A circuit simulation package was developed for REE extraction and purification using a spreadsheet-based platform (Microsoft Excel). The REESim circuit simulation package is configured to track the mass and volume flows of components passing through a series of unit operations specified and configured by the user. The mass rates can then be utilized by the user to determine important performance indicators such as product mass yields, concentrate purity levels, element-by-element recoveries, and so forth. The techno-economic analysis showed that the roasting and leaching operations were the most expensive capital items, each contributing approximately 30% to the total capital cost. One notable contributor to the high production costs was the low REE recovery observed in the pilot scale trials. The product basket price was shown to have a strong influence on the economic viability of the scenarios, with the scandium price being the most significant influencer. Operating cost was shown to be extremely sensitive to REE recovery, REE feed grade, and leaching acid consumption. An analysis of ten different scenarios for a 500 t/h commercial operation revealed that three were economically favorable, producing internal rates of return varying from 27.7% to 33.1% and payback periods of 4 to 5 years. The project successfully developed and demonstrated a process to recover REEs from coal and coal byproducts in a pilot-plant operation which consistently produced over 90% grade REO mix from varies types of feedstocks. Commercialization analysis showed that the technology readiness level successfully achieved TRL 6 at the end of the project and demonstrated the need and the potential for scaling the process to further advance the technologies toward the goal of providing a domestic supply of REEs at a commercial scale.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2016

Grizzly: A Multi-scale and Multi-Physics Tool to Model Aging of Nuclear Power Plant Components

Grizzly is a multi-scale multi-physics tool being developed at Idaho National Laboratory (INL) as part of the US Department of Energy's Light Water Reactor Sustainability program to provide improved safety assessments of systems, components, and structures in nuclear power plants subjected to age-related degradation. Its goal is to provide an improved scientific basis for decisions surrounding license renewal, which would permit operation of commercial nuclear power plants beyond 60 years. Grizzly is based on INL's MOOSE framework, which enables multi-physics simulations in a massively parallel computing environment. The reactor pressure vessel (RPV) has been chosen as the initial application for Grizzly. The primary concern for the RPV is that it becomes increasingly susceptible to fracture initiation at the sites of pre-existing flaws over time. This is because exposure to irradiation and high temperatures embrittles the steel. Assessing the safety of an RPV under transient conditions involves simulations at multiple scales. At the engineering scale, strongly coupled equations of heat conduction and solid mechanics are solved to simulate the global response of the RPV to accident conditions. Sub-models are used to represent regions with pre-existing flaws and the stress intensity factors are evaluated from the domain integrals. Crack propagation is then determined using a model that provides the temperature and neutron fluence dependent fracture toughness of the RPV steel. The ability of Grizzly to provide accurate engineering scale estimates of RPV safety beyond 60 years of operation strongly relies on the predictability of the embrittlement model. Multiple models for embrittlement have been proposed, and the physically based empirical model is implemented in Grizzly. However, the formation of late blooming phases greatly reduces the confidence on the extrapolated values obtained from embrittlement model, which is further exacerbated by the absence of real field data for validation. Experimentally validated lower length scale models can significantly improve the predictability of the embrittlement model beyond the calibrated range and the development of such models is another ongoing effort in Grizzly. The modeling activity in Grizzly spans across different length scales, and can be divided into: (i) Irradiation induced microstructure evolution at the atomistic and single crystal length scale using molecular dynamics, atomistic kinetic Monte-Carlo, phase-field and rate theory; (ii) Effect of the irradiation induced defects on the flow stress and fracture processes at the polycrystalline scale using crystal plasticity. Subsequently, the prediction from these models can be used to obtain irradiation induced ductile-to-brittle transition temperature (DBTT) shifts from specimen-scale simulations; (iii) Global thermo-mechanical response of the RPV under normal and accident scenarios, and probabilistic risk analyses at the engineering scale using the DBTT shifts obtained from the lower length scale models. The item (i) has been discussed in a previous report. In this extended abstract an overview of items (ii) and (iii) is provided. The extended abstract summarizes the current status of the modeling effort under Grizzly. At the engineering scale, reduced order model development is underway using the statistical tool RAVEN. This will allow probabilistic analysis considering a wide range of flaw topologies in a computationally efficient manner. Coupling the irradiation induced defect growth predictions to crystal plasticity is also in progress. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Pore-Scale Simulation of Intragranular Diffusion: Effects of Incomplete Mixing on Macroscopic Manifestations

Diffusive mass transfer into and out of intragranular micropores ("intragranular diffusion") plays an important role in the large-scale transport of some groundwater contaminants. We are interested in understanding the combined effect of pore-scale advection and intragranular diffusion on solute transport at the effective porous medium scale. To study this problem, we have developed a three-dimensional pore-scale model of fluid flow and solute transport that incorporates diffusion into and out of intragranular pore spaces. Our model is based on the Smoothed Particle Hydrodynamics (SPH) simulation method, which represents fluid and solid phases by a mesh-free particle discretization. In the pore spaces, fluid flow is simulated by discretizing the Navier-Stokes equations using the SPH approach. Solute transport is represented by advection, diffusion within the fluid phase, and diffusion between the fluid and solid phases. Our model is implemented on large-scale parallel computing hardware, allowing us to simulate millions of computational particles and represent fully three-dimensional systems of pores and grains with arbitrarily complex physical geometry. We have used this model system to perform numerical experiments using various model porous media systems, which allows us to draw comparisons between macroscopic measures computed from the pore-scale simulations (such as breakthrough curves) and those predicted by macroscopic formulations that assume complete mixing over the representative volume. In this paper we present results of 3D simulations of pore-scale flow and transport, including cases with and without intragranular diffusion, in two model porous media, one with randomly-packed uniform spherical grains and a second with randomly-packed spheres drawn from a binary grain size distribution. Breakthrough curves were computed from the 3D simulations at various transport distances. Comparable breakthrough curves were computed using 1D macroscopic models with parameters determined independently to the degree possible. Based on comparisons of the pore-scale and macroscopic model results, we draw two primary conclusions. First, non-Fickian behavior is persistent and ubiquitous at the scales considered, and most cases are better represented by a multi-rate mass transfer model even when there is no distinct secondary porosity (i.e., no intragranular diffusion). This suggests that diffusive mass transfer processes between preferential flow paths and relatively immobile zones within the primary porosity may have significant impact on transport, particular in low-concentration tails. Second, the application of mass transfer rate parameters based on an assumption of well-mixed concentrations at the pore-scale tend to overestimate the amount of mass transfer that occurs in heterogeneous pore geometries in which preferential flow leads to incomplete pore-scale lateral mixing.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2025

Scaling Arctic landscape and permafrost features improves active layer depth modeling

Tundra ecosystems in the Arctic store up to 40% of global below-ground organic carbon but are exposed to the fastest climate warming on Earth. However, accurately monitoring landscape changes in the Arctic is challenging due to the complex interactions among permafrost, micro-topography, climate, vegetation, and disturbance. This complexity results in high spatiotemporal variability in permafrost distribution and active layer depth (ALD). Moreover, these key tundra processes interact at different scales, and an observational mismatch can limit our understanding of intrinsic connections and dynamics between above and below-ground processes. Consequently, this could limit our ability to model and anticipate how ALD will respond to climate change and disturbances across tundra ecosystems. In this paper, we studied the fine-scale heterogeneity of ALD and its connections with land surface characteristics across spatial and spectral scales using a combination of ground, unoccupied aerial system, airborne, and satellite observations. We showed that airborne sensors such as AVIRIS-NG and medium-resolution satellite Earth observation systems like Sentinel-2 can capture the average ALD at the landscape scale. We found that the best observational scale for ALD modeling is heavily influenced by the vegetation and landform patterns occurring on the landscape. Landscapes characterized by small-scale permafrost features such as polygon tussock tundra require high-resolution observations to capture the intrinsic connections between permafrost and small-scale land surface and disturbance patterns. Conversely, in landscapes dominated by water tracks and shrubs, permafrost features manifest at a larger scale and our model results indicate the best performance at medium resolution (5 m), outperforming both higher (0.4 m) and lower resolution (10 m) models. This transcends our study to show that permafrost response to climate change may vary across dominant ecosystem types, driven by different above- and below-ground connections and the scales at which these connections are happening. We thus recommend tailoring observational scales based on landforms and characteristics for modeling permafrost distribution, thereby mitigating the influences of spatial-scale mismatches and improving the understanding of vegetation and permafrost changes for the Arctic region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Spatially Local Surrogate Modeling of Subgrid-Scale Effects in Idealized Atmospheric Flows: A Deep Learned Approach Using High-Resolution Simulation Data

Abstract We introduce a machine learned surrogate model from high-resolution simulation data to capture the subgrid-scale effects in dry, stratified atmospheric flows. We use deep neural networks (NNs) to model the spatially local state differences between a coarse-resolution simulation and a high-resolution simulation. The setup enables the capture of both dissipative and antidissipative effects in the state differences. The NN model is able to accurately capture the state differences in offline tests outside the training regime. In online tests intended for production use, the NN-coupled coarse simulation has higher accuracy over a significant period of time compared to the coarse-resolution simulation without any correction. We provide evidence of the capability of the NN model to accurately capture high-gradient regions in the flow field. With the accumulation of the errors, the NN-coupled simulation becomes computationally unstable after approximately 90 coarse simulation time steps. Insights gained from these surrogate models further pave the way for formulating stable, complex, physics-based spatially local NN models which are driven by traditional subgrid-scale turbulence closure models. Significance Statement Flows in the atmosphere are highly chaotic and turbulent, comprising flow structures of broad scales. For effective computational modeling of atmospheric flows, the effects of the small- and large-scale structures need to be captured by the simulations. Capturing the small-scale structures requires fine-resolution simulations. Even with the current state-of-the-art supercomputers, it can be prohibitively expensive to simulate these flows when computed for the entire earth over climate time scales. Thus, it is necessary to focus on the larger-scale structures using a coarse-resolution simulation while capturing the effects of the smaller-scale structures using some parameterization (approximation) scheme and incorporating it into the coarse-resolution simulation. We use machine learning to model the effects of the small-scale structures (subgrid-scale effects) in atmospheric flows. Data from a fine-resolution simulation is used to compute the missing subgrid-scale effects in coarse-resolution simulations. We then use machine learning models to approximate these differences between the coarse- and fine-resolution simulations. We see improved accuracy for the coarse-resolution simulations when corrected using these machine learned models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Impact of scaling on the nitric-glycolic acid flowsheet

Savannah River Remediation (SRR) is considering using glycolic acid as a replacement for formic acid in Sludge Receipt and Adjustment Tank (SRAT) processing in the Defense Waste Processing Facility (DWPF). Catalytic decomposition of formic acid is responsible for the generation of hydrogen, a potentially flammable gas, during processing. To prevent the formation of a flammable mixture in the offgas, an air purge is used to dilute the hydrogen concentration below the 60% of the Composite Lower Flammability Limit (CLFL). The offgas is continuously monitored for hydrogen using Gas Chromatographs (GCs). Since formic acid is much more volatile and toxic than glycolic acid, a formic acid spill would lead to the release of much larger quantities to the environment. Switching from formic acid to glycolic acid is expected to eliminate the hydrogen flammability hazard leading to lower air purges, thus downgrading of Safety Significant GCs to Process Support GCs, and minimizing the consequence of a glycolic acid tank leak in DWPF. Overall this leads to a reduction in process operation costs and an increase in safety margin. Experiments were completed at three different scales to demonstrate that the nitric-glycolic acid flowsheet scales from the 4-L lab scale to the 22-L bench scale and 220-L engineering scale. Ten process demonstrations of the sludge-only flowsheet for SRAT and Slurry Mix Evaporator (SME) cycles were performed using Sludge Batch 8 (SB8)-Tank 40 simulant. No Actinide Removal Process (ARP) product or strip effluent was added during the runs. Six experiments were completed at the 4-L scale, two experiments were completed at the 22-L scale, and two experiments were completed at the 220-L scale. Experiments completed at the 4-L scale (100 and 110% acid stoichiometry) were repeated at the 22-L and 220-L scale for scale comparisons.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗