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 685 records · Page 38

A new determination of the halo luminosity density of the galaxy

The distribution of tangential velocities for arbitrary halo velocity dispersions and rotation rates are explicitly computed, assuming a particular form of the halo velocity distribution function near the sun. It is found that for a specific model of the halo velocity dispersion tensor that Schmidt's (1975) estimate of the mass carried by halo stars of 0.25-0.70 solar masses should be revised by a factor of between 1.0 and 2.0, depending on the rotation rate of the halo. Although obtaining a mass density in all stars from the stellar density in that mass range is extremely uncertain, it is found, in agreement with earlier work, that even under rather extreme assumptions about the remnants the observed stars and expected remnants fail to provide the theoretically expected mass. Comparison of the luminosity to the theoretical mass yields a mass to light ratio consistent with that of rich clusters of galaxies.

Richstone, D. O.↗

Multi-Body Analysis of the 1/5 Scale Wind Tunnel Model of the V-22 Tiltrotor

The paper presents a multi-body analysis of the 1/5 scale wind tunnel model of the V-22 tiltrotor, the Wing and Rotor Aeroelastic Testing System (WRATS), currently tested at NASA Langley Research Center. An original multi-body formulation has been developed at the Dipartimento di Ingegneria Aerospaziale of the Politecnico di Milano, Italy. It is based on the direct writing of the equilibrium equations of independent rigid bodies, connected by kinematic constraints that result in the addition of algebraic constraint equations, and by dynamic constraints, that directly contribute to the equilibrium equations. The formulation has been extended to the simultaneous solution of interdisciplinary problems by modeling electric and hydraulic networks, for aeroservoelastic problems. The code has been tailored to the modeling of rotorcrafts while preserving a complete generality. A family of aerodynamic elements has been introduced to model high aspect aerodynamic surfaces, based on the strip theory, with quasi-steady aerodynamic coefficients, compressibility, post-stall interpolation of experimental data, dynamic stall modeling, and radial flow drag. Different models for the induced velocity of the rotor can be used, from uniform velocity to dynamic in flow. A complete dynamic and aeroelastic analysis of the model of the V-22 tiltrotor has been performed, to assess the validity of the formulation and to exploit the unique features of multi-body analysis with respect to conventional comprehensive rotorcraft codes; These are the ability to model the exact kinematics of mechanical systems, and the possibility to simulate unusual maneuvers and unusual flight conditions, that are particular to the tiltrotor, e.g. the conversion maneuver. A complete modal validation of the analytical model has been performed, to assess the ability to reproduce the correct dynamics of the system with a relatively coarse beam model of the semispan wing, pylon and rotor. Particular care has been used to model the kinematics of the gimbal joint, that characterizes the rotor hub, and of the control system, consisting in the entire swashplate mechanism. The kinematics of the fixed and the rotating plates have been modeled, with variable length control links used to input the controls, the rotating flexible links, the pitch horns and the pitch bearings. The investigations took advantage of concurring wind tunnel test runs, that were performed in August 1998, and allowed the acquisition of data specific to the multi-body analysis.

Ghiringhelli, G. L.↗

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗

Numerical integrators for Lagrangian oceanography

Abstract. A common task in Lagrangian oceanography is to calculate a large number of drifter trajectories from a velocity field precalculated with an ocean model. Mathematically, this is simply numerical integration of an ordinary differential equation (ODE), for which a wide range of different methods exist. However, the discrete nature of the modelled ocean currents requires interpolation of the velocity field in both space and time, and the choice of interpolation scheme has implications for the accuracy and efficiency of the different numerical ODE methods. We investigate trajectory calculation in modelled ocean currents with 800 m, 4 km, and 20 km horizontal resolution, in combination with linear, cubic and quintic spline interpolation. We use fixed-step Runge–Kutta integrators of orders 1–4, as well as three variable-step Runge–Kutta methods (Bogacki–Shampine 3(2), Dormand–Prince 5(4) and 8(7)). Additionally, we design and test modified special-purpose variants of the three variable-step integrators, which are better able to handle discontinuous derivatives in an interpolated velocity field. Our results show that the optimal choice of ODE integrator depends on the resolution of the ocean model, the degree of interpolation, and the desired accuracy. For cubic interpolation, the commonly used Dormand–Prince 5(4) is rarely the most efficient choice. We find that in many cases, our special-purpose integrators can improve accuracy by many orders of magnitude over their standard counterparts, with no increase in computational effort. Equivalently, the special-purpose integrators can provide the same accuracy as standard methods at a reduced computational cost. The best results are seen for coarser resolutions (4 and 20 km), thus the special-purpose integrators are particularly advantageous for research using regional to global ocean models to compute large numbers of trajectories. Our results are also applicable to trajectory computations on data from atmospheric models.

58 GEOSCIENCES↗

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗

Ion heating mechanism in a modified Penning discharge

Ions with Maxwellian energy distributions and kinetic temperatures ranging from 20 eV to 7 keV have been observed in a modified Penning discharge operating in the steady state. Investigation of the plasma revealed two distinct spoke-like concentrations of charge rotating with different velocities in the sheath between the plasma and the anode ring. The faster spoke consists of electrons rotating with the E/B drift velocity, where E is the electric field and B is the magnetic field strength. The slow spoke consists of ions, the thermal velocity of which is observed to be proportional to the spoke velocity. The experimental data are consistent with a model whereby the ion drift velocity in this spoke, corresponding to kilovolt ion energies, is Maxwellianized by strong electrostatic turbulence in the sheath. Theoretical expressions are derived for the frequency of the electron and ion spoke rotation, for the ion kinetic temperature, and for the ion heating efficiency as functions of the discharge parameters. These expressions are shown to be consistent with extensive experimental data.

Roth, J. R.↗

Observations of downward moving O(+) in the polar topside ionosphere

New investigations of the core ion motion within high-latitude topside ionosphere and near-Earth magnetosphere, using data from the Dynamics Explorer (DE) retarding ion mass spectrometer (RIMS), reveal the existence of significant regions of downward moving O(+). The occurences of downgoing versus upgoing O(+) are not clearly separable in terms of either polar zenith angle or Kp but are well distinquished by the direction of the z component of the interplanetary magnetic field (IMF). On the average, down flow dominates when IMFB2 less than O, while upflow dominates for IMFB2 greater than 0. Combining cross-field convection velocities derived from an empirical convection electric field model with the observed parallel velocities yields a two dimensional ion velocity field. This velocity field is consistent with a senario which has O(+) of cusp/cleft and auroral zone origin concvecting into the polar cap and, because of the dominance of gravitational energy over the upward kinetic energy, falling back into the inonsphere. This provides additional confirmation of the results of studies of the cleft ion fountain. Estimates of the flux of O(+) in the upflowing and downflowing regions for Lambda greater than or equal to 60 deg give a total upflow of approximately 6 X 10(exp 25) ions/sec for IMFB2 greater than 0 and total upflow and downflow of approximately 4 X 10 (exp 25) ions/sec and 1 x 10(exp 25) ions/sec, respectively, for IMFB less than 0. In all cases the magnitude of the dayside outflow is consistent with previous work on upwelling ions. While the magnitudes vary for high and low Kp the ratios of upward to downward flow are roughly the same at approximately 1.7. The downflowing O(+) shows a correlation with the magnitude of the outflow of light ions in the same region but the cause and effect of this relationship is not distinquishable. hable.

Chandler, M. O.↗

Velocity Dependence of Rate‐And‐State Friction in Granular Fault Gouge and Implications for Slow‐Slip Events

Abstract The standard rate‐and‐state friction (RSF) has extensively captured frictional behaviors, but it fails to explain the velocity dependence of frictional stability transition and widespread slow‐slip events (SSEs) in experiments and nature adequately. An alternative microphysical Chen‐Niemeijer‐Spiers (CNS) model can well describe the velocity dependence of frictional behaviors of granular gouges. Using the original CNS model, standard RSF parameters can be quantified microphysically. However, some micro‐parameters are not easy to estimate quantitatively, making it difficult to extrapolate to natural and experimental conditions. Here, we simplify the microphysically‐derived RSF parameters including direct effect a , evolution effect b , and critical slip distance D c , as well as equivalent values ( a eq , b eq , and D eq ). The simplified friction parameters directly illustrate their velocity dependence, namely the essentially constant a , a eq , and D c , negatively velocity‐dependent b and b eq , as well as varying D eq for different laws. They are roughly consistent with experimental results in various fault gouges. A modified CNS model is further derived from the original CNS model, establishing a direct link between the standard RSF and CNS models. The modified CNS model exhibits virtually identical frictional behaviors to the original CNS, but differs from the standard RSF at large velocity perturbations. Moreover, the linearized stability analysis indicates that the critical stiffness for the modified CNS model is velocity‐dependent. Compared with the standard RSF, the modified CNS model not only explains the velocity dependence of frictional stability transition, but also exhibits a more gradual transition for SSEs with a broader range of stiffness ratios.

Geochemistry & Geophysics↗

Physical and mathematical modelling of ladle metallurgy operations

Experimental measurements are reported, on the velocity fields and turbulence parameters on a water model of an argon stirred ladle. These velocity measurements are complemented by direct heat transfer measurements, obtained by studying the rate at which ice rods immersed into the system melt, at various locations. The theoretical work undertaken involved the use of the turbulence Navier-Stokes equations in conjunction with the kappa-epsilon model to predict the local velocity fields and the maps of the turbulence parameters. Theoretical predictions were in reasonably good agreement with the experimentally measured velocity fields; the agreement between the predicted and the measured turbulence parameters was less perfect, but still satisfactory. The implications of these findings to the modelling of ladle metallurgical operations are discussed.

El-Kaddah, N.↗

A numerical study of thunderstorm electrification - Model development and case study

A numerical model was developed for examining the thunderstorm electrification process in which it is assumed that the electrification is entirely due to noninductive charge transfer between colliding ice crystals and hail. Since this ice-hail charge mechanism is very dependent on particle sizes and distributions, an explicit microphysical framework is used. To maintain simplicity, the electrification model is kinematic; thus the temperature and velocity fields are input into the electrification model. The cloud model of Taylor (1989) was used to generate the temperature and velocity fields to examine the July 19, 1981, Cooperative Convective Precipitation Experiment thundercloud. Using these fields, the electrification model produced time-dependent ice particle concentrations, radar reflectivities, charge, and vertical electric field distributions in good general agreement with those observed. The model produced a maximum electric field strength of 1.27 kV/cm, which is on the order of that needed for lightning initiation, and this maximum occurred very close to the time of the observed discharge. Thus, the ice-hail charge mechanism appears to have played an important role in the electrical development of the July 19 cloud.

Norville, K.↗

Recovering Photospheric Velocities from Vector Magnetograms by using a Three-Dimensional, Fully Magnetohydrodynamic (MHD) Model

We introduce a numerical simulation method for recovering the photospheric velocity field from the vector magnetograms. The traditional method is local correlation tracking (LCT) which is based on measuring the relative displacements of features in blocks of pixels between successive white-light images or magnetograms. Within this method, there are a variety of implementations. One of recently developed implementations is induction local correlation tracking (ILCT) as described by Welsch et al. (2004). They employ the normal component of magnetic induction equation as a constraint to assure consistent solutions. Our numerical method uses the fully three-dimensional MHD equations to recover the photospheric velocity field with individual vector magnetograms. We compare our method to the ILCT method using NOAA AR8210 as an example. The differences and similarities are discussed in detail.

Wang, A. H.↗

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗

Stellar model chromospheres. XIII - M dwarf stars

Single-component, homogeneous model chromospheres that are consistent with high-resolution profiles of the Ca II K line calibrated in surface flux units for three dMe and 2 dM stars observed at quiescent times are constructed. The models reveal several systematic trends. Large values of the ratio of T(min) to T(eff) are derived, indicating a large amount of nonradiative heating present in the upper photospheres of M dwarf stars. It is also found that the lower chromospheric temperature gradient is similar for all the M dwarf stars. Since for the models here the chromospheric K line emission strength is most sensitive to the total amount of chromospheric material present within the approximate temperature range T(min)-6000 K, increasing the emission strength is not simply due to increasing chromospheric temperature gradients. It is also found that both the electron density and electron temperature at one thermalization length in the K line below the top of the chromospheres are greater in the dMe stars than in the dM stars. The M dwarf models here have microturbulent velocities between 1 and 2 km/sec, which are much smaller than for solar chromosphere models.

Giampapa, M. S.↗

Effects of AP particle size on combustion response to crossflow

An analytical model is developed for the linearized velocity-coupled combustion response function. The model treats elements of response to perturbations in pressure, in composition due to the heterogeneity of composite propellants, and in crossflow velocity. The effects of AP particle size are accounted for in terms of effects on controlling ballistics properties and in terms of fluctuations in propellant composition. There are two facets of the crossflow problem: the effect of crossflow velocity on the various response elements, and the response to velocity perturbations. Both are dealt with in this paper. Important trends derived from series of parametric computations are described.

Cohen, N. S.↗

Simulation of Tropical Biomass Burning

The research objectives that were achieved during the course of our studies include the following: (1) Over the last few years, a model has been developed in the Atmospheric Chemistry and Dynamics Branch at Ames Research Center in collaboration with the Physics Department at San Jose State University. It is referred to as the Global/Regional Atmospheric Chemistry Event Simulator (GRACES). Currently, the GRACES model system combines an atmospheric chemistry and transport model, and a regional mesoscale meteorological model. Therefore this system is suitable for simulating the conditions observed by the tropical observation missions, such as the Pacific Exploratory Mission in the 'Tropics (PEM-Tropics), Study of Ozone and Nitrogen oxides Experiment (SONEX), and other periods. Specifically, the research carried out included the evaluation of the behavior of several components of the MM5 (I.e., Meteorological Model 5, version 2) and the GRACES combined modeling system. We initiated research on (a) the ability of the MM5 model to assimilate downward vertical velocities at least as high as the analyses, (b) the ability of the Graces model to incorporate the vertical velocities from MM5, and (c) other factors related to transport patterns required to transport CO in the observed manner. We carried out improved calculations of the transport of tracers for both NASA airborne missions, SONEX and PEM-Tropics. We also made improved source strength estimates fopr isoprene dust, and similar emissions from the Earths surface. This required the use of newly available databases on the Earth's surface and vegetation. We completed atmospheric chemistry simulations of radicals and nitrogen oxide species. We have greatly improved the handling of cumulonimnbus convection by modifing an existing scheme.

Hamill, Patrick↗

Jet Noise Reduction Potential from Emerging Variable Cycle Technologies

Acoustic and flow-field experiments were conducted on exhaust concepts for the next generation supersonic, commercial aircraft. The concepts were developed by Lockheed Martin (LM), Rolls-Royce Liberty Works (RRLW), and General Electric Global Research (GEGR) as part of an N+2 (next generation forward) aircraft system study initiated by the Supersonics Project in NASA s Fundamental Aeronautics Program. The experiments were conducted in the Aero-Acoustic Propulsion Laboratory at the NASA Glenn Research Center. The exhaust concepts utilized ejectors, inverted velocity profiles, and fluidic shields. One of the ejector concepts was found to produce stagnant flow within the ejector and the other ejector concept produced discrete-frequency tones that degraded the acoustic performance of the model. The concept incorporating an inverted velocity profile and fluid shield produced overall-sound-pressure-level reductions of 6 dB relative to a single stream nozzle at the peak jet noise angle for some nozzle pressure ratios. Flow separations in the nozzle degraded the acoustic performance of the inverted velocity profile model at low nozzle pressure ratios.

Henderson, Brenda↗

Jet Noise Reduction Potential From Emerging Variable Cycle Technologies

Acoustic and flow-field experiments were conducted on exhaust concepts for the next generation supersonic, commercial aircraft. The concepts were developed by Lockheed Martin (LM), Rolls-Royce Liberty Works (RRLW), and General Electric Global Research (GEGR) as part of an N+2 (next generation forward) aircraft system study initiated by the Supersonics Project in NASA s Fundamental Aeronautics Program. The experiments were conducted in the Aero-Acoustic Propulsion Laboratory at the NASA Glenn Research Center. The exhaust concepts utilized ejectors, inverted velocity profiles, and fluidic shields. One of the ejector concepts was found to produce stagnant flow within the ejector and the other ejector concept produced discrete-frequency tones that degraded the acoustic performance of the model. The concept incorporating an inverted velocity profile and fluid shield produced overall-sound-pressure-level reductions of 6 dB relative to a single stream nozzle at the peak jet noise angle for some nozzle pressure ratios. Flow separations in the nozzle degraded the acoustic performance of the inverted velocity profile model at low nozzle pressure ratios.

Source record↗

Coefficient of restitution for low velocity transverse impact of thin graphite-epoxy laminates

A simple model is presented for the low velocity impact behavior of thin plates composed of T300/5208 graphite/epoxy. The model, which is constructed on the assumption that the plate is inextensible in the fiber direction and that the material is incompressible in the z direction, is able to predict the energy absorbed by the target during an impact event, or the 'coefficient of restitution'. The plate essentially deforms by shear, allowing the model to neglect bending deformations. The coefficient of restitution is predicted to increase with large interlaminar shear strength and low transverse shear modulus of the laminate. The coefficient's experimentally measured values agree with model-predicted values within a reasonable error.

Lal, K. M.↗