Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical transport”

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 325 records · Page 18

The simulated features of heliospheric cosmic-ray modulation with a time-dependent drift model. III - General energy dependence

The time-dependent cosmic-ray transport equation is solved numerically in an axially symmetric heliosphere. Gradient and curvature drifts are incorporated, together with an emulated wavy neutral sheet. This model is used to simulate heliospheric cosmic-ray modulation for the period 1985-1989 during which drifts are considered to be important. The general energy dependence of the modulation of Galactic protons is studied as predicted by the model for the energy range 1 MeV to 10 GeV. The corresponding instantaneous radial and latitudinal gradients are calculated, and it is found that, whereas the latitudinal gradients follow the trends in the waviness of the neutral sheet to a large extent for all energies, the radial gradients below about 200 MeV deviate from this general pattern. In particular, these gradients increase when the waviness decreases for the simulated period 1985-1987.3, after which they again follow the neutral sheet by increasing rapidly.

Potgieter, M. S.↗

Recent advances in MPD thruster research at Princeton

A summary of last years anode, plasma, and cathode findings is presented. A summary of this years activities and findings is also presented. A brief discussion of previous and current understanding is given and covers the following topics: existence of microinstabilities; the scaling of Va with the Hall parameter; the scaling of anomalous resistivity with the Hall parameter; the relation between anomalous resistivity and the anode drop; the presence of turbulence in the anode region; numerical simulation with anomalous transport; the use of magnets to decrease dissipation; performance testing with the new anode; the mechanisms behind the ionization sink; and lithiated cathode research.

Jahn, Robert G.↗

Compression strength of composite primary structural components

The linear elastic response is determined for an internally pressurized, long circular cylindrical shell stiffened on the inside by a regular arrangement of identical stringers and identical rings. Periodicity of this configuration permits the analysis of a portion of the shell wall centered over a generic stringer-ring joint; i.e., a unit cell model. The stiffeners are modeled as discrete beams, and the stringer is assumed to have a symmetrical cross section and the ring an asymmetrical section. Asymmetery causes out-of-plane bending and torsion of the ring. Displacements are assumed as truncated double Fourier series plus simple terms in the axial coordinate to account for the closed and pressure vessel effect (a non-periodic effect). The interacting line loads between the stiffeners and the inside shell wall are Lagrange multipliers in the formulation, and they are also assumed as truncated Fourier series. Displacement continuity constraints between the stiffeners and shell along the contact lines are satisfied point-wise. Equilibrium is imposed by the principle of virtual work. A composite material crown panel from the fuselage of a large transport aircraft is the numerical example. The distributions of the interacting line loads, and the out-of-plane bending moment and torque in the ring, are strongly dependent on modeling the deformations due to transverse shear and cross-sectional warping of the ring in torsion. This paper contains the results from the semiannual report on research on 'Pressure Pillowing of an Orthogonally Stiffened Cylindrical Shell'. The results of the new work are illustrated in the included appendix.

Johnson, Eric R.↗

CALIPSO Satellite Lidar Identification Of Elevated Dust Over Australia Compared With Air Quality Model PM60 Forecasts

Global measurements of the vertical distribution of clouds and aerosols have been recorded by the lidar on board the CALIPSO (Cloud Aerosol Lidar Infrared Pathfinder Satellite Observations) satellite since June 2006. Such extensive, height-resolved measurements provide a rare and valuable opportunity for developing, testing and validating various atmospheric models, including global climate, numerical weather prediction, chemical transport and air quality models. Here we report on the initial results of an investigation into the performance of the Australian Air Quality Forecast System (AAQFS) model in forecasting the distribution of elevated dust over the Australian region. The model forecasts of PM60 dust distribution are compared with the CALIPSO lidar Vertical Feature Mask (VFM) data product. The VFM classifies contiguous atmospheric regions of enhanced backscatter as either cloud or aerosols. Aerosols are further classified into six subtypes. By comparing forecast PM60 concentration profiles to the spatial distribution of dust reported in the CALIPSO VFM, we can assess the model s ability to predict the occurrence and the vertical and horizontal extents of dust events within the study area.

Young, Stuart A.↗

Two-Dimensional Quantum Model of a Nanotransistor

A mathematical model, and software to implement the model, have been devised to enable numerical simulation of the transport of electric charge in, and the resulting electrical performance characteristics of, a nanotransistor [in particular, a metal oxide/semiconductor field-effect transistor (MOSFET) having a channel length of the order of tens of nanometers] in which the overall device geometry, including the doping profiles and the injection of charge from the source, gate, and drain contacts, are approximated as being two-dimensional. The model and software constitute a computational framework for quantitatively exploring such device-physics issues as those of source-drain and gate leakage currents, drain-induced barrier lowering, and threshold voltage shift due to quantization. The model and software can also be used as means of studying the accuracy of quantum corrections to other semiclassical models.

Govindan, T. R.↗

Molecular Contamination Investigation Facility (MCIF) Capabilities

This facility was used to guide the development of ASTM E 1559 center dot Multiple Quartz Crystal Microbalances (QCMs), large sample and spectral effects capability center dot Several instrumented, high vacuum chamber systems are used to evaluate the molecular outgassing characteristics of materials, flight components and other sensitive surfaces. Test materials for spacecraft/instrument selection center.Test flight components for acceptable molecular outgas levels center dot Determine time/temperature vacuum bake-out requirements center. Data used to set limits for use of materials and specific components center. Provide Input Data to Contamination Transport Models -Applied to numerous flight projects over the past 20 years.

Soules, David M.↗

Barium Depletion in the NSTAR Discharge Cathode After 30,000 Hours of Operation

Dispenser hollow cathodes rely on a consumable supply of barium released by impregnant materials in the pores of a tungsten matrix to maintain a low work function surface. Examinations of cathode inserts from long duration ion engine tests show deposits of tungsten at the downstream end that appear to block the flow of barium from the interior. In addition, a numerical model of barium transport in the insert plasma indicates that the barium partial pressure in the insert may exceed the equilibrium vapor pressure of the dominant barium-producing reaction, and it was postulated previously that this would suppress barium loss in the upstream part of the insert. New measurements of the depth of barium depletion from a cathode insert operated for 30,352 hours reveal that barium loss is confined to a narrow region near the downstream end, confirming this hypothesis.

plasma transport↗

Numerical Behaviour of a Smooth Local Correlation-based Transition Model in a Newton-Krylov Flow Solver

The numerical behaviour of transport-equation-based transition models, including both iterative and grid convergence, is influenced by the source terms. Transition models contain source terms that are large and highly nonlinear, and can be destabilizing in a strong implicit solver. Linearization strategies with varying levels of coupling are evaluated in conjunction with a source-term time step restriction to determine best-practices for solving the SA-sLM2015smooth local correlation-based transition model in an implicit Newton-Krylov flow solver. Achieving deep iterative convergence facilitates a detailed investigation of the grid convergence of these free-transition simulations, which are evaluated relative to fully-turbulent simulations performed using the Spalart-Allmaras turbulence model. Simulations of the NLF0416 general aviation airfoil, VA-2 supercritical airfoil, and NASA CRM-NLF wing-body geometry are performed over a range of grid levels. The results demonstrate that both a fully-coupled linearization strategy and a source-term time step restriction improve nonlinear convergence as the complexity of the free-transition simulations increases. In general, additional grid resolution is required for free-transition simulations relative to fully-turbulent simulations in order to achieve a similar level of accuracy, with the grid convergence of free-transition simulations sensitive to the streamwise grid spacings in the transition regions.

AATT↗

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗

Numerical Study on the Effect of Methane Doping in Hydrogen-Air Rotating Detonation Engines for Various Temperatures and Pressures

Rotating detonation engines (RDEs) have gained attention as a promising technology for future aviation engines. However, the numerical studies of these systems pose severe challenges due to the broad range of spatial and temporal scales. In this study, we use an adaptive mesh refinement based compressible, reactive solver PeleC to resolve the broad range of scales and accurately capture shock and detonation waves using high-resolution numerical schemes. Multi- species transport along with compressible Navier-Stokes equations are solved in the model along with a finite-rate based chemistry model. Embedded boundary method is used to model the complex geometry consisting of discrete fuel nozzles and the combustion chamber. The fuel consisting predominantly of hydrogen is doped with varying levels of methane while air is used as the oxidizer. For a specified total pressure and temperature, the number of stable detonation waves is found to decrease with increasing methane concentration in the fuel mixture. Additionally, no stable detonation solutions are observed for methane composition higher than 20% by volume for the range of operating conditions studied (300-900K, 10 Atm). The increased presence of high temperature zones is also indicative of higher thermal NOx emissions at low methane concentrations. The effect of fuel-air mixture composition and temperature on the detonability, detonation wave structure, mode transitions and their stability are analyzed in this study in addition to their implication on NOx emission.

automatic mesh refinement↗

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗

A viscous flow study of shock-boundary layer interaction, radial transport, and wake development in a transonic compressor

A numerical study based on the 3D Reynolds-averaged Navier-Stokes equation has been conducted to investigate the detailed flow physics inside a transonic compressor. 3D shock structure, shock-boundary layer interaction, flow separation, radial mixing, and wake development are all investigated at design and off-design conditions. Experimental data based on laser anemometer measurements are used to assess the overall quality of the numerical solution. An additional experimental study to investigate end-wall flow with a hot-film was conducted, and these results are compared with the numerical results. Detailed comparison with experimental data indicates that the overall features of the 3D shock structure, the shock-boundary layer interaction, and the wake development are all calculated very well in the numerical solution. The numerical results are further analyzed to examine the radial mixing phenomena in the transonic compressor. A thin sheet of particles is injected in the numerical solution upstream of the compressor. The movement of particles is traced with a 3D plotting package. This numerical survey of tracer concentration reveals the fundamental mechanisms of radial transport in this transonic compressor.

Hah, Chunill↗

An Improved Neutron Transport Algorithm for HZETRN

Long term human presence in space requires the inclusion of radiation constraints in mission planning and the design of shielding materials, structures, and vehicles. In this paper, the numerical error associated with energy discretization in HZETRN is addressed. An inadequate numerical integration scheme in the transport algorithm is shown to produce large errors in the low energy portion of the neutron and light ion fluence spectra. It is further shown that the errors result from the narrow energy domain of the neutron elastic cross section spectral distributions, and that an extremely fine energy grid is required to resolve the problem under the current formulation. Two numerical methods are developed to provide adequate resolution in the energy domain and more accurately resolve the neutron elastic interactions. Convergence testing is completed by running the code for various environments and shielding materials with various energy grids to ensure stability of the newly implemented method.

Slaba, Tony C.↗

Modeling transient edge plasma transport with dynamic recycling

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Space Transport and Fatigue

Space travel presents numerous challenges to fatigue management. Fatigue in space results from many of the same challenges that shiftworkers face on Earth, including inadequate sleep, circadian misalignment, and elevated workload. However, spaceflight also presents unique challenges including circadian misalignment due to non-24- hour light-dark cycles, erratic work schedules due to operational constraints, a microgravity sleep environment, constraints to habitat design, and an isolated and confined living environment. These challenges can be mitigated with fatigue management strategies and countermeasures specifically targeted for space travel. Those strategies can include light interventions, advanced scheduling tools, improved habitat design, and pharmacological countermeasures. What is needed is a future-facing approach to most effectively manage the challenges that may arise as advancements in space travel continue to develop at a rapid pace.

sleep deprivation and performance↗