Search NASA⌕ Search

SEARCH · Search NASA

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 415 records · Page 23

Supersonic flow calculation using a Reynolds-stress and an eddy thermal diffusivity turbulence model

A second-order model for the velocity field and a two-equation model for the temperature field are used to calculate supersonic boundary layers assuming negligible real gas effects. The modeled equations are formulated on the basis of an incompressible assumption and then extended to supersonic flows by invoking Morkovin's hypothesis, which proposes that compressibility effects are completely accounted for by mean density variations alone. In order to calculate the near-wall flow accurately, correction functions are proposed to render the modeled equations asymptotically consistent with the behavior of the exact equations near a wall and, at the same time, display the proper dependence on the molecular Prandtl number. Thus formulated, the near-wall second order turbulence model for heat transfer is applicable to supersonic flows with different Prandtl numbers. The model is validated against flows with different Prandtl numbers and supersonic flows with free-stream Mach numbers as high as 10 and wall temperature ratios as low as 0.3. Among the flow cases considered, the momentum thickness Reynolds number varies from approximately 4,000 to approximately 21,000. Good correlation with measurements of mean velocity, temperature, and its variance is obtained. Discernible improvements in the law-of-the-wall are observed, especially in the range where the big-law applies.

Sommer, T. P.↗

Evaluation of the maximum cross-correlation method of estimating sea surface velocities from sequential satellite images

The maximum cross correlation (MCC) method of Leese et al. (1971) for estimating sea surface velocities from sequences of satellite images is evaluated by comparing the MCC fields obtained from sequences of AVHRR and CZCS with in situ data and velocity fields calculated with a high-resolution quasi-geostrophic model. The rms differences and vector correlations between the velocity field produced by the MCC method and the model's field are presented. It is shown that much of the difference between the MCC fields and either the in situ data or the model velocity fields can be accounted for by considering physical and biological processes not included in the MCC method. The conditions under which the method is likely to be most successful are discussed.

Tokmakian, Robin↗

Towards LES Models of Jets and Plumes

As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Webb, A. T.↗

Phenomena Simulation for Heavy Doping and Surface Recombination Velocity

The theoretical models now available that characterize heavily doped (highly conducting) regions in silicon are survyed. Analytical and numerical approaches that determine the influence of such regions on the conversion efficiency of solar cells are examined. Although dilutely doped silicon is well characterized except for some disagreement about optical absorption coefficients, what exists now for heavily doped silicon and its interplay with adjoining regions is an incomplete theory in which not all contributers to transport, recombination, generation, and trapping are defined. Further, the parameters relating to these mechanisms and their values as determined by experiment are subject to various interpretations. The characterization of heavily doped silicon is treated not as a theory but rather as an imperfectly articulated and incompletely formalized body of experience. This view is intended to help point the way toward the attainment of a more complete of heavily doped silicon and thereby toward more informed designs of solar cells. Because computer programs constitute tools both for design and for estimating performance limits, the review includes some remarks pertinent to existing and developing programs.

Lindholm, F. A.↗

Deterrent Concentration Measurement with FTIR and Subsequent Ballistic Performance in Medium Caliber Propellant

To meet ballistic requirements, medium and small caliber propellants use deterrent coatings to obtain burn rate progressivity. The required amount and distribution of deterrent varies between gun systems, propellant types, and often between lots. Micro Fourier Transform Infrared (FTIR) spectroscopy was used to measure deterrent gradients in RP36 propellants coated with methyl centralite (MC) at different deterrent levels and different processing conditions. The aromatic C-C bonds at 1496 cm(exp -1) wavenumber were used to monitor the deterrent profiles through the grain. Deterrent gradients measured with FTIR spectroscopy were then used to estimate burn rate gradients in the deterred grains. Burn rates were calculated from literature models and from closed bomb data of RP36 containing uniform deterrent concentration. Finally, the burn rate gradients were input into an IBHFG2 model of a 200 cc-closed bomb. The early flame spreading portion of the closed bomb ballistic cycle (0 to 0.2 P/Pmax) was roughly modeled by dividing the charge up into five propellant decks and igniting them at different times in the ballistic cycle. Pressure traces and vivacity curves from closed bomb shots were compared to predictions. In addition to the burn rate gradient, the closed bomb pressure trace was heavily dependent on ignition and flame spread. These two phenomena were not readily distinguishable from one another in deterred grains. The same RP-36 propellant was shot in a 25 mm M793TP round which was again modeled with IBHVG2. Peak pressure and muzzle velocity were accurately modeled when erosive burning effects were empirically factored into the model.

Furrow, Keith W.↗

A Dissection of Spatially Resolved AGN Feedback across the Electromagnetic Spectrum

We present optical SuperNova Integral Field Spectrograph integral field spectroscopy, Hubble Space Telescope optical imaging, Chandra X-ray imaging, and Very Large Array radio interferometry of the merging galaxy 2MASX J04234080+0408017, which hosts a Seyfert 2 active galactic nucleus (AGN) at z = 0.046. Our observations reveal that radiatively driven, ionized gas outflows are successful to distances >10 kpc due to the low mass of the host system, encompassing the entirety of the observed optical emission. We also find that at large radii, where observed velocities cannot be reproduced by radiative driving models, high-velocity kinematics are likely due to mechanical driving from AGN winds impacting high-density host material. This impacting deposits sufficient energy to shock the host material, producing thermal X-ray emission and cosmic rays, which in turn promote the formation of in situ radio structure in a pseudo-jet morphology along the high-density lanes.

Travis Fischer↗

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗

On the choice of boundary conditions in continuum models of continental deformation

Recent studies of continental deformation have treated the lithosphere as a viscous medium and investigated the time evolution of the deformation caused by tectonic and buoyancy forces. This paper examines the differences between (1) continuum models that keep velocity boundary conditions constant with time and (2) models that keep stress boundary conditions constant with time. These differences are demonstrated by using a simple example of a continental lithosphere that is subjected to horizontal compression. The results show that in (2) the indentation velocity decreases with time, while in (1) the indentation velocity remains constant with time.

Wdowinski, Shimon↗

Determining the ballistic threshold velocity for a composite structure using multiple failure models

A composite structure, such as a laminated composite panel, for example, comprises one or more layers or “plies” embedded in a matrix material or otherwise fixed together in an arrangement, commonly referred to as a “stack up.” Each material in the structure has a corresponding material failure model (MFM) defining the physical characteristics of that material. A ballistic threshold velocity computing device obtains the MFMs for each material in the composite structure, generates a predicted ballistic velocity threshold velocity for each MFM, and then generates a parametric model to compute a composite ballistic velocity threshold velocity for the composite structure.

36 MATERIALS SCIENCE↗

A Tidal Hydrodynamic Model for Cook Inlet, Alaska, to Support Tidal Energy Resource Characterization

Cook Inlet in Alaska has been identified as a prime site in the U.S. for potential tidal energy development, because of its enormous tidal power potential that accounts for nearly one-third of the national total. As one important step to facilitate tidal energy development, a tidal hydrodynamic model based on the unstructured-grid, finite-volume community ocean model (FVCOM) was developed for Cook Inlet to characterize the tidal stream energy resource. The model has a grid resolution that varies from about 1000 m at the open boundary to 100–300 m inside the Inlet. Extensive model validation was achieved by comparing model predictions with field observations for tidal elevation and velocity at various locations in Cook Inlet. The error statistics confirmed the model performs reasonably well in capturing the tidal dynamics in the system, e.g., R2 > 0.98 for tidal elevation and generally > 0.9 for velocity. Model results suggest that tides in Cook Inlet evolve from progressive waves at the entrance to standing waves at the upper Inlet, and that semi-diurnal tidal constituents are amplified more rapidly than diurnal constituents. The model output was used to identify hotspots that have high energy potential and warrant additional velocity and turbulence measurements such as East Foreland, where averaged power density exceeds 5 kw/m2. Lastly, a tidal energy extraction simulation was conducted for a hypothetical turbine farm configuration at the Forelands cross section to evaluate tidal energy extraction and resulting changes in far-field hydrodynamics.

16 TIDAL AND WAVE POWER↗

A Comparison of Quasi-Static Indentation Testing to Low Velocity Impact Testing

The need for a static test method for modeling low-velocity foreign object impact events to composites would prove to be very beneficial to researchers since much more data can be obtained from a static test than from an impact test. In order to examine if this is feasible, a series of static indentation and low velocity impact tests were carried out and compared. Square specimens of many sizes and thickness were utilized to cover the array of types of low velocity impact events. Laminates with a n/4 stacking sequence were employed since this is by the most common type of engineering laminate. Three distinct flexural rigidities under two different boundary conditions were tested in order to obtain damage due to large deflections, contact stresses and both to examine if the static indentation-impact comparisons are valid under the spectrum of damage modes that can be experienced. Comparisons between static indentation and low velocity impact tests were based on the maximum applied transverse load. The dependent parameters examined included dent depth, back surface crack length, delamination area and to a limited extent, load-deflection behavior. Results showed that no distinct differences could be seen between the static indentation tests and the low velocity impact tests, indicating that static indentation can be used to represent a low velocity impact event.

Nettles, Alan T.↗

GRAIL Refinements to Lunar Seismic Structure

The present ]day internal structure of the Moon provides insight not only into its own formation and evolution, but also that of all rocky planetary bodies. The most direct way to probe a planet fs interior structure is through seismology. As part of the Apollo lunar missions, four seismometers were deployed on the nearside surface of the Moon between the years 1969 and 1972. These instruments operated continuously until 1977, forming the only substantial extraterrestrial seismic data set in existence. These data have been used to constrain various aspects of the seismic velocity and density structure of the Moon. Typical 1-D models recognize a 30-60 km thick crust overlying a nearly constant ]velocity mantle, and extend to a depth of approximately 1000 km, below which the lack of penetrating moonquake ray ]paths precludes the seismic determination of deeper structure. Previously, the lack of observed moonquakes from the far side of the Moon has been used to infer the presence of a highly attenuating (possibly molten) core. Indirect geophysical measurements such as moment of inertia, magnetic induction, lunar laser ranging, and elemental abundances of mare basalts also place varying constraints on core size and state. In combination with seismic studies, these indirect measurements have been used to arrive at a commonly accepted model of the Moon's deepest interior that includes a solid inner and fluid outer core, overlain by a partial melt boundary layer. We recently applied modern array seismology techniques to the Apollo data and revealed detailed core structure, including the first direct confirmation of the presence of a solid inner core. Our study focused on the identification of core ]reflected phases in deep moonquake seismograms. The resulting model of the Moon fs innermost structure was found to be consistent with the commonly accepted model. However, the modeled layer radii may vary by tens of kilometers, as is expected when accounting for uncertainties such as moonquake location, timing errors, and potential seismic heterogeneities. In addition, the modeled velocities may vary with a 1-to-1 trade ]off with the modeled reflector depth. The GRAIL (Gravity Recovery and Interior Laboratory) mission, launched in Sept. 2011, placed two nearly identical spacecraft in lunar orbit. The two satellites make extremely high-resolution measurements of the lunar gravity field, which can be used to constrain the interior structure of the Moon using a "crust to core" approach. GRAIL fs constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten) will complement seismic investigations in several ways. Here we present a progress report on our efforts to advance our knowledge of the Moon fs internal structure using joint gravity and seismic analyses. We will focus on methodology, including 1) refinements to the seismic core constraint accomplished through array processing of Apollo seismic data, made by applying a set of travel time corrections based on GRAIL structure estimates local to each Apollo seismic station; 2) modeling deep lunar structure through synthetic seismograms, to test whether the seismic core model can reproduce the core reflections observed in the Apollo seismograms; and 3) a joint seismic and gravity inversion in which we attempt to fit a family of seismic structure models with the gravity constraints from GRAIL, resulting in maps of seismic velocities and densities that vary from a nominal model both laterally and with depth.

Weber, Renee↗

Experimental stress analysis of large plastic deformations in a hollow sphere deformed by impact against a concrete block

An experimental plastic strain measurement system is presented for use on the surface of high velocity impact test models. The system was used on a hollow sphere tested in impact against a reinforced concrete block. True strains, deviatoric stresses, and true stresses were calculated from experimental measurements. The maximum strain measured in the model was small compared to the true failure strain obtained from static tensile tests of model material. This fact suggests that a much greater impact velocity would be required to cause failure of the model shell structure.

Morris, R. E.↗

Velocity Fourier transform solution of a model collision operator

A simpler technique than those introduced by Lenard and Bernstein (1958), and Dougherty (1964) is employed to obtain the perturbed species density from a specified kinetic equation for a plasma in a given uniform magnetic field. The technique is a generalization of the velocity-Fourier transform method employed by Karpman (1967) for B sub 0 identical to zero, and relies on the fact that in transform space the model collision operator, used to obtain the kinetic equation for waves in a magnetized plasma, contains only first derivatives. The technique is illustrated by evaluating the perturbed density of an arbitrary species.

Catto, P. J.↗

New Hubble Space Telescope Observations of High-Velocity Ly(alpha) and H(alpha) in SNR 1987A

We describe and model high-velocity (approximately 15,000 km S(exp -1)) Ly Alpha and H Alpha emission from the supernova remnant SNR 1987A seen in 1997 September and October with the Space Telescope Imaging Spectrograph. Part of this emission comes from a reverse shock located at approximately 75% of the radius of the inner boundary ofthe innercircumstellar ring and confined within + or - 30 degrees of the equatorial plane. Departure from axisymmetry in the Ly Alpha and H Alpha emission correlates with that seen in nonthermal radio emission and reveals an asymmetry in the circumstellar gas distribution. We also see diffuse high-velocity Ly-Alpha emission from supernova debris inside the reverse shock that may be due to excitation by nonthermal particles accelerated by the shock.

Michael, Eli↗

A transport model of the turbulent scalar-velocity

Performance tests of the third-order turbulence closure for predictions of separating and recirculating flows in backward-facing steps were studied. Computations of the momentum and temperature fields in the flow domain being considered entail the solution of time-averaged transport equations containing the second-order turbulent fluctuating products. The triple products, which are responsible for the diffusive transport of the second-order products, attain greater significance in separating and reattaching flows. The computations are compared with several algebraic models and with the experimental data. The prediction was improved considerably, particularly in the separated shear layer. Computations are further made for the temperature-velocity double products and triple products. Finally, several advantages were observed in the usage of the transport equations for the evaluation of the turbulence triple products; one of the most important features is that the transport model can always take the effects of convection and diffusion into account in strong convective shear flows such as reattaching separated layers while conventional algebraic models cannot account for these effects in the evaluation of turbulence variables.

Amano, R. S.↗