Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “drag”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hyperdiffusion of dust particles in a turbulent tokamak plasma

The effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad-hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number $St\propto a_d/n_0$, with $a_d$ the dust radius and $n_0$ the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories $\Delta r^2$ is shown to depend on the ratio $Ku/St$, with $Ku\propto u_{rms}$ the Kubo number and $u_{rms}$ the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting $\Delta r^2\propto (Ku/St)^2t^3$, or $\Delta r^2\propto (u_{rms}n_0/a_d)^2t^3$. As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible, and drives the transition from a hyperdiffusive regime towards a ballistic one.

dust↗
DOE OSTI2026

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation) [SWR-26-095]

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation): Multifidelity aerodynamic polar data generation for hydrofoil/tidal-turbine airfoil sections. Foilpolars ties together three pieces: *AeroSandbox supplies the baseline airfoil coordinates (UIUC database). *G2Aero parameterizes those shapes on a Grassmannian manifold (Karcher mean + PGA basis) and samples new perturbed shapes around that basis. *XFoil (panel method) and NeuralFoil (neural-network surrogate, shipped with AeroSandbox) each solve the resulting shapes for lift, drag, moment, and pressure at the swept angles of attack, Reynolds numbers, and n_crit values. Design optimization of foil shapes in a computationally efficient way requires polars data across many candidate shapes, not just a handful of baseline foils. However, high-fidelity CFD at that scale is too costly, and naive shape perturbation strays from realistic geometries. FOILPOLARS addresses this by loading baseline airfoils (via AeroSandbox) and mapping them onto a Grassmannian manifold (via G2Aero), computing a Karcher mean and principal geodesic analysis (PGA) basis. New shapes are sampled by perturbing PGA coefficients, keeping them close to the manifold of realistic foils. Each sampled shape is evaluated across a configurable sweep of angle of attack, Reynolds number, and critical amplification factor using two solvers: XFoil (panel method) and NeuralFoil (neural-network surrogate), producing a paired dataset of lift, drag, moment, pressure, convergence, and confidence, indexed alongside each shape's PGA coefficients and shared Grassmannian basis in a single xarray dataset. From this, FOILPOLARS produces convergence summaries and comparison plots per shape, Reynolds number, and n_crit. A command-line interface exposes each pipeline stage independently, supporting data-driven design, optimization, and machine-learning workflows for foils.

Sandhu, Rimple [National Laboratory of the Rockies↗
DOE OSTI2019

Improving Wind Energy Forecasting through Numerical Weather Prediction Model Development

The primary goal of the Second Wind Forecast Improvement Project (WFIP2) is to advance the state-of-the-art of wind energy forecasting in complex terrain. To achieve this goal, a comprehensive 18-month field measurement campaign was conducted in the region of the Columbia River basin. The observations were used to diagnose and quantify systematic forecast errors in the operational High-Resolution Rapid Refresh (HRRR) model during weather events of particular concern to wind energy forecasting. Examples of such events are cold pools, gap flows, thermal troughs/marine pushes, mountain waves, and topographic wakes. WFIP2 model development has focused on the boundary layer and surface-layer schemes, cloud–radiation interaction, the representation of drag associated with subgrid-scale topography, and the representation of wind farms in the HRRR. Additionally, refinements to numerical methods have helped to improve some of the common forecast error modes, especially the high wind speed biases associated with early erosion of mountain–valley cold pools. This study describes the model development and testing undertaken during WFIP2 and demonstrates forecast improvements. Specifically, WFIP2 found that mean absolute errors in rotor-layer wind speed forecasts could be reduced by 5%–20% in winter by improving the turbulent mixing lengths, horizontal diffusion, and gravity wave drag. The model improvements made in WFIP2 are also shown to be applicable to regions outside of complex terrain. Ongoing and future challenges in model development will also be discussed.

17 WIND ENERGY↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
DOE OSTI2021

How Snow Aggregate Ellipsoid Shape and Orientation Variability Affects Fall Speed and Self-Aggregation Rates

Snow aggregate shapes and orientations have long been known to exhibit substantial variability. Despite this observed variability, most weather and climate prediction models use fixed power-law functions that deterministically map particle size to mass and fall speed. As such, integrated quantities like precipitation and self-aggregation rates currently ignore nonlinear effects resulting from variation in shape and orientation for aggregates of the same size. This study therefore develops an analytic framework that couples an empirically based bivariate distribution of ellipsoid shapes to classical hydrodynamic theory so as to capture an appropriate dispersion of masses, projected areas, and fall speeds for an assumed size distribution. For a fixed aggregate size, shape variations produce approximately ±0.13 m s -1 standard deviation of fall speed which increases the mass flux fall speed dispersion by more than 100% over traditional microphysics models. This increased fall speed dispersion results predominantly from shape-induced mass dispersion whereas orientation and drag dispersion play a lesser role. Shape variations can increase mass- and reflectivity-weighted fall speeds by up to 60% of traditional models whereas self-aggregation rates can increase by a factor of 100 for very small slope parameters. This implies that aggregate shape variations effectively forestall the theorized onset of fall speed distribution narrowing and subsequent quenching of the aggregation process. As a result, it is likely that secondary ice formation is necessary to prevent an ever decreasing slope parameter. The mathematical theory presented in this study is used to develop simple correction factors for snow forecast and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

The Unsteady Aerodynamics Module For FAST8

The new modularization framework of FAST v.8 (Jonkman 2013) was accompanied by a complete overhaul of the aerodynamics routines. AeroDyn is an aerodynamics module comprised of four submodels: rotor wake/induction, blade airfoil aerodynamics, tower influence on the blade nodes, and tower drag. Throughout the software overhaul, several improvements to the original theoretical treatments were achieved, including more accurate skewed-wake and unsteady airfoil aerodynamics, and the possibility of modeling highly flexible and curved blades. Under asymmetric conditions, such as wind shear, yawed, and tilted flow, the individual blade elements undergo variations in angle of attack that lead to unsteady aerodynamics phenomena, which can no longer be captured through the static airfoil lift and drag look-up tables. This document covers the main theory and the organization of the modularization framework of the new unsteady aerodynamics submodule (UAM)1, which includes unsteady aerodynamics under attached flow conditions and dynamic stall. The UAM can be called by either blade element momentum theory (only interface discussed in this report), or other future wake models (e.g., dynamic blade element momentum theory, or generalized dynamic wake).

17 WIND ENERGY↗
DOE OSTI2020

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Experimental and Simulation Analysis of Binary Mixtures of Biomass and Inert Material

A literature review was conducted to identify experiments on various scales of biomass combustion specific to circulating fluidized bed (CFB) combustors. Typical feedstock and inert bed material characteristics for biomass combustion in a CFB were gathered. This information was used to design laboratory-scale, non-reacting experiments on biomass-inert bed material mixtures in a 1-in. cylindrical test rig in the Multiphase Flow Analysis Laboratory (MFAL) in Morgantown, WV. The biomass chosen was a hardwood. Pellets from the hardwood were milled and sieved to obtain three size ranges—small, medium, and large. Three mass loadings of each of the biomass samples were fluidized with glass beads as the inert bed material to study the mixing characteristics under varying fluidizing gas velocities. Hydrodynamic data including pressure drop measurements at select vertical locations and video recordings of the mixture behavior were generated in this effort.

09 BIOMASS FUELS↗
DOE OSTI2023

Unsupervised Learning Based Interaction Force Model for Nonspherical Particles in Incompressible Flows

This project provides a neural network-based interaction force model for gas-solid flows from low to intermediate Reynolds numbers and concentration, which can be linked to MFiX-DEM. We have constructed a database of the interaction force between the irregular-shaped particles using a spherical harmonic method and the fluid phase based on the particle-resolved direct numerical simulation (PR-DNS) with immersed boundary-based gas kinetic scheme. Unsupervised learning method, i.e., variational auto-encoder (VAE) has been applied to extract the primitive shape factors determining the drag force, lifting forces, and torque. The interaction force model has been trained and validated with a simple but effective multi-layer feed-forward neural network: multi-layer perceptron (MLP), which will be concatenated after the encoder of the previously trained VAE for geometry feature extraction for single, irregular particles. We have trained transpose convolutional neural networks with the PR-DNS data to predict the velocity and pressure gradient of the single particle systems and utilized them to calculate drag force of multi-particle systems. This model can provide high computational efficiency because it does not require collecting multiparticle system data from PR-DNS.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

43-0434 Bridge Inspection Report 2021

The documentation includes the following: 1. 2021 NMDOT Bridge Inspection Report including Element Level Data Collection (in digital formats prepared by NMSU)--conforms to the National Bridge Inspection Standards and AASHTO Manual for Bridge Element Inspection; 2. 2021 Supplemental Report (in digital formats prepared by NMSU)--provides detailed information related to current condition of major bridge components; 3. 2021 Inspection Pictures (in digital format prepared by NMSU); 4. 2021 Delamination Map (in digital format prepared by LANL). During the 2021 inspection, three critical findings were reported to our LANL contact, Mr. Jonathan Stein, by text and/or email on June 26, 2021. 1. The south joint of the bridge had a modular section of the joint that could potentially come loose during the passing of vehicles. This was reported as a critical finding for safety. If the modular section becomes dislodged or deformed it could pose a serious hazard that could result in a punctured tire to vehicles, motorcyclists and/or bicyclist causing drivers/riders to lose control; 2. The bracket plate located at the north end of the bridge on the pedestrian walkway was corroded through and provides little protection to pedestrians and bicyclists. This was reported as a critical finding for safety to prevent injuries to people; 3. The north approach rail has three missing posts. These posts help ensure that traffic is redirected and the energy is absorbed by the rail. This was reported as a critical finding for safety. Additionally, the approach rail is on a curve with a nearby drop off. Immediate repair is recommended. Based on the 2021 inspection, the bridge deck is rated in "fair" condition. The chain drag performed on the deck identified several areas with delamination that are concentrated near the expansion joints, in the closure joint of the deck near the bridge centerline, and at the south end of the northbound lanes. The chain drag performed during the 2021 inspection revealed 215,333 sq. in. (1495 sq. ft) of delaminations and patched areas (not including the sidewalk). This is approximately a 600% increase from 2019. It is recommended that the delaminations and spalls with exposed rebar be repaired. The "delaminated area" map for the 2021 inspection is provided in the supplemental report.

42 ENGINEERING↗
DOE OSTI2024

History and State of the Art in Advanced Thermal Protection Systems

This report documents the current state of the art of thermal protection systems (TPS) relevant to missile application materials and MTCR control language. The information is intended to provide data on current TPS for determining the current state of development and identifying potential research directions. Various vehicle configurations and types of TPS are described, including the evolution of materials and systems. Vehicle configurations such as high drag rocket and capsule systems and low-drag air-breathers are discussed. TPS are complex, transitioning from often very low temperatures to temperatures exceeding 3000°C and heat fluxes of up to 30,000 W/cm 2 . They can take the form of blankets, sheets, tiles and fabricated rigid structures, and can be made from high temperature capable alloys, refractory metals, ceramics, ceramic matrix composites, ultra-high temperature ceramics, carbon-carbon composites or ablative carbonaceous materials. Current state of the art in TPS is based on materials any systems developed for the initial space program. TPS was upgraded to the point of the Space Shuttle development and has stagnated since them.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design of a Preliminary Family of Airfoils for High Reynolds Number Wind Turbine Applications

For the past 30 years, offshore wind turbines exhibited a continual pattern of growth that is expected to continue as the industry pushes for higher efficiency. Current designs for the next generation of wind turbines are so large that the chordwise Reynolds number of the blades is well beyond the design range of existing open-source airfoil families. This paper presents a preliminary family of new airfoils designed specifically for the needs of these next-generation offshore turbines, ranging from 21% thick to 30% thick with operating Reynolds numbers between 12 million and 18 million. These airfoils are intended to be alternative to the FFA airfoils that are commonly used on reference turbines such as the IEA 15MW and 22MW designs. In this work, airfoil performance metrics and design targets are developed, the design process is outlined, an optimization scheme is presented, and finally the airfoils and their simulated performance are compared to existing baselines. Lift to drag ratios in a clean condition were improved by up to 49.3% from the baseline FFA airfoil, and rough condition lift to drag ratio was improved by up to 9.3%. It is estimated that the cumulative improvements provided by this airfoil family would result in an approximately 1% increase of Annual Expected Power (AEP) for the 22 MW turbine compared to the current baseline.

airfoils↗
DOE OSTI2021

NUMERICAL CALCULATION OF LOSSES OF TRAPPED VORTICES UNDER STRONG RF MEISSNER CURRENT AND DC SUPERHEATING FIELD IN TYPE II SUPERCONDUCTORS

Research on the vortex dynamics and enhancing of superheating field in superconductors has attracted much attention in accelerator physics community to develop next-generation high-performance accelerator cavities. However, the extreme dynamics of curvilinear elastic vortices driven by very strong currents close to the depairing limit or superheating field of a superconductor with a nanostructured surface has not been well understood. We calculated the superheating field Hsh and critical momentum kc characterizing the wavelength of the instability ?m of the Meissner state to flux penetration by solving numerically the Ginzburg?Landau equations. A bulk superconductor, superconductor with the inhomogeneous surface disorder (S-S), and multilayered surface (S-I-S) have been thoroughly investigated in this work. Our result showed that S-S and S-I-S structures can enhance the superheating field well above their clean limit. In this work extensive numerical simulation of the power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was performed. Our simulations take into account the nonlinear vortex line tension, vortex mass, Bardeen?Stephen viscous vortex drag applicable at low fields, and nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient ?(v) at high fields and pinning force. We showed that the LO decrease of ?(v) with the vortex velocity v could radically change the field dependence of the surface resistance Ri(H) caused by trapped vortices. At low frequencies Ri(H) exhibits a conventional increase with H. However, as frequency increases, the surface resistance becomes a nonmonotonic function of H which decreases with H at higher fields irrespective of the pinning distribution. Overheating can mask the descending field dependence of Ri(H) as frequency increases. Our numerical simulations also show that the LO effect can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along with the vortex when a vortex is pinned strongly to one end. Nonlinear losses of trapped vortices in thick films under high-amplitude RF fields as functions of frequency, mean free path and pinning characteristics have been calculated.

Pathirana, Walive↗
DOE OSTI2021

Fluctuations and Pinning for Individually Manipulated Skyrmions

We numerically examine the dynamics of individually dragged skyrmions interacting simultaneously with an array of other skyrmions and quenched disorder. For drives just above depinning, we observe a broadband noise signal with a 1/ f characteristic, while at higher drives, narrowband or white noise appears. Even in the absence of quenched disorder, the threshold force that must be applied to translate the driven skyrmion is finite due to elastic interactions with other skyrmions. The depinning threshold increases as the strength of the quenched disorder is raised. Above the depinning force, the skyrmion moves faster in the presence of quenched disorder than in a disorder-free system since the pinning sites prevent other skyrmions from being dragged along with the driven skyrmion. For strong pinning, we find a stick-slip motion of the driven skyrmion which produces a telegraph noise signature. The depinning threshold increases monotonically with skyrmion density in the absence of quenched disorder, but when pinning is present, the depinning threshold changes nonmonotonically with skyrmion density, and there are reentrant pinned phases due to a competition between pinning induced by the quenched disorder and that produced by the elastic interactions of the skyrmion lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Giant Angular Nernst Effect in the Organic Metal α-(BEDT-TTF) 2 KHg(SCN) 4

We have detected a large Nernst effect in the charge density wave state of the multiband organic metal α-(BEDT-TTF) 2 KHg(SCN) 4 . We find that apart from the phonon drag effect, the energy relaxation processes that govern the electron–phonon interactions and the momentum relaxation processes that determine the mobility of the q1D charge carriers have a significant role in observing the large Nernst signal in the CDW state in this organic metal. The emphasised momentum relaxation dynamics in the low field CDW state (CDW 0 ) is a clear indicator of the presence of a significant carrier mobility that might be the main source for observation of the largest Nernst signal. The momentum relaxation is absent with increasing angle and magnetic field, i.e., in the high-field CDW state (CDW x ) as evident from the much smaller Nernst effect amplitude in this state. In this case, only the phonon drag effect and electron–phonon interactions are contributing to the transverse thermoelectric signal. Our findings advance and change previous observations on the complex properties of this organic metal.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Applying Physics-Informed Neural Networks to Solve Navier–Stokes Equations for Laminar Flow around a Particle

In recent years, Physics-Informed Neural Networks (PINNs) have drawn great interest among researchers as a tool to solve computational physics problems. Unlike conventional neural networks, which are black-box models that “blindly” establish a correlation between input and output variables using a large quantity of labeled data, PINNs directly embed physical laws (primarily partial differential equations) within the loss function of neural networks. By minimizing the loss function, this approach allows the output variables to automatically satisfy physical equations without the need for labeled data. The Navier–Stokes equation is one of the most classic governing equations in thermal fluid engineering. This study constructs a PINN to solve the Navier–Stokes equations for a 2D incompressible laminar flow problem. Flows passing around a 2D circular particle are chosen as the benchmark case, and an elliptical particle is also examined to enrich the research. The velocity and pressure fields are predicted by the PINNs, and the results are compared with those derived from Computational Fluid Dynamics (CFD). Additionally, the particle drag force coefficient is calculated to quantify the discrepancy in the results of the PINNs as compared to CFD outcomes. The drag coefficient maintained an error within 10% across all test scenarios.

Hu, Beichao (ORCID:0009000163215151)↗
DOE OSTI2021

Spitzer Phase-curve Observations and Circulation Models of the Inflated Ultrahot Jupiter WASP-76b

The large radii of many hot Jupiters can only be matched by models that have hot interior adiabats, and recent theoretical work has shown that the interior evolution of hot Jupiters has a significant impact on their atmospheric structure. Due to its inflated radius, low gravity, and ultrahot equilibrium temperature, WASP-76b is an ideal case study for the impact of internal evolution on observable properties. Hot interiors should most strongly affect the nonirradiated side of the planet, and thus full phase-curve observations are critical to ascertain the effect of the interior on the atmospheres of hot Jupiters. In this work, we present the first Spitzer phase-curve observations of WASP-76b. We find that WASP-76b has an ultrahot dayside and relatively cold nightside with brightness temperatures of 2471 ± 27 K/1518 ± 61 K at 3.6 μm and 2699 ± 32 K/1259 ± 44 K at 4.5 μm, respectively. These results provide evidence for a dayside thermal inversion. Both channels exhibit small phase offsets of 0.68 ± 0.°48 at 3.6 μm and 0.67 ± 0.°2 at 4.5 μm. We compare our observations to a suite of general circulation models (GCMs) that consider two endmembers of interior temperature along with a broad range of frictional drag strengths. Strong frictional drag is necessary to match the small phase offsets and cold nightside temperatures observed. From our suite of cloud-free GCMs, we find that only cases with a cold interior can reproduce the cold nightsides and large phase-curve amplitude at 4.5 μm, hinting that the hot interior adiabat of WASP-76b does not significantly impact its atmospheric dynamics or that clouds blanket its nightside.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hydrodynamics and Survivability during Post-main-sequence Planetary Engulfment

The engulfment of substellar bodies (SBs), such as brown dwarfs and planets, by giant stars is a possible explanation for rapidly rotating giants, lithium-rich giants, and the presence of SBs in close orbits around subdwarfs and white dwarfs. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an engulfed SB. We model the SB as a rigid body with a reflective surface because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram-pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory inside the star. We find that engulfment can increase the luminosity of a 1 M ⊙ star by up to a few orders of magnitude. The time for the star to return to its original luminosity is up to a few thousand years when the star has evolved to ≈10 R ⊙ and up to a few decades at the tip of the red giant branch (RGB). No SBs can eject the envelope of a 1 M ⊙ star before it evolves to ≈10 R ⊙ if the orbit of the SB is the only energy source contributing to the ejection. In contrast, SBs as small as ≈10 MJup can eject the envelope at the tip of the RGB. The numerical framework we introduce here can be used to study planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB.

79 ASTRONOMY AND ASTROPHYSICS↗