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 127 records · Page 7

Numerical prediction of wakes in cascades and compressor rotors including the effects of mixing. II - Rotor passage flow and wakes including the effects of spanwise mixing

The results of a numerical investigation to predict the flow field including wakes and mixing in axial-flow compressor rotors are presented. The wake behavior in a moderately loaded compressor rotor is studied numerically using a 3D incompressible Navier-Stokes solver with a high Reynolds number form of a turbulence model. The equations are solved using a time dependent implicit technique. The agreement between the measured data and the predictions is good; including the blade boundary-layer profiles, wake mean-velocity profiles, and decay. The ability of the pseudocompressibility scheme to predict the entire flow field including the near and far wake profiles and its decay characteristics, effect of loading, and the viscous losses of a 3D rotor flow field are demonstrated. The mixing in the downstream regions away from the hub and annulus walls is dominated by wake diffusion. In regions away from the walls the radial mixing is predominantly caused by the transport of mass, momentum, and energy by the radial component of velocity in the wake.

Lakshminarayana, B.↗

Characteristics of 3D turbulent jets in crossflow

Three dimensional turbulent jets in crossflow at low to medium jet-to-crossflow velocity ratios are computed with a finite volume numerical procedure which utilizes a second-moment closure model to approximate the Reynolds stresses. A multigrid method is used to accelerate the convergence rate of the procedure. Comparison of the computations to measured data show good qualitative agreement. All trends are correctly predicted, though there is some uncertainty on the height of penetration of the jet. The evolution of the vorticity field is used to explore the jet-crossflow interaction.

Demuren, A. O.↗

Characteristics of 3D turbulent jets in crossflow

Three dimensional turbulent jets in crossflow at low to medium jet-to-crossflow velocity ratios are computed with a finite volume numerical procedure which utilizes a second-moment closure model to approximate the Reynolds stresses. A multigrid method is used to accelerate the convergence rate of the procedure. Comparison of the computations to measured data show good qualitative agreement. All trends are correctly predicted, though there is some uncertainty on the height of penetration of the jet. The evolution of the vorticity field is used to explore the jet-crossflow interaction.

Demuren, A. O.↗

The GMRT archive atomic gas survey – II. Mass modelling and dark matter halo properties across late-type spirals

Studying the kinematics and mass modelling of galaxies from H i 21 cm data provides valuable insights into the properties of both the baryonic components and the dark matter halo in nearby galaxies. Despite many observational studies, mass modelling of galaxies remains challenging due to different limitations. For example, most of the previous studies involving mass modelling are based on rotation curves derived from 2D velocity fields from H i or H α spectroscopic observation which are often affected by beam smearing and projection effect. However, kinematic modelling done by fitting the ‘Tilted ring model’ to 3D data cube is not affected by these issues. In this study, we present and compare 3D kinematic modelling of a pilot sample of 11 galaxies from the GMRT archive atomic gas survey (GARCIA) using two different publicly available pipelines. We model the observed H i rotation curve using 3.6-μm infrared data and SDSS r-band data for stellar contribution, H i surface density profile for gas, and Navarro–Frenk–White profile for dark matter halo; and employ the Markov chain Monte Carlo optimization method for parameter estimation. Further, to validate our analysis, we revisit important scaling relations, e.g. the M gas –M star relation, M star –M halo relation, M gas –M halo relation and Baryonic Tully–Fisher relation. Here, the scaling relations from our analysis are broadly consistent with that reported in the literature. A larger sample of galaxies from GARCIA in the near future will allow studying these scaling relations in greater details.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating Isotropic Rayleigh Wave Horizontal to Vertical Ratios at Rock Valley

Elastic anisotropy in the subsurface can manifest from pre-existing structure or stress-induced cracking – both phenomena are of interest at the Rock Valley site in Nevada. The azimuthal dependence of Rayleigh wave ellipticity, also known as H/V, can reveal elastic anisotropy, but the effect of isotropic variability in H/V must first be accounted for. To that end, we simulate synthetic seismograms from a ring of sources in a fully 3-D isotropic velocity model that captures the geology at Rock Valley. We measure H/V from these waveforms and compare them with H/V predictions obtained using 1D models extracted at specific locations from the 3D model. We find that an average of the azimuthal H/V measurements from the 3-D simulations and 1-D predictions agree, but that the azimuthal variability of H/V is a large fraction of the mean value. These results estimate the noise character in azimuthal H/V so it can be factored out of future anisotropic simulations.

58 GEOSCIENCES↗

Effects of stratospheric lapse rate on numerically simulated thunderstorm cloud top structure

The study of overshooting midlatitude severe thunderstorm tops has led to the detection of some intriguing patterns. This paper provides a summary of an investigation of the effects of stratospheric stability on thunderstorm cloud top structure, taking into account the utilization of a fully three-dimensional (3D) anelastic model discussed by Schlesinger (1975, 1978, 1984). The layout of model experiments is discussed, giving attention to initial environmental soundings, and initial forcing. The time variation of the vertical velocity extreme is examined, and a description of the mature storm structure is presented. Constant-altitude plane views throughout storm depth are considered along with cloud top height fields. With the stratospheric lapse, the warm region is small and breaks down into separate spots with secondary cold spots downshear of them.

Schlesinger, R. E.↗

Upgrades and Modifications of the NASA Ames HFFAF Ballistic Range

The NASA Ames Hypervelocity Free Flight Aerodynamics Facility ballistic range is described. The various configurations of the shadowgraph stations are presented. This includes the original stations with film and configurations with two different types of digital cameras. Resolution tests for the 3 shadowgraph station configurations are described. The advantages of the digital cameras are discussed, including the immediate availability of the shadowgraphs. The final shadowgraph station configuration is a mix of 26 Nikon cameras and 6 PI-MAX2 cameras. Two types of trigger light sheet stations are described visible and IR. The two gunpowders used for the NASA Ames 6.251.50 light gas guns are presented. These are the Hercules HC-33-FS powder (no longer available) and the St. Marks Powder WC 886 powder. The results from eight proof shots for the two powders are presented. Both muzzle velocities and piston velocities are 5 9 lower for the new St. Marks WC 886 powder than for the old Hercules HC-33-FS powder (no longer available). The experimental and CFD (computational) piston and muzzle velocities are in good agreement. Shadowgraph-reading software that employs template-matching pattern recognition to locate the ballistic-range model is described. Templates are generated from a 3D solid model of the ballistic-range model. The accuracy of the approach is assessed using a set of computer-generated test images.

Ballistic Range↗

Upgrades and Modifications of the NASA Ames HFFAF Ballistic Range

The NASA Ames Hypervelocity Free Flight Aerodynamics Facility ballistic range is described. The various configurations of the shadowgraph stations are presented. This includes the original stations with film and configurations with two different types of digital cameras. Resolution tests for the 3 shadowgraph station configurations are described. The advantages of the digital cameras are discussed, including the immediate availability of the shadowgraphs. The final shadowgraph station configuration is a mix of 26 Nikon cameras and 6 PI-MAX2 cameras. Two types of trigger light sheet stations are described visible and IR. The two gunpowders used for the NASA Ames 6.251.50 light gas guns are presented. These are the Hercules HC-33-FS powder (no longer available) and the St. Marks Powder WC 886 powder. The results from eight proof shots for the two powders are presented. Both muzzle velocities and piston velocities are 5 9 lower for the new St. Marks WC 886 powder than for the old Hercules HC-33-FS powder (no longer available). The experimental and CFD (computational) piston and muzzle velocities are in good agreement. Shadowgraph-reading software that employs template-matching pattern recognition to locate the ballistic-range model is described. Templates are generated from a 3D solid model of the ballistic-range model. The accuracy of the approach is assessed using a set of computer-generated test images.

Ballistic Range↗

Determination of 3D coronal structures from UVCS/SOHO synoptic observations

Daily coronal synoptic observations with the ultraviolet coronagraph spectrometer (UVCS) on SOHO are used to create Carrington maps of line of sight intensities and spectral line widths for H I Lyman alpha and the O 6 1032/1037 doublet. These 2D maps can then be used as inputs to a tomographic inversion routine to produce 3D models of the corona out to the 2.5 solar radius in the polar regions (3.0 solar radius in equatorial regions). The initial results for the determination of O(5+) outflow velocities as a function of heliographic latitude is presented.

Strachan, L.↗

The Causes of Propeller Pitching Moment and the Conditions for its Significance

Recent development of vertical takeoff and landing (VTOL) aircraft has renewed interest in the study of propellers. One metric in particular, the propeller pitching moment, has been observed to be important to VTOL aircraft stability and control in the past. Propellers at angles of attack could not be accurately modeled in generations past due to a lack of computational power, but even with advances in computer technology, modern designers seem to possess insufficient knowledge in this area. In this dissertation, we study the physics behind propeller pitching moment in the context of an isolated propeller and a propeller upstream of a wing. An unsteady 3D vortex lattice method is developed specifically to model propellers at angles of attack and is validated by comparing to high-fidelity CFD analyses. We then use the model to isolate velocity influences to show that the propeller pitching moment is largely caused by two effects: a skewed wake and the presence of wing circulation. Generated maps of propeller pitching moment over a range of operational parameters corresponding to VTOL transition show that the low flight speeds and high angles of attack encountered during transition lead to significant magnitudes of propeller pitching moment that would be difficult to trim using passive methods. Also, derivation of a generalizable metric of significance shows that the peak contribution of propeller pitching moment to aircraft stability is comparable to a longitudinal displacement of the center of gravity by several percent of the wing chord. Finally, we give a concluding discussion on the impact of propeller pitching moment on VTOL aircraft design.

Xiaofan Fei↗

3D Cell-centered hydrodynamics with subscale closure model and multi-material remap

We extend a higher-order finite volume cell-centered hydrodynamic (CCH) formulation to include an interface-aware subscale closure model and a multi-material remap for simulating 3D compressible hydrodynamic problems within an arbitrary Lagrangian-Eulerian (ALE) framework. This CCH formulation involves a multidirectional approximate Riemann solution using quadratic polynomial reconstructions of the stress tensor and the velocity. At the subscale level, we determine pair-wise material interactions by solving a distinct approximate Riemann problem at the common interface, using the volume of fluids (VOF) method to find the interface. Material interactions are constrained to ensure smooth pressure equilibration among materials. The accuracy and robustness of the ALE method is demonstrated by simulating a suite of 3D Cartesian multi-material problems covering both gas and solid dynamics, where each test case has two or more materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonresonant particle acceleration in strong turbulence: Comparison to kinetic and MHD simulations

Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. Additionally we also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling Stellar Jitter and its Effects on the Stellar Spectrum for More Precise Radial Velocity Measurements

Detection of Earth-size exoplanets is a technological and data analysis challenge. Future progress in Earth-mass exoplanet detection is expected from the development of extreme precision radial velocity (EPRV) measurements. The required precision of these measurements necessitates developing a new physics-based data analysis methodology to discriminate planetary signals from host-star-related effects, taking stellar variability and instrumental uncertainties into account. In this work, we investigate and quantify stellar disturbances of the planet-hosting solar-type star HD121504 from 3D radiative modeling obtained with the StellarBox code. The model was used to determine the statistical properties of the turbulent plasma and to obtain synthetic spectroscopic observations for several Fe I lines at different locations on the stellar disk to mimic high-resolution spectroscopic observations. Presented are the initial results of our effort to characterize stellar jitter properties as well as future objectives.

SMD↗

Development of a coupled experimental–computational approach for engineering optimization of spout-fluidized bed particle coating systems

The design of spout-fluidized bed (SFB) coating systems for nuclear particle fuels typically relies on trial-and-error processes, comprising iterative and time-consuming coating deposition experiments and post-deposition characterization. At an engineering scale, this approach to guided SFB system design is inefficient, highlighting the need for streamlined experimental methodologies which can correlate fluidization conditions to downstream coating outcomes. In this study, we combine time-resolved particle image velocimetry (PIV) with CFD–DEM simulations to benchmark hydrodynamic behavior in a 3D spout-fluidized bed. By exploiting easily accessible optical measurements of particle motion at the bed wall and within the spouting region, we obtain quantitative velocity fields that can be directly compared with model predictions of the occluded bed region, without resorting to complex imaging and characterization techniques such as X-ray or magnetic resonance tomography. Experimental benchmarking reveals strong agreement between CFD–DEM and PIV in the spout and annulus regions, while discrepancies near the wall highlight areas for future model development. Here, the proposed integrated experimental–numerical framework will enable a direct connection between measured variables and numerically predicted fluidization performance of dense, surrogate nuclear particle fuel feedstock such that experimental SFB component design can be rapidly evaluated, informing design decisions for nozzle geometry and operating conditions. Future work will extend this framework by correlating quantified fluidization metrics across nozzle geometries and operating conditions with the resulting coating morphology, microstructure, and uniformity. Establishing these correlations will enable predictive links between hydrodynamic performance and coating quality, providing a rational, scalable basis for optimizing SFB design prior to coating deposition.

CFD/DEM↗

Revisiting the Hail Radar Reflectivity–Kinetic Energy Flux Relation by Combining T-Matrix and Discrete Dipole Approximation Calculations to Size Distribution Observations

The retrieval of hail kinetic energy with weather radars or its simulation in numerical models is challenging because of the shape complexity and variable density of hailstones. We combine 3D scans of individual hailstones with measurements of the particle size distributions (PSD) and T-matrix calculations to understand how hail reflectivity Z changes when approximating hailstones as spheroids, as compared to the realistic shapes obtained by 3D scanning technology. Additionally, recent terminal velocity relations are used to compare Z to the hail kinetic energy flux E ˙ . We parameterize the hail backscattering cross sections at L, S, C, and X bands as a function of size between 0.5 and 5.0 cm, matching the range of the observed PSDs. The scattering calculations use the T-matrix method for size parameters below 1.0 and the discrete dipole approximation (DDA) method otherwise. The DDA calculations are done for 48 digital models of realistic hailstones of sizes between 1 and 5 cm. The DDA cross sections are calculated for multiple orientations and averaged assuming a fully random orientation distribution to provide a single value per hailstone. The T-matrix reflectivity assuming solid ice spheres presents negligible differences to DDA results for size parameters below 1.0. Therefore, T matrix was used to fill in the gaps left by the DDA calculations. The results are mapped to the same size bins of the observed PSDs, allowing the calculation of the radar reflectivity. This is then correlated to E ˙ , allowing a potential improvement of past retrieval methods of E ˙ from Z in multiple wavelengths.

54 ENVIRONMENTAL SCIENCES↗

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗

Saturn's Magnetosphere and Properties of Upstream Flow at Titan: Preliminary Results

Using Cassini Plasma Spectrometer (CAPS) Ion Mass Spectrometer (IMS) measurements, we present the ion fluid properties and its ion composition of the upstream flow for Titan's interaction with Saturn's magnetosphere. A 3D ion moments algorithm is used which is essentially model independent with only requirement is that ion flow is within the CAPS IMS 2(pi) steradian field-of-view (FOV) and that the ion 'velocity distribution function (VDF) be gyrotropic. These results cover the period from TA flyby (2004 day 300) to T22 flyby (2006 363). Cassini's in situ measurements of Saturn's magnetic field show it is stretched out into a magnetodisc configuration for Saturn Local Times (SLT) centered about midnight local time. Under those circumstances the field is confined near the equatorial plane with Titan either above or below the magnetosphere current sheet. Similar to Jupiter's outer magnetosphere where a magnetodisc configuration applies, one expects the heavy ions within Saturn's outer magnetosphere to be confined within a few degrees of the current sheet while at higher magnetic latitudes protons should dominate. We show that when Cassini is between dusk-midnight-dawn local time and spacecraft is not within the current sheet that light ions (H, 142) tend to dominate the ion composition for the upstream flow. If true, one may expect the interaction between Saturn's magnetosphere, locally devoid of heavy ions and Titan's upper atmosphere and exosphere to be significantly different from that for Voyager 1, TA and TB when heavy ions were present in the upstream flow. We also present observational evidence for Saturn's magnetosphere interaction with Titan's extended H and H2 corona which can extend approx. 1 Rs from Titan.

Sittler, E. C., Jr.↗

Lag Model Predictions for UFAST SBLI Flowfield

Presentation for Shock Boundary Layer Interaction Workshop. Shows results for Lag turbulence model on one of the international workshop test cases the UFAST 8 degree test case. Comparison with PIV velocity measurements as well as computed tunnel wall flowfields are shown, emphasizing that the interaction is a 3D phenomena, and is reasonably well predicted.

Olsen, Mike↗