Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity steps”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Magnetic flux rope versus the spheromak as models for interplanetary magnetic clouds

Magnetic clouds form a subset of interplanetary ejecta with well-defined magnetic and thermodynamic properties. Observationally, it is well established that magnetic clouds expand as they propagate antisunward. The aim of this paper is to compare and contrast two models which have been proposed for the global magnetic field line topology of magnetic clouds: a magnetic flux tube geometry, on the one hand, and a spheromak geometry (including possible higher multiples), on the other. Traditionally, the magnetic structure of magnetic clouds has been modeled by force-free configurations. In a first step, we therefore analyze the ability of static force-free models to account for the asymmetries observed in the magnetic field profiles of magnetic clouds. For a cylindrical flux tube the magnetic field remains symmetric about closest approach to the magnetic axis on all spacecraft orbits intersecting it, whereas in a spheromak geometry one can have asymmetries in the magnetic field signatures along some spacecraft trajectories. The duration of typical magnetic cloud encounters at 1 AU (1 to 2 days) is comparable to their travel time from the Sun to 1 AU and thus magnetic clouds should be treated as strongly nonstationary objects. In a second step, therefore, we abandon the static approach and model magnetic clouds as self-similarly evolving MHD configurations. In our theory, the interaction of the expanding magnetic cloud with the ambient plasma is taken into account by a drag force proportional to the density and the velocity of expansion. Solving rigorously the full set of MHD equations, we demonstrate that the asymmetry in the magnetic signature may arise solely as a result of expansion. Using asymptotic solutions of the MHD equations, we least squares fit both theoretical models to interplanetary data. We find that while the central part of the magnetic cloud is adequately described by both models, the 'edges' of the cloud data are modeled better by the magnetic flux tube. Further comparisons of the two models necessarily involve thermodynamic properties, since real magnetic configurations are never exactly force-free and gas pressure plays an essential role. We consider a polytropic gas. Our theoretical analysis shows that the self-similar expansion of a magnetic flux tube requires the polytropic index gamma to be less than unity. For the spheromak, however, self-similar, radially expanding solutions are known only for gamma equal to 4/3. This difference, therefore, yields a good way of distinguishing between the two geometries. It has been shown recently that the polytropic relationship is applicable to magnetic clouds and that the corresponding polytropic index is approximately 0.5. This observational result is consistent with the self-similar model of the magnetic flux rope but is in conflict with the self-similar spheromak model.

Farrugia, C. J.↗
DOE OSTI2023

Adjoint DSMC for nonlinear spatially-homogeneous Boltzmann equation with a general collision model

We derive an adjoint method for the Direct Simulation Monte Carlo (DSMC) method for the spatially homogeneous Boltzmann equation with a general collision law. This generalizes our previous results in Caflisch et al., which was restricted to the case of Maxwell molecules, for which the collision rate is constant. The main difficulty in generalizing the previous results is that a rejection sampling step is required in the DSMC algorithm in order to handle the variable collision rate. We find a new term corresponding to the so-called score function in the adjoint equation and a new adjoint Jacobian matrix capturing the dependence of the collision parameter on the velocities. The new formula works for a much more general class of collision models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Sensitivity of optimized control systems to minor structural modifications

A procedure for checking whether small changes in a structure have the potential for significant enhancements of its optimized vibration control system is described. The first step in the procedure consists of the calculation of the sensitivity of the parameters of the optimized control system to small changes in the structural parameters. Theh second step consists of the optimization of the structural parameters to produce maximal increase in the performance of the control system with minimal change in the structure. The procedure has been demonstrated for a flexible laboratory structure controlled by several rate-feedback colocated force-actuator velocity-sensor pairs. Significant improvements in the performance of the control system were obtained with small structural modifications. Analytical predictions of such effects have also been validated experimentally.

Haftka, R. T.↗
NASA NTRS1991

Software manual for operating particle displacement tracking data acquisition and reduction system

The software manual is presented. The necessary steps required to record, analyze, and reduce Particle Image Velocimetry (PIV) data using the Particle Displacement Tracking (PDT) technique are described. The new PDT system is an all electronic technique employing a CCD video camera and a large memory buffer frame-grabber board to record low velocity (less than or equal to 20 cm/s) flows. Using a simple encoding scheme, a time sequence of single exposure images are time coded into a single image and then processed to track particle displacements and determine 2-D velocity vectors. All the PDT data acquisition, analysis, and data reduction software is written to run on an 80386 PC.

Wernet, Mark P.↗
NASA NTRS1994

Formation of iron metal and grain coagulation in the solar nebula

The interstellar grain population in the giant molecular cloud from which the sun formed contained little or no iron metal. However, thermal processing of individual interstellar silicates in the solar nebula is likely to result in the formation of a population of very small iron metal grains. If such grains are exposed to even transient magnetic fields, each will become a tiny dipole magnet capable of interacting with other such dipoles over spatial scale orders of magnitude larger than the radii of individual grains. Such interactions will greatly increase the coagulation cross-section for this grain population. Furthermore, the magnetic attraction between two iron dipoles will significantly increase both the collisional sticking coefficient and the strength of the interparticle binding energy for iron aggregates. Formation of iron metal may therefore be a key step in the aggregation of planetesimals in a protoplanetary nebula. Such aggregates may have already been observed in protoplanetary systems. The enhancement in the effective interaction distance between two magnetic dipoles is directly proportional to the strength of the magnetic dipoles and inversely proportional to the relative velocity. It is less sensitive to the reduced mass of the interacting particles (alpha M(exp -1/2)) and almost insensitive to the initial number density of magnetic dipoles (alpha n(sub o)(exp 1/6)). We are in the process of measuring the degree of coagulation in our condensation flow apparatus as a function of applied magnetic field and correlating these results by means of magnetic remanance acquisition measurements on our iron grains with the strength of the magnetic field to which the grains are exposed. Results of our magnetic remanance acquisition measurements and the magnetic-induced coagulation study will be presented as well as an estimate of the importance of such processes near the nebular midplane.

Nuth, Joseph A., III↗
DOE OSTI2024

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
NASA NTRS1980

Interaction of oscillations of channel flow and flow separation at duct discontinuities

A two dimensional experiment was used to study the flow fluctuations in the separated region resulting from a wall slot step (convergent step). This geometry corresponds to the geometry present in a segmented solid rocket motor in the region of the ends of adjoining segments. A modulated tunnel flow of controlled oscillation frequency was used, and the effect on flow fluctuations in the separated region was examined by hot wire velocity measurements. The transfer function between free stream oscillations and separated flow fluctuations at the driven frequency was determined over a range of frequencies and main stream flow velocities. The transfer function exhibited a frequency dependence with maximum in the 10 - 25 Hz range, with the maximum occurring at higher frequency when the mean flow velocity was higher.

Price, E. W.↗
NASA NTRS1985

Application of radiative boundary conditions to nonhydrostatic primitive equation models

A procedure is presented for applying a radiative boundary condition to the lateral boundary of a limited domain model. The model considered, the Terminal Area Simulation System (TASS), has nonhydrostatic compressible and unsteady governing equations. The primary applications to date have been studies of downbursts, wind shear, thunderstorms and mass fires. A time-splitting integration is used with TASS, with the acoustic terms being integrated separately with a smaller time step. A first-order closure approximation is employed to ensure subgrid turbulence closure. TASS applications for two-dimensional axisymmetric and three-dimensional simulations are described. The radiative boundary conditions were applied in the two-dimensional case to the temperature, pressure and non-normal velocity components. Solutions in both cases avoided runaway circulations and experienced minimal domain-wide mass loss. Also, distorsion was minimized at the boundaries.

Proctor, F. H.↗
NASA NTRS1994

A New Reynolds Stress Algebraic Equation Model

A general turbulent constitutive relation is directly applied to propose a new Reynolds stress algebraic equation model. In the development of this model, the constraints based on rapid distortion theory and realizability (i.e. the positivity of the normal Reynolds stresses and the Schwarz' inequality between turbulent velocity correlations) are imposed. Model coefficients are calibrated using well-studied basic flows such as homogeneous shear flow and the surface flow in the inertial sublayer. The performance of this model is then tested in complex turbulent flows including the separated flow over a backward-facing step and the flow in a confined jet. The calculation results are encouraging and point to the success of the present model in modeling turbulent flows with complex geometries.

Shih, Tsan-Hsing↗
NASA NTRS1995

Exhausted Plume Flow Field Prediction Near the Afterbody of Hypersonic Flight Vehicles in High Altitudes

A two-dimensional computer code to solve the Burnett equations has been developed which computes the flow interaction between an exhausted plume and hypersonic external flow near the afterbody of a flight vehicle. This Burnett-2D code extends the capability of Navier-Stokes solver (RPLUS2D code) to include high-order Burnett source terms and slip-wall conditions for velocity and temperature. Higher-order Burnett viscous stress and heat flux terms are discretized using central-differencing and treated as source terms. Blocking logic is adopted in order to overcome the difficulty of grid generation. The computation of exhaust plume flow field is divided into two steps. In the first step, the thruster nozzle exit conditions are computed which generates inflow conditions in the base area near the afterbody. Results demonstrated that at high altitudes, the computations of nozzle exit conditions must include the effects of base flow since significant expansion exists in the base region. In the second step, Burnett equations were solved for exhaust plume flow field near the afterbody. The free stream conditions are set at an altitude equal to 80km and the Mach number is equal to 5.0. The preliminary results show that the plume expansion, as altitude increases, will eventually cause upstream flow separation.

Chou, Lynn Chen↗
DOE OSTI2022

A correction procedure for self-induced velocity of a finite-sized particle in two-way coupled Euler–Lagrange simulations

The importance of incorporating a correction to undo the self-induced perturbation velocity of a particle, when its size becomes comparable to the Eulerian grid, in a two-way coupled Euler–Lagrange (EL) simulation is now well appreciated. The present work improves upon the prior correction procedures in a few important ways. First, the past correction procedures have been scalar-based with the assumption that the quasi-steady force is the source of self-induced velocity perturbation. Here we generalize to a vector correction procedure and thereby the directions of feedback force and relative velocity can be different. This allows the correction procedure to be used even in the presence of added-mass, history, and lift forces. Second, the effect of a nearby wall has been systematically included in the correction procedure. The correction procedure depends on fundamental Oseen solutions of streamwise and transverse regularized feedback forces. We present a Fourier transform-based analytical approach to obtaining these regularized Oseen solutions. We also present a step-by-step numerical procedure for obtaining the Oseen solutions in any EL code. With the analytical or numerical Oseen functions, the correction procedure can be easily implemented in any EL code. Iterations are required in solving the implicit correction equations and it is demonstrated that the correction procedure converges rapidly within three or four iterations. In conclusion, a simple empirical approach is also presented to account for unsteady effects in the correction procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

SUGAR: An improved empirical model of Type Ia supernovae based on spectral features

Type Ia supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature. This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements. This model was constructed from SNe Ia spectral properties and spectrophotometric data from the Nearby Supernova Factory collaboration. In a first step, a principal component analysis-like method was used on spectral features measured at maximum light, which allowed us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties were used to extract the average extinction curve. Third, an interpolation using Gaussian processes facilitated using data taken at different epochs during the lifetime of an SN Ia and then projecting the data on a fixed time grid. Finally, the three steps were combined to build the SED model as a function of time and wavelength. This is the SUGAR model. Results. The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity and the second correlates with their calcium lines. The addition of these parameters, as well as the high quality of the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade. The performance of this model makes it an excellent SED model for experiments like the Zwicky Transient Facility, the Large Synoptic Survey Telescope, or the Wide Field Infrared Survey Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Influence of Cs Loading on Pt/m-ZrO 2 Water–Gas Shift Catalysts

Certain alkali metals (Na, K) at targeted loadings have been shown in recent decades to significantly promote the LT-WGS reaction. This occurs at alkali doping levels where a redshift in the C-H band of formate occurs, indicating electronic weakening of the bond. The C-H bond breaking of formate is the proposed rate-limiting step of the formate associative mechanism, lending support to the occurrence of this mechanism in H2-rich environments of the LT-WGS stage of fuel processors. Continuing in this vein of research, 2%Pt/m-ZrO 2 was promoted with various levels of Cs in order to explore its influence on the rate of formate intermediate decomposition, as well as that of LT-WGS in a fixed bed reactor. In situ DRIFTS experiments revealed that Cs promoter loadings of 3.87% to 7.22% resulted in significant acceleration of the forward formate decomposition in steam at 130 °C. Of all of the alkali metals tested to date, the redshift in the formate ν(CH) band with the incorporation of Cs was the greatest. XANES difference experiments at the Pt L 2 and L 3 edges indicated that the electronic effect was not likely due to an enrichment of electronic density on Pt. CO 2 TPD experiments revealed that, unlike Na and K promoters, Cs behaves more like Rb in that the decomposition of the second intermediate in LT-WGS, carbonate species, is hindered due to (1) increased basicity of Cs, (2) the tendency of Cs to cover Pt sites that facilitate CO 2 decomposition, and (3) the tendency of Cs to increase Pt particle size as shown by EXAFS results, resulting in fewer Pt sites that facilitate CO2 decomposition. As such, the LT-WGS rate was hindered overall and the rate-limiting step shifted to carbonate decomposition (CO 2 removal). Like its Rb counterpart, low levels of added Cs (e.g., 0.72%Cs) were found to improve the stability of the catalyst relative to the unpromoted catalyst; the stability comparison was made at similar CO conversion level as well as similar space velocity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Diffusion flame extinction in slow convenctive flow under microgravity environment

A theoretical analysis is presented to study the extinction characteristics of a diffusion flame near the leading edge of a thin fuel plate in slow, forced convective flows in a microgravity environment. The mathematical model includes two-dimensional Navier-Stokes momentum, energy and species equations with one-step overall chemical reaction using second-order finite rate Arrhenius kinetics. Radiant heat loss on the fuel plate is applied in the model as it is the dominant mechanism for flame extinguishment in the small convective flow regime. A parametric study based on the variation of convective flow velocity, which varies the Damkohler number (Da), and the surface radiant heat loss parameter (S) simultaneously, is given. An extinction limit is found in the regime of slow convective flow when the rate of radiant heat loss from fuel surface outweighs the rate of heat generation due to combustion. The transition from existent envelope flame to extinguishment consists of gradual flame contraction in the opposed flow direction together with flame temperature reduction as the convective flow velocity decreases continuously until the extinction limit is reached. A case of flame structure subjected to surface radiant heat loss is also presented and discussed.

Chen, C. H.↗
NASA NTRS1986

Diffusion flame extinction in slow convection flow under microgravity environment

A theoretical analysis is presented to study the extinction characteristics of a diffusion flame near the leading edge of a thin fuel plate in slow, forced convective flows in a microgravity environment. The mathematical model includes two-dimensional Navier-Stokes momentum, energy and species equations with one-step overall chemical reaction using second-order finite rate Arrhenius kinetics. Radiant heat loss on the fuel plate is applied in the model as it is the dominant mechanism for flame extinguishment in the small convective flow regime. A parametric study based on the variation of convective flow velocity, which varies the Damkchler number (Da), and the surface radiant heat loss parameter (S) simultaneously, is given. An extinction limit is found in the regime of slow convective flow when the rate of radiant heat loss from fuel surface outweighs the rate of heat generation due to combustion. The transition from existent envelope flame to extinguishment consists of gradual flame contraction in the opposed flow direction together with flame temperature reduction as the convective flow velocity decreases continuously until the extinction limit is reached. A case of flame structure subjected to surface radiant heat loss is also presented and discussed.

Chen, Chiun-Hsun↗
NASA NTRS1991

Three-component laser velocimeter surveys of the flow over a backward-facing step

A three-component laser velocimeter is used to investigate the flow over a backward-facing step. The backward-facing step had an expansion ratio of 2, a boundary layer height to step height ratio of 0.34 and a Reynolds number based on step height of 19,000. Results from three-component velocimeter surveys of the flow over the backward-facing step are presented with comparisons of the current experiment with previous experiments and computational results. The present results compared well with previous experiments with the exception of the reattachment length. The short reattachment length was due to the short length of the channel downstream. The measurement of the lateral velocity component showed that there is a mean flow in and out of the centerline plane as high as 7 percent of the freestream velocity. However, the shear stresses show no correlation between the lateral fluctuations and the longitudinal and vertical fluctuations, indicating that the flow is 2D in terms of the turbulence quantities.

Kjelgaard, Scott O.↗
DOE OSTI2023

Plasma surface ionization wave interactions with single channels

Abstract The study of plasma surface ionization waves (SIWs) in recent years has primarily focused on planar surfaces and periodic two dimensional structures. In application, substrates are likely to have non-planar morphology such as cracks, pores, and steps. Additionally, targets for the applications of medicine or catalysis may have targets with heterogeneous composition. This classification of targets are brought under the umbrella of complex interfaces. In this work, plasma SIWs were incident on a complex target consisting of a single channel cut into glass slides. The SIW velocities for the in-channel portion of the wave and radially propagating portion of the wave were tracked. It was found that surface wave velocities are not significantly affected by channel geometry, but primarily increase with pulse energy. A third propagation direction for the SIW is characterized in the azimuthal direction relative to the radial portion of the wave. Channel geometry is found to greatly effect the area treated by the plasma but not the propagation velocity of the surface wave. Surface wave morphology and the impact on application is also discussed. A simple model was introduced to understand the mechanisms behind SIW escape from a channel. It was found that the ratio of pulse energy to a geometry dependent minimum energy can predict the escape angle of a SIW from single channels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗