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 721 records · Page 40

Detailed Velocity, Temperature, and Heat Flux Measurements on a Large Scale Film Cooling Model

To investigate the flow physics of turbine film cooling, detailed mean and fluctuating velocity and temperature surveys were made in the turbulent jet flow field behind a row of large scale film cooling holes. Additionally, surface heat transfer, film effectiveness, and turbulent heat flux values were determined. Measurements were made in the NASA Engine Research Building (ERB) SW-6 test facility, which consists of an 8.15-in2 open-inlet wind tunnel with a temperature-controlled secondary flow system. The film-cooling test plate fit into the wind tunnel such that the heat transfer surface formed the wind tunnel floor. This plate consisted of a three-hole array of film cooling holes that were fed from a plenum of different temperature.

Fluid Mechanics↗

Light curves for bump Cepheids computed with a dynamically zoned pulsation code

The dynamically zoned pulsation code developed by Castor, Davis, and Davison was used to recalculate the Goddard model and to calculate three other Cepheid models with the same period (9.8 days). This family of models shows how the bumps and other features of the light and velocity curves change as the mass is varied at constant period. The use of a code that is capable of producing reliable light curves demonstrates that the light and velocity curves for 9.8 day Cepheid models with standard homogeneous compositions do not show bumps like those that are observed unless the mass is significantly lower than the 'evolutionary mass.' The light and velocity curves for the Goddard model presented here are similar to those computed independently by Fischel, Sparks, and Karp. They should be useful as standards for future investigators.

Adams, T. F.↗

A non-gaussian model of continuous atmospheric turbulence for use in aircraft design

A non-Gaussian model of atmospheric turbulence is presented and analyzed. The model is restricted to the regions of the atmosphere where the turbulence is steady or continuous, and the assumptions of homogeneity and stationarity are justified. Also spatial distribution of turbulence is neglected, so the model consists of three independent, stationary stochastic processes which represent the vertical, lateral, and longitudinal gust components. The non-Gaussian and Gaussian models are compared with experimental data, and it is shown that the Gaussian model underestimates the number of high velocity gusts which occur in the atmosphere, while the non-Gaussian model can be adjusted to match the observed high velocity gusts more satisfactorily. Application of the proposed model to aircraft response is investigated, with particular attention to the response power spectral density, the probability distribution, and the level crossing frequency. A numerical example is presented which illustrates the application of the non-Gaussian model to the study of an aircraft autopilot system. Listings and sample results of a number of computer programs used in working with the model are included.

Reeves, P. M.↗

Theoretical study of hull-rotor aerodynamic interference on semibuoyant vehicles

Analytical methods are developed to predict the pressure distribution and overall loads on the hulls of airships which have close coupled, relatively large and/or high disk loading propulsors for attitude control, station keeping, and partial support of total weight as well as provision of thrust in cruise. The methods comprise a surface-singularity, potential-flow model for the hull and lifting surfaces (such as tails) and a rotor model which calculates the velocity induced by the rotor and its wake at points adjacent to the wake. Use of these two models provides an inviscid pressure distribution on the hull with rotor interference. A boundary layer separation prediction method is used to locate separation on the hull, and a wake pressure is imposed on the separated region for purposes of calculating hull loads. Results of calculations are shown to illustrate various cases of rotor-hull interference and comparisons with small scale data are made to evaluate the method.

Spangler, S. B.↗

Localized and distributed boundary-layer receptivity to convected unsteady wake in free stream

Receptivity to a model convected disturbance in the presence of localized and distributed variations in wall geometry and wall-suction velocity is examined. The model free-stream disturbance corresponds to the time-harmonic wake of a vibrating ribbon that is placed at a suitable distance above the surface of a thin airfoil. The advantages of using this disturbance for experiments on receptivity to convected disturbances are outlined. A brief parametric study is presented for a flat-plate boundary layer. The study quantifies the effect of wake position as well as wake width; in addition, it should be helpful in the choice of an optimal setting for a controlled experiment of the above type, which the above parametric study shows as feasible.

Choudhari, Meelan↗

Convective and Turbulent Motions in Nonprecipitating Cu. Part II: LES Simulated Cloud Represented by a Starting Plume

The dynamic structure of a small trade wind cumulus (Cu) is analyzed using a novel approach. Cu developing in a shear-free environment is simulated by 10-m-resolution LES model with spectral bin microphysics. The aim is to clarify the dynamical nature of cloud updraft zone (CUZ) including entrainment and mixing in growing Cu. The validity of concept stating that a cloud at developing state can be represented by a parcel or a jet is tested. To investigate dynamical entrainment in CUZ performed by motions with scales larger than the turbulence scales, the modeled fields of air velocity were filtered by wavelet filter that separated convective motions from turbulent ones. Two types of objects in developing cloud were investigated: small volume ascending at maximal velocity (point parcel) and CUZ. It was found that the point parcel representing the upper part of cloud core is adiabatic. The motion of the air in this parcel ascending from cloud base determines cloud-top height. The top-hat (i.e., averaged) values of updraft velocity and adiabatic fraction in CUZ are substantially lower than those in the point parcel. Evaluation of the terms in the dynamical equation typically used in 1D cloud parcel models show that this equation can be applied for calculation of vertical velocities at the developing stage of small Cu, at least up to the heights of the inversion layer. Dynamically, the CUZ of developing cloud resembles the starting plume with the tail of nonstationary jet. Both the top-hat vertical velocity and buoyancy acceleration linearly increase with the height, at least up to the inversion layer. An important finding is that lateral entrainment of convective (nonturbulent) nature has a little effect on the top-hat CUZ velocity and cannot explain the vertical changes of conservative variables qt and θl. In contrast, entrained air lifting inside CUZ substantially decreases top-hat liquid water content and its adiabatic fraction. Possible reasons of these effects are discussed.

54 ENVIRONMENTAL SCIENCES↗

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗

Velocity correlations of galaxy clusters

We determine the velocity correlation function, pairwise peculiar velocity difference, and rms pairwise peculiar velocity dispersion of rich clusters of galaxies, as a function of pair separation, for three cosmological models: Omega = 1 and Omega = 0.3 cold dark matter (CDM), and Omega = 0.3 primeval baryonic isocurvature (PBI) models (all flat and Cosmic Background Explorer (COBE)-normalized). We find that close cluster pairs, with separation r is less than or equal to 10/h Mpc, exhibit strong attractive peculiar velocities in all models; the cluster pairwise velocities depend sensitively on the model. The mean pairwise attractive velocity of clusters on 5/h Mpc scale ranges from approximately 1700 km/s for Omega = 1 CDM to approximately 1000 km/s for PBI to approximately 700 km/s for Omega = 0.3 CDM. The small-scale pairwise velocities depend also on cluster mass: richer, more massive clusters exhibit stronger attractive velocities than less massive clusters. On large scales, from approximately 20 to 200/h Mpc, the cluster peculiar velocities are increasingly dominated by bulk and random motions; they are independent of cluster mass. The cluster velocity correlation function is negative on small scales for Omega = 1 and Omega = 0.3 CDM, indicating strong pairwise motion relative to bulk motion on small scales; PBI exhibits relatively larger bulk motions. The cluster velocity correlation function is positive on very large scales, from r approximately 10/h Mpc to r approximately 200/h Mpc, for all models. These positive correlations, which decrease monotonically with scale, indicate significant bulk motions of clusters up to approximately 200/h Mpc. The strong dependence of the cluster velocity functions on models, especially at small separations, makes them useful tools in constraining cosmological models when compared with observations.

Cen, Renyue↗

Modeling of Cube Array Roughness: RANS, Large Eddy Simulation, and Direct Numerical Simulation

Abstract Flow over arrays of cubes is an extensively studied model problem for rough wall turbulent boundary layers. While considerable research has been performed in computationally investigating these topologies using direct numerical simulation (DNS) and large eddy simulation (LES), the ability of sublayer-resolved Reynolds-averaged Navier–Stokes (RANS) to predict the bulk flow phenomena of these systems is relatively unexplored, especially at low and high packing densities. Here, RANS simulations are conducted on six different packing densities of cubes in aligned and staggered configurations. The packing densities investigated span from what would classically be defined as isolated, up to those in the d-type roughness regime, filling in the gap in the present literature. Three different sublayer-resolved turbulence closure models were tested for each case: a low Reynolds number k–ϵ model, the Menter k–ω SST model, and a full Reynolds stress model. Comparisons of the velocity fields, secondary flow features, and drag coefficients are made between the RANS results and existing LES and DNS results. There is a significant degree of variability in the performance of the various RANS models across all comparison metrics. However, the Reynolds stress model demonstrated the best accuracy in terms of the mean velocity profile as well as drag partition across the range of packing densities.

Engineering↗

Stellar winds driven by multi-line scattering

This paper presents a model of a radiation-driven stellar wind with overlapping spectral lines. It is based on the Castor, Abbott, and Klein (CAK) theory. The presence of overlapping lines allows a photon to be scattered many times in different lines. The properties of the wind at any point depend on the wavelength-averaged intensity, which in turn depends on the structure of the wind. A self-consistent wind model is found. The mass loss rate does not saturate as line overlap becomes more pronounced, but continues to increase. The terminal velocity is much larger than in the CAK model, while the velocity law is shallower. This model might help explain the massive winds from Wolf-Rayet stars.

Friend, D. B.↗

Ionospheric Signature of Surface Mine Blasts from Global Positioning System Measurements

Sources such as atmospheric or buried explosions and shallow earthquakes are known to produce infrasonic pressure waves in the atmosphere. Because of the coupling between neutral particles and electrons at ionospheric altitudes, these acoustic and gravity waves induce variations of the ionospheric electron density. The Global Positioning System (GPS) provides a way of directly measuring the total electron content in the ionosphere and, therefore, of detecting such perturbations in the upper atmosphere. In July and August 1996, three large surface mine blasts (1.5 Kt each) were detonated at the Black Thunder coal mine in eastern Wyoming. As part of a seismic and acoustic monitoring- experiment, we deployed five dual-frequency GPS receivers at distances ranging from 50 to 200 km from the mine and were able to detect the ionospheric perturbation caused by the blasts. The perturbation starts 10 to 15 min after the blast, lasts for about 30 min, and propagates with an apparent horizontal velocity of 1200 meters per second. Its amplitude reaches 3 x 10 (exp 14) el per square meters in the 7-3 min period band, a value close to the ionospheric perturbation caused by the M = 6.7 Northridge earthquake. The small signal-to-noise ratio of the perturbation can be improved by slant-stacking the electron content time-series recorded by the different GPS receivers taking into account the horizontal propagation of the perturbation. The energy of the perturbation is concentrated in the 200 to 300 second period band, a result consistent with previous observations and numerical model predictions. The 300 second band probably corresponds to gravity modes and shorter periods to acoustic modes, respectively. Using a 1-D stratified velocity model of the atmosphere we show that linear acoustic ray tracing fits arrival times at all GPS receivers. We interpret the perturbation as a direct acoustic wave caused by the explosion itself. This study shows that even relatively small subsurface events can produce ionospheric perturbations that are above the detection threshold of the GPS technique. By sensing derivative signals, which can be detected over a relatively broad region, it appears that GPS might be particularly useful for source characterization within the first acoustic quiet zone where infrasound would probably be ineffective. This suggests that dual-frequency GPS monitoring could contribute to Comprehensive Nuclear Test Ban Treaty verification.

Calais, Eric↗

A novel transformation of the ice sheet Stokes equations and some of its properties and applications

We introduce a novel transformation of the Stokes equations into a form closely resembling the shallow Blatter–Pattyn equations. The two forms differ by only a few additional terms, while their variational formulations differ only by a single term in each horizontal direction. Specifically, the variational formulation of the Blatter–Pattyn model drops the vertical velocity in the second invariant of the strain rate tensor. Here we make use of the new transformation in two ways. First, we consider incorporating the transformed equations into a code that can be very easily converted from a Stokes to a Blatter–Pattyn model, and vice versa, by switching these terms on or off. This may be generalized so that the Stokes model is switched on adaptively only where the Blatter–Pattyn model loses accuracy. Second, the key role played by the vertical velocity in the Blatter–Pattyn approximation motivates new approximations. Two examples are presented. These require a mesh that enables the discrete continuity equation to be invertible for the vertical velocity in terms of the horizontal velocity components. Examples of such meshes, such as the first-order P1–E0 mesh and the second-order P2–E1 mesh, are given in both 2D and 3D. However, the transformed Stokes model has the same type of gravity forcing as the Blatter–Pattyn model, determined by the ice surface slope, thereby forgoing some of the mesh generality of the traditional formulation of the Stokes model.

58 GEOSCIENCES↗

On the local standard of rest

Under the influence of a spiral gravitational field, there should be differences among the mean motions of different types of objects with different dispersion velocities in a sipral galaxy. The old stars with high dispersion velocity should have essentially no mean motion normal to the galactic rotation. On the other hand, young objects and interstellar gas may be moving relative to the old stars at a velocity of a few kilometer per second in both the radial (galacto-centric), and circular directions, depending on the spiral model adopted. Such a velocity is usually referred as the systematic motion or the streaming motion. The conventionally adopted local standard of rest is indeed co-moving with the young objects of the solar vicinity. Therefore, it has a net systematic motion with respect to the circular motion of an equilibrium galactic model, defined by the old stars.

Yuan, C.↗

On the local standard of rest

Under the influence of a spiral gravitational field, there should be differences among the mean motions of different types of objects with different dispersion velocities in a spiral galaxy. The old stars with high dispersion velocity should have essentially no mean motion normal to the galactic rotation. On the other hand, young objects and interstellar gas may be moving relative to the old stars at a velocity of a few kilometer per second in both the radial (galacto-centric), and circular directions, depending on the spiral model adopted. Such a velocity is usually referred as the systematic motion or the streaming motion. The conventionally adopted local standard of rest is indeed co-moving with the young objects of the solar vicinity. Therefore, it has a net systematic motion with respect to the circular motion of an equilibrium galactic model, defined by the old stars. Previously announced in STAR as N83-24443

Yuan, C.↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Effects of spoilers and gear on B-747 wake vortex velocities

Vortex velocities were measured in the wakes of four configurations of a 0.61-m span model of a B-747 aircraft. The wakes were generated by towing the model underwater in a ship model basin. Tangential and axial velocity profiles were obtained with a scanning laser velocimeter as the wakes aged to 35 span lengths behind the model. A 45 deg deflection of two outboard flight spoilers with the model in the landing configuration resulted in a 36 percent reduction in wake maximum tangential velocity, altered velocity profiles, and erratic vortex trajectories. Deployment of the landing gear with the inboard flaps in the landing position and outboard flaps retracted had little effect on the flap vortices to 35 spans, but caused the wing tip vortices to have: (1) more diffuse velocity profiles; (2) a 27 percent reduction in maximum tangential velocity; and (3) a more rapid merger with the flap vortices.

Luebs, A. B.↗

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗

Inverse problems in the modeling of vibrations of flexible beams

The formulation and solution of inverse problems for the estimation of parameters which describe damping and other dynamic properties in distributed models for the vibration of flexible structures is considered. Motivated by a slewing beam experiment, the identification of a nonlinear velocity dependent term which models air drag damping in the Euler-Bernoulli equation is investigated. Galerkin techniques are used to generate finite dimensional approximations. Convergence estimates and numerical results are given. The modeling of, and related inverse problems for the dynamics of a high pressure hose line feeding a gas thruster actuator at the tip of a cantilevered beam are then considered. Approximation and convergence are discussed and numerical results involving experimental data are presented.

Banks, H. T.↗