Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “advection”

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

Nonsteady Load Responses to Mountain-Generated Turbulence Eddies on the DOE 1.5 MW Wind Turbine at the National Wind Technology Center

Field data collected from the NREL/GE 1.5MW wind turbine and met tower at the NREL Wind Technology Center near Boulder, Colorado from June-October 2018 were analyzed to quantify the impacts of turbulence eddies on the load responses measured from sensors on the main shaft, blade and tower. The passage of individual mountain-generated eddies from the met tower to the wind turbine were critically determined by correlating the optimal time shifts in signal between met tower and nacelle anemometers with mean advection time. Loading responses from mountain eddy passage were compared with atmospheric eddies from the north/south, unimpeded by the mountains, and found to be similar. Whereas time variations in torque were highly correlated with time changes in horizontal eddy velocity, the out-of-plane bending moments on the main shaft (directly forcing the main bearing) were uncorrelated with horizontal eddy velocity. This result is consistent with a previous LES study indicating that the main bearing is forced by asymmetrical interactions between the rotor and turbulence eddies, while power fluctuations respond primarily to advective eddy velocity. Surprisingly, the nacelle anemometer produced statistics very similar to the met tower.

ENGINEERING,WIND ENERGY↗
DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

The Influence of Aqueous Dispersivity on Transport of Radiological Material at LANL MDA G - 20450

Aqueous dispersivity is a phenomenon by which a solute moving in groundwater mixes with water that does not contain the solute and becomes diluted. Solutes take paths of varying lengths through pore networks or experience slower velocities due to interactions with pore walls. Aqueous dispersivity leads to breakthrough curves with earlier arrival than expected in purely advective systems and persistent 'tails,' or late-arriving, dilute concentrations of the solute. Dispersivity can be longitudinal (in the direction of flow) or transverse (orthogonal to flow). LANL has disposed of radioactive material at Material Disposal Area (MDA) G since the late 1950's. Most of the disposed materials at MDA G are low-level waste (LLW), but other types of radioactive wastes are disposed or stored at the site. For example, there is a collection of remote-handled transuranic waste (RH-TRU) waste in 33 shafts, which were intended to be stored retrievably. Aqueous-phase radionuclides included in models of contaminant transport at MDA G are subject to dispersive processes in porous media of the vadose zone and saturated zone. To model dispersivity at MDA G probabilistically, we created site-specific distributions of aqueous dispersivity. First, we performed a comprehensive review of peer-reviewed publications, reports from other performance assessments of MDA G, and documentation from sites that were considered similar based on expert opinion, such as Yucca Mountain. Data were included from a wide range of sources, so that distributions developed based on the data reflect the current state of knowledge in properties of the subsurface environment. Data were grouped based on value type (e.g., determined via experiment or used in a different model), rock type, zone type (saturated or unsaturated), and site, and were weighted in the distributions according to their relevance to the site and zone type. We combined Monte Carlo simulation and bootstrapping techniques to produce approximate distributions of mean aqueous dispersivity. Distributions of dispersivity for saturated and unsaturated zone transport are differentiated by data weight and the length scales over which transport occurs. We find that underestimating or disregarding transverse dispersivity may lead to overstating simulated well capture in the saturated zone at MDA G. Neglecting transverse dispersivity is a common assumption due in part to use of one-dimensional transport models that do not consider solute spreading orthogonal to flow. Also, we find that variability in longitudinal dispersivity is associated with substantial travel time alterations in the MDA G unsaturated zone and in simulated contaminant breakthrough times at a downstream observation plane in the MDA G saturated zone. Finally, at large scales (>1000 m), a linear trend between the dispersivity coefficient and length over which transport occurs may overestimate the dispersivity coefficient. Future modeling work with application to solute spreading could include interacting advective and diffusive zones in the highly fractured basalt below MDA G. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

A remapping scheme for history-dependent material state variables

Here, we present an interpolation scheme for use with material models involving history-dependent state variables. The interpolation scheme is based on the solution of a constrained optimization problem, to identify the smoothest monotone curve which produces an admissible solution. Unlike classical schemes used for transport equations, the proposed methodology accounts for the interdependence of the history variables and produces admissible interpolations of the material's state variables. After introducing a general solution scheme for the problem, a simplified and more computationally efficient method is derived. Finally, as a validation problem for the proposed methods, we consider the advection of internal state variables belonging to material models with highly nonlinear yield surfaces.

42 ENGINEERING↗
DOE OSTI2020

Matrix-free subcell residual distribution for Bernstein finite elements: Monolithic limiting

This paper is focused on the aspects of limiting in residual distribution (RD) schemes for high-order finite element approximations to advection problems. Both continuous and discontinuous Galerkin methods are considered in this work. Discrete maximum principles are enforced using algebraic manipulations of element contributions to the global nonlinear system. The required modifications can be carried out without calculating the element matrices and assembling their global counterparts. The components of element vectors associated with the standard Galerkin discretization are manipulated directly using localized subcell weights to achieve optimal accuracy. Low-order nonlinear RD schemes of this kind were originally developed to calculate local extremum diminishing predictors for flux-corrected transport (FCT) algorithms. In the present paper, we incorporate limiters directly into the residual distribution procedure, which makes it applicable to stationary problems and leads to well-posed nonlinear discrete problems. To circumvent the second-order accuracy barrier, the correction factors of monolithic limiting approaches and FCT schemes are adjusted using smoothness sensors based on second derivatives. Finally, the convergence behavior of presented methods is illustrated by numerical studies for two-dimensional test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

1-D coupled surface flow and transport equations revisited via the physics-informed neural network approach

The de Saint-Venant equation (SVE) and advection–diffusion equation (ADE) are commonly employed to solve solute transport problems in surface water. In this work, we propose a mesh-free method based on the physics-informed neural network (PINN) to solve the one dimensional (1-D) SVE, ADE, and the coupled SVE and ADE (SVE-ADE) under various initial and boundary conditions. The PINN model extends the architecture of deep neural networks (DNNs) with implementation of loss function, which are additionally subject to constraints imposed by the physical laws of SVE and ADE, along with their initial and boundary conditions. In such a manner, PINNs can be quickly steered to the true solution while obeying the physical laws. The results of PINN model are compared with the analytical and/or numerical solutions under various conditions to investigate its accuracy and efficiency in solving the SVE, ADE, and SVE-ADE. Our results indicate PINN can accurately simulate the shock wave morphology and avoid numerical dissipation in unsteady flow condition. The PINN method outweighed traditional numerical methods in several aspects, including its ability to function with small amounts of data, no grid discretization, and random selection of sampling points, etc. Additionally, the PINN method is also suitable for solving inverse problems with sparse and noisy data. With 1% noise and 2000 initial and boundary condition points (N u ), the errors of the estimated flow rate (v) and diffusion coefficient (D) are 0.003% and 0.105%, respectively, which indicate the accuracy and robustness of the proposed method. Finally, our results indicate the capability and robustness of the proposed PINN methodology for solving multi-physics problems, irrespective of the presence of sparse and noisy data in the training dataset.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coupled neutronics and species transport simulation of the Molten Salt Reactor Experiment

This paper presents the development of coupling between the molten salt reactor species transport code Mole and the reactor physics code Griffin. Here, in this study, tracking of delayed neutron precursors was investigated in the Molten-Salt Reactor Experiment (MSRE), accounting for changes in fuel flow velocity as a function of position in the primary loop. The neutron transport calculations in Griffin were performed using 11 energy groups, and the species advection calculations in Mole used 6 delayed neutron precursor groups to predict spatial distribution of the neutron flux and neutron precursors in the MSRE. Mole–Griffin was used to calculate $k$ eff and $β$ eff in the reactor as a function of different volumetric flow rates.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

A comparison of Eulerian and semi-Lagrangian approaches for modeling stream water quality

This paper describes and compares some of the advantages and limitations of Eulerian and Lagrangian approaches to water quality modeling and introduces a mixed Eulerian-Lagrangian (or semi-Lagrangian) methodology that captures the strengths of both approaches. The semi-Lagrangian modeling approach is applied to advection-dominated rivers, and flexibly ensures unconditional stability for all time step durations and grid segmentations. The semi-Lagrangian modeling approach is demonstrated by applying it to estimate the dissolved oxygen concentrations in the Sava River in Slovenia, focusing on aspects of the methodology and findings that would be of broad interest to managers of water quality in fluvial water bodies. Results of comparisons of the semi-Lagrangian model with the Eulerian-based QUAL2K model in steady and non-steady scenarios demonstrate that while both models are fully capable of producing satisfactory results when optimally configured, the semi-Lagrangian approach offers accuracy and stability without sensitivity to the interaction of time step size and computational grid segmentation scheme.

Sava River↗
DOE OSTI2021

Hyperbolic Reformulation Approach to Enable Efficient Simulation of Groundwater Flow and Reactive Transport

We apply Cattaneo's relaxation approach to the one-dimensional coupled Boussinesq groundwater flow and advection-diffusion-reaction equations, commonly used in engineering applications to simulate contaminant transport in the subsurface. The diffusion-type governing equations are reformulated as a hyperbolic system, augmented by an equation that can be interpreted as a momentum balance. The hyperbolization enables an efficient unified computation of the primary variable and its gradients, for example piezometric head and unit discharge in the Boussinesq equation. An augmented Roe scheme is used to solve the hyperbolic system. The hyperbolized system of equations is studied in a set of steady state and transient test cases with idealized geometry. These test cases confirm the equivalence of the hyperbolic system to its original formulation. The larger time step size of the hyperbolic equation is verified theoretically by means of a stability analysis and numerically in the test cases. Finally, a reach-scale application of flow and transport across a river meander is considered. This application case shows that the performance of the hyperbolic relaxation approach holds for more realistic groundwater flow and transport problems, relevant to water resources management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Utah FORGE: Geochemical Data for Cold Groundwaters and Produced Geothermal Fluids

Geochemical data for cold groundwaters and produced geothermal fluids around the Utah FORGE site. The data is compiled into four tables in the attached Excel File. Table 1 is a compilation of compositions (anions, cations, weak acids, oxygen, hydrogen, and carbon isotopes) for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 2 is a compilation of noble gas (He, Ne, Ar) and He and Ne isotopic compositions for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 3 provides values for calculated advective and diffusive fluxes of helium. Table 4 provides values of calculated subsurface stored heat between the Opal Mound fault and the Utah FORGE site, which are related to volumes of recently solidified magmatic heat sources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Manifold Learning-Based Polynomial Chaos Expansions for High-Dimensional Surrogate Models

In this work we introduce a manifold learning-based method for uncertainty quantification (UQ) in systems describing complex spatiotemporal processes. Our first objective is to identify the embedding of a set of high-dimensional data representing quantities of interest of the computational or analytical model. For this purpose, we employ Grassmannian diffusion maps, a two-step nonlinear dimension reduction technique which allows us to reduce the dimensionality of the data and identify meaningful geometric descriptions in a parsimonious and inexpensive manner. Polynomial chaos expansion is then used to construct a mapping between the stochastic input parameters and the diffusion coordinates of the reduced space. An adaptive clustering technique is proposed to identify an optimal number of clusters of points in the latent space. The similarity of points allows us to construct a number of geometric harmonic emulators which are finally utilized as a set of inexpensive pretrained models to perform an inverse map of realizations of latent features to the ambient space and thus perform accurate out-of-sample predictions. Thus, the proposed method acts as an encoder-decoder system which is able to automatically handle very high-dimensional data while simultaneously operating successfully in the small-data regime. The method is demonstrated on two benchmark problems and on a system of advection-diffusion-reaction equations which model a first-order chemical reaction between two species. In all test cases, the proposed method is able to achieve highly accurate approximations which ultimately lead to the significant acceleration of UQ tasks.

42 ENGINEERING↗
DOE OSTI2024

Characterizing Land-Atmosphere Interactions During the Afternoon-to-Evening Transition Using ARM SGP Observations

The land-atmosphere (LA) system consists of the soil, vegetation, and lowest layer of the atmosphere. The energy and water transport processes in this system contribute to the structure and the daily cycle of the atmospheric boundary layer (ABL), which is the layer between the land surface and the free troposphere that controls the exchange of momentum, energy, and water. The state of the ABL is the result of complex two-way interactions between local processes such as surface fluxes and entrainment at the ABL top as well as large-scale and mesoscale processes such as advection and subsidence.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

A Linearized Prognostic Cloud Scheme in NASAs Goddard Earth Observing System Data Assimilation Tools

A linearized prognostic cloud scheme has been developed to accompany the linearized convection scheme recently implemented in NASA's Goddard Earth Observing System data assimilation tools. The linearization, developed from the nonlinear cloud scheme, treats cloud variables prognostically so they are subject to linearized advection, diffusion, generation, and evaporation. Four linearized cloud variables are modeled, the ice and water phases of clouds generated by large-scale condensation and, separately, by detraining convection. For each species the scheme models their sources, sublimation, evaporation, and autoconversion. Large-scale, anvil and convective species of precipitation are modeled and evaporated. The cloud scheme exhibits linearity and realistic perturbation growth, except around the generation of clouds through large-scale condensation. Discontinuities and steep gradients are widely used here and severe problems occur in the calculation of cloud fraction. For data assimilation applications this poor behavior is controlled by replacing this part of the scheme with a perturbation model. For observation impacts, where efficiency is less of a concern, a filtering is developed that examines the Jacobian. The replacement scheme is only invoked if Jacobian elements or eigenvalues violate a series of tuned constants. The linearized prognostic cloud scheme is tested by comparing the linear and nonlinear perturbation trajectories for 6-, 12-, and 24-h forecast times. The tangent linear model performs well and perturbations of clouds are well captured for the lead times of interest.

Precipitation↗
NASA NTRS2024

Increasing Contribution of the Atmospheric Vertical Motion to Precipitation in a Warming Climate

Global warming already influences precipitation, with more intense precipitation in many locations. Although the ‘wet-get-wetter, dry-get-drier’ tendency in mean precipitation holds in many locations, the situations for precipitation extremes are more complex, due to changes in dynamic and thermodynamic influences on atmospheric moisture distributions. Here, we build a dynamically interactive atmospheric moisture model for the present (2006–2025) and the future climate (2081–2100), using outputs from coupled ocean-atmosphere general circulation models. We find that the dynamic process of vertical advection of moisture dominates the same-day precipitation, while the smaller impact of the thermodynamic process provides available moisture for several days. As climate warms, we find that the dynamical-induced precipitation more completely exhausts the vertically-integrated moisture and the distribution of the dynamic process’s impact on precipitation exhibits a greater spread in the warmer future. The dynamical process is primarily responsible for more extreme heavy precipitation as climate warms, at all latitudes.

Climate and Earth system modelling↗
DOE OSTI2020

Vertical Transport of Water/Contaminant Due to Intermittent Water Influx on the Ground Surface of Vadose Zone - 20088

For movement of moisture and transport of contaminants in the vertical direction, two equations, the Richards equation and the advection-dispersion equation, are solved in a vadose zone for prescribed water and contaminant influx conditions at the ground surface. Two different boundary conditions are considered, one with a single peak of sinusoidal variation and four peaks of finite duration, each of which is sinusoidal in shape. The moisture migration is solved and then the contaminant transport, which is coupled to moisture movement, is solved. In the case of a single influx on the ground surface, at first the hydraulic head and moisture content show a peak. Then above the peak and below the ground surface, the peaks tend to die out, showing a tendency to return to the initial state. The vertical profile of concentration also shows a single peak, which soon becomes stagnant because downward water migration is faster than concentration profile movement. In the case of four intermittent array of peaks, the hydraulic head and moisture content profiles tend to merge and combine together. The concentration profiles show appearance of four peaks clearly. However, although the peak value of concentration remains quite constant, the values may vary with the sorption coefficient because of repeated influx of water and contaminant through the top boundary, the ground surface. The penetration depth decreases with an increasing sorption coefficient. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Fundamental Study of Migration Behavior of Cesium under the Unsaturated Condition through One-dimensional Column Simulating Soil around the Disposal Site - 20262

Using a packed column simulating unsaturated soil around the disposal site, flow experiments of cesium ions were conducted in order to examine the difference of retardation effect between saturated and unsaturated conditions. The column was packed with silica sand particles of 350 μm in diameter. In this study, it was confirmed that the X-ray CT images of the column showed that the water saturation was mostly uniform within the column under unsaturated conditions. The values of water saturation in the flow experiments of Cs ions were in a range from 0.85 to 0.93. The retardation coefficients(Rd) were estimated by applying the one-dimensional advection-dispersion equation to the experimental results. The comparison of the Rd-values estimated experimentally and theoretically showed that the experimental value was estimated to be larger than the theoretical value. This means that the theoretical value evaluated through a conventional method would bring the conservative value as an apparent retardation coefficient, at least when the flow direction is from top to bottom as shown in this study. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Choosing the Best Modeling Platform for Radiological Risk Assessment Models - 20468

Radiological risk assessments, in the form of performance or safety assessments, are often required under regulations or guidance for remediation of contaminated land, decommissioning of contaminated buildings or structures, and radioactive waste disposal. These risk assessments are usually supported by fate and transport models that address decay and ingrowth of radionuclides, as well as their movement through engineered systems and the natural environment. These models are often projected thousands, or more, years into the future, largely because the radioactive species change through decay and ingrowth, and hence the magnitude of the radioactive effect changes with time. There are many computer codes that are available to address this type of modeling. They range from addressing specific pathways or processes such as infiltration of water, groundwater, surface water, air, biota, diffusion and advection of water and gases, to those that try to couple all processes together to evaluate the impact of fate and transport through the entire system to places in space and time to which access is assumed. These different types of codes are sometimes separated with the monikers process-level and systems-level codes, although it is often not clear that this separation does justice to the capabilities of the many codes that are available to evaluate fate and transport of radionuclides. The focus of this paper is the latter group of modeling codes. Several systems level modeling codes exist and are used. There are differences between these codes in terms of utility, flexibility, complexity and cost. The purpose of this paper is to compare a few of these codes in the context of work currently being performed by the International Atomic Energy Agency (IAEA) Modeling and Data for Radiological Impact Assessments (MODARIA) II Working Group 1 (WG1). The MODARIA II WG1's main focus is how stakeholder engaged decision analysis can, or should, be applied to radiological contamination problems so that better, longstanding, sustainable, solutions are reached. However, the WG1 also recognizes the potential impact of the modeling tools that are chosen to address radiological risk, which is often a primary objective of decision making for radiological problems. Other objectives might also be important, such as constraining costs, obtaining financing, minimizing impact on ecosystems, saving cultural resources, saving jobs, farmland, environmental justice, etc., in a full decision analysis for a given radiological contamination problem, but none of these other objectives have the same types of complex modeling needs as minimize radiological dose. Consequently, a further focus of WG1 is to evaluate the potential impacts on decision making of the choice of fate and transport, and risk assessment, modeling codes that are used to support decision making. The WG1 will produce a report at the end of 2020 that will focus on an approach to effective decision making and stakeholder engagement. The report will also consider the role that performance assessment modeling should play in the decision-making process, including the impact of the choice of modeling tools or computer codes on risk-informed decision making. Several sites around the World have been made available by Member States for these model comparisons, and several modeling tools have been considered. However, the focus of this paper is on two of the sites, one in Belgium and one in Ukraine, and on three of the tools: NORMALYSA (NORM And Legacy Site Assessment); GoldSim{sup C}, and AMBER{sup C}. The final report from this working group will also cover other modeling tools, including RESRAD, and PC-Cream{sup R}. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗