Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Velocity Model”

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 1,783 records · Page 99

Noise from Supersonic Coaxial Jets: Mean Flow Predictions - Part 1

Recent theories for supersonic jet noise have used an instability wave noise generation model to predict radiated noise. This model requires a known mean flow that has typically been described by simple analytic functions for single jet mean flows. The mean flow of supersonic coaxial jets is not described easily in terms of analytic functions. To provide these profiles at all axial locations, a numerical scheme is developed to calculate the mean flow properties of a coaxial jet. The Reynolds-averaged, compressible, parabolic boundary layer equations are solved using a mixing length turbulence model. Empirical correlations are developed to account for the effects of velocity and temperature ratios and Mach number on the shear layer spreading. Both normal velocity profile and inverted velocity profile coaxial jets are considered. The mixing length model is modified in each case to obtain reasonable results when the two stream jet merges into a single fully developed jet. The mean flow calculations show both good qualitative and quantitative agreement with measurements in single and coaxial jet flows.

Dahl, Milo D.↗

PREDICTING CME EJECTA AND SHEATH FRONT ARRIVAL AT L1 WITH A DATA-CONSTRAINED PHYSICAL MODEL

We present a method for predicting the arrival of a coronal mass ejection (CME) flux rope in situ, as well as the sheath of solar wind plasma accumulated ahead of the driver. For faster CMEs, the front of this sheath will be a shock. The method is based upon geometrical separate measurement of the CME ejecta and sheath. These measurements are used to constrain a drag-based model, improved by including both a height dependence and accurate de-projected velocities. We also constrain the geometry of the model to determine the error introduced as a function of the deviation of the CME nose from the Sun–Earth line. The CME standoff-distance in the heliosphere fit is also calculated, fit, and combined with the ejecta model to determine sheath arrival. Combining these factors allows us to create predictions for both fronts at the L1 point and compare them against observations. We demonstrate an ability to predict the sheath arrival with an average error of under 3.5 hr, with an rms error of about 1.58 hr. For the ejecta the error is less than 1.5 hr, with an rms error within 0.76 hr. We also discuss the physical implications of our model for CME expansion and density evolution. We show the power of our method with ideal data and demonstrate the practical implications of having a permanent L5 observer with space weather forecasting capabilities, while also discussing the limitations of the method that will have to be addressed in order to create a real-time forecasting tool.

79 ASTRONOMY AND ASTROPHYSICS↗

Multifrequency spectra of BL Lac objects

A program to obtain simultaneous multifrequency spectra of BL Lacertae objects that are known X-ray sources is discussed. The IUE spectra are generally featureless and well-fitted by power law models. For the faintest exposures, Gaussian extraction of the spectrum can greatly impove the signal-to-noise. Most program objects vary in the ultraviolet, although the time scales are not known because of limited observing time. The broadband spectra of BL Lacs exhibit a range of characteristics but the curvature is always downward and the shape is generally smooth. This can be interpreted as synchrotron emission from a relativistic jet; different jet models are possible, and each allows a range of values for the bulk velocity, magnetic field strength, and electron density. Synchrotron models are not required, however an accretion disk model also gives a good fit to the ultraviolet-through-X-ray continuum.

Urry, C. M.↗

Model Refinement Using Eigensystem Assignment

This paper presented a novel approach for the refinement of finite-element-based analytical models of flexible structures is presented. The proposed approach models the possible refinements in the mass, damping, and stiffness matrices of the finite element model in the form of a constant gain feedback with acceleration, velocity, and displacement measurements, respectively. Once, the free elements of the structural matrices have been defined, the problem of model refinement reduces to obtaining position, velocity, and acceleration gain matrices, which reassign a desired subset of the eigenvalues of the model, along with partial mode shapes, from their baseline values to those obtained from system identification test data. A sequential procedure is used to assign one self-conjugate pair of closed-loop eigenvalues at each step using symmetric output feedback gain matrices, and the closed-loop eigenvectors are partially assigned, while ensuring that the eigenvalues assigned in the previous steps are not disturbed. The procedure can also impose that gain matrices be dissipative in order to guarantee the stability of the refined model. A numerical example, involving finite element model refinement for a structural testbed at NASA Langley (CSI Evolutionary Model) is presented to demonstrate the feasibility of the proposed approach.

Maghami, Pieman G.↗

Phase space distribution of halo particles and detection of WIMPS

We calculate the velocity distribution of halo particles using realistic models for the density profiles of the Galactic disk, spheroid, and halo components. We get the resulting velocity distribution, which is different from the Maxwellian distribution of isothermal halo. We investigate the consequences concerning the time averaged count rate and the amplitude of the annual modulation in WIMPS detection experiments. Compared to the previous estimate based on the Maxwellian distribution, the time averaged count rate is up to 40 percent larger. On the other hand, the amplitude of the annual modulation is up to a factor of two larger.

Ryu, Dongsu↗

Spatiotemporal changes of seismic attenuation caused by injected CO 2 at the Frio-II pilot site, Dayton, TX, USA

A continuous active source seismic monitoring data set was collected with crosswell geometry during CO 2 injection at the Frio-II brine pilot, near Liberty, TX. Previous studies have shown that spatiotemporal changes in the P wave first arrival time reveal the movement of the injected CO 2 plume in the storage zone. To further constrain the CO 2 saturation, particularly at higher saturation levels, we investigate spatial-temporal changes in the seismic attenuation of the first arrivals. The attenuation changes over the injection period are estimated by the amount of the centroid frequency shift computed by local time-frequency analysis. We observe that (1) at receivers above the injection zone seismic attenuation does not change in a physical trend; (2) at receivers in the injection zone attenuation sharply increases following injection and peaks at specific points varying with distributed receivers, which is consistent with observations from time delays of first arrivals; then, (3) attenuation decreases over the injection time. The attenuation change exhibits a bell-shaped pattern during CO 2 injection. Under Frio-II field reservoir conditions, White's patchy saturation model can quantitatively explain both the P wave velocity and attenuation response observed. We have combined the velocity and attenuation change data in a crossplot format that is useful for model-data comparison and determining patch size. Our analysis suggests that spatial-temporal attenuation change is not only an indicator of the movement and saturation of CO 2 plumes, even at large saturations, but also can quantitatively constrain CO 2 plume saturation when used jointly with seismic velocity.

58 GEOSCIENCES↗

Steady detonation propagation in a circular arc: a Detonation Shock Dynamics model

We study the physics of steady detonation wave propagation in a two-dimensional circular arc via a Detonation Shock Dynamics (DSD) surface evolution model. The dependence of the surface angular speed and surface spatial structure on the inner arc radius ( $R_{i}$ ), the arc thickness ( $R_{e}-R_{i}$ , where $R_{e}$ is the outer arc radius) and the degree of confinement on the inner and outer arc is examined. We first analyse the results for a linear $D_{n}$ – $\unicode[STIX]{x1D705}$ model, in which the normal surface velocity $D_{n}=D_{CJ}(1-B\unicode[STIX]{x1D705})$ , where $D_{CJ}$ is the planar Chapman–Jouguet velocity, $\unicode[STIX]{x1D705}$ is the total surface curvature and $B$ is a length scale representative of a reaction zone thickness. An asymptotic analysis assuming the ratio $B/R_{i}\ll 1$ is conducted for this model and reveals a complex surface structure as a function of the radial variation from the inner to the outer arc. For sufficiently thin arcs, where $(R_{e}-R_{i})/R_{i}=O(B/R_{i})$ , the angular speed of the surface depends on the inner arc radius, the arc thickness and the inner and outer arc confinement. For thicker arcs, where $(R_{e}-R_{i})/R_{i}=O(1)$ , the angular speed does not depend on the outer arc radius or the outer arc confinement to the order calculated. It is found that the leading-order angular speed depends only on $D_{CJ}$ and $R_{i}$ , and corresponds to a Huygens limit (zero curvature) propagation model where $D_{n}=D_{CJ}$ , assuming a constant angular speed and perfect confinement on the inner arc surface. Having the normal surface speed depend on curvature requires the insertion of a boundary layer structure near the inner arc surface. This is driven by an increase in the magnitude of the surface wave curvature as the inner arc surface is approached that is needed to meet the confinement condition on the inner arc surface. For weak inner arc confinement, the surface wave spatial variation with the radial coordinate is described by a triple-deck structure. The first-order correction to the angular speed brings in a dependence on the surface curvature through the parameter $B$ , while the influence of the inner arc confinement on the angular velocity only appears in the second-order correction. For stronger inner arc confinement, the surface wave structure is described by a two-layer solution, where the effect of the confinement on the angular speed is promoted to the first-order correction. We also compare the steady-state arc solution for a PBX 9502 DSD model to an experimental two-dimensional arc geometry validation test.

Short, Mark↗

Intercomparison of 3D pore-scale flow and solute transport simulation methods

Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include methods that 1) explicitly model the three-dimensional geometry of pore spaces and 2) those that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of class 1, based on direct numerical simulation using computational fluid dynamics (CFD) codes, against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of class 1 based on the immersed-boundary method (IMB), lattice Boltzmann method (LBM), smoothed particle hydrodynamics (SPH), as well as a model of class 2 (a pore-network model or PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and nonreactive solute transport, and intercompare the model results with previously reported experimental observations. Experimental observations are limited to measured pore-scale velocities, so solute transport comparisons are made only among the various models. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations).

computational fluid dynamics↗

Planck intermediate results: XXXVII. Evidence of unbound gas from the kinetic Sunyaev-Zeldovich effect

By looking at the kinetic Sunyaev-Zeldovich effect (kSZ) in Planck nominal mission data, we present in this paper a significant detection of baryons participating in large-scale bulk flows around central galaxies (CGs) at redshift z ≈ 0.1. We estimate the pairwise momentum of the kSZ temperature fluctuations at the positions of the Central Galaxy Catalogue (CGC) samples extracted from Sloan Digital Sky Survey (SDSS-DR7) data. For the foreground-cleaned SEVEM, SMICA, NILC, and COMMANDER maps, we find 1.8–2.5σ detections of the kSZ signal, which are consistent with the kSZ evidence found in individualPlanck raw frequency maps, although lower than found in the WMAP-9yr W-band (3.3σ). We further reconstruct the peculiar velocity field from the CG density field, and compute for the first time the cross-correlation function between kSZ temperature fluctuations and estimates of CG radial peculiar velocities. This correlation function yields a 3.0–3.7σ detection of the peculiar motion of extended gas on Mpc scales in flows correlated up to distances of 80–100 h -1 Mpc. Both the pairwise momentum estimates and the kSZ temperature-velocity field correlation find evidence for kSZ signatures out to apertures of 8 arcmin and beyond, corresponding to a physical radius of >1 Mpc, more than twice the mean virial radius of halos. This is consistent with the predictions from hydrodynamical simulations that most of the baryons are outside the virialized halos. We fit a simple model, in which the temperature-velocity cross-correlation is proportional to the signal seen in a semi-analytic model built upon N-body simulations, and interpret the proportionality constant as an effective optical depth to Thomson scattering. Finally, we find τ T = (1.4 ± 0.5) × 10 -4 ; the simplest interpretation of this measurement is that much of the gas is in a diffuse phase, which contributes little signal to X-ray or thermal Sunyaev-Zeldovich observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Dust particles in the anomalous tail of Comet Kohoutek (1973 XII)

Intermediate results from theoretical models of the anomalous spikelike dust distribution observed in the tail of Comet Kohoutek using the Skylab white-light coronograph on December 28.996, 1973, are reported. The time evolution is studied using the model of Finson and Probstein (1968) and assuming power-law variation of the dust-production rate and the particle-size distribution. Brightness profiles are calculated and presented graphically for comparison with those derived from the observation photograph: the model predicts a longer duration of the spike phenomenon than was observed. It is shown qualitatively that this discrepancy can be removed by introducing a finite dust-emission velocity (e.g., 1 m/s) into the model.

Keller, H. U.↗

Laser transit anemometer measurements on a slender cone in the Langley unitary plan wind tunnel

A laser transit anemometer (LTA) system was used to probe the boundary layer on a slender (5 degree half angle) cone model in the Langley unitary plan wind tunnel. The anemometer system utilized a pair of laser beams with a diameter of 40 micrometers spaced 1230 micrometers apart to measure the transit times of ensembles of seeding particles using a cross-correlation technique. From these measurements, boundary layer profiles around the model were constructed and compared with CFD calculations. The measured boundary layer profiles representing the boundary layer velocity normalized to the edge velocity as a function of height above the model surface were collected with the model at zero angle of attack for four different flow conditions, and were collected in a vertical plane that bisected the model's longitudinal center line at a location 635 mm from the tip of the forebody cone. The results indicate an excellent ability of the LTA system to make velocity measurements deep into the boundary layer. However, because of disturbances in the flow field caused by onboard seeding, premature transition occurred implying that upstream seeding is mandatory if model flow field integrity is to be maintained. A description and results of the flow field surveys are presented.

Humphreys, William M., Jr.↗

Hybrid (Vlasov-Fluid) simulation of ion-acoustic solitons chain formation including trapped electrons

Disintegration of a Gaussian profile into ion-acoustic solitons in the presence of trapped electrons [H. Hakimi Pajouh and H. Abbasi, Phys. Plasmas 15, 082105 (2008)] is revisited. Through a hybrid (Vlasov-Fluid) model, the restrictions associated with the simple modified Korteweg de-Vries (mKdV) model are studied. For instance, the lack of vital information in the phase space associated with the evolution of electron velocity distribution, the perturbative nature of mKdV model which limits it to the weak nonlinear cases, and the special spatio-temporal scaling based on which the mKdV is derived. Remarkable differences between the results of the two models lead us to conclude that the mKdV model can only monitor the general aspects of the dynamics, and the precise picture including the correct spatio-temporal scales and the properties of solitons should be studied within the framework of hybrid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The characteristics of a lift cruise fan V/STOL configuration in near proximity to a small deck with finite edge positions

Two V/STOL model configurations utilizing lift-cruise-fans for the propulsion system have been tested in the presence of a small landing platform with finite deck edge positions. Results and analysis of these tests are presented. The model was tested at both static (zero velocity) and forward speed conditions with several heights and deck edge positions. Configurations utilizing three fans and four fans were investigated. Incremental forces and moments induced by the various deck edge positions are analyzed as relating to the model height above the ground and the thrust velocity relationships. The test and analysis show the induced forces to be relatively unaffected by deck edge.

Stewart, V. R.↗

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps to reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co‐located at one corner. Due to this limited geographical extent, near‐surface ray coverage from moonquakes is low, but increases with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefore well suited to joint inversion, since the complementary information reduces inherent model ambiguity. Previous joint inversions of the Apollo seismic data (seismic phase arrival times) and Clementine‐ or Lunar Prospector‐derived gravity data (mass and moment of inertia) attempted to recover the subsurface structure of the Moon by focusing on hypothetical lunar compositions that explore the density/velocity relationship. These efforts typically search for the best fitting thermodynamically calculated velocity/density model, allowing variables like core size, velocity, and/or composition to vary freely. Seismic velocity profiles previously derived from the Apollo seismic data through inversion of travel times vary both in the depth of the crust and mantle layers, and the seismic velocities and densities assigned to those layers. The lunar mass and moment of inertia likewise only constrain gross variations in the density profile beyond that of a uniform density sphere. As a result, composition and structure models previously obtained by jointly inverting these data retain the original uncertainties inherent in the input data sets. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by the GRAIL lunar gravity mission, in order to recover seismic velocities and density as a function of latitude, longitude, and depth within the Moon. We will relate density to seismic velocity using a linear relationship that is allowed to be depth‐dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of density, velocity, and B‐coefficient perturbations that minimize (in a least‐squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗

Diffusion and chaos from near AdS 2 horizons

We calculate the thermal diffusivity D = $\kappa/c_\rho$ and butterfy velocity $\upsilon_\beta$ in holographic models that flow to $AdS_2$ x $R^d$ fixed points in the infra-red. We show that both these quantities are governed by the same irrelevant deformation of $AdS_2$ and hence establish a simple relationship between them. When this deformation corresponds to a universal dilaton mode of dimension $\Delta$ = 2 then this relationship is always given by D</> = $\upsilon_B^2$/(2$\pi$T).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Building blocks for a digital twin of additive manufacturing

Properties and serviceability of additively manufactured components are affected by their geometry, microstructure and defects. These important attributes are now optimized by trial and error because the essential process variables cannot currently be selected from scientific principles. A recourse is to build and rigorously validate a digital twin of the additive manufacturing process that can provide accurate predictions of the spatial and temporal variations of metallurgical parameters that affect the structure and properties of components. Key building blocks of a computationally efficient first-generation digital twin of laser-based directed energy deposition additive manufacturing utilize a transient, three-dimensional model that calculates temperature and velocity fields, cooling rates, solidification parameters and deposit geometry. The measured profiles of stainless steel 316L and Alloy 800H deposits as well as the secondary dendrite arm spacing (SDAS) and Vickers hardness measurements are used to validate the proposed digital twin. The predicted cooling rates, temperature gradients, solidification rates, SDAS and micro-hardness values are shown to be more accurate than those obtained from a commonly used heat conduction calculation. Finally, these metallurgical building blocks serve as a phenomenological framework for the development of a digital twin that will make the expanding knowledge base of additive manufacturing usable in a practical way for all scientists and engineers.

36 MATERIALS SCIENCE↗

Natural attenuation of uranium in a fluvial Wetland: Importance of hydrology and speciation

A nuclear fuel fabrication facility released 43,500 kg of uranium into a riparian wetland located on the Savannah River Site between 1955 and 1988. Studies were undertaken to evaluate hydrological and geochemical processes influencing uranium accumulation in the wetland. Gamma-radiation-mapping surveys were conducted by systematically walking over the contaminated wetland with backpacks equipped with global positioning systems and NaI gamma detectors. Based on maps compiled from >700,000 gamma spectra and eight sediment uranium depth profiles, it was determined that 94% of the released uranium remained in the wetland. The uranium in the wetland is concentrated in five multi-hectare areas along the stream, accounting for ~11% of the land area adjacent to the stream. While land type (upland or wetland) and topography provided a reasonable first approximation of where much of the uranium was deposited, hydrological watershed modeling revealed that the stream velocity was especially slow through many of the hot spots. Here, using autoradiography combined with SEM/EDX measurements of contaminated sediments, surprisingly few hot particles were detected. Instead, uranium was evenly distributed throughout the sampled sediment, suggesting that dissolved uranium had bound to sediment particles that became suspended and later deposited in low energy (low flow velocity) portions of the stream. EXAFS suggested that U atoms were present as individual ions in disordered complexes within the sediment. Furthermore, linear combination analyses suggested that the predominant component of the U(VI) was adsorbed to sediment minerals (~70%) and a minor component (~30%) was associated with organic matter phases. Furthermore, these studies show that wetlands can be extraordinarily effective at binding and retaining uranium, thereby providing a natural barrier to the transport of uranium out of a watershed. However, significant anthropogenic or climatic changes to wetlands, such as those associated with flooding, forest fires, or land use, may disrupt the complex hydrological and biogeochemical balance necessary to maintain long-term immobilization of uranium.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗