Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity model evaluation”

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 487 records · Page 27

Evaluating the carbon footprint of the integrated DBD‐plasma bi‐reforming unit via laboratory scale experiments and scaled‐up process modeling

Catalytic dielectric barrier discharge (DBD) plasma reactor experiments were performed in a tubular glass reactor with a 2 mm gap at 550°C to facilitate the reaction kinetics of steam added dry reforming or bireforming. The best specific energy input obtained was 11.2 eV/molecule feed at CO 2 :CH 4 :H 2 O of 4.5:1:4.5 ratio and gas hour space velocity (GHSV) = 432 h −1 . This value was used to design a conceptual process and assess the environmental impact of methane steam reforming-based H 2 production 18.4 kmol/h CO 2 emission processing into H 2 :CO = 2 syngas, with an emphasis on the carbon footprint.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of a locally homogeneous model of spray evaporation

A model of spray evaporation which employs a second-order turbulence model in conjunction with the locally homogeneous flow approximation, which implies infinitely fast interphase transport rates is presented. Measurements to test the model were completed for single phase constant and variable density jets, as well as an evaporating spray in stagnant air. Profiles of mean velocity, composition, temperature and drop size distribution as well as velocity fluctuations and Reynolds stress, were measured within the spray. Predictions were in agreement with measurements in single phase flows and also with many characteristics of the spray, e.g. flow width, radial profiles of mean and turbulent quantities, and the axial rate of decay of mean velocity and mixture fraction.

Shearer, A. J.↗

Use of Remote Sensing and In-Situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the ACME Model

The overachieving goal of the whole CMDV-MCS project is to improve understanding of warm season continental convection and to develop treatments of convection and microphysics capable of representing mesoscale convective systems (MCSs) features in large-scale models. Our tasks for this project contributing to the overachieving goal include: (1) Improve the ice nucleation formulation for MG2 and P3 cloud microphysics schemes; (2) Improve the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM; and (3) Test the performance of improved ice microphysics in E3SM with observation data. In this project, we have (1) Improved the ice nucleation parameterization for MG2 and P3 in E3SM by implementing two advanced empirical parameterizations with connection to aerosols. The two deterministic heterogeneous ice nucleation parameterizations (i.e., DeMott et al., 2015; Niemand et al., 2012) were merged with the MG2 and P3 cloud microphysics schemes in E3SM. Long-term simulations were conducted to examine the impacts of these new parameterizations on simulated cloud properties; (2) Improved the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM. We evaluated the double Gaussian PDF of vertical velocity simulated by the Cloud Layers Unified By Binormals (CLUBB) and the sub-column vertical velocity sampled from the Subgrid Importance Latin Hypercube Sampler (SILHS) in E3SM. We introduced the vertical velocity variance induced by topographic gravity waves for ice nucleation and droplet activation; and (3) Tested the performance of improved ice microphysics in E3SM with observation data. We tested the new treatments of ice nucleation in the single column model (SCM) mode for the stratiform mixed-phase clouds observed during 9-10 October 2004 in the DOE ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) and for the convective clouds observed on 20 May 2011 in the Midlatitude Continental Convective Clouds Experiment (MC3E). Modeled ice nucleating particles (INPs) concentrations were compared against observations collected around the globe.

54 ENVIRONMENTAL SCIENCES↗

Effects of the low Earth orbital environment on spacecraft materials

It is evident from space flights during the last three years that the low Earth orbital (LEO) environment interacts with spacecraft surfaces in significant ways. One manifestation of these interactions is recession of, in particular, organic-polymer-based surfaces presumably due to oxidation by atomic oxygen, the major component of the LEO environment. Three experiments have been conducted on Space Shuttle flights 5, 8 and 41-G to measure reaction rates and the effects of various parameters on reaction rates. Surface recession on these flights indicates reaction efficiencies approximately 3 x 10(-24) cu cm/atoms for unfilled organic polymers. Of the metals, silver and osmium are very reactive. Effects on spacecraft or experiment surfaces can be evaluated using the derived reaction efficiencies and a definition of the total exposure to atomic oxygen. This exposure is obtained using an ambient density model, solar activity data and spacecraft parameters of altitude, attitude and operational date. Oxygen flux on a given surface is obtained from the ambient density and spacecraft velocity and can then be integrated to provide the total exposure or fluence. Such information can be generated using simple computational programs and can be converted to various formats. Overall, the extent of damage is strongly dependent on the type of surface and total exposure time.

Leger, L. J.↗

One-equation near-wall turbulence modeling with the aid of direct simulation data

The length scales appearing in the relations for the eddy viscosity and dissipation rate in one-equation models were evaluated from direct numerical simulation data for developed channel and boundary-layer flow at two Reynolds numbers each. To prepare the ground for the evaluation, the distribution of the most relevant mean-flow and turbulence quantities is presented and discussed with respect to Reynolds-number influence and to differences between channel and boundary-layer flow. An alternative model is also examined in which bar-(v'(exp 2))(exp 1/2) is used as velocity scale instead of k(exp 1/2). With this velocity scale, the length scales now appearing in the model follow very closely a linear relationship near the wall so that no damping is necessary. For the determination of bar-v'(exp 2) in the context of a one-equation model, a correlation is provided between bar-(v'(exp 2))/k and bar-(u'v')/k.

Rodi, W.↗

Lyapunov Orbits in the Jupiter System Using Electrodynamic Tethers

Various researchers have proposed the use of electrodynamic tethers for power generation and capture from interplanetary transfers. The effect of tether forces on periodic orbits in Jupiter-satellite systems are investigated. A perturbation force is added to the restricted three-body problem model and a series of simplifications allows development of a conservative system that retains the Jacobi integral. Expressions are developed to find modified locations of equilibrium positions. Modified families of Lyapunov orbits are generated as functions of tether size and Jacobi integral. Zero velocity curves and stability analyses are used to evaluate the dynamical properties of tether-modified orbits.

Cassini↗

Prediction of blade stresses due to gust loading

An analysis is developed for investigating the response of a rotor-fuselage system in a three-dimensional gust field wherein the gust velocity conponents can have arbitrary variation in space and time. Each rotor blade undergoes flap bending, lag bending and torsional deflections. The blades are divided into beam elements and each element consists of fifteen nodal degrees of freedom. Quasi-steady strip theory is used to obtain the aerodynamic loads. Unsteady aerodynamic effects are introduced through dynamic inflow modeling. Dynamic stall and reverse flow effects are also included. The fuselage is allowed five degrees of freedom: vertical, longitudinal, lateral, pitch and roll motions. The gust response equations are linearized about the vehicle trim state and the blade steady-state deflected position, and then solved by time integration. The blade bending moments, which determine blade stresses, are evaluated using the force summation technique. Systematic studies are made to identify the importance of several parameters including dynamic stall, forward speed, lag stiffness, gust profile, gust penetration rate and gust velocity direction.

Bir, Gunjit↗

An analytical and experimental study of sound propagation and attenuation in variable-area ducts

The performance of sound suppression techniques in ducts that produce refraction effects due to axial velocity gradients was evaluated. A computer code based on the method of multiple scales was used to calculate the influence of axial variations due to slow changes in the cross-sectional area as well as transverse gradients due to the wall boundary layers. An attempt was made to verify the analytical model through direct comparison of experimental and computational results and the analytical determination of the influence of axial gradients on optimum liner properties. However, the analytical studies were unable to examine the influence of non-parallel ducts on the optimum linear conditions. For liner properties not close to optimum, the analytical predictions and the experimental measurements were compared. The circumferential variations of pressure amplitudes and phases at several axial positions were examined in straight and variable-area ducts, hard-wall and lined sections with and without a mean flow. Reasonable agreement between the theoretical and experimental results was obtained.

Nayfeh, A. H.↗

An analysis of temperature effect in a finite journal bearing with spatial tilt and viscous dissipation

The analysis presented herein deals with the evaluation of the pressure, velocity, and temperature profiles in a finite-length plane journal bearing. The geometry of the case under study consists of a spatially tilted shaft. The two-dimensional Reynolds equation accounts for the variation of the clearance gap h with x and z and is used to model the pressure field. The latter is solved for a variety of shaft tilt angles and then used to calculate the two-dimensional flow field. Finally, the flow field is used in the energy equation to solve for the film temperature profile, when the effect of viscous dissipation is taken into account.

Braun, M. J.↗

Liquid droplet radiator thermal characteristics

The effects of droplet size, emissivity, and spacing on the liquid droplet radiator (LDR) performance space heat rejection system are investigated. The procedures for predicting surface emissivity of semitransparent droplets as a function of absorptivity and radius are examined. The relation between absorptivity, reflectivity, transmissivity, and emissivity is studied. It is observed that for droplets in the 25-200 micron range the radius affects emissivity. The thermal performance of an array of droplet streams is described using a thermal radiation model. Droplet temperatures are calculated; it is shown that outside droplet streams are more effective than interior droplets for arrays in the range of radius to spacing ratio from 0.05-0.2. Temperature distributions as a function of droplet size, emissivity, and velocity are presented. The heat transfer performance of a LDR is evaluated in terms of fin effectiveness. The equivalent sheet emissivities of droplet streams are examined.

Edgerton, R. H.↗

Tail lobe ion composition at energies of 0.1 to 16 keV/e: Evidence for mass-dependent density gradients

A large set of energetic (0.1 to 16 keV/e) ion composition data from the central magnetotail, obtained by the ISEE 1 spacecraft between 10 and 23R(sub E) from Earth, is sorted according to measured total ion beta value, in order to investigate whether bulk properties of different ions vary in different ways between plasma sheets and lobes, as suggested by a comparison of certain experimental and theoretical results in the literature. Despite inevitable difficulties with extracting statistically valid data at very low beta (10(exp -2) or lower), the results seem to support a recent theoretical model suggesting that lighter ions have a steeper density gradient than heavier ions, especially when comparing H(+) ions with O(+) ions. The results also indicate that ion velocity distributions are fairly isotropic even at low beta, at least those of H(+) ions, although field-aligned flows are common. The results are evaluated in the context of plasma transport and are found to lend some support to the notion that tail lobe convection may be directed inward from the dawn and dusk flanks.

Lennartsson, O. W.↗

Numerical simulation of steady and unsteady viscous flow in turbomachinery using pressure based algorithm

The objective of this research is to simulate steady and unsteady viscous flows, including rotor/stator interaction and tip clearance effects in turbomachinery. The numerical formulation for steady flow developed here includes an efficient grid generation scheme, particularly suited to computational grids for the analysis of turbulent turbomachinery flows and tip clearance flows, and a semi-implicit, pressure-based computational fluid dynamics scheme that directly includes artificial dissipation, and is applicable to both viscous and inviscid flows. The values of these artificial dissipation is optimized to achieve accuracy and convergency in the solution. The numerical model is used to investigate the structure of tip clearance flows in a turbine nozzle. The structure of leakage flow is captured accurately, including blade-to-blade variation of all three velocity components, pitch and yaw angles, losses and blade static pressures in the tip clearance region. The simulation also includes evaluation of such quantities of leakage mass flow, vortex strength, losses, dominant leakage flow regions and the spanwise extent affected by the leakage flow. It is demonstrated, through optimization of grid size and artificial dissipation, that the tip clearance flow field can be captured accurately. The above numerical formulation was modified to incorporate time accurate solutions. An inner loop iteration scheme is used at each time step to account for the non-linear effects. The computation of unsteady flow through a flat plate cascade subjected to a transverse gust reveals that the choice of grid spacing and the amount of artificial dissipation is critical for accurate prediction of unsteady phenomena. The rotor-stator interaction problem is simulated by starting the computation upstream of the stator, and the upstream rotor wake is specified from the experimental data. The results show that the stator potential effects have appreciable influence on the upstream rotor wake. The predicted unsteady wake profiles are compared with the available experimental data and the agreement is good. The numerical results are interpreted to draw conclusions on the unsteady wake transport mechanism in the blade passage.

Lakshminarayana, B.↗

Tail lobe ion composition at energies of 0.1 to 16 keV/e: Evidence for mass-dependent density gradients

A large set of energetic (0.1- to 16-keV/e) ion composition data from the central magnetotail, obtained by the International Sun Earth Explorer (ISEE) 1 spacecraft between 10 and 23 R(sub E) from Earth, is sorted according to measured total ion beta value, in order to investigate whether bulk properties of different ions vary in different ways between plasma sheet and lobes, as suggested by a comparison of certain experimental and theoretical results in the literature. Despite inevitable difficulties with extracting statistically valid data at very low beta (10(exp -2) or lower), the results seem to support a recent theoretical model suggesting that lighter ions have a steeper density gradient than heavier ions, especially when comparing H(+) ions with O(+) ions. The results also indicate that ion velocity distributions are fairly isotropic even at low beta, at least those of H(+) ions, although field-aligned flows are common. The results are evaluated in the context of plasma transport and are found to lend some support to the notion that tail lobe convection may be directed inward from the dawn and dusk flanks.

Lennartsson, O. W.↗

Revision of Paschen's Law Relating to the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a version of Paschens law that takes into account the flow of ambient gas past electrode surfaces. Paschens law does not consider the flow of gas past an aerospace vehicle whose surfaces may be triboelectrically charged by dust or ice crystal impingement while traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance between electrode surfaces is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised theoretical model must be a function of the mean velocity vxm of the ambient gas and reduce to Paschens law when the mean velocity is zero. A new theoretical formulation of Paschens law, taking into account the Mach number and compressible dynamic pressure, derived by the authors, will be discussed. This equation has been evaluated by wind tunnel experimentation. Initial data of the baseline wind tunnel experiments show results consistent with the hypothesis. This work may enhance the safety of aerospace vehicles through a redefinition of electrostatic launch commit criteria. It is also possible for new products, such as antistatic coatings, to be formulated based on this data.

discharge↗

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

A comparison of washout filters using a human dynamic orientation model

The Ormsby model of human dynamic orientation, a discrete time computer program, was used to provide a vestibular explanation for observed differences between two washout schemes. These washout schemes, a linear washout and a nonlinear washout, were subjectively evaluated. It was found that the linear washout presented false rate cues, causing pilots to rate the simulation fidelity of the linear scheme much lower than the nonlinear scheme. By inputting these motion histories into the Ormsby model, it was shown that the linear filter causes discontinuities in the pilot's perceived angular velocity, resulting in the sensation of an anomalous rate cue. This phenomenon does not occur with the use of the nonlinear filter.

Riedel, S. A.↗

Does lower Omega allow a resolution of the large-scale structure problem?

The intermediate angular scale anisotropy of the cosmic microwave background, peculiar velocities, density correlations, and mass fluctuations for both neutrino and baryon-dominated universes with Omega less than one are evaluated. The large coherence length associated with a low-Omega, hot dark matter-dominated universe provides substantial density fluctuations on scales up to 100 Mpc: there is a range of acceptable models that are capable of producing large voids and superclusters of galaxies and the clustering of galaxy clusters, with Omega roughly 0.3, without violating any observational constraint. Low-Omega, cold dark matter-dominated cosmologies are also examined. All of these models may be reconciled with the inflationary requirement of a flat universe by introducing a cosmological constant 1-Omega.

Silk, Joseph↗