Search NASA⌕ Search

SEARCH · Search NASA

Results for “lateral force”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamics of streamflow permanence in a headwater network: Insights from catchment-scale model simulations

The hillslope and channel dynamics that govern streamflow permanence in headwater systems have important implications for ecosystem functioning and downstream water quality. Recent advancements in process-based, semi-distributed hydrologic models that build upon empirical studies of streamflow permanence in well-monitored headwater catchments show promise for characterizing the dynamics of streamflow permanence in headwater systems. However, few process-based models consider the continuum of hillslope-stream network connectivity as a control on streamflow permanence in headwater systems. The objective of this study was to expand a process-based, catchment-scale hydrologic model to better understand the spatiotemporal dynamics of headwater streamflow permanence and to identify controls of streamflow expansion and contraction in a headwater network. Further, we aimed to develop an approach that enhanced the fidelity of model simulations, yet required little additional data, with the intent that the model might be later transferred to catchments with limited long-term and spatially explicit measurements. This approach facilitated network-scale estimates of the controls of streamflow expansion and contraction, albeit with higher degrees of uncertainty in individual reaches due to data constraints. Our model simulated that streamflow permanence was highly dynamic in first-order reaches with steep slopes and variable contributing areas. The simulated stream network length ranged from nearly 98±2% of the geomorphic channel extent during wet periods to nearly 50±10% during dry periods. The model identified a discharge threshold of approximately 1 mm d –1 , above which the rate of streamflow expansion decreases by nearly an order of magnitude, indicating a lack of sensitivity of streamflow expansion to hydrologic forcing during high-flow periods. Overall, we demonstrate that process-based, catchment-scale models offer important insights on the controls of streamflow permanence, despite uncertainties and limitations of the model. We encourage researchers to increase data collection efforts and develop benchmarks to better evaluate such models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Deterministic Fabrication of Large-Area, High-Crystallinity Oxide Moiré Superlattices

Oxide twistronics extends moiré engineering beyond van der Waals materials, offering a promising platform for accessing emergent interfacial phenomena arising from the strong coupling of lattice, charge, and orbital degrees of freedom in complex oxides. However, deterministic fabrication of high-crystallinity oxide moiré superlattices over large lateral dimensions remains challenging due to the three-dimensional bonding network of oxides. Here, we demonstrate a scalable, generalized fabrication strategy that enables the formation of high-crystallinity oxide moiré superlattices with clean, chemically bonded interfaces and precisely controlled twist angles down to nominal values of 0.1°, achieving subdegree twist-angle accuracy across large contiguous lateral dimensions approaching the millimeter scale. Using NaNbO3 as a model system, we show that the resulting interlayer coupling drives pronounced structural reconstruction that modifies both the phase structure and ferroelectric domain configuration. Synchrotron-based X-ray 3D reciprocal space mapping reveals the emergence of a single-phase state in twisted bilayers, in contrast to the mixed-phase structure observed in single-layer membranes prior to twist assembly. The structural signatures are further consistent with gradual lattice rotation distributed along the thickness direction that may accommodate interfacial shear strain, distinct from reconstruction observed in van der Waals moiré systems which primarily occurs through in-plane stacking rearrangement. This collective lattice response is correlated with twist-dependent nanoscale electromechanical modulations observed by piezoresponse force microscopy. These results establish a scalable materials platform for oxide twistronics and support the implementation of twist-engineered functionalities in practical, macroscale device architectures.

Ghanbari, Reza [North Carolina State University (N↗
DOE OSTI2025

Unprecedented Beaufort Sea ice loss in late summer 2021 and its relationship to an extended period of unusually stormy weather

Previous case studies have linked cyclone-induced atmospheric forcing and/or upper-ocean processes to notable Arctic sea ice loss events in the summers of 2012 and 2016. This study examines a more recent and noteworthy case in late summer 2021 in which substantial sea ice loss followed a period of surface meteorological extremes in the Beaufort Sea region of the Arctic. We focus on the period from mid-August to mid-September 2021 that coincided with the Office of Naval Research THINICE Pilot Field Campaign and investigate stormy and windy conditions with respect to air-sea processes impacting sea ice conditions. We find that during the stormy first half of the campaign, cyclone-induced energy fluxes into the marginal ice zone and surrounding waters preconditioned the ice pack for more rapid melt later in the campaign. The second half of the campaign, in contrast, was marked by non-cyclone wind events that enhanced turbulent (namely sensible) heat fluxes into the ice and upper ocean that increased melt. Moreover, this latter period had enhanced advection of the Beaufort Sea ice pack into above-freezing waters, increasing bottom melt to >1 cm d −1 over the remainder of the campaign. While findings are shown to vary by surface type and at relatively small (i.e. ice-floe) scales, insights are offered on the roles of late summer coupled processes on rapid ice loss events in today’s Arctic environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effects of 5-Ion Beam Irradiation and Hindlimb Unloading on Metabolic Pathways in Plasma and Brain of Behaviorally Tested WAG/Rij Rats

A limitation of simulated space radiation studies is that radiation exposure is not the only environmental challenge astronauts face during missions. Therefore, we characterized behavioral and cognitive performance of male WAG/Rij rats 3 months after sham-irradiation or total body irradiation with a simplified 5-ion mixed beam exposure in the absence or presence of simulated weightlessness using hindlimb unloading (HU) alone. Six months following behavioral and cognitive testing or 9 months following sham-irradiation or total body irradiation, plasma and brain tissues (hippocampus and cortex) were processed to determine whether the behavioral and cognitive effects were associated with long-term alterations in metabolic pathways in plasma and brain. Sham HU, but not irradiated HU, rats were impaired in spatial habituation learning. Rats irradiated with 1.5 Gy showed increased depressive-like behaviors. This was seen in the absence but not presence of HU. Thus, HU has differential effects in sham-irradiated and irradiated animals and specific behavioral measures are associated with plasma levels of distinct metabolites 6 months later. The combined effects of HU and radiation on metabolic pathways in plasma and brain illustrate the complex interaction of environmental stressors and highlights the importance of assessing these interactions.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2026

Characterization of lateral amorphous selenium photodetectors for low-photon and VUV detection at cryogenic temperatures

The performance of amorphous selenium (a-Se) as a cryogenic photodetector material is evaluated through a series of experiments using laterally structured devices operated in a custom optical test stand. These studies investigate the response of a-Se detectors to low-photon fluxes at high electric fields near avalanche conditions, the linearity of the photoconductive response over a wide dynamic range and the direct detection of narrowband 130 nm vacuum ultraviolet (VUV) illumination. At 87 K, matched-filter analysis shows reliable single-shot detection with efficiencies ≥80% and area under the curve (AUC) ≥ 0.85 using as few as ∼ 6800 incident 401 nm photons, corresponding to ∼ 3400 photons within field-active regions after accounting for geometric constraints. Measurements are performed at cryogenic temperatures using calibrated photon fluxes derived from a silicon photomultiplier reference and a characterized optical filter stack. Additional experiments using a tellurium-doped a-Se (a-SeTe) device explore the material's behavior under identical test conditions and demonstrate that avalanche is achievable in a-SeTe at cryogenic temperatures. The results demonstrate reproducible low-noise operation, VUV sensitivity and field-dependent gain behavior in a lateral a-Se architecture, representing the first reported observation of avalanche multiplication in laterally structured a-Se and a-SeTe devices at cryogenic temperatures. These findings support the potential integration of laterally structured a-Se devices into next-generation pixelated liquid-argon time projection chambers (TPCs) requiring scalable, high-field-compatible photon detection systems.

Amorphous selenium↗
DOE OSTI2024

Permafrost Region Greenhouse Gas Budgets Suggest a Weak CO 2 Sink and CH 4 and N 2 O Sources, But Magnitudes Differ Between Top-Down and Bottom-Up Methods

Large stocks of soil carbon (C) and nitrogen (N) in northern permafrost soils are vulnerable to remobilization under climate change. However, there are large uncertainties in present-day greenhouse gas (GHG) budgets. We compare bottom-up (data-driven upscaling and process-based models) and top-down (atmospheric inversion models) budgets of carbon dioxide (CO 2 ), methane (CH 4 ) and nitrous oxide (N 2 O) as well as lateral fluxes of C and N across the region over 2000–2020. Bottom-up approaches estimate higher land-to-atmosphere fluxes for all GHGs. Both bottom-up and top-down approaches show a sink of CO 2 in natural ecosystems (bottom-up: -29 (-709, 455), top-down: -587 (-862, -312) Tg CO 2 -C yr -1 ) and sources of CH 4 (bottom-up: 38 (22, 53), top-down: 15 (11, 18) Tg CH 4 -C y -1 ) and N 2 O (bottom-up: 0.7 (0.1, 1.3), top-down: 0.09 (-0.19, 0.37) Tg N 2 O-N yr -1 ). The combined global warming potential of all three gases (GWP-100) cannot be distinguished from neutral. Over shorter timescales (GWP-20), the region is a net GHG source because CH 4 dominates the total forcing. The net CO 2 sink in Boreal forests and wetlands is largely offset by fires and inland water CO 2 emissions as well as CH 4 emissions from wetlands and inland waters, with a smaller contribution from N 2 O emissions. Priorities for future research include the representation of inland waters in process-based models and the compilation of process-model ensembles for CH 4 and N 2 O. Discrepancies between bottom-up and top-down methods call for analyses of how prior flux ensembles impact inversion budgets, more and well-distributed in situ GHG measurements and improved resolution in upscaling techniques.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Tidal disruption discs formed and fed by stream–stream and stream–disc interactions in global GRHD simulations

When a star passes close to a supermassive black hole (BH), the BH’s tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code h-amr. We carry out the first grid-based simulation of a deep-penetration TDE (β = 7) with realistic system parameters: a black hole-to-star mass ratio of 10 6 , a parabolic stellar trajectory, and a non-zero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disc forms due to the dissipation of orbital energy with ~20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream–disc interactions near the pericentre. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 h and a flaring in the accretion rate. We also find that the disc is eccentric with mean eccentricity e ≈ 0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericentre. Although these partial intersections eject gas out of the orbital plane, an accretion disc still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disc formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Chemistry and Interfacial Structure Promoting Quasi-van der Waals Epitaxial Growth of WS 2 Nanosheets on Sapphire for Prospective Application in Field-Effect Transistors

How do chemical and structural modifications to the supporting crystal surface affect the subsequent van der Waals (vdW) or quasi(Q)-vdW epitaxial growth of 2D nanocrystals? Developing an atomic-scale picture of such an interfacial system is crucial for understanding its impact on the physical and chemical properties of the supported 2D materials. The elucidation of the interfacial structure and chemistry needed to promote the Q-vdW epitaxial growth of 2D tungsten disulfide (WS 2 ) nanocrystals contributes to the growth mechanism understanding, thus pushing forward the integration of such atomically thin semiconductors toward real field-effect transistor applications. In addition to an atomic-force microscopy top view, we showcase a combination of X-ray techniques for a top-to-bottom investigation of the complexities of the buried interface structures. Furthermore, this approach uses X-ray photoelectron spectroscopy, X-ray standing wave excited X-ray fluorescence, and crystal truncation rod scattering to produce a highly resolved chemical-state-specific 3D atomic map for the extended interface structure of WS 2 /α-Al 2 O 3 (001). Employing these detailed analysis methods, along with density functional theory to further refine the picoscale structure, we demonstrate how two different types of interface engineering during the pregrowth stage lead to significant differences in the chemical and structural modifications to the terminal surface of c-face sapphire, which in turn leads to substantial differences in the submonolayer growth of supported WS 2 2D nanocrystals in terms of lateral domain sizes, epitaxial registry, vdW gaps, and stability.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Operation Upshot-Knothole

A thundering boom and sheet of fire announced a dramatic, once in a lifetime event – the firing of a nuclear projectile from a 280 mm cannon. Travelling seven miles downrange, the projectile detonated 500 feet over the Nevada Test Site’s Frenchman Flat with a force of fifteen kilotons. Although the test was successful, the massive cannon, designed for tactical battlefield use, had a critical flaw. It was too large and cumbersome for mobile combat. The cannon test, codenamed Grable, was part of Operation Upshot-Knothole, the fourth test series conducted in Nevada. Beginning with the Annie test on March 17, 1953, and concluding with the Climax test on June 4th, the Operation consisted of three airdrops, seven tower detonations, and the one artillery shot. As with previous Nevada tests, the demarcation between the weapon development series, Upshot, and the DOD effects tests, Knothole, was blurred. All eleven tests generated data in support of weapon development as well as tactical battlefield applications. During Upshot-Knothole, the University of California Radiation Laboratory (UCRL), conducted its first ever tests, Ruth and Ray. Although neither test was successful, these failures, along with another a year later at Operation Castle, provided the new laboratory with design data leading to future successful tests. The Federal Civil Defense Administration (FCDA) participated in a test series with the general goal of preparing for societal recovery from a nuclear attack. Known as Operation Doorstep, the FCDA’s experiments evaluated the effects of a nuclear attack on food supplies, housing and urban structures, and even a copse of trees transplanted from nearby Mt. Charleston. A contingent of the press witnessed test Annie from an outcropping of rocks, subsequently named News Nob. The presence of the news media, along with the FCDA’s participation, marked a fundamental change in the Nevada test program - a change that that now included the general populace, whose survival was, of course, the rationale for testing. The Department of Defense continued and expanded its Desert Rock exercises of troop exercises, helicopter operations, and evaluation of blast damage to military aircraft, tanks, and guns. Annie, Nancy, Badger, Simon, Encore, and Grable involved larger numbers of personnel than the other five shots. Below is an edited and condensed account of Desert Rock activities compiled by the Defense Threat Reduction Agency.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Exploring the gravitational structure of the proton with the dilepton final state using the CLAS12 detector at Jefferson Lab

The Gravitational Form Factors (GFFs) give access to the internal distributions of mass, pressure and shear force inside the proton. They were considered experimentally unmeasurable for decades due to the very weak gravitational interaction [1]. However, the Generalized Parton Distributions (GPDs), which describe the correlations between the longitudinal momentum and the transverse position of the partons inside the nucleon, have lately been related to the GFFs. For the first time, this relation gives the opportunity to extract GFFs experimentally. In this talk, I will present two ways to access GFFs using data taken in 2018 by the CLAS12 detector with a 10.6 GeV electron beam impinging on a liquid-hydrogen target. First, I will present the first measurement of the Timelike Compton Scattering reaction (the hard photoproduction of a lepton pair), that gives access to the quark GFFs via the angular asymmetry of the electron/positron pair [2]. I will then present the current effort to extract the near-threshold J/? photoproduction cross section using the same dataset. This later measurement is expected to provide direct insight on the gluons GFFs of the proton. [1] H. PAGELS. Energy-Momentum Structure Form Factors of Particles. Phys.Rev. 144 (1966) 1250-1260 [2] P. Chatagnon. First Measurement of Timelike Compton Scattering. Phys.Rev.Lett. 127 (2021) 26, 262501

Chatagnon, Pierre↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2022

Magnetic cores with high reluctance differences in flux paths

Embodiment of the present invention includes a magnetic structure and a magnetic structure used in a direct current (DC) to DC energy converter. The magnetic structure has an E-core and a plate, with the plate positioned in contact or in near contact with the post surfaces of the E-core. The E-core has a base, a no-winding leg, a transformer leg, and an inductor leg. The no-winding leg, the transformer leg, and the inductor leg are perpendicular and magnetically in contact with the base. The plate is a flat slab with lateral dimensions generally larger than its thickness. The plate has a plate nose that overlaps a top no-winding leg surface of the no-winding leg with a no-winding gap area to form a no-winding gap with a no-winding gap reluctance. The plate also has a plate end that overlaps a top inductor leg surface of the inductor leg with an inductor gap area to form an inductor gap with an inductor gap reluctance. In some embodiments, e.g., where the duty cycle is less than 50 percent, the inductor gap reluctance will be designed to be less than the no-winding gap reluctance. In these cases, the majority of the magnetic flux that passes through the transformer leg will return through the inductor leg, instead of through the no-winding leg. The inductor and no-winding gap reluctances can he adjusted, so that the electromotive force applied to a charge passing through the inductor will partially cancel the electromotive force applied by the transformer secondary. The gap reluctance ratio can be defined, so that the difference in secondary and inductor electromotive forces is equal to the output voltage defined by an optimal no-ripple duty cycle. In this way no changing current is required through the inductor to create a dI/dt inductive voltage drop across the output inductor. Zero output current ripple is achieved.Various embodiments of the plate, plate shape, and no-winding leg are disclosed. These embodiments allow achieving a high ratio of no-winding gap reluctance to inductor gap reluctance, for practical, affordable magnetic material structures and aspect ratios. A high gap reluctance ratio enables zero output current ripple for the high transformer turns ratios that are needed to achieve high input to output voltage ratios. The embodiments therefore allow achieving low output current ripple for 48 V or higher input voltages, 1 V or lower output voltages, and high output currents.

Yao, Yuan↗
DOE OSTI2021

Holocene thinning of Darwin and Hatherton glaciers, Antarctica, and implications for grounding-line retreat in the Ross Sea

Chronologies of glacier deposits in the Transantarctic Mountains provide important constraints on grounding-line retreat during the last deglaciation in the Ross Sea. However, between Beardmore Glacier and Ross Island – a distance of some 600 km – the existing chronologies are generally sparse and far from the modern grounding line, leaving the past dynamics of this vast region largely unconstrained. We present exposure ages of glacial deposits at three locations alongside the Darwin–Hatherton Glacier System – including within 10 km of the modern grounding line – that record several hundred meters of Late Pleistocene to Early Holocene thickening relative to present. As the ice sheet grounding line in the Ross Sea retreated, Hatherton Glacier thinned steadily from about 9 until about 3 ka. Our data are equivocal about the maximum thickness and Mid-Holocene to Early Holocene history at the mouth of Darwin Glacier, allowing for two conflicting deglaciation scenarios: (1) ~500 m of thinning from 9 to 3 ka, similar to Hatherton Glacier, or (2) ~950 m of thinning, with a rapid pulse of ~600 m thinning at around 5 ka. We test these two scenarios using a 1.5-dimensional flowband model, forced by ice thickness changes at the mouth of Darwin Glacier and evaluated by fit to the chronology of deposits at Hatherton Glacier. The constraints from Hatherton Glacier are consistent with the interpretation that the mouth of Darwin Glacier thinned steadily by ~500 m from 9 to 3 ka. Rapid pulses of thinning at the mouth of Darwin Glacier are ruled out by the data at Hatherton Glacier. This contrasts with some of the available records from the mouths of other outlet glaciers in the Transantarctic Mountains, many of which thinned by hundreds of meters over roughly a 1000-year period in the Early Holocene. The deglaciation histories of Darwin and Hatherton glaciers are best matched by a steady decrease in catchment area through the Holocene, suggesting that Byrd and/or Mulock glaciers may have captured roughly half of the catchment area of Darwin and Hatherton glaciers during the last deglaciation. An ensemble of three-dimensional ice sheet model simulations suggest that Darwin and Hatherton glaciers are strongly buttressed by convergent flow with ice from neighboring Byrd and Mulock glaciers, and by lateral drag past Minna Bluff, which could have led to a pattern of retreat distinct from other glaciers throughout the Transantarctic Mountains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Sockeye Code Validation Against UNIST Heat Pipe - Poster - Internship 2025

Sockeye is an engineering level code being developed under the MOOSE framework for modeling heat pipes. A heat pipe is a sealed tube with a wick structure filled with a working fluid that transfers heat efficiently and passively via phase change of the working fluid. Heat applied to the evaporator end creates a pressure gradient which causes vapor to migrate to the condenser end where the vapor deposits its energy and condenses. Capillary force generated by the wick structure then draws the liquid back to the evaporator. A validation study was performed to compare the Sockeye code against data produced at the Ulsan National Institute of Standards and Technology for a slightly overfilled sodium heat pipe. The heat pipe was operated under natural convection cooling and a decreasing inactive length was observed. The experiment was modeled in Sockeye using the condenser pool model, the front non-condensable gas (NCG) model, and the mixture NCG model separately to reproduce the inactive length phenomenon. Also, a new capability was implemented in MOOSE to allow for conjugate heat transfer from the condenser based on the Churchill-Chu correlation for natural convection. The condenser pool model showed an increasing inactive length, demonstrating that the behavior was likely not caused by a pool of excess liquid. The front NCG model was able to show good agreement with the experiment, but it suffered convergence issues at the front. Finally, the mixture NCG model gave good results when the axial mesh was sufficiently refined. This work culminated in additions to the Sockeye documentation and a contribution to a journal article that will be published later. This poster is a summary of my work which I can take back to my university for presentation.

42 - ENGINEERING↗
DOE OSTI2023

Analysis of an Arctic cold air outbreak during autumn and related air mass transformations forced by surface changes and advection in higher altitudes

During marine cold air outbreaks (MCAOs), cold and dry Arctic air masses are transported from the central Arctic southward across the closed sea ice and much warmer open oceans. They experience significant transformations including a rapid heating and moistening, often leading to cloud formation. While intense wintertime MCAOs have been analyzed widely, the air mass transformations during other seasons have been studied sparsely. We address this gap by investigating an MCAO case observed in September 2020. To study the transformation processes, we combine the fifth generation of atmospheric reanalyses of the global climate (ERA5), trajectory calculations, as well as shipborne and airborne measurements. In the central Arctic, observations acquired from aboard the research vessel (RV) Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition characterized the initial state of the air mass over closed sea ice. Trajectories indicated the pathway the air mass took from RV Polarstern southward to the Fram Strait. For the first 24 h of the southbound drift, the air masses remained quasi-stationary. Then, still 15 h ahead of the marginal sea ice zone, differential advection across the boundary layer flow introduced humidity and clouds at higher altitudes between 1.5 and 2.5 km. ERA5-derived temperature and humidity tendencies indicated complex vertical interactions. Radiative cloud-top cooling, entrainment, and turbulence were significantly reduced in the lower and enhanced in the upper advected cloud layer. Eventually, the lower cloud deck dissipated. After this confluence of 2 different air masses, observations gathered by Polar 5 in Fram Strait as part of the MOSAiC Airborne observations in the Central Arctic campaign revealed cloudy, moist layers throughout the lowest 3.5 km and an increasing boundary layer height. Comparing the initial with the final state 48 h later, the largest net heating of +8 K was found close to the surface, yet the largest net moistening of +2.5 g kg -1 at an altitude of 1 km, as the initial profile was exceptionally dry here. We conclude that the observed air mass transformations were driven by the surface changes from sea ice to open ocean but additionally strongly impacted by the differential advection of clouds and moisture across the near-surface MCAO flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Curled-Skewed Wakes behind Yawed Wind Turbines Subject to Veered Inflow

This work presents a new engineering analytical model that predicts the effect of both the turbine yaw misalignment and the inflow wind veer on the wake flow distribution downwind of a wind turbine. To consider the veered inflow, two methods were examined. In the first method, the curled shape of the wake due to the yaw offset is initially modelled. The wake shape is then laterally skewed at each height due to the wind veer based on the assumption that the turbine wake is transported downstream by the incoming flow. The second method is a more realistic approach that accounts for the effect of wind veer on the wind velocity direction and the yaw angle seen by the wind turbine. This models the wake region in a local coordinate system defined based on the wind direction at each height. A coordinate transformation is then performed to represent the wake flow distribution in the global coordinate system attached to the ground. The results show that while the two methods provide similar outputs for small variations in the wind direction across the rotor, the difference becomes more evident with an increase in wind veer. High-fidelity simulations for a turbine subject to a neutral atmospheric boundary layer were employed to validate model predictions for different operating conditions.

17 WIND ENERGY↗
DOE OSTI2025

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗