SEARCH · Engineering Papers

Results for “spatial discretization”

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 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1990

Clutter filter design considerations for Airborne Doppler radar detection of windshear

The problem of clutter rejection when processing down-looking Doppler radar returns from a low altitude airborne platform is a paramount problem. With radar as a remote sensor for detecting and predicting windshear in the vicinity of an urban airport, dynamic range requirements can exceed 50 dB because of high clutter to signal ratios. This presentation describes signal processing considerations in the presence of distributed and/or discrete clutter interference. Previous analyses have considered conventional range cell processing of radar returns from a rigidly mounted radar platform using either the Fourier or the pulse-pair method to estimate average windspeed and windspeed variation within a cell. Clutter rejection has been based largely upon analyzing a particular environment in the vicinity of the radar and employing a variety of techniques to reduce interference effects including notch filtering, Fourier domain line editing, and use of clutter maps. For the airborne environment the clutter characteristics may be somewhat different. Conventional clutter rejection methods may have to be changed and new methods will probably be required to provide useful signal to noise ratios. Various considerations are described. A major thrust has been to evaluate the effect of clutter rejection filtering upon the ability to derive useful information from the post-filter radar data. This analysis software is briefly described. Finally, some ideas for future analysis are considered including the use of adaptive filtering for clutter rejection and the estimation of windspeed spatial gradient directly from radar returns as a means of reducing the effects of clutter on the determination of a windshear hazard.

Baxa, Ernest G., Jr.↗

DOE OSTI2020

Devising novel methods for the controlled synthesis with morphology and size control of scintillator materials

Nanoscale scintillators represent an important class of materials underlying the functioning of high resolution, high ‘light yield’ radiation detection equipment, applicable for medicine in addition to security and environmental monitoring. This current study has focused on the controlled synthesis of Cs 2 ZnX 4 (X = Cl, Br) not only due to its very fast and promising radioluminescence decay time but also because no reliable synthesis protocols have been reported for producing nanoscale motifs of this material, to date. Solution-based methods can be used to synthesize products, possessing different, discrete morphologies with the ability to dictate particle size and shape. Specifically, we have demonstrated that a hot injection technique can be used to fabricate nanoscale rods and plates. Using this method, we systematically studied the effect of varying reaction variables, such as (i) reaction temperature, (ii) surfactants, and (iii) reaction time. As synthesized nanorods of Cs 2 ZnCl 4 possessed lengths of 101.3 ± 24.4 nm with widths of 17.3 ± 4.5 nm, whereas Cs 2 ZnBr 4 , created under identical conditions, were characterized by lengths of 206 ± 50 nm and widths of 26.7 ± 12.6 nm. Ligand-assisted re-precipitation (LARP) was used to produce the corresponding micron-sized particles. Under these conditions, we successfully synthesized spindles and spherical particles of Cs 2 ZnCl 4 with sizes of 1.58 ± 0.16 μm and 1.03 ± 0.34 μm, respectively, whereas plates of Cs 2 ZnBr 4 were prepared with average sizes of 0.76 ± 0.21 μm. High resolution synchrotron mapping, using the hard X-ray nanoprobe at NSLS II, was used to confirm the expected spatial distribution of Cs, Zn, and Cl within the as-produced crystals. Nanorods of Cs 2 ZnCl 4 and Cs 2 ZnBr 4 were subsequently tested for scintillation light yield and decay times. Preliminary results showed that Cs 2 ZnCl 4 maintained a light yield of 100–300 photons per MeV with a fast decay component of 2.6 ns. Finally, the Cs 2 ZnBr 4 analogue did not give rise to any perceptible scintillation emission.

36 MATERIALS SCIENCE↗

NASA NTRS1994

Wavelet compression techniques for hyperspectral data

Hyperspectral sensors are electro-optic sensors which typically operate in visible and near infrared bands. Their characteristic property is the ability to resolve a relatively large number (i.e., tens to hundreds) of contiguous spectral bands to produce a detailed profile of the electromagnetic spectrum. In contrast, multispectral sensors measure relatively few non-contiguous spectral bands. Like multispectral sensors, hyperspectral sensors are often also imaging sensors, measuring spectra over an array of spatial resolution cells. The data produced may thus be viewed as a three dimensional array of samples in which two dimensions correspond to spatial position and the third to wavelength. Because they multiply the already large storage/transmission bandwidth requirements of conventional digital images, hyperspectral sensors generate formidable torrents of data. Their fine spectral resolution typically results in high redundancy in the spectral dimension, so that hyperspectral data sets are excellent candidates for compression. Although there have been a number of studies of compression algorithms for multispectral data, we are not aware of any published results for hyperspectral data. Three algorithms for hyperspectral data compression are compared. They were selected as representatives of three major approaches for extending conventional lossy image compression techniques to hyperspectral data. The simplest approach treats the data as an ensemble of images and compresses each image independently, ignoring the correlation between spectral bands. The second approach transforms the data to decorrelate the spectral bands, and then compresses the transformed data as a set of independent images. The third approach directly generalizes two-dimensional transform coding by applying a three-dimensional transform as part of the usual transform-quantize-entropy code procedure. The algorithms studied all use the discrete wavelet transform. In the first two cases, a wavelet transform coder was used for the two-dimensional compression. The third case used a three dimensional extension of this same algorithm.

Evans, Bruce↗

NASA NTRS1997

Spatial layout affects speed discrimination

We address a surprising result in a previous study of speed discrimination with multiple moving gratings: discrimination thresholds decreased when the number of stimuli was increased, but remained unchanged when the area of a single stimulus was increased [Verghese & Stone (1995). Vision Research, 35, 2811-2823]. In this study, we manipulated the spatial- and phase relationship between multiple grating patches to determine their effect on speed discrimination thresholds. In a fusion experiment, we merged multiple stimulus patches, in stages, into a single patch. Thresholds increased as the patches were brought closer and their phase relationship was adjusted to be consistent with a single patch. Thresholds increased further still as these patches were fused into a single patch. In a fission experiment, we divided a single large patch into multiple patches by superimposing a cross with luminance equal to that of the background. Thresholds decreased as the large patch was divided into quadrants and decreased further as the quadrants were maximally separated. However, when the cross luminance was darker than the background, it was perceived as an occluder and thresholds, on average, were unchanged from that for the single large patch. A control experiment shows that the observed trend in discrimination thresholds is not due to the differences in perceived speed of the stimuli. These results suggest that the parsing of the visual image into entities affects the combination of speed information across space, and that each discrete entity effectively provides a single independent estimate of speed.

NASA Discipline Neuroscience↗

DOE OSTI2021

Visualization of Plasmonic Couplings Using Ultrafast Electron Microscopy

Hybrids of graphene and metal plasmonic nanostructures are promising building blocks for applications in optoelectronics, surface-enhanced scattering, biosensing, and quantum information. An understanding of the coupling mechanism in these hybrid systems is of vital importance to its applications. Previous efforts in this field mainly focused on spectroscopic studies of strong coupling within the hybrids with no spatial resolution. In this study, we report direct imaging of the local plasmonic coupling between single Au nanocapsules and graphene step edges at the nanometer scale by photon-induced near-field electron microscopy in an ultrafast electron microscope for the first time. The proximity of a step in the graphene to the nanocapsule causes asymmetric surface charge density at the ends of the nanocapsules. Computational electromagnetic simulations confirm the experimental observations. The results reported here indicate that this hybrid system could be used to manipulate the localized electromagnetic field on the nanoscale, enabling promising future plasmonic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2022

Neutron diffusion calculation in heterogeneous geometry based on local/global iteration using proper orthogonal decomposition

This study newly proposes a heterogeneous core calculation method based on local/global iteration using proper orthogonal decomposition (POD). By using the singular value decomposition (SVD) and the low-rank approximation, appropriate POD bases for expanding the neutron flux can be obtained from snapshot data of the neutron flux obtained by fine mesh calculations. By projection using the POD bases, the dimension of the target equation (e.g., discretized neutron diffusion equation) can be dramatically reduced. In the proposed method, POD is effectively applied to each single assembly calculation (local calculation). Furthermore, using the local/global iteration, the effective neutron multiplication factor and the neutron flux distribution in the whole core geometry can be obtained by combining the numerical results of the local calculation for each fuel assembly and the global calculation for the whole core. As a feasibility study, the proposed method is applied to a one-dimensional heterogeneous core analysis, and the accuracy is investigated by changing the total number of POD bases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2010

NASA's Advancements in Space-Based Spectrometry Lead to Improvements in Weather Prediction and Understanding of Climate Processes

AIRS (Atmospheric Infra-Red Sounder), was launched, in conjunction with AMSU-A (Advanced Microwave Sounding Unit-A) on the NASA polar orbiting research satellite EOS (Earth Observing System) Aqua satellite in May 2002 as a next generation atmospheric sounding system. Atmospheric sounders provide information primarily about the vertical distribution of atmospheric temperature and water vapor distribution. This is achieved by measuring outgoing radiation in discrete channels (spectral intervals) which are sensitive primarily to variations of these geophysical parameters. The primary objectives of AIRS/AMSU were to utilize such information in order to improve the skill of numerical weather prediction as well as to measure climate variability and trends. AIRS is a multi-detector array grating spectrometer with 2378 channels covering the spectral range 650/cm (15 microns) to 2660/cm (3.6 microns) with a resolving power (i/a i) of roughly 1200 where a i is the spectral channel bandpass. Atmospheric temperature profile can be determined from channel observations taken within the 15 micron (the long-wave CO2 absorption band) and within the 4.2 micron (the short-wave CO2 absorption band). Radiances in these (and all other) spectral intervals in the infrared are also sensitive to the presence of clouds in the instrument?s field of view (FOV), which are present about 95% of the time. AIRS was designed so as to allow for the ability to produce accurate Quality Controlled atmospheric soundings under most cloud conditions. This was achieved by having 1) extremely low channel noise values in the shortwave portion of the spectrum and 2) a very flat spatial response function within a channel?s FOV. IASI, the high spectral resolution IR interferometer flying on the European METOP satellite, does not contain either of these important characteristics. The AIRS instrument was also designed to be extremely stabile with regard to its spectral radiometric characteristics, which is critical with regard to the ability to measure accurate long term trends.

Susskind, Joel↗

DOE OSTI2019

Experimental and Computational Analysis of NEAMS Pebble Bed Reactors

The Advanced High-Temperature Pebble Bed Reactor concept leverages a particle-based fuel format consisting of discrete spherical graphite pebbles arrayed in a packed bed architecture. Thermal regulation achieved via flow of gas (e.g., helium) or liquid (e.g., molten salt) coolants through the void spaces between pebbles in the bed (characteristic pebble diameters: ~ 6 cm, gas cooled; ~ 3 cm, liquid cooled). This design is of interest owing to inherent advantages including passive safety features and highly efficient heat transfer characteristics that enable high power densities to be achieved. Fully leveraging these advantages requires a fundamental understanding of the complex coolant flow structure within the pebble-shaped fuel bed. Unfortunately, complexities associated with non-invasively probing the micro-scale flow phenomena within the complex network of void spaces between pebbles have hindered efforts to characterize the underlying transport phenomena. Without this fundamental foundation, it becomes challenging to rationally construct macroscopic transport models based upon volume averages of micro-scale parameters that accurately capture the observed flow behavior. Consequently, this lack of fundamental knowledge has limited previous work to coarse comparisons between theory and extremely limited experimental data. Moreover, these comparisons have been made only in terms of integral or macroscopic results. Therefore, multipoint measurements with a high level of spatial and temporal resolution are critically needed to fully map the complex flow patterns and to provide data at high spatial density to permit accurate volume averaging in the pebble bed. For this project, we performed a coordinated experimental and computational effort to quantitatively map the full-field velocity and temperature fields in the interstitial spaces within the pebble bed. The major accomplishments of this project included the important knowledge acquired for the flow mixing and heat transfer in the interstitial pore scales, the effects of flow mixing in the near-wall and far-wall regions, and the by-pass flows of pebble bed systems. The produced high-resolution experimental data (not broadly available in existing literature) were then used for the numerical studies enabling the refinement and validation of standard simulation tools (Nek5000).

42 ENGINEERING↗

DOE OSTI2025

Bottomonium properties in the quark-gluon plasma from a lattice-QCD informed 𝑇-matrix approach

Recent computations of bottomonium correlation functions with extended sources in lattice-discretized quantum chromodynamics (lQCD) provide new insights into heavy-quark dynamics at distance scales which are of the order of the inverse temperature. We analyze these results employing the thermodynamic 𝑇-matrix framework, in a continued effort to interpret lQCD data for quarkonium correlation functions in a nonperturbative and self-consistently solved quantum-many-body approach to a strongly coupled quark-gluon plasma (QGP). Its key inputs are the in-medium driving kernel (potential) of the scattering equation and an interference function which implements 3-body effects in the quarkonium coupling to the thermal medium. A simultaneous description of lQCD results for the bottomonium correlators with extended operators and the previously analyzed Wilson line correlators only requires minor refinements of the potential but calls for stronger interference effects at larger separation of the bottom quark and antiquark. We then analyze the poles of the self-consistent 𝑇 matrices on the real axis to assess the survival of the various bound states. Here, we estimate the pertinent temperatures where the poles disappear for the various bottomonium states and discuss the relation to the corresponding peaks in the bottomonium spectral functions. We also recalculate the spatial diffusion coefficient of the QGP and find it to be similar to that in our previous study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRSDate not supplied

Joint Acceptance Attenuation Factor of Integrated Pressure with Unsteady Pressure-Sensitive Paint Measurements

The Unsteady Pressure-Sensitive Paint (uPSP) is widely used to measure the surface pressure of scaled models in wind tunnel tests. Compared to the conventional pressure transducers, uPSP has the advantage of high spatial resolution. With multiple high-speed Complementary Metal Oxide Superconductor (CMOS) cameras, the uPSP data collected with the camera pixels are mapped to the surface grid of the scaled model and converted to pressure. The shot noise is the dominant component of the noise in the uPSP measurement. The integrated pressure is usually computed on the grid nodes of a user-defined patch. The effect of shot noise is reduced in the integrated pressure on the patch; however, the measurement of the aerodynamic pressure may also be attenuated by the decorrelation of the flow pressure field being measured. This paper discusses the Joint Acceptance Attenuation Factor (JAAF) of the integrated pressure with uPSP measurements. The JAAF, a function of frequency, is defined as the ratio of the Power Spectral Density (PSD) of the integrated aerodynamic pressure on the patch to the average PSD of the aerodynamic pressure on the grid nodes of the patch. In this paper, the JAAF is investigated for the integrated pressure on rectangular patches, whose edges are defined in the direction along the streamline or across the streamline. Based on the assumption that the surface pressure field can be described by the Corcos model, the closed-form formulas to compute the JAAFs of the integrated pressure on the discrete grid nodes of a rectangular patch and over the continuous area of a rectangular patch are derived respectively. It is shown that the JAAF of the integrated pressure over the continuous area of a rectangular patch is the limit of that on discrete nodes of the rectangular patch when the number of nodes in each row or column goes to infinity. The closed-form formulas of the JAAF derived in this paper, with estimated parameters of the model, are verified with the measurements of the uPSP and the conventional pressure transducer collected in the Space Launch System Ascent Unsteady Aerodynamics Test at NASA Ames Research Center in November 2017. The closed-form formulas of the JAAF of the integrated pressure on the rectangular patches, based on the Corcos model, provide an efficient method to estimate the attenuation of integration by the decorrelation of the flow pressure field and set references for the comparison of the spectrum of the integrated uPSP measurements and the conventional pressure transducer measurements. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding was provided by the NASA Aerosciences Evaluation and Test Capabilities Portfolio Office.

acoustics↗

NASA NTRSDate not supplied

Joint Acceptance Attenuation Factor of Integrated Pressure with Unsteady Pressure-Sensitive Paint Measurements

acoustics↗

NASA NTRSDate not supplied

Joint Acceptance Attenuation Factor of Integrated Pressure with Unsteady Pressure-Sensitive Paint Measurements

The Unsteady Pressure-Sensitive Paint (uPSP) is widely used to measure the surface pressure of the scaled model in the wind tunnel tests. Compared to the conventional pressure transducers, uPSP has the advantage of high spatial resolution. With multiple high-speed cameras, the uPSP data collected with the camera pixels are mapped to the surface grid of the scaled model and converted to pressure. The shot noise, which is associated with the quantum processes both in the generation of the photons by the luminescent uPSP and in the conversion of the photons into electrons within the camera, is the dominant component of the noise in the uPSP measurement. The integrated pressure is usually computed as an area-weighted average of the uPSP measured pressures on the grid nodes of a user-defined patch. The effect of the shot noise is reduced in the integrated pressure; however, the measurement of the aerodynamic pressure is also attenuated. This paper discusses the Joint Acceptance Attenuation Factor (JAAF) of the integrated pressure with uPSP measurements. The JAAF is defined as the ratio of the auto or cross Power Spectral Density (PSD) of the integrated aerodynamic pressure on the patch to the average auto PSD of the aerodynamic pressure on the grid nodes of the patch. Under the assumption that the cross spectrum of the aerodynamic pressure can be described with the Corcos model, the closed-form formula of the JAAF is derived for the integrated pressure on rectangular patches via the summation of uPSP measurements on discrete grid nodes. The closed-form formula is verified via the integration of the pressure field of Corcos model over continuous areas. The properties of the JAAF with different parameters are presented and the applications in the comparison of the spectrum of the pressure transducers and uPSP measurements, which were collected in the NASA Space Launch System (SLS) tests with the 11-by-11-foot transonic wind tunnel at NASA Ames Research Center in December 2017 and September 2019 respectively, are also discussed. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding was provided by the NASA Aerosciences Evaluation and Test Capabilities(AETC) Portfolio Office.

Unsteady Pressure-Sensitive Paint↗

NASA NTRS2011

Stereoscopic Planar Laser-Induced Fluorescence Imaging at 500 kHz

A new measurement technique for obtaining time- and spatially-resolved image sequences in hypersonic flows is developed. Nitric-oxide planar laser-induced fluorescence (NO PLIF) has previously been used to investigate transition from laminar to turbulent flow in hypersonic boundary layers using both planar and volumetric imaging capabilities. Low flow rates of NO were typically seeded into the flow, minimally perturbing the flow. The volumetric imaging was performed at a measurement rate of 10 Hz using a thick planar laser sheet that excited NO fluorescence. The fluorescence was captured by a pair of cameras having slightly different views of the flow. Subsequent stereoscopic reconstruction of these images allowed the three-dimensional flow structures to be viewed. In the current paper, this approach has been extended to 50,000 times higher repetition rates. A laser operating at 500 kHz excites the seeded NO molecules, and a camera, synchronized with the laser and fitted with a beam-splitting assembly, acquires two separate images of the flow. The resulting stereoscopic images provide three-dimensional flow visualizations at 500 kHz for the first time. The 200 ns exposure time in each frame is fast enough to freeze the flow while the 500 kHz repetition rate is fast enough to time-resolve changes in the flow being studied. This method is applied to visualize the evolving hypersonic flow structures that propagate downstream of a discrete protuberance attached to a flat plate. The technique was demonstrated in the NASA Langley Research Center s 31-Inch Mach 10 Air Tunnel facility. Different tunnel Reynolds number conditions, NO flow rates and two different cylindrical protuberance heights were investigated. The location of the onset of flow unsteadiness, an indicator of transition, was observed to move downstream during the tunnel runs, coinciding with an increase in the model temperature.

Medford, Taylor L.↗

DOE OSTI2024

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗

DOE OSTI2022

Development of a high fidelity model of the CROCUS experimental reactor

Measurements of scalar flux distributions with fine spatial and energy resolutions are needed to remedy one of the validation shortcomings of the novel neutronics full core solvers, such as MPACT and nTRACER. While a very detailed resolution of the flux can be calculated with such codes, only a limited experimental data set is available to check their accuracy. Such type of measurements are on-going at the zero power reactor CROCUS, operated at the Laboratory for Reactor Physics and System Behaviours of the EPFL, thanks to the development of advanced miniature neutron detection systems. This kind of experimental data would provide the community with a suitable benchmark for the validation of high fidelity neutronics solvers. In parallel, a multi-physics solver for steady-state and transient analysis of nuclear reactors, named GeN-Foam, has been developed. Based on the finite-volume OpenFOAM library, GeN-Foam provides us with enough flexibility to analyze non-conventional reactor geometries such as that of CROCUS. While CROCUS heterogeneities cannot be modeled by MPACT and nTRACER for the moment, GeN-Foam offers a unique opportunity to build a high fidelity model which mimics these codes' method to reach sub-pin simulation resolution. This document aims at describing the work achieved to get from the existing GeN-Foam model of the CROCUS reactor based on a structured mesh and using the neutron diffusion, the first high-fidelity model using discrete ordinates method as an approximation to neutron transport and an unstructured mesh for inter lattice water gap description and sub-pin heterogeneous modeling. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS2010

Hypersonic Laminar Boundary Layer Velocimetry with Discrete Roughness on a Flat Plate

Laminar boundary layer velocity measurements are made on a 10-degree half-angle wedge in a Mach 10 flow. Two types of discrete boundary layer trips were used to perturb the boundary layer gas. The first was a 2-mm tall, 4-mm diameter cylindrical trip. The second was a scaled version of the Orbiter Boundary Layer Transition (BLT) Detailed Test Objective (DTO) trip. Both 1-mm and 2.5-mm tall BLT DTO trips were tested. Additionally, side-view and plan-view axial boundary layer velocity measurements were made in the absence of these tripping devices. The free-stream unit Reynolds numbers tested for the cylindrical trips were 1.7x10(exp 6)/m and 3.3x10(exp 6)/m. The free-stream unit Reynolds number tested for the BLT DTO trips was 1.7x10(exp 6)/m. The angle of attack was kept at approximately 5-degrees for most of the tests resulting in a Mach number of approximately 8.3. These combinations of unit Reynolds numbers and angle of attack resulted in laminar flowfields. To study the precision of the measurement technique, the angle of attack was varied during one run. Nitric-oxide (NO) molecular tagging velocimetry (MTV) was used to obtain averaged axial velocity values and associated uncertainties. These uncertainties are as low as 20 m/s. An interline, progressive scan CCD camera was used to obtain separate images of the initial reference and shifted NO molecules that had been tagged by the laser. The CCD configuration allowed for sub-microsecond sequential acquisition of both images. The maximum planar spatial resolution achieved for the side-view velocity measurements was 0.07-mm in the wall-normal direction by 1.45-mm in the streamwise direction with a spatial depth of 0.5-mm. For the plan-view measurements, the maximum planar spatial resolution in the spanwise and streamwise directions was 0.69-mm by 1.28-mm, respectively, with a spatial depth of 0.5-mm. Temperature sensitive paint (TSP) measurements are provided to compliment the velocity data and to provide further insight into the behavior of the boundary layers. The experiments were performed at the NASA Langley Research Center 31-Inch Mach 10 Air tunnel.

Bathel, Brett↗