Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D velocity models”

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 271 records · Page 15

Planet Formation by Gas-assisted Accretion of Small Solids

We compute the accretion efficiency of small solids, with radii 1 cm ≤ R s ≤ 10 m, on planets embedded in gaseous disks. Planets have masses 3 ≤ M p ≤ 20 Earth masses (M ⊕ ) and orbit within 10 au of a solar mass star. Disk thermodynamics is modeled via 3D radiation-hydrodynamics calculations that typically resolve the planetary envelopes. Both icy and rocky solids are considered, explicitly modeling their thermodynamic evolution. The maximum efficiencies of 1 ≤ R s ≤ 100 cm particles are generally ≲10%, whereas 10 m solids tend to accrete efficiently or be segregated beyond the planet’s orbit. A simplified approach is applied to compute the accretion efficiency of small cores, with masses M p ≤ 1 M ⊕ and without envelopes, for which efficiencies are approximately proportional to $M^{2/3}_{p}$. The mass flux of solids, estimated from unperturbed drag-induced drift velocities, provides typical accretion rates dM p /dt ≲ 10 -5 M ⊕ yr -1 . In representative disk models with an initial gas-to-dust mass ratio of 70–100 and total mass of 0.05–0.06 M ⊙ , the solids’ accretion falls below 10 -6 M ⊕ yr -1 after 1–1.5 Myr. The derived accretion rates, as functions of time and planet mass, are applied to formation calculations that compute dust opacity self-consistently with the delivery of solids to the envelope. Assuming dust-to-solid coagulation times of ≈0.3 Myr and disk lifetimes of ≈3.5 Myr, heavy-element inventories in the range 3–7 M ⊕ require that ≈90–150 M ⊕ of solids cross the planet’s orbit. The formation calculations encompass a variety of outcomes, from planets a few times M ⊕ , predominantly composed of heavy elements, to giant planets. The peak luminosities during the epoch of the solids’ accretion range from ≈10 -7 to ≈10 -6 L ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗

Kinematic Evidence for an Embedded Planet in the IM Lupi Disk

Abstract We test the hypothesis that an embedded giant planet in the IM Lupi protostellar disk can produce velocity kinks seen in CO line observations as well as the spiral arms seen in scattered light and continuum emission. We inject planets into 3D hydrodynamics simulations of IM Lupi, generating synthetic observations using Monte Carlo radiative transfer. We find that an embedded planet of 2–3 M Jup can reproduce non-Keplerian velocity perturbations, or “kinks”, in the 12 CO J = 2–1 channel maps. Such a planet can also explain the spiral arms seen in 1.25 mm dust continuum emission and 1.6 μ m scattered-light images. We show that the wake of the planet can be traced in the observed peak velocity map, which appears to closely follow the morphology expected from our simulations and from analytic models of planet–disk interaction.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmark microgravity experiments and computations for 3D dendritic-array stability in directional solidification

In this study, we present a comprehensive quantitative analysis of stability bands for dendritic arrays during directional solidification of a transparent succinonitrile-0.46 wt % camphor alloy, spanning a broad range of pulling velocities. Taking advantage of the microgravity environment aboard the International Space Station where most convection effects are suppressed, we obtain unique measurements that quantify the stable primary spacing range of spatially extended three-dimensional dendritic array structures under purely diffusive growth conditions. Through carefully designed velocity jump experiments and detailed examination of sub-grain boundary dynamics, we characterize key instabilities, including elimination and tertiary branching, shedding new light on the mechanisms governing dynamic dendritic spacing selection in extended 3D arrays. Phase field simulations are performed to characterize the stability limits of dendritic array structures for quantitative comparison with the flight experiments. Although the simulations capture general trends, significant deviations are noted at the upper stability boundary, indicating the influence of additional, unexplored factors. These findings contribute to a deeper understanding of dendritic growth dynamics and offer valuable benchmark data that could aid in refining predictive models and improving control of dendritic microstructures in metallurgical applications.

36 MATERIALS SCIENCE↗

Multigrid acceleration and turbulence models for computations of 3D turbulent jets in crossflow

A multigrid method is presented for the calculation of three-dimensional turbulent jets in crossflow. Turbulence closure is achieved with either the standard k-epsilon model or a Reynolds Stress Model (RSM). Multigrid acceleration enables convergence rates which are far superior to that for a single grid method. With the k-epsilon model the rate approaches that for laminar flow, but with RSM it is somewhat slower. The increased stiffness of the system of equations in the latter may be responsible. Computed results with both turbulence models are compared with experimental data for a pair of opposed jets in crossflow. Both models yield reasonable agreement with mean flow velocity but RSM yields better prediction of the Reynolds stresses.

Demuren, A. O.↗

Analysis of core asymmetries in inertial confinement fusion implosions using three-dimensional hot-spot reconstruction

Three-dimensional effects play a crucial role during the hot-spot formation in inertial confinement fusion (ICF) implosions. A data analysis technique for 3D hot-spot reconstruction from experimental observables has been developed to characterize the effects of low modes on 3D hot-spot formations. In nuclear measurements, the effective flow direction, governed by the maximum eigenvalue in the velocity variance of apparent ion temperatures, has been found to agree with the measured hot-spot flows for implosions dominated by mode ℓ = 1. Asymmetries in areal-density (ρR) measurements were found to be characterized by a unique cosine variation along the hot-spot flow axis. In x-ray images, a 3D hot-spot x-ray emission tomography method was developed to reconstruct the 3D hot-spot plasma emissivity using a generalized spherical-harmonic Gaussian function. The gradient-descent algorithm was used to optimize the mapping between the projections from the 3D hot-spot emission model and the measured x-ray images along multiple views. Furthermore, this work establishes a platform to analyze 3D low-mode core asymmetries in ICF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hot electron preheat effects in two-photon polymerization plastic lattice laser-driven shock-tube experiments on the OMEGA laser

We present simulations of Omega shock tube experiments designed to investigate hot electron preheat effects in 3D-printed, two-photon polymerization (2PP) plastic lattices. Preheat is inferred in the experiments from the expansion of a plastic witness disk embedded in the lattice. Using the Eulerian radiation-hydrodynamics code xRAGE, we model shock propagation and preheat from both radiative and hot electron energy sources to evaluate their relative impact. To simulate the transport of laser-generated hot electrons, the nonlocal electron heat transport model proposed by Schurtz, Nicolaï, and Busquet (SNB) is extended with a hot electron source term and an energy cascade algorithm. We explore how variations in ablator, lattice geometry, and laser drive affect the shock velocity and witness disk expansion. Simulations show that the inclusion of a 5 μm gold layer reduces shock pressure by 60% and shock speeds by 30%–40% but does not significantly reduce the hot electron preheat, and that different lattice geometries lead to enhanced shock velocities—up to 40% faster than in homogeneous foams. However, radiative and conductive preheat from classical mechanisms alone fail to match experiment. By including a hot electron source term, we reproduce experimental observables such as disk expansion rates and spatial radiographic features. We find that a hot electron population corresponding to 4%–8% of the incident laser energy with T hot = 50 keV produces expansion which agrees with the experimental data, suggesting hot electron preheat is the most plausible explanation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-dimensional diagnostics and measurements of inertial confinement fusion plasmas

Recent inertial confinement fusion measurements have highlighted the importance of 3D asymmetry effects on implosion performance. One prominent example is the bulk drift velocity of the deuterium–tritium plasma undergoing fusion (“hotspot”), v HS . Upgrades to the National Ignition Facility neutron time-of-flight diagnostics now provide v HS to better than 1 part in 10 4 and enable cross correlations with other measurements. Here we present the impact of v HS on the neutron yield, downscatter ratio, apparent ion temperature, electron temperature, and 2D x-ray emission. The necessary improvements to diagnostic suites to take these measurements are also detailed. The benefits of using cross-diagnostic analysis to test hotspot models and theory are discussed, and cross-shot trends are shown.

47 OTHER INSTRUMENTATION↗

Comparing Methods to Extract Sectional Angle of Attack From Blade-Resolved Rotor Simulations

The angle of attack is widely used to describe the aerodynamics of rotor blades. Models based on lifting lines, actuator lines, and actuator disks use it as a primary variable, assuming the aerodynamics can be represented by the flow around a collection of equivalent 2D airfoils. This paper evaluates the accuracy of existing methods that extract the angle of attack from blade-resolved CFD simulations, towards enabling a direct comparison of angle of attack between blade-resolved and blade-modeled simulations. Results are compared for a 2D airfoil for which the angle of attack is known by definition. The hovering rotor case is also considered, where we look at two different rotor blades. We discuss the factors that impact the accuracy of the 2D predictions. For the 3D rotors, we show that three-dimensional effects in the root and tip regions lead to large discrepancies between methods. This suggests that the very definition of angle of attack is unable to accurately represent the flow in these regions. Finally, based on guidelines that we establish, we propose an alternate method that uses local velocities sampled on a contour located 0.25 chord, or less, away from the airfoil surface and generally yields improved accuracy compared to most existing techniques.

ARMD↗

Kinetic modeling of solitary wave dynamics in a neutralizing ion beam

In this work, we characterize the formation and evolution of electrostatic solitary waves (ESWs) in the space-charge neutralization of ion beams using particle-in-cell simulations. These waves become excited when the electrons emitted from an external filament source initiate a two-stream instability in the beam. We show that such electrostatic waves become excited in both two-dimensional (2D) and three-dimensional (3D) beams with different shapes and sizes. Through a 1D Bernstein–Greene–Kruskal (BGK) analysis of the 2D beam, we find that the non-Maxwellian nature of the beam electrons gives rise to large-sized ESWs that are not predicted by BGK theory since it assumes a Maxwellian electron velocity distribution in the beam. Finally, we show that a 1D BGK theory is inadequate to describe ESWs in 3D beams because of complex electron trajectories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental investigation of the relationship between impact crater morphology and impacting particle velocity and direction

Interpretation of the wealth of impact data available from the Long Duration Exposure Facility, in terms of the absolute and relative populations of space debris and natural micrometeoroids, requires three dimensional models of the distribution of impact directions, velocities and masses of such particles, as well as understanding of the impact processes. Although the stabilized orbit of LDEF provides limited directional information, it is possible to determine more accurate impact directions from detailed crater morphology. The applicability of this technique has already been demonstrated but the relationship between crater shape and impactor direction and velocity has not been derived in detail. We present the results of impact experiments and simulations: (1) impacts at micron dimensions using the Unit's 2MV Van de Graaff accelerator; (2) impacts at mm dimensions using a Light Gas Gun; and (3) computer simulations using AUTODYN-3D from which an empirical relationship between crater shape and impactor velocity, direction and particle properties we aim to derive. Such a relationship can be applied to any surface exposed to space debris or micrometeoroid particles for which a detailed pointing history is available.

Mackay, N. G.↗

Multigrid acceleration and turbulence models for computations of 3D turbulent jets in crossflow

A multigrid method is presented for the calculation of three-dimensional turbulent jets in crossflow. Turbulence closure is achieved with either the standard k-epsilon model or a Reynolds stress model (RSM). Multigrid acceleration enables convergence rates which are far superior to that for a single grid method to be obtained with both turbulence models. With the k-epsilon model the rate approaches that for laminar flow, but with RSM it is somewhat slower. The increased stiffness of the system of equation in the latter may be responsible. Computed results with both turbulence models are compared to experimental data for a pair of opposed jets in crossflow. Both models yield reasonable agreement for the mean flow velocity, but RSM yields better predictions of the Reynolds stresses.

Demuren, A. O.↗

Towards a quantitative understanding of bonding in supersonic single particle impacts: A three-dimensional FIB-SEM exploration

Particle bonding is crucial to coating quality in cold spray, but it has been a challenge to accurately quantify bonding even in single particle impacts. This paper uses FIB-SEM to explicitly map the particle-substrate interface for Cu-on-Cu single microparticle impacts in a full 3D rendering that spans a wide range of impact velocities. This approach permits a detailed quantification of the total bonding area and all of its associated components. In addition to revealing why prior 2D characterization efforts have missed important details about impact bonding, these data quantitatively reveal the evolution of bonding from its onset at the “critical velocity” V cr (where bonding is generally poor, ∼6 %) to its peak at around 1.3‧V cr (where almost 90 % of the particle bonds). Further increase in the velocity to 1.5‧V cr and beyond finds the onset of hydrodynamic penetration and a decrease in bonding. These data then support the development of a simple analytical model based on oxide rarefication and extrusion of bare metal through gaps in the oxide layer as driving the development of bonding. As a result, the model reproduces the experiments and provides guidance on optimization of bonding as a function of material and process parameters.

Analytical methods↗

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗

Intrinsic Elastic Anisotropy of Westerly Granite Observed by Ultrasound Measurements, Microstructural Investigations, and Neutron Diffraction

Westerly granite (WG) has been generally accepted as an isotropic homogeneous rock. Here, we return to WG and observe significant elastic anisotropy using multidirectional ultrasonic sounding on spherical samples at pressures up to 400 MPa. Thermal treatment of WG leads to formation of microcracks that reduce elastic wave velocities and increase its elastic anisotropy. The 3D distribution of P-wave velocities at low pressure is close to orthorhombic symmetry. Application of hydrostatic pressure closes most of thermally induced microcracks and decreases elastic anisotropy of WG, but at high pressure the anisotropy is practically reversed compared to low pressure: maximum P-wave velocity direction at low pressures is near minimum velocity direction at high pressure and vice versa. To understand this effect, microstructures of the rock were investigated by optical and scanning electron microscopy. Preferred orientations of four major rock-forming minerals—quartz, orthoclase, plagioclase, and biotite—were measured by time-of-flight neutron diffraction, which confirms significant crystal alignment. All these data were used to numerically model anisotropic elastic properties of WG. It is shown that WG possesses weak intrinsic elastic anisotropy related mainly to the preferred orientation of feldspars formed during igneous crystallization. Observed microcracks are mostly related to the cleavage planes of feldspars and biotite, and thus also demonstrate preferred orientation. Higher preheating temperatures produce larger quantity of longer microcracks. A numerical model shows that these microcracks act against the weak intrinsic elastic anisotropy of WG, and define the elastic anisotropy at low pressures.

58 GEOSCIENCES↗

Controls and process planning strategies for 5-axis laser directed energy deposition of Ti-6Al-4V using an 8-axis industrial robot and rotary motion

Interest in multi-axis additive manufacturing has grown significantly in recent years due to associated benefits that range from improved quality to support for greater geometric complexity. A method of 5-axis Laser Hot-Wire Directed Energy Deposition is developed and demonstrated in this study using an 8-axis industrial robot, comprised of a 6-axis serial link robot and a 2 degree-of-freedom positioner. Coordinated 5-axis motion is used to print Ti-6Al-4V objects predominantly through build plate rotation. Methods of slicing and toolpath generation are presented, including the use of an alternating toolpath direction philosophy and the manipulation of the CAD model surface normal vectors for determining the required build plate tilt angles for maintaining consistent deposition orientations relative to gravity. Another objective of the study was to determine the limits of conventional slicing and layer height control techniques for printing geometries with significant unsupported overhangs. A comparison is made between the behavior of normal and conformal bead alignment in this regard, and the performance of the layer height controller is reported as it relates to process stability and the limits of printing overhang geometries. Unsupported overhangs were successfully printed at 35° without build plate tilt and 45° with build plate tilt. Toolpath velocity and net-shape printing performance are reported as well, with 3D scans of the printed surfaces used for net-shape characterization. In all, the printing of three geometries is presented, including the overhang geometry, a 5-axis toolpath geometry, and a large-scale, multi-bead demonstrator part, which exceeded 500 mm in height and weighed 44 kg (97 lbs).

36 MATERIALS SCIENCE↗

The simulation of 3D hypervelocity spallation using a hydrocode PAGOSA with FLIP+MPM

n this work, a hydrocode PAGOSA with FLIP+MPM is presented and exercised to investigate the fracture in ductile material. The merit of PAGOSA with FLIP+MPM to solve the advection problem is first illustrated by a solid piston periodically moving in a sealed tube with air. Furthermore, the ability of PAGOSA with FLIP+MPM to capture the fracture in material is shown by a simple stretching fracture in ductile material. In both of two benchmark problems, the PAGOSA results and analytical solutions are also presented for comparison. Then PAGOSA with FLIP+MPM is used to model complex spall in ductile material, which is a challenging problem in engineering applications. The convergences of PAGOSA with FLIP+MPM—based on both the mesh size and the marker density—are investigated by monitoring free surface velocity. To further show the ability of PAGOSA with FLIP+MPM to predict fracture in ductile material, the numerical results are compared with the experimental results and other published numerical results. Moreover, the effect of the spall parameter in PAGOSA with FLIP+MPM on the numerical simulation is also analyzed by investigating the free surface velocity. Finally, the effects of the peak compressive stress, the tensile strain rate, and the loading rate on the spallation are further investigated. The numerical results show that PAGOSA with FLIP+MPM can improve PAGOSA's performance when applied to predicting fracture in ductile materials, and is robust enough to accurately predict spall fracture in ductile materials.

42 ENGINEERING↗

Solar Wind Strahl Observations and Their Implication to the Core-Halo Formation due to Scattering

A study of the kinetic properties of the strahl electron velocity distribution functions (VDF?s) in the solar wind is presented. This study focuses on the mechanisms that control and regulate the electron VDF?s and the stability of the strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. These parameters are used to investigate the stability of the strahl population. The analysis check for whether the strahl electrons are constrained by some instability (e.g., the whistler or KAW instabilities), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling of the 3D-VDF?s from which the moments properties of the various populations are obtained. The results of this study have profound implication on the current hypothesis about the probable formation of the solar wind halo electrons produced from the scattering of the strahl population. This hypothesis is strengthened by direct observations of the strahl electrons being scattered into the core-halo in an isolated event. The observation implies that the scattering of the strahl is not a continuous process but occurs in bursts in regions where conditions for wave growth providing the scattering are optimum. Sometimes, observations indicate that the strahl component is anisotropic (Tper/Tpal approx. 2). This provides a possible free energy source for the excitation of whistler waves as a possible scattering mechanism, however this condition is not always observed. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer.

Vinas, Adolfo F.↗