SEARCH · Engineering Papers

Results for “smoothing”

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.

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Diffuse X-ray scattering from polished silicon: application of the distorted wave Born approximation

Measured diffuse X-ray scattering data for a `smooth' as well as for a `rough' silicon sample were fit to theoretical expressions within the distorted wave Born approximation (DWBA). Data for the power spectral density (PSD) for both samples were also obtained by means of atomic force microscopy and optical interferometry. The Fourier transforms of trial correlation functions were fit to the PSD data and then applied to the DWBA formalism. The net correlation functions needed to fit the PSD data for each sample comprised the sum of two terms with different cutoff lengths and different self-affine fractal exponents. At zero distance these correlation functions added up to yield net values of σ 2 = (2) 2 and (71) 2 Å 2 for the smooth and rough samples, respectively. X-ray scattering data were obtained at beamline 1-BM of the Advanced Photon Source. Data and fits at values of q z = 0.05 and 0.10 Å −1 for the smooth sample are reported. Good fits for the smooth sample were obtained at both q z values simultaneously, that is, identical fitting parameters were applied at both values of q z . The smooth sample also exhibited weak Yoneda wings and a clear distinction between the strong specular scattering and the weak diffuse scattering. Data for the rough sample were qualitatively different and exhibited very weak scattering at the specular condition in contrast to extremely large Yoneda wings. Fits for the rough sample are reported for q z = 0.04, 0.05, and 0.06 Å −1 . Although the large Yoneda wings could be fit quite well in both position and amplitude, scattering near the specular condition could not be equally well fit by applying the same fitting parameters at all values of q z . Albeit imperfect, best-fitting results at the specular condition were obtained by invoking only diffuse scattering, that is, without including a separate theoretical expression for specular scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

A Sparse and Low Rank Penalized Signal Decomposition Model with Constraints: Anomaly Detection in PV Systems

Recently, robust PCA has seen its wide application in various industries for its ability to perform the task of anomaly detection. The essence of robust PCA approach is to break down the signal into a low rank component and sparse component. In many applications, a simple breakdown of the signal without accounting for the signs of low rank components and sparse components would violate the physical constraints of the decomposed signal. In addition, often times, the signals in the real world collected for a long duration has smooth changes within a day and between days. As an example, the power signals collected in a photovoltaic (PV) system are cyclostationary, exhibiting these characteristics. Neglecting the smoothness of signals would result in miss detection of anomalous signals which are smooth within a day but non-smooth between days and vice versa. In this paper, we developed a signal decomposition approach for the purpose of anomaly detection based on the idea of low rank and sparse decomposition taking into consideration the signs of the decomposed low rank and sparse components and the within-day and between-day smooth changes in the original signals. The proposed unsupervised approach for fault detection eliminates the need for faulty samples required by other machine learning methods. It does not require the full I-V characteristics to work. Furthermore, there is no need for complex modelling of PV systems as in the case of power loss analysis. Using Monte Carlo simulations, we demonstrate the ability of our proposed approach for detecting anomalies of different duration and severity in PV systems.

14 SOLAR ENERGY↗

NASA NTRS1989

The area of isodensity contours in cosmological models and galaxy surveys

The contour crossing statistic, defined as the mean number of times per unit length that a straight line drawn through the field crosses a given contour, is applied to model density fields and to smoothed samples of galaxies. Models in which the matter is in a bubble structure, in a filamentary net, or in clusters can be distinguished from Gaussian density distributions. The shape of the contour crossing curve in the initially Gaussian fields considered remains Gaussian after gravitational evolution and biasing, as long as the smoothing length is longer than the mass correlation length. With a smoothing length of 5/h Mpc, models containing cosmic strings are indistinguishable from Gaussian distributions. Cosmic explosion models are significantly non-Gaussian, having a bubbly structure. Samples from the CfA survey and the Haynes and Giovanelli (1986) survey are more strongly non-Gaussian at a smoothing length of 6/h Mpc than any of the models examined. At a smoothing length of 12/h Mpc, the Haynes and Giovanelli sample appears Gaussian.

Ryden, Barbara S.↗

NASA NTRS1993

A variational method for finite element stress recovery and error estimation

A variational method for obtaining smoothed stresses from a finite element derived nonsmooth stress field is presented. The method is based on minimizing a functional involving discrete least-squares error plus a penalty constraint that ensures smoothness of the stress field. An equivalent accuracy criterion is developed for the smoothing analysis which results in a C sup 1-continuous smoothed stress field possessing the same order of accuracy as that found at the superconvergent optimal stress points of the original finite element analysis. Application of the smoothing analysis to residual error estimation is also demonstrated.

Tessler, A.↗

NASA NTRS1993

The giant impact produced a precipitated Moon

The author's current simulations of Giant Impacts on the protoearth show the development of large hot rock vapor atmospheres. The Balbus-Hawley mechanism will pump mass and angular momentum outwards in the equatorial plane; upon cooling and expansion the rock vapor will condense refractory material beyond the Roche distance, where it is available for lunar formation. During the last seven years, the author together with several colleagues has carried out a series of numerical investigations of the Giant Impact theory for the origin of the Moon. These involved three-dimensional simulations of the impact and its aftermath using Smooth Particle Hydrodynamics (SPH), in which the matter in the system is divided into discrete particles whose motions and internal energies are determined as a result of the imposed initial conditions. Densities and pressures are determined from the combined overlaps of the particles, which have a bell-shaped density distribution characterized by a smoothing length. In the original series of runs all particle masses and smoothing lengths had the same values; the matter in the colliding bodies consisted of initial iron cores and rock (dunite) mantles. Each of 41 runs used 3,008 particles, took several weeks of continuous computation, and gave fairly good representations of the ultimate state of the post-collision body or bodies but at best crude and qualitative information about individual particles in orbit. During the last two years an improved SPH program was used in which the masses and smoothing lengths of the particles are variable, and the intent of the current series of computations is to investigate the behavior of the matter exterior to the main parts of the body or bodies subsequent to the collisions. These runs are taking times comparable to a year of continuous computation in each case; they use 10,000 particles with 5,000 particles in the target and 5,000 in the impactor, and the particles thus have variable masses and smoothing lengths (the latter are dynamically adjusted so that a particle typically overlaps a few tens of its neighbors). Since the matter in the impactor provides the majority of the mass left in orbit after the collision, and since the masses of the particles that originated in the impactor are smaller than those in the target, the mass resolution in the exterior parts of the problem is greatly improved and the exterior particles properly simulate atmospheres in hydrostatic equilibrium.

Cameron, A. G. W.↗

NASA NTRS1993

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

NASA NTRS1995

The effect of adding roughness and thickness to a transonic axial compressor rotor

The performance deterioration of a high speed axial compressor rotor due to surface roughness and airfoil thickness variations is reported. A 0.025 mm (0.001 in.) thick rough coating with a surface finish of 2.54-3.18 RMS microns (100-125 RMS microinches) is applied to the pressure and suction surface of the rotor blades. Coating both surfaces increases the leading edge thickness by 10% at the hub and 20% at the tip. Application of this coating results in a loss in efficiency of 6 points and a 9% reduction in the pressure ratio across the rotor at an operating condition near the design point. To separate the effect of thickness and roughness, a smooth coating of equal thickness is also applied to the blade. The smooth coating surface finish is 0.254-0.508 RMS microns (10-20 RMS microinches), compared to the bare metal blade surface finish of 0.508 RMS microns (20 RMS microinches). The smooth coating results in approximately half of the performance deterioration found from the rough coating. Both coatings are then applied to different portions of the blade surface to determine which portions of the airfoil are most sensitive to thickness/roughness variations. Aerodynamic performance measurements are presented for a number of coating configurations at 60%, 80%, and 100% of design speed. The results indicate that thickness/roughness over the first 10% of blade chord accounts for virtually all of the observed performance degradation for the smooth coating, compared to about 70% of the observed performance degradation for the rough coating. The performance deterioration is investigated in more detail at design speed using laser anemometer measurements as well as predictions generated by a quasi-3D Navier-Stokes flow solver which includes a surface roughness model. Measurements and analysis are performed on the baseline blade and the full-coverage smooth and rough coatings. The results indicate that coating the blade causes a thickening of the blade boundary layers. The interaction between the rotor passage shock and the thickened suction surface boundary layer then results in an increase in blockage which reduces the diffusion level in the rear half of the blade passage, thus reducing the aerodynamic performance of the rotor.

Suder, Kenneth L.↗

NASA NTRS2003

A Wind Tunnel Study of Icing Effects on a Business Jet Airfoil

Aerodynamic wind tunnel tests were conducted to study the effects of various ice accretions on the aerodynamic performance of a 36-inch chord, two-dimensional business jet airfoil. Eight different ice shape configurations were tested. Four were castings made from molds of ice shapes accreted in an icing wind tunnel. Two were made using computationally smoothed tracings of two of the ice shapes accreted in the icing tunnel. These smoothed profiles were then extended in the spanwise direction to form a two-dimensional ice shape. The final two configurations were formed by applying grit to the smoothed ice shapes. The ice shapes resulted in as much as 48% reduction in maximum lift coefficient from that of the clean airfoil. Large increases in drag and changes in pitching moment were also observed. The castings and their corresponding smoothed counterparts yielded similar results. Little change in performance was observed with the addition of grit to the smoothed ice shapes. Changes in the Reynolds number (from 3 x 10(exp 6) to 10.5 x 10(exp 6) and Mach number (from 0.12 to 0.28) did not significantly affect the iced-airfoil performance coefficients.

Addy, Harold E., Jr.↗

NASA NTRS2004

Kernel PLS Estimation of Single-trial Event-related Potentials

Nonlinear kernel partial least squaes (KPLS) regressior, is a novel smoothing approach to nonparametric regression curve fitting. We have developed a KPLS approach to the estimation of single-trial event related potentials (ERPs). For improved accuracy of estimation, we also developed a local KPLS method for situations in which there exists prior knowledge about the approximate latency of individual ERP components. To assess the utility of the KPLS approach, we compared non-local KPLS and local KPLS smoothing with other nonparametric signal processing and smoothing methods. In particular, we examined wavelet denoising, smoothing splines, and localized smoothing splines. We applied these methods to the estimation of simulated mixtures of human ERPs and ongoing electroencephalogram (EEG) activity using a dipole simulator (BESA). In this scenario we considered ongoing EEG to represent spatially and temporally correlated noise added to the ERPs. This simulation provided a reasonable but simplified model of real-world ERP measurements. For estimation of the simulated single-trial ERPs, local KPLS provided a level of accuracy that was comparable with or better than the other methods. We also applied the local KPLS method to the estimation of human ERPs recorded in an experiment on co,onitive fatigue. For these data, the local KPLS method provided a clear improvement in visualization of single-trial ERPs as well as their averages. The local KPLS method may serve as a new alternative to the estimation of single-trial ERPs and improvement of ERP averages.

Rosipal, Roman↗

NASA NTRS2004

Power Spectrum of Atmospheric Scintillation for the Deep Space Network Goldstone Ka-Band Downlink

Dynamic signal fluctuations due to atmospheric scintillations may impair the Ka-band (around 32-GHz) link sensitivities for a low-margin Deep Space Network (DSN) receiving system. The ranges of frequency and power of the fast fluctuating signals (time scale less than 1 min) are theoretically investigated using the spatial covariance and turbulence theory. Scintillation power spectrum solutions are derived for both a point receiver and a finite-aperture receiver. The aperture-smoothing frequency ((omega(sub s)), corner frequency ((omega(sub c)), and damping rate are introduced to define the shape of the spectrum for a finite-aperture antenna. The emphasis is put on quantitatively describing the aperture-smoothing effects and graphically estimating the corner frequency for a large aperture receiver. Power spectral shapes are analyzed parametrically in detail through both low- and high-frequency approximations. It is found that aperture-averaging effects become significant when the transverse correlation length of the scintillation is smaller than the antenna radius. The upper frequency or corner frequency for a finite-aperture receiver is controlled by both the Fresnel frequency and aperture-smoothing frequency. Above the aperture-smoothing frequency, the spectrum rolls off at a much faster rate of exp (-omega(sup 2)/omega(sup 2, sub s), rather than omega(sup -8/3), which is customary for a point receiver. However, a relatively higher receiver noise level can mask the fast falling-off shape and make it hard to be identified. We also predict that when the effective antenna radius a(sub r) less than or = 6 m, the corner frequency of its power spectrum becomes the same as that for a point receiver. The aperture-smoothing effects are not obvious. We have applied these solutions to the scenario of a DSN Goldstone 34-m-diameter antenna and predicted the power spectrum shape for the receiving station. The maximum corner frequency for the receiver (with omega(sub s) = 0.79 omega(sub 0) is found to be 0.44 Hz (or 1.0 omega(sub 0), while the fading rate (or fading slope) is about 0.06 dB/s.

Ho, C.↗

NASA NTRS2001

Use of a genetic algorithm for the analysis of eye movements from the linear vestibulo-ocular reflex

It is common in vestibular and oculomotor testing to use a single-frequency (sine) or combination of frequencies [sum-of-sines (SOS)] stimulus for head or target motion. The resulting eye movements typically contain a smooth tracking component, which follows the stimulus, in which are interspersed rapid eye movements (saccades or fast phases). The parameters of the smooth tracking--the amplitude and phase of each component frequency--are of interest; many methods have been devised that attempt to identify and remove the fast eye movements from the smooth. We describe a new approach to this problem, tailored to both single-frequency and sum-of-sines stimulation of the human linear vestibulo-ocular reflex. An approximate derivative is used to identify fast movements, which are then omitted from further analysis. The remaining points form a series of smooth tracking segments. A genetic algorithm is used to fit these segments together to form a smooth (but disconnected) wave form, by iteratively removing biases due to the missing fast phases. A genetic algorithm is an iterative optimization procedure; it provides a basis for extending this approach to more complex stimulus-response situations. In the SOS case, the genetic algorithm estimates the amplitude and phase values of the component frequencies as well as removing biases.

Non-NASA Center↗

NASA NTRS2007

Composite Elastic Skins for Shape-Changing Structures

Composite elastic skins having tailorable mechanical properties have been invented for covering shape-changing ( morphable ) structures. These skins are intended especially for use on advanced aircraft that change shapes in order to assume different aerodynamic properties. Many of the proposals for aircraft that could perform large aerodynamic shape changes require flexible skins that could follow shape changes of internal structures driven by actuators. Examples of such shape changes can include growth or shrinkage of bumps, conformal changes in wing planforms, cambers, twists, and bending of integrated leading- and trailing-edge flaps. Prior to this invention, there was no way of providing smooth aerodynamic surfaces capable of large deflections while maintaining smoothness and sufficient rigidity. Although latex rubber, silicone rubber, and similar conventional materials can be made into smooth coverings, they are not suitable for this purpose because, in order to impart required stiffness against out-of-plane bending, it would be necessary to make the coverings excessively thick, thereby necessitating the use of impractically large actuation forces. The basic idea of the invention is that of smoothly wrapping an underlying variable structure with a smooth skin that can be stretched or otherwise warped with low actuation force in one or both in-plane direction(s) and is relatively stiff against out-of-plane bending. It is envisioned that a skin according to the invention could be stretched as much as 20 percent in a desired direction. Because this basic idea admits of numerous variations, the following description is necessarily oversimplified for the sake of brevity.

Cagle, Christopher M.↗

NASA NTRS2004

Alternative Attitude Commanding and Control for Precise Spacecraft Landing

A report proposes an alternative method of control for precision landing on a remote planet. In the traditional method, the attitude of a spacecraft is required to track a commanded translational acceleration vector, which is generated at each time step by solving a two-point boundary value problem. No requirement of continuity is imposed on the acceleration. The translational acceleration does not necessarily vary smoothly. Tracking of a non-smooth acceleration causes the vehicle attitude to exhibit undesirable transients and poor pointing stability behavior. In the alternative method, the two-point boundary value problem is not solved at each time step. A smooth reference position profile is computed. The profile is recomputed only when the control errors get sufficiently large. The nominal attitude is still required to track the smooth reference acceleration command. A steering logic is proposed that controls the position and velocity errors about the reference profile by perturbing the attitude slightly about the nominal attitude. The overall pointing behavior is therefore smooth, greatly reducing the degree of pointing instability.

Singh, Gurkirpal↗

NASA NTRS2012

Prediction of Liquid Slosh Damping Using a High Resolution CFD Tool

Propellant slosh is a potential source of disturbance critical to the stability of space vehicles. The slosh dynamics are typically represented by a mechanical model of a spring mass damper. This mechanical model is then included in the equation of motion of the entire vehicle for Guidance, Navigation and Control analysis. Our previous effort has demonstrated the soundness of a CFD approach in modeling the detailed fluid dynamics of tank slosh and the excellent accuracy in extracting mechanical properties (slosh natural frequency, slosh mass, and slosh mass center coordinates). For a practical partially-filled smooth wall propellant tank with a diameter of 1 meter, the damping ratio is as low as 0.0005 (or 0.05%). To accurately predict this very low damping value is a challenge for any CFD tool, as one must resolve a thin boundary layer near the wall and must minimize numerical damping. This work extends our previous effort to extract this challenging parameter from first principles: slosh damping for smooth wall and for ring baffle. First the experimental data correlated into the industry standard for smooth wall were used as the baseline validation. It is demonstrated that with proper grid resolution, CFD can indeed accurately predict low damping values from smooth walls for different tank sizes. The damping due to ring baffles at different depths from the free surface and for different sizes of tank was then simulated, and fairly good agreement with experimental correlation was observed. The study demonstrates that CFD technology can be applied to the design of future propellant tanks with complex configurations and with smooth walls or multiple baffles, where previous experimental data is not available.

Yang, H. Q.↗

NASA NTRS2022

Global geologic map of asteroid (101955) Bennu indicates heterogeneous resurfacing in the past 500,000 years

Global geologic maps are useful tools for efficient interpretation of a planetary body, and they provide global context for the diversity and evolution of the surface. We used data acquired by the OSIRIS-REx spacecraft to create the first global geologic map of the near-Earth asteroid (101955) Bennu. As this is the first geologic map of a small, non-spherical, rubble-pile asteroid, we discuss the distinctive mapping challenges and best practices that may be useful for future exploration of similar asteroids, such as those to be visited with the Hera and Janus missions. By mapping on two centimeter-scale global image mosaics (2D projected space) and a centimeter-scale global shape model (3D space), we generated three input maps respectively describing Bennu’s shape features, geologic features, and surface texture. Based on these input maps, we defined two geologic units: the Smooth Unit and the Rugged Unit. The units are differentiated primarily on the basis of surface texture, concentrations of boulders, and the distributions of lineaments, mass movement features, and craters. They are bounded by several scarps. The Rugged Unit contains abundant boulders and signs of recent mass movement. It also has fewer small(<20 m), putatively fresh craters than the Smooth Unit, suggesting that such craters have been erased in the former. Based on these geologic indicators, we infer that the Rugged Unit has the younger surface of the two. Differential crater size-frequency distributions and the distribution of the freshest craters suggest that both unit surfaces formed ~10–65 million years ago, when Bennu was located in the Main Asteroid Belt, and the Smooth Unit has not been significantly resurfaced in the past 2 million years. Meanwhile, the Rugged Unit has experienced resurfacing within the past ~500,000years during Bennu’s lifetime as a near-Earth asteroid. The geologic units are consistent with global diversity in slope, surface roughness, normal albedo, and thermal emission spectral characteristics. The site on Bennu where the OSIRIS-REx mission collected a regolith sample is located in the Smooth Unit, in a small crater nested within a larger one. So although the Smooth Unit is an older surface than the Rugged Unit, the impact-crater settingindicates that the material sampled was recently exposed. Several similarities are apparent between Bennu and asteroid (162173) Ryugu from a global geologic perspective, including two geologic units distinguishable by variations in the number density of boulders, as well as in other datasets such as brightness.

E R Jawin↗

DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Heat Transfer Coefficients of Additively Manufactured Tubes with Internal Pin Fins for Supercritical Carbon Dioxide Cycle Recuperators

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in additively manufactured tubes with internal pin fins at the US DOE’s National Energy Technology Laboratory in Morgantown, WV. The measurement procedures were validated by conducting benchmark tests with smooth stainless-steel tube and comparing the results with published correlations for Nusselt number (Nu) and friction factor. Over Reynolds numbers (ReD) ranging from 5 × 10 4 to 2.5 × 10 5 , measured Nu was within 5% of the Dittus-Boelter correlation and measured friction factors were within 5% of the McAdams correlation for smooth tube flow. The candidate pin fin patterned pipes were additively manufactured (AM) at the Oak Ridge National Laboratory. The pins were circular or elliptical in cross-section and were printed at a 30⁰ angle relative to the inner wall (to meet AM constraints). The pin arrangement was helical to promote enhanced heat transfer due to swirl flow. Pin length to diameter aspect ratio was 1.33, 2, and 8, while the pin diameter to tube diameter ratio was 0.188, 0.125, and 0.063. Tests were performed for ReD varying from 6.9 × 10 4 to 2.2 × 10 5 and at conditions equivalent to the low pressure (LP) outlet (8.69 MPa, 361 K) and the high pressure (HP) inlet (20.7 MPa, 350 K) of the low temperature recuperator (LTR) in an indirect sCO 2 cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients. For the best performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.19), the local heat transfer coefficient increased by 136% relative to the Dittus-Boelter correlation at the LTR low pressure outlet and 194% at the LTR high pressure inlet. These correspond to a 282% and a 271% increase in the product of the heat transfer coefficient and surface area (adjusted for fin efficiency) product, respectively. Large pressure drops across the test articles were observed. For Tube Design A, the average friction factor, across the range of ReD considered, was significantly larger than the McAdams correlation at both the LTR LP outlet and the LTR HP inlet. A thermal performance factor was utilized to express the ratio of material required to build a finned heat exchanger relative to a finless heat exchanger with the same heat duty and pumping power. Tube Design A was estimated to decrease the required heat exchanger material by 13%.

Searle, Matthew↗

DOE OSTI2024

Simulation of Microcapsule Transport in Fractured Media Using Coupled CFD‐DEM

Geothermal energy is sustainable and gaining momentum as a solution to energy crises and environmental issues. However, challenges like production temperature and thermal breakthrough can impact geothermal project efficiency. One innovative solution to alleviate the thermal breakthrough is to inject polymer-based materials that are encapsulated in microcapsules into fractures to modify fracture permeability and prevent preferential flow. In our study, we utilized a coupled computational fluid dynamics and discrete element method to simulate the transport of microcapsules under various scenarios controlled by microcapsule size, microcapsule concentration, and fracture roughness. For a smooth fracture, the results indicate that small microcapsules can travel through a smooth fracture regardless of their concentrations. Large microcapsules can transport through a smooth fracture when present in lower concentrations. However, medium and mixed-size microcapsules tend to cause the sealing of a smooth fracture, irrespective of their concentrations. For a rough fracture, the transport of microcapsules is complicated by their interactions with the rough fracture walls. The presence of two sealing positions in a rough fracture adds further complexity to this transport phenomenon. The size and concentration of microcapsules control one sealing location, while the rough fracture walls determine the other sealing location. The rough walls substantially affect microcapsule transport, rendering the role of microcapsule size and concentration less significant. The simulation results suggest that complex fracture surfaces significantly elevate the occurrence of sealing behavior. To mitigate sealing behavior within more complex fractures, it would be beneficial to use smaller and lower concentrations of microcapsules.

15 GEOTHERMAL ENERGY↗