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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Investigation of Biomass Fouling on Screw Feeder in Preconversion of Pyrolysis

Biomass is a source of renewable energy and can undergo thermochemical processes to be converted into gaseous and liquid fuels. Thermal decomposition-induced biomass fouling on the screw feeder is a major challenge in preconversion of pyrolysis, leading to plugging of the feed line. Here, this work intends to gain fundamental understanding of the biomass fouling phenomena across the heat gradient on the screw feeder via a combination of materials characterization and thermal simulation. An actual screw feeder that failed because of biomass fouling was examined to reveal the deposit’s morphology, composition, and mechanical properties. Thermal simulations were developed for both the screw feeder and the wood particles being transported using a combined analytical and numerical approach. The simulated temperature profiles correlated well with the deposit observations, based on which potential mitigations were proposed from both the operation and screw design perspectives and supported by preliminary experimental validation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Progress Report on Model Development for the Transport of Aerosol through Microchannels

This report summarizes the current progress in the development of a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose is to introduce a generic, reliable numerical model for the prediction of aerosol transport, deposition, and plugging in leak paths while accounting for potential plugging formation, to a user community involving researchers, regulators, and industry. In that regard, a Graphical User Interface (GUI) was generated by integrating the individual MATLAB scripts that make up the model and adding additional features to aid the general understanding of the problem to a user. This report focuses on the model development and features included. Further, predictions from the GUI are compared with experimental data for validation. The strength of the GUI is the ability of a user to plug in basic parameters such as the initial pressure conditions and canister model specifications in order to obtain a first principal approximation of vital information such as the blowdown pressure differential, aerosol penetration, and deposition as a function of time. The user can obtain this without the know-how of the underlying MATLAB scripts and thus enables the model to be readily applied by regulators, industry, and shareholders to reduce the uncertainty in the off-site radiological consequences. The report also lists future tasks under this work scope. This includes ongoing efforts to include additional crack geometries (divergent and divergent-convergent slots) to the model; GUI development and improvement based on feedback from users, and integration of aerosol source term data from Sibling Pin tests into the model as we approach realistic canister stress corrosion cracking–induced aerosol release scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development of the Table of Initial Isolation and Protective Action Distances for the 2024 Emergency Response Guidebook

The transportation of hazardous materials creates numerous opportunities for the release of toxic substances into the environment, whether caused by traffic accidents, train derailments, equipment failures, or human error. Such releases can pose acute hazards to the general public and to emergency response personnel who are the first to arrive at the scene. To help first responders determine whether a shipment is potentially hazardous and decide what actions should be taken if a toxic spill does occur, the Emergency Response Guidebook (ERG) is published by the U.S. Department of Transportation (DOT), Transport Canada, and the Secretariat of Transport and Communications of Mexico; with contributions from Centro de Informaciòn Quìmica para Emergencias of Argentina. The most recent version is the 2024 edition of the ERG (ERG 2024), titled 2024 Emergency Response Guidebook (ERG2024). The ERG provides essential information about firefighting, spill response, and potential public health effects. For chemicals that are toxic by inhalation (TIH) and chemicals that produce TIH gases upon reaction with water (TIH by water reactivity or TIHWR), the ERG provides initial isolation distances (IIDs) and protective action distances (PADs). The IID defines the radius of the zone around the spill that should be accessed solely by people who are directly involved in emergency response. The PAD is the distance downwind of the source of the release within which persons should be either evacuated or sheltered in place, depending on the severity of the incident and the nature of the population (e.g., density, age, health).

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
NASA NTRS1993

Development of deterministic transport methods for low energy neutrons for shielding in space

Transport of low energy neutrons associated with the galactic cosmic ray cascade is analyzed in this dissertation. A benchmark quality analytical algorithm is demonstrated for use with BRYNTRN, a computer program written by the High Energy Physics Division of NASA Langley Research Center, which is used to design and analyze shielding against the radiation created by the cascade. BRYNTRN uses numerical methods to solve the integral transport equations for baryons with the straight-ahead approximation, and numerical and empirical methods to generate the interaction probabilities. The straight-ahead approximation is adequate for charged particles, but not for neutrons. As NASA Langley improves BRYNTRN to include low energy neutrons, a benchmark quality solution is needed for comparison. The neutron transport algorithm demonstrated in this dissertation uses the closed-form Green's function solution to the galactic cosmic ray cascade transport equations to generate a source of neutrons. A basis function expansion for finite heterogeneous and semi-infinite homogeneous slabs with multiple energy groups and isotropic scattering is used to generate neutron fluxes resulting from the cascade. This method, called the FN method, is used to solve the neutral particle linear Boltzmann transport equation. As a demonstration of the algorithm coded in the programs MGSLAB and MGSEMI, neutron and ion fluxes are shown for a beam of fluorine ions at 1000 MeV per nucleon incident on semi-infinite and finite aluminum slabs. Also, to demonstrate that the shielding effectiveness against the radiation from the galactic cosmic ray cascade is not directly proportional to shield thickness, a graph of transmitted total neutron scalar flux versus slab thickness is shown. A simple model based on the nuclear liquid drop assumption is used to generate cross sections for the galactic cosmic ray cascade. The ENDF/B V database is used to generate the total and scattering cross sections for neutrons in aluminum. As an external verification, the results from MGSLAB and MGSEMI were compared to ANISN/PC, a routinely used neutron transport code, showing excellent agreement. In an application to an aluminum shield, the FN method seems to generate reasonable results.

Ganapol, Barry↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A low-rank power iteration scheme for neutron transport criticality problems

Computing effective eigenvalues for neutron transport often requires a fine numerical resolution. Here, the main challenge of such computations is the high memory effort of classical solvers, which limits the accuracy of chosen discretizations. In this work, we derive a method for the computation of effective eigenvalues when the underlying solution has a low-rank structure. This is accomplished by utilizing dynamical low-rank approximation (DLRA), which is an efficient strategy to derive time evolution equations for low-rank solution representations. The main idea is to interpret the iterates of the classical inverse power iteration as pseudo-time steps and apply the DLRA concepts in this framework. In our numerical experiment, we demonstrate that our method significantly reduces memory requirements while achieving the desired accuracy. Analytic investigations show that the proposed iteration scheme inherits the convergence speed of the inverse power iteration, at least for a simplified setting.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1980

On the calculation of turbulent heat transport downstream from an abrupt pipe expansion

A numerical study of flow and heat transfer in the separated flow region produced by an abrupt pipe explosion is reported, with emphasis on the region in the immediate vicinity of the wall where turbulent transport gives way to molecular conduction and diffusion. The analysis is based on a modified TEACH-2E program with the standard k-epsilon model of turbulence. Predictions of the experimental data of Zemanick and Dougall (1970) for a diameter ratio of 0.54 show generally encouraging agreement with experiment. At a diameter ratio of 0.43 different trends are discernable between measurement and calculation, though this appears to be due to effects unconnected with the wall region studied here.

Chieng, C. C.↗
NASA NTRS1983

A coupled ice-ocean model of upwelling in the marginal ice zone

A dynamical coupled ice-ocean numerical model for the marginal ice zone (MIZ) is suggested and used to study upwelling dynamics in the MIZ. The nonlinear sea ice model has a variable ice concentration and includes internal ice stress. The model is forced by stresses on the air/ocean and air/ice surfaces. The main coupling between the ice and the ocean is in the form of an interfacial stress on the ice/ocean interface. The ocean model is a linear reduced gravity model. The wind stress exerted by the atmosphere on the ocean is proportional to the fraction of open water, while the interfacial stress ice/ocean is proportional to the concentration of ice. A new mechanism for ice edge upwelling is suggested based on a geostrophic equilibrium solution for the sea ice medium. The upwelling reported in previous models invoking a stationary ice cover is shown to be replaced by a weak downwelling due to the ice motion. Most of the upwelling dynamics can be understood by analysis of the divergence of the across ice edge upper ocean transport. On the basis of numerical model, an analytical model is suggested that reproduces most of the upwelling dynamics of the more complex numerical model.

Roed, L. P.↗
NASA NTRS1988

Transient natural convection heat and mass transfer in a rectangular enclosure - A numerical analysis

A numerical analysis of transient heat and solute transport across a rectangular cavity with combined horizontal temperature and concentration gradients is performed by a numerical method based on the SIMPLE. Numerical results show that the average Nusselt and Sherwood numbers both decrease markedly when the solutal and thermal buoyancy forces act in the opposite directions. When the solutal and thermal buoyancy forces act in the same directions, however, the average Sherwood number increases significantly and yet the average Nusselt number decreases slightly.

Han, Samuel S.↗
NASA NTRS1995

The mesoscale forcing of a midlatitude upper-tropospheric jet streak by a simulated convective system. 2: Kinetic energy and resolution analysis

A kinetic energy (KE) analysis of the forcing of a mesoscale upper-tropospheric jet streak by organized diabatic processes within the simulated convective system (SCS) that was discussed in Part 1 is presented in this study. The relative contributions of the ageostrophic components of motion to the generation of KE of the convectively generated jet streak are compared, along with the KE generation by the rotational (nondivergent) and irrotational (divergent) mass transport. The sensitivity of the numerical simulations of SCS development to resolution is also briefly examined. Analysis within isentropic coordinates provides for an explicit determination of the influence of the diabatic processes on the generation of KE. The upper-level production of specific KE is due predominatly to the inertial advective ageostrophic component (IAD), and as such represents the primary process through which the KE of the convectively generated jet streak is realized. A secondary contribution by the inertial diabatic (IDI) term is observed. Partitioning the KE generation into its rotational and irrotational components reveals that the latter, which is directly linked to the diabatic heating within the SCS through isentropic continuity requirements, is the ultimate source of KE generation as the global area integral of generation by the rotational component vanishes. Comparison with an identical dry simulation reveals that the net generation of KE must be attributed to latent heating. Both the IAD and IDI ageostrophic components play important roles in this regard. Examination of results from simulations conducted at several resolutions supports the previous findings in that the effects of diabatic processes and ageostrophic motion on KE generation remain consistent. Resolution does impact the location and timing of SCS development, a result that has important implications in forecasting the onset of convection that develops from evolution of the large-scale flow and moisture transport. Marked differences are observed in the momentum field aloft subsequent to the life cycle of the SCS in the 1 deg, 30-level base case (MP130) simulation discussed in Part 1 versus its 2 deg counterparts in that the MP130 simulation with higher spatial resolution contains from 14% to 30% more total KE.

Wolf, Bart J.↗
NASA NTRS1991

Suppression of turbulent transport by a weak magnetic field

Two-dimensional numerical simulations with high spatial resolution are used to study the effects of a large-scale magnetic field on its turbulent transport. It is commonly believed that the nonlinear back-reaction of the magnetic field on the turbulence becomes important when the field strength is close to equipartition. However, turbulent diffusion is effectively reduced even when the large-scale field is much weaker than equipartition.

Cattaneo, Fausto↗
DOE OSTI2020

Vortex guidance and transport in channeled pinning arrays

Here, we numerically examine vortices interacting with pinning arrays where a portion of the pinning sites have been removed in order to create coexisting regions of strong and weak pinning. The region without pinning sites acts as an easy-flow channel. For driving in different directions with respect to the channel, we observe distinct types of vortex flow. When the drive is parallel to the channel, the flow first occurs in the pin free region followed by a secondary depinning transition in the pinned region. At high vortex densities there is also an intermediate plastic flow phase due to the coupling between the weak and strong pinning regions. For driving applied perpendicular to the channel, we observe a jammed phase in which vortices accumulate on the boundary of the pinned region due to the vortex-vortex repulsion, while at higher drives the vortices begin to flow through the pinning array. For driving at an angle to the channel, depending on the filling we observe a drive-induced reentrant pinning effect as well as negative differential mobility which occurs when vortices move from the unpinned to the pinned portion of the sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Current collection through the transport of electrons across magnetic field lines

Numerical simulations are used to investigate the collection of electrons by a long conducting cylinder in a flowing plasma. The plasma flow simulates the relative motion between a spacecraft and the plasma. The sheath structures and the levels of electron current collections for the cases with and without an ambient magnetic field are examined. It is found that for the flow perpendicular to the magnetic field, the current is considerably enhanced depending on the relative drift velocity. In the case of a nonzero magnetic field perpendicular to the cylinder axis, the potential structure is a 2D double layer with dimensions L-perpendicular much less than L-parallel. L-perpendicular is found to be the current limiting radius given by the Parker-Murphy model (1967). For the flow along the ambient magnetic field, the electron current is found to be smaller than that for the flow perpendicular to this field. This is explained in terms of the potential structures.

Singh, N.↗
DOE OSTI2021

Plant SWEETs: from sugar transport to plant–pathogen interaction and more unexpected physiological roles

Abstract Sugars Will Eventually be Exported Transporters (SWEETs) have important roles in numerous physiological mechanisms where sugar efflux is critical, including phloem loading, nectar secretion, seed nutrient filling, among other less expected functions. They mediate low affinity and high capacity transport, and in angiosperms this family is composed by 20 paralogs on average. As SWEETs facilitate the efflux of sugars, they are highly susceptible to hijacking by pathogens, making them central players in plant–pathogen interaction. For instance, several species from the Xanthomonas genus are able to upregulate the transcription of SWEET transporters in rice (Oryza sativa), upon the secretion of transcription-activator-like effectors. Other pathogens, such as Botrytis cinerea or Erysiphe necator, are also capable of increasing SWEET expression. However, the opposite behavior has been observed in some cases, as overexpression of the tonoplast AtSWEET2 during Pythium irregulare infection restricted sugar availability to the pathogen, rendering plants more resistant. Therefore, a clear-cut role for SWEET transporters during plant–pathogen interactions has so far been difficult to define, as the metabolic signatures and their regulatory nodes, which decide the susceptibility or resistance responses, remain poorly understood. This fuels the still ongoing scientific question: what roles can SWEETs play during plant–pathogen interaction? Likewise, the roles of SWEET transporters in response to abiotic stresses are little understood. Here, in addition to their relevance in biotic stress, we also provide a small glimpse of SWEETs importance during plant abiotic stress, and briefly debate their importance in the particular case of grapevine (Vitis vinifera) due to its socioeconomic impact.

Breia, Richard↗
DOE OSTI2023

Model data for numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

This model-data archive pertains to a study aimed at the numerical evaluation of the photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively). Separating effects of SSZ and HSZ are important for accurately representing stream function as HSZs and SSZs expose solutes to significantly different biogeochemical conditions. We perform numerical experiments using a multiscale reactive transport model for stream corridors implemented in ATS code, which allows for representing multiple storage zones with their respective travel time distributions and biogeochemistry. For each of the numerical experiment, we provide in this model-data archive python wrapper to drive ATS (forward_model.py), synthetic observation (synthetic_btc.py, BTC_observed.csv), forward model files including ATS input (multiscale_transport.tpl), PFLOTRAN inputs (reactions_channel.tpl, reactions_hz.tpl, reactions_sz.tpl), MCMC (mcmc_run.py), predictive uncertainty (pred_uncert.py, BTCs_simulated_day.csv, BTCs_simulated_night.csv), MCMC outputs (tracer_test_GR.npy, tracer_test_logps.npy, tracer_test_parameters.npy) and Jupyter notebook for post-processing and visualization (post-processing.ipynb). For the denitrification application, ATS input file (denitrification_multisubgrid.tpl), PFLOTRAN input files (denitrification_channel.in, denitrification_hz.in), predictive uncertainty (pred_uncert_denitrification.py, BTCs_simulated_DO.csv, BTCs_simulated_DOC.csv, BTCs_simulated_Nitrate.csv).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Future of Land-Based Wind Turbine Rotor Technology: The Perspective from NREL

The Big Adaptive Rotor (BAR) Project looks at potential innovative pathways for future land-based wind turbine technology. BAR is led by the National Renewable Energy Laboratory and is sponsored by the Wind Energy Technology Office of the US Department of Energy. This talk will describe innovations and challenges faced by the teams working on BAR. The innovations include numerical models predicting the stability and performance of highly flexible blades, downwind rotors to maximize wind farm power, and controlled bending of 100-meter-long blades during rail transport. The challenges of developing predictive numerical models and formulating successful value propositions supporting new technologies such as distributed aero control devices will also be discussed.

BAR↗
DOE OSTI2022

Deep-learning-based workflow for boundary and small target segmentation in digital rock images using UNet++ and IK-EBM

We report three-dimensional (3D) X-ray micro-computed tomography (μCT) has been widely used in petroleum engineering because it can provide detailed pore structural information for a reservoir rock, which can be imported into a pore-scale numerical model to simulate the transport and distribution of multiple fluids in the pore space. The partial volume blurring (PVB) problem is a major challenge in segmenting raw μCT images of rock samples, which impacts boundaries and small targets near the resolution limit. We developed a deep-learning (DL)-based workflow for accurate and fast partial volume segmentation. The DL model's performance depends primarily on the training data quality and model architecture. This study employed the entropy-based-masking indicator kriging (IK-EBM) to segment 3D Berea sandstone images as training datasets. The comparison between IK-EBM and manual segmentation using a 3D synthetic sphere pack, which had a known ground truth, showed that IK-EBM had higher accuracy on partial volume segmentation. We then trained and tested the UNet++ model, a state-of-the-art supervised encoder-decoder model, for binary (i.e., void and solid) and four-class segmentation. We compared the UNet++ with the commonly used U-Net and wide U-Net models and showed that the UNet++ had the best performance in terms of pixel-wise and physics-based evaluation metrics. Specifically, boundary-scaled accuracy demonstrated that the UNet++ architecture outperformed the regular U-Net architecture in the segmentation of pixels near boundaries and small targets, which were subjected to the PVB effect. Feature map visualization illustrated that the UNet++ bridged the semantic gaps between the feature maps extracted at different depths of the network, thereby enabling faster convergence and more accurate extraction of fine-scale features. The developed workflow significantly enhances the performance of supervised encoder-decoder models in partial volume segmentation, which has extensive applications in fundamental studies of subsurface energy, water, and environmental systems.

02 PETROLEUM↗
DOE OSTI2021

Direct Observation of the Depth of Active Groundwater Circulation in an Alpine Watershed

Abstract The depth of active groundwater circulation is a fundamental control on stream flows and chemistry in mountain watersheds, yet it remains challenging to characterize and is rarely well constrained. We collected hydraulic conductivity, hydraulic head, temperature, chemical, noble gas, and 3 H/ 3 He groundwater age data from discrete levels in two boreholes 46 and 81 m deep in an alpine watershed, in combination with chemical and age data from shallow groundwater discharge, to discern groundwater flow rates at different depths and directly observe active and inactive groundwater. Vertical head gradients are steep (average of 0.4) and thermal profiles are consistent with typical linear conductive continental geotherms. Groundwater deeper than ∼20 m is distinct from shallow groundwater and creek water in its chemistry, noble gas signature, and age (dominantly >65 years compared to <9 years). Together these results suggest low vertical groundwater flow velocities and a relatively shallow active circulation depth of ∼20 m. This hypothesis is tested with a simple 2‐D numerical fluid flow and heat transport model representing a hillslope transect through the two boreholes. The modeling indicates that the subhorizontally bedded sedimentary rocks underlying the basin are highly anisotropic with low vertical hydraulic conductivity, and at most ∼10% of bedrock recharge (equivalent to <2% of stream baseflow) flows below a depth of 20 m. The study demonstrates the considerable value of discrete‐depth hydrogeologic, chemical, and age data for determining active circulation depth, and illustrates an approach for maximizing the utility of individual boreholes drilled for mountain bedrock aquifer characterization.

54 ENVIRONMENTAL SCIENCES↗