SEARCH · Engineering Papers

Results for “numeric profiling”

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 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Development of an accelerator-based neutron source to prototype Mo-99 production, part I: A liquid LBE windowless target

In this article, Molybdenum-99 (Mo-99)’s decay product, technetium-99 (Tc-99 m), is one of the most critical isotopes for medical diagnostics. To provide U.S. domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. This study discusses the development of the accelerator-based neutron source. The high-energy electrons from the accelerator irradiate a liquid lead-bismuth eutectic (LBE) target to produce neutrons. Part I of this work focuses on numerical and experimental analysis towards the development of a liquid LBE windowless target. Unlike the existing windowless targets in literature, the current design creates a vertical free surface for a beam to irradiate. First, a hydrodynamic analysis of the LBE windowless target is performed. Simplified analytical calculations are assisted by 2D computational fluid dynamics (CFD) simulations to design the target, with the focus on eliminating recirculation zones and avoiding cavitation. With the optimized geometry, the experimental study is performed to investigate the flow hydrodynamics using liquid LBE. The experiments (1) compare pressure drop in the system to correlation predictions; (2) visualize the free surface liquid LBE flow from the beam view; (3) validate the LBE flow profile using temperature sensitive paint from the side view; and (4) validate the liquid LBE film thickness using gamma densitometer measurements. Second, the power handling capability of the designed windowless target is investigated. The divider plate in the current design is susceptible to overheating due to the thin LBE film in front. As LBE erosion and corrosion is likely to occur at an LBE velocity of 2.0 m/s and temperature above 500 °C, a power limit of 10 kW of beam power was established to prevent this corrosion from occurring, which is calculated by a Nusselt number correlation. The divider plate surface temperature at 10 kW agrees well with the 3D CFD simulation results. Part I demonstrates the fundamental physics in liquid LBE windowless target design and associated testing. A companion paper, Part II will demonstrate how to couple this windowless target into the Mo-99 production system, including an accelerator system operating under an ultra-high vacuum and the UTA cooled by water at room temperature.

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

Low-frequency shear Alfvén waves at DIII-D: Theoretical interpretation of experimental observations

The linear properties of the low-frequency shear Alfvén waves such as those associated with the beta-induced Alfvén eigenmodes (BAEs) and the low-frequency modes observed in reversed-magnetic-shear DIII-D discharges [W. Heidbrink et al., Nucl. Fusion 61, 066031 (2021)] are theoretically investigated and delineated based on the theoretical framework of the general fishbone-like dispersion relation (GFLDR). By adopting representative experimental equilibrium profiles, it is found that, even though both modes are predominantly of Alfvénic polarization, the low-frequency mode is a reactive unstable mode with weak coupling to the energetic particles, while the BAE involves a dissipative instability due to resonant excitation by the energetic ions. Thus, the low-frequency mode is more appropriately called a low-frequency Alfvén mode (LFAM). Moreover, the ascending frequency spectrum patterns of the experimentally observed BAEs and LFAMs can be theoretically reproduced by varying qmin and also be well interpreted based on the GFLDR. The present analysis illustrates the solid predictive capability of the GFLDR and its practical usefulness in enhancing the interpretative capability of both experimental and numerical simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Various facets of magnetic charge correlation: Micromagnetic and distorted-wave Born approximation simulations study

The emergent concept of the magnetic charge quasiparticle provides a new realm to study the evolution of magnetic properties in two-dimensional artificially frustrated magnets. Here we report on the exploration of magnetic phases due to various magnetic charge correlation using the complementary numerical techniques of micromagnetic and distorted-wave Born approximation simulations in artificial permalloy honeycomb lattice. The honeycomb element length varies between 10 nm and 100 nm, while the width and thickness are kept within the single domain limit. In addition to the charge ordered loop state, we observe disordered charge arrangement, characterized by the random distribution of ±Q charges, in single domain size honeycomb lattice. As the length of the honeycomb element increases, low multiplicity magnetic charges tend to form contiguous bands in thinner lattice. Thin honeycomb lattice with 100 nm element length exhibits a perfect spin ice pattern, which remains unaffected by the modest increase in the width of element size. We simulate scattering profiles under the pretext of distorted-wave Born approximation formalism for the micromagnetic phases. The results are expected to provide useful guidance in the experimental investigation of magnetic phases in an artificial honeycomb magnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Turbulence-Resolving Simulations of Hurricane Laura (2020): Insights Into Extreme Winds and Eyewall Turbulence

Extreme weather events pose significant risks to coastal and offshore energy infrastructure. In this work, we investigate the structure of mean winds and turbulence near the surface ( z < 300 m) that occurred during Hurricane Laura in 2020 on the US Gulf Coast. To this end, we perform turbulence-resolving simulations (..delta..x = 33.33 m) that encompass the entire eyewall of Hurricane Laura by introducing a novel method that we will refer to as Moving-Mesoscale to Static-LES Integrated Coupling (MOSAIC). The simulation results are validated against surface, boundary-layer, and mid-tropospheric observations. Our analysis of the extreme wind conditions near the surface indicates that the mean wind and turbulence profiles vary greatly at and near the eyewall, resulting in extreme values of mean wind speed (U > 50 m * s-1) and turbulence kinetic energy (k ~ 50 m2 * s-2) at altitudes above 50 m . To investigate in detail the nature of turbulence within the eyewall of the storm, we conduct higher-resolution simulations (..delta..x = 11.11 m) of a portion of the hurricane. We provide a comprehensive statistical description of turbulence in the eyewall region, highlighting the need to employ ..delta..x ~ 10 m grid spacing to capture the integral length-scales near the surface, second- and higher-order central moments (i.e., skewness and kurtosis), and spectral coherence in the flow accurately. These numerical simulations provide the most detailed description of mean winds and turbulent conditions within the eyewall of a historical tropical cyclone to date, illustrating how high-resolution simulations can enhance limited turbulence measurements in hurricanes and giving insight into the extreme conditions threatening offshore and coastal infrastructure.

17 WIND ENERGY↗

DOE OSTI2021

Numerical study to predict the effect of surface roughness on the thermal and hydraulic performance of additively manufactured heat exchangers

Heat exchangers are used in many industries to transfer thermal energy between two or more objects with an associated thermal and hydraulic performance indicated by thermal effectiveness and pressure drop respectively. Many investigations have aimed to enhance the thermal effectiveness of heat exchangers by introducing turbulent flow inducers that increase the heat transfer rate at the expense of increasing pressure drop. Flow inducers like surface roughness can be controlled with advanced manufacturing methods like additive manufacturing for heat exchanger production. Selecting and manufacturing an optimal surface roughness for heat exchangers requires a more comprehensive approach that accounts for the surface roughness effect on the thermal and hydraulic performance of a sample. Since fabricating samples with additive manufacturing to a desired surface roughness can be an extensive experimental process, this study will fundamentally investigate the effect of surface roughness on the thermal and velocity boundary layers for various surface roughness profiles with physics-based computational modeling. To accurately model the surface roughness effect on thermohydraulic performance, correlations like roughness ratio, friction factor, and Nusselt number are investigated. This study expects to indicate correlations for different flow regimes which can successfully predict the effect of surface roughness on thermohydraulic performance. Results generated from modeling surface roughness effect with unique correlations are then compared to experimental data obtained from literature for validation. The results of this study will provide insight to understand the thermal fluid physics associated with surface roughness and design guidelines for heat exchangers produced with additive manufacturing.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Prediction and uncertainty quantification of shale well performance using multifidelity Monte Carlo

Uncertainty quantification is an integral component of reservoir management, especially considering the inherent uncertainty in subsurface systems. While a standard practice to estimate the uncertainty, Monte Carlo (MC) simulation is computationally intense when the sampling population comprises high-fidelity simulations. Alternatively, the Multi-fidelity Monte Carlo (MFMC) simulation overcomes this computational intensity by integrating low- and high-fidelity simulations. Our goal is to minimize the number of expensive high-fidelity simulations while maintaining accuracy and using numerous fast and cheap low-fidelity simulations to efficiently sample to input parameter space of interest. We selected gas production from unconventional wells to demonstrate the potential speedups and accuracy of the MFMC approach. The model fidelity usually determines the trade-off between accuracy and efficiency. While the high-fidelity model is more accurate, the low-fidelity model is more efficient. Our high-fidelity simulation includes reservoir simulations of a hydraulically fractured well. On the other hand, our low-fidelity model comprises the parallel-plate flow model. We used differential programming to efficiently solve the 1D flow model, where automatic differentiation is used to efficiently compute the gradients. We matched the production profile of high-fidelity simulations with our low-fidelity simulations. Then, we used a support vector regression to map the high- and low-fidelity input parameters. The mapping function is essential to tune the low-dimensional parameter space of the low-fidelity model to the high-dimensional parameter space of the high-fidelity model. We found that we can use a combination of 9 high fidelity and 10,000 low fidelity simulations to efficiently and accurately simulate pressure management. This method is at least two orders of magnitude faster than only using high-fidelity simulations. Finally, from a broader perspective, MFMC could efficiently estimate the uncertainty of various systems and models, integrating low- and high-fidelity models.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2022

Ion dynamics driven by a strongly nonlinear plasma wake

Abstract In plasma wakefield accelerators, the wave excited in the plasma eventually breaks and leaves behind slowly changing fields and currents that perturb the ion density background. We study this process numerically using the example of a Facility for Advanced aCcelerator Experimental Tests (FACET) experiment where the wave is excited by an electron bunch in the bubble regime in a radially bounded plasma. Four physical effects underlie the dynamics of ions: (1) attraction of ions toward the axis by the fields of the driver and the wave, resulting in formation of a density peak, (2) generation of ion-acoustic solitons following the decay of the density peak, (3) positive plasma charging after wave breaking, leading to acceleration of some ions in the radial direction, and (4) plasma pinching by the current generated during the wave-breaking. The interplay of these effects results in the formation of various radial density profiles, which are difficult to produce in any other way.

Physics↗

DOE OSTI2020

Probing the nature of dark matter with accreted globular cluster streams

ABSTRACT The steepness of the central density profiles of dark matter (DM) in low-mass galaxy haloes (e.g. dwarf galaxies) is a powerful probe of the nature of DM. We propose a novel scheme to probe the inner profiles of galaxy subhaloes using stellar streams. We show that the present-day morphological and dynamical properties of accreted globular cluster (GC) streams – those produced from tidal stripping of GCs that initially evolved within satellite galaxies and later merged with the Milky Way (MW) – are sensitive to the central DM density profile and mass of their parent satellites. GCs that accrete within cuspy cold dark matter (CDM) subhaloes produce streams that are physically wider and dynamically hotter than streams that accrete inside cored subhaloes. A first comparison of MW streams ‘GD-1’ and ‘Jhelum’ (likely of accreted GC origin) with our simulations indicates a preference for cored subhaloes. If these results hold up in future data, the implication is that either the DM cusps were erased by baryonic feedback, or their subhaloes naturally possessed cored density profiles implying particle physics models beyond CDM. Moreover, accreted GC streams are highly structured and exhibit complex morphological features (e.g. parallel structures and ‘spurs’). This implies that the accretion scenario can naturally explain the recently observed peculiarities in some of the MW streams. We also propose a novel mechanism for forming ‘gaps’ in stellar streams when the remnant of the parent subhalo (which hosted the GC) later passes through the GC stream. This encounter can last a longer time (and have more of an impact) than the random encounters with DM subhaloes previously considered, because the GC stream and its parent subhalo are on similar orbits with small relative velocities. Current and future surveys of the MW halo will uncover numerous faint stellar streams and provide the data needed to substantiate our preliminary tests with this new probe of DM.

(Galaxy:) globular clusters: individual↗

DOE OSTI2021

A novel physics-based and data-supported microstructure model for part-scale simulation of laser powder bed fusion of Ti-6Al-4V

The elasto-plastic material behavior, material strength and failure modes of metals fabricated by additive manufacturing technologies are significantly determined by the underlying process-specific microstructure evolution. In this work a novel physics-based and data-supported phenomenological microstructure model for Ti-6Al-4V is proposed that is suitable for the part-scale simulation of laser powder bed fusion processes. The model predicts spatially homogenized phase fractions of the most relevant microstructural species, namely the stable β-phase, the stable α s -phase as well as the metastable Martensite α m -phase, in a physically consistent manner. In particular, the modeled microstructure evolution, in form of diffusion-based and non-diffusional transformations, is a pure consequence of energy and mobility competitions among the different species, without the need for heuristic transformation criteria as often applied in existing models. The mathematically consistent formulation of the evolution equations in rate form renders the model suitable for the practically relevant scenario of temperature- or time-dependent diffusion coefficients, arbitrary temperature profiles, and multiple coexisting phases. Due to its physically motivated foundation, the proposed model requires only a minimal number of free parameters, which are determined in an inverse identification process considering a broad experimental data basis in form of time-temperature transformation diagrams. Subsequently, the predictive ability of the model is demonstrated by means of continuous cooling transformation diagrams, showing that experimentally observed characteristics such as critical cooling rates emerge naturally from the proposed microstructure model, instead of being enforced as heuristic transformation criteria. Eventually, the proposed model is exploited to predict the microstructure evolution for a realistic selective laser melting application scenario and for the cooling/quenching process of a Ti-6Al-4V cube of practically relevant size. Numerical results confirm experimental observations that Martensite is the dominating microstructure species in regimes of high cooling rates, e.g., due to highly localized heat sources or in near-surface domains, while a proper manipulation of the temperature field, e.g., by preheating the base-plate in selective laser melting, can suppress the formation of this metastable phase.

Inverse parameter identification↗

DOE OSTI2026

Benchmark for two-dimensional large scale coherent structures in partially magnetized E × B plasmas—community collaboration & lessons learned

Low-temperature plasmas (LTPs) are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of LTP codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E × B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development.

benchmarking↗

DOE OSTI2025

Effect of pattern transfer process on roughness of block copolymer patterns from directed self-assembly

Block copolymer-directed self-assembly (DSA) remains promising for improving pattern quality and reducing the stochastic variations that challenge high numerical aperture extreme ultraviolet lithography. Equally critical is refining pattern transfer methods for accurately transferring the rectified DSA features to the underlying substrate. We compare two atomic layer deposition (ALD)-based techniques: sequential infiltration synthesis (SIS) and dry liftoff, applied to polystyrene-block-poly(methyl methacrylate) (PS-b-PMMA) DSA patterns. Both methods utilize aluminum oxide hard masks, with one synthesized through infiltration into the PMMA domains and the other through conformal ALD coating. High-resolution scanning electron micrographs were analyzed to measure line edge, width, and placement roughness for both the line (PMMA) and space (PS) features. Although both methods yielded similar overall 3σ rms roughness, they differed significantly in the frequency-dependent power spectral density (PSD) profiles. SIS reduced line placement roughness at length scales associated with the polymer pitch, but increased space width roughness at low frequencies, whereas dry liftoff mimicked the frequency content of the original guiding pattern. This study underscores the importance of PSD evaluation in selecting optimal pattern transfer strategies for specific applications.

Block copolymers↗

DOE OSTI2021

Simulation-aided characterization of a versatile water-based condensation particle counter for atmospheric airborne research

Capturing the vertical profiles and horizontal variations of atmospheric aerosols often requires accurate airborne measurements. With the advantage of avoiding health and safety concerns related to the use of butanol or other chemicals, water-based condensation particle counters have emerged to provide measurements under various environments. However, airborne deployments are relatively rare due to the lack of instrument characterization under reduced pressure at flight altitudes. This study investigates the performance of a commercial “versatile” water-based condensation particle counter (vWCPC, model 3789, TSI, Shoreview, MN, USA) under various ambient pressure conditions (500–920 hPa) with a wide range of particle total number concentrations (1500–70 000 cm -3 ). The effect of conditioner temperature on vWCPC 3789 performance at low pressure is examined through numerical simulation and laboratory experiments. We show that the default instrument temperature setting of 30 °C for the conditioner is not suitable for airborne measurement and that the optimal conditioner temperature for low-pressure operation is 27°. Under the optimal conditioner temperature (27°), the 7 nm cut-off size is also maintained. Additionally, we show that insufficient droplet growth becomes more significant under the low-pressure operation. The counting efficiency of the vWCPC 3789 can vary up to 20 % for particles of different chemical compositions (e.g., ammonium sulfate and sucrose particles). However, such variation is independent of pressure.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Estimates of Southern Hemispheric Gravity Wave Momentum Fluxes across Observations, Reanalyses, and Kilometer-Scale Numerical Weather Prediction Model

Abstract Gravity waves (GWs) are among the key drivers of the meridional overturning circulation in the mesosphere and upper stratosphere. Their representation in climate models suffers from insufficient resolution and limited observational constraints on their parameterizations. This obscures assessments of middle atmospheric circulation changes in a changing climate. This study presents a comprehensive analysis of stratospheric GW activity above and downstream of the Andes from 1 to 15 August 2019, with special focus on GW representation ranging from an unprecedented kilometer-scale global forecast model (1.4 km ECMWF IFS), ground-based Rayleigh lidar (CORAL) observations, modern reanalysis (ERA5), to a coarse-resolution climate model (EMAC). Resolved vertical flux of zonal GW momentum (GWMF) is found to be stronger by a factor of at least 2–2.5 in IFS compared to ERA5. Compared to resolved GWMF in IFS, parameterizations in ERA5 and EMAC continue to inaccurately generate excessive GWMF poleward of 60°S, yielding prominent differences between resolved and parameterized GWMFs. A like-to-like validation of GW profiles in IFS and ERA5 reveals similar wave structures. Still, even at ∼1 km resolution, the resolved waves in IFS are weaker than those observed by lidar. Further, GWMF estimates across datasets reveal that temperature-based proxies, based on midfrequency approximations for linear GWs, overestimate GWMF due to simplifications and uncertainties in GW wavelength estimation from data. Overall, the analysis provides GWMF benchmarks for parameterization validation and calls for three-dimensional GW parameterizations, better upper-boundary treatment, and vertical resolution increases commensurate with increases in horizontal resolution in models, for a more realistic GW analysis. Significance Statement Gravity wave–induced momentum forcing forms a key component of the middle atmospheric circulation. However, complete knowledge of gravity waves, their atmospheric effects, and their long-term trends are obscured due to limited global observations, and the inability of current climate models to fully resolve them. This study combines a kilometer-scale forecast model, modern reanalysis, and a coarse-resolution climate model to first compare the resolved and parameterized momentum fluxes by gravity waves generated over the Andes, and then evaluate the fluxes using a state-of-the-art ground-based Rayleigh lidar. Our analysis reveals shortcomings in current model parameterizations of gravity waves in the middle atmosphere and highlights the sensitivity of the estimated flux to the formulation used.

Meteorology & Atmospheric Sciences↗

DOE OSTI2025

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Single-cell quantification of ribosome occupancy in early mouse development

Translation regulation is critical for early mammalian embryonic development. However, previous studies had been restricted to bulk measurements, precluding precise determination of translation regulation including allele-specific analyses. Here, to address this challenge, we developed a novel microfluidic isotachophoresis (ITP) approach, named RIBOsome profiling via ITP (Ribo-ITP), and characterized translation in single oocytes and embryos during early mouse development. We identified differential translation efficiency as a key mechanism regulating genes involved in centrosome organization and N 6 -methyladenosine modification of RNAs. Our high-coverage measurements enabled, to our knowledge, the first analysis of allele-specific ribosome engagement in early development. These led to the discovery of stage-specific differential engagement of zygotic RNAs with ribosomes and reduced translation efficiency of transcripts exhibiting allele-biased expression. By integrating our measurements with proteomics data, we discovered that ribosome occupancy in germinal vesicle-stage oocytes is the predominant determinant of protein abundance in the zygote. The Ribo-ITP approach will enable numerous applications by providing high-coverage and high-resolution ribosome occupancy measurements from ultra-low input samples including single cells.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Numerical assessment of the new V-shape small-angle slot divertor on DIII-D

The small-angle slot (SAS) divertor of the DIII-D tokamak, and its upcoming upgrade, the V-shape small-angle slot (SAS-V) divertor, are numerically investigated using the SOLPS-ITER code package, including the effect of particle drifts, for a range of plasma density, heating power, strike point position in the slot, and for both magnetic field directions. The simulations show that the electron temperature near the strike point is reduced in SAS-V compared to SAS, for both magnetic field directions, such that SAS-V achieves divertor detachment at a lower value of the outboard mid-plane separatrix electron density. The detachment threshold is lower because the V-shape focuses recycling neutrals on the V-end, densifying and cooling the plasma in the slot. At sufficiently high density, the V-shape also reduces the radial gradient of the temperature profile at the target, which in turns reduces the radial electric field and the E x B drift velocities, further densifying and cooling the plasma in the slot and leading to detachment. The V-shape effect, however, is reduced for higher heating power. With more heating power, the detachment density increases, reducing the ionization mean free path of recycled neutrals, which therefore become less sensitive to target shape changes. Furthermore, this suggests that in a fusion reactor, where the heating power is high, optimization of the divertor target shape needs to be combined with other strategies to lower the detachment density, such as in-slot injection of low-Z impurities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Structure of Offshore Low-Level Jet Turbulence and Implications to Mesoscale-to-Microscale Coupling

This paper explores realistic nonstationary atmospheric boundary layer (ABL) turbulence arising from nonstationarity at the mesoscale, particularly within offshore low-level jets with implications to offshore wind farms, using high-fidelity multiscale large-eddy simulations (LES). To this end, we analyzed the single-point turbulence statistical structure of a North-Atlantic offshore LLJ event simulated using high-resolution LES (AMR-Wind). The nonstationary LLJ is simulated using a mesoscale-to-microscale coupled (MMC) simulation procedure involving data assimilation of mesoscale velocity and temperature data from the Weather Research and Forecasting (WRF) model. Unlike the assimilation of mesoscale velocity data into the LES, the direct assimilation of temperature profiles had a strong impact on turbulence stratification, thereby causing erroneous predictions of turbulence both above and within the jet layer. Various approaches to mitigate this effect have resulted in multiple (four) variants of this MMC strategy. Outcomes from this work clearly show that the turbulence within the low-level jet is a strong function of the MMC approach as the turbulence structure within the low-level jet is dependent on the flux of residual turbulence from outside the jet, which in turn depends on the temperature forcing history. Additionally, the turbulence predicted by all these different methods (as well as the observation data) show similar deviations from equilibrium as evidenced by comparisons with idealized atmospheric turbulence structure obtained using the same numerical method. In general, we observe that the predicted LLJ turbulence tends to differ from canonical ABL turbulence with comparable shear. Particularly, the combination of shear and turbulence observed in such nonstationary low-level turbulence cannot be matched using equilibrium settings and therefore, represents a critical use-case for both testing and leveraging meso–micro coupling strategies.

17 WIND ENERGY↗

DOE OSTI2023

COUPLING SMOOTHED PARTICLE HYDRODYNAMICS WITH FINITE ELEMENT METHOD TO SIMULATE RESIDUAL STRESSES FROM FRICTION STIR PROCESSING

Friction stir processing (FSP) is a solid-state material processing technique that locally modifies the microstructure but also induces undesirable residual stresses. A robust numerical model for the FSP can help in mitigating these residual stresses. Heat source models within a finite element method (FEM) framework suffer from inaccuracies. In contrast, smoothed particle hydrodynamics (SPH) model that explicitly captures the material flow near the tool and the associated heat generation are accurate. However, the computational expense of SPH simulations can be prohibitive. In this work, we propose a coupled SPH-FEM framework. SPH is used to model the heat generation accurately near the tool and which is then inserted into to the FEM model as a heat source. To verify this proposed coupling approach, a test case is set up with typical FSP conditions and it is modeled in both SPH and SPH-FEM. The temperatures profiles were compared after the simulations have reached steady-state temperatures. The similarity of the temperature profiles from SPH-FEA and SPH validated the proposed coupling approach. This proposed approach achieves the accuracy of the SPH method while potentially retaining the low computational expense of FEM.

smoothed particle hydrodynamics, Finite Element Me↗