Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model”

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 685 records · Page 38

Multiscale modeling and nonlinear model predictive control for flue gas desulfurization

The primary source of sulfur dioxide (SO 2 ) emissions is flue gas from fossil fuels-based power plants. SO 2 emissions are known to not only cause health issues, but also have an adverse effect on the environment in various ways. Several Flue Gas Desulfurization (FGD) technologies have been incorporated in power plants. The most popular technology is Wet FGD, where a limestone slurry is used to absorb SO 2 from the flue gas. A detailed droplet scale model describing the instantaneous and finite rate chemistry is developed. The ill-posed Differential-Algebraic Equation (DAE) droplet model is reformulated to a well-posed index-1 DAE through index reduction. The droplet model is integrated with the bulk phase by incorporating gas-liquid mass transfer, and an oxidation reactor model to simulate the dynamic operation of the counter-current spray scrubber. As a result, the model has a well-conditioned Jacobian and overcomes the modeling challenges of previous works and enables numerical solution without requiring carefully selected initialization or specialized solution procedures. The model is successfully validated using power plant measurements, and nonlinear model predictive control (NMPC) studies are demonstrated to optimize recycle stream flowrates to minimize pumping costs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of Cd adsorption to goethite-bacteria composites

The accurate modeling of heavy metal adsorption in complex systems is fundamental for risk assess-ments in soils and associated environments. Bacteria-iron (hydr)oxide associations in soils and sedi-ments play a critical role in heavy metal immobilization. The reduced adsorption of heavy metals on these composites have been widely reported using the component additivity (CA) method. However, there is a lack of a mechanism model to account for these deviations. In this study, we established models for Cd adsorption on goethite-Pseudomonas putida composites at 1:1 and 5:1 mass ratios. Cadmium adsorption on the 5:1 composite was consistent with the additivity method. However, the CA method over predicted Cd adsorption by approximately 8% on the 1:1 composite at high Cd concentration. The deviation was corrected by adding the site blockage reactions between P. putida and goethite. Both CA and “CA-site masking” models for Cd adsorption onto the composites were in line with the ITC data. These results indicate that CA method in simulating Cd adsorption on bacteria-iron oxides composites is limited to low bacterial and Cd concentrations. Therefore the interfacial complexation reactions that occur between iron (hydr)oxides and bacteria should be taken into account when high concentrations of bacteria and heavy metals are present.

Qu, Chenchen↗

Effective moisture penetration depth model for residential buildings: Sensitivity analysis and guidance on model inputs

Moisture buffering of building materials has a significant impact on the building's indoor humidity, and building energy simulations need to model this buffering to accurately predict the humidity. Researchers requiring a simple moisture-buffering approach typically rely on the effective-capacitance model, which has been shown to be a poor predictor of actual indoor humidity. This paper describes an alternative two-layer effective moisture penetration depth (EMPD) model and its inputs. While this model has been used previously, there is a need to understand the sensitivity of this model to uncertain inputs. In this paper, we use the moisture-adsorbent materials exposed to the interior air: drywall, wood, and carpet. We use a global sensitivity analysis to determine which inputs are most influential and how the model's prediction capability degrades due to uncertainty in these inputs. We then compare the model's humidity prediction with measured data from five houses, which shows that this model, and a set of simple inputs, can give reasonable prediction of the indoor humidity.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Equation-based object-oriented modeling and simulation of data center cooling systems

In this paper, we introduce a newly developed open source data center package in the Modelica Buildings library to support modeling and simulation of cooling and control systems of data centers. The data center package contains major thermal and control component models, such as Computer Room Air Handler, Computer Room Air Conditioner, models of different subsystem configurations such as chillers with differently configured waterside economizers, as well as templates for different systems. Two case studies based on this package are performed to investigate the performances of the cooling and electrical system under normal conditions and emergency situations such as a blackout: one is for a data center powered by conventional energy, and the other is for a data center powered by renewable energy. Simulation results show that the dynamic modeling and multi-domain simulation in the Modelica-based tool make it convenient for users to investigate both normal and emergent operations for conventional and renewable data centers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Conditional distribution estimation of building characteristics with diffusion models for urban energy modeling

Understanding current energy consumption behavior in communities is critical for informing future energy use decisions and enabling efficient energy management. Urban energy models, which are used to simulate these energy use patterns, require large datasets with detailed building characteristics for accurate outcomes. However, such detailed characteristics at the individual building level are often unknown and costly to acquire, or unavailable. Through this work, we propose using a generative modeling approach to generate realistic building attributes to fill in the data gaps and finally provide complete characteristics as inputs to energy models. Our model learns complex, building-level patterns from training on a large-scale residential building stock model containing 2.2 million buildings. We employ a tabular diffusion-based framework that is designed to handle heterogeneous (discrete and continuous) features in tabular building data, such as occupancy, floor area, heating, cooling, and other equipment details. We develop a capability for conditional diffusion, enabling the imputation of missing building characteristics conditioned on known attributes. We conduct a comprehensive validation of our conditional diffusion model, firstly by comparing the generated conditional distributions against the underlying data distribution, and secondly, by performing a case study for a Baltimore residential region, showing the practical utility of our approach. Our work is one of the first to demonstrate the potential of generative modeling to accelerate building energy modeling workflows.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Modelling soil erosion in a Mediterranean watershed: Comparison between SWAT and AnnAGNPS models

Highlights: • Outputs, input data and time requirements are compared for SWAT and AnnAGNPS. • Models are assessed also in wet and dry seasons. • Most of the basin is under moderate or high risk of erosion. • SWAT and AnnAGNPS are compared with RUSLE2015 and PESERA. In this study, the simulations generated by two of the most widely used hydrological basin-scale models, the Annualized Agricultural Non-Point Source (AnnAGNPS) and the Soil and Water Assessment Tool (SWAT), were compared in a Mediterranean watershed, the Carapelle (Apulia, Southern Italy). Input data requirements, time and efforts needed for input preparation, strength and weakness points of each model, ease of use and limitations were evaluated in order to give information to users. Models were calibrated and validated at monthly time scale for hydrology and sediment load using a four year period of observations (streamflow and suspended sediment concentrations). In the driest year, the specific sediment load measured at the outlet was 0.89 t ha{sup −1} yr{sup −1}, while the simulated values were 0.83 t ha{sup −1} yr{sup −1} and 1.99 t ha{sup −1} yr{sup −1} for SWAT and AnnAGNPS, respectively. In the wettest year, the specific measured sediment load was 7.45 t ha{sup −1} yr{sup −1}, and the simulated values were 8.27 t ha{sup −1} yr{sup −1} and 6.23 t ha{sup −1} yr{sup −1} for SWAT and AnnAGNPS, respectively. Both models showed from fair to a very good correlation between observed and simulated streamflow and satisfactory for sediment load. Results showed that most of the basin is under moderate (1.4–10 t ha{sup −1} yr{sup −1}) and high-risk erosion (> 10 t ha{sup −1} yr{sup −1}). The sediment yield predicted by the SWAT and AnnAGNPS models were compared with estimates of soil erosion simulated by models for Europe (PESERA and RUSLE2015). The average gross erosion estimated by the RUSLE2015 model (12.5 t ha{sup −1} yr{sup −1}) resulted comparable with the average specific sediment yield estimated by SWAT (8.8 t ha{sup −1} yr{sup −1}) and AnnAGNPS (5.6 t ha{sup −1} yr{sup −1}), while it was found that the average soil erosion estimated by PESERA is lower than the other estimates (1.2 t ha{sup −1} yr{sup −1}).

54 ENVIRONMENTAL SCIENCES↗

What are the effects of Agro-Ecological Zones and land use region boundaries on land resource projection using the Global Change Assessment Model?

Understanding the potential impacts of climate change is complicated by mismatched spatial representations between gridded Earth System Models (ESMs) and Integrated Assessment Models (IAMs), whose regions are typically larger and defined by geopolitical and biophysical criteria. In this study we address uncertainty stemming from the construction of land use regions in an IAM, the Global Change Assessment Model (GCAM), whose regions are currently based on historical climatic conditions (1961-1990). We re-define GCAM’s regions according to projected climatic conditions (2070-2099), and investigate how this changes model outcomes for land use, agriculture, and forestry. By 2100, we find potentially large differences in projected global and regional area of biomass energy crops, fodder crops, harvested forest, and intensive pasture. These land area differences correspond with changes in agricultural commodity prices and production. These results have broader implications for understanding policy scenarios and potential impacts, and for evaluating and comparing IAM and ESM simulations.

54 ENVIRONMENTAL SCIENCES↗

Molecular thermodynamic modeling of aqueous Na + -K + -Mg 2+ -Ca 2+ -SO 4 2- quinary system

As part of the effort to develop a comprehensive molecular thermodynamic model for high salinity produced water based on the electrolyte nonrandom two-liquid theory, this work presents the model results for the aqueous Na + -K + -Mg 2+ -Ca 2+ -SO 4 2- quinary system. With the binary interaction parameters fully parameterized for the water-electrolyte pairs and the electrolyte-electrolyte pairs, the model satisfactorily describes the thermodynamic properties for the aqueous quinary system and its binary, ternary, and quaternary subsystems for the temperature range 273.15–473.15 K and the salt concentrations up to saturation. The model results are compared with the salt solubility and thermodynamic data available for the subsystems with special attention to systems involving (Ca 2+ SO 4 2- ).

Chemistry↗

A unified non-equilibrium phase change model for injection flow modeling

The homogenous relaxation model (HRM) is one of the most widely used models to describe the liquid- gas phase transition. However, in its original formulation, it is unable to handle multispecies vapor-liquid equilibrium (VLE), which limits its applicability to single-component fluids. In this work, a unified non-equilibrium phase change model that considers the VLE of multicomponent mixtures is proposed building upon the HRM's structure. A time factor is introduced to mimic the effect of different phase change timescales due to different mechanisms, e.g., cavitation, flash-boiling, and evaporation. Here to assess the model's performance, computational fluid dynamics simulations of the internal and near-nozzle injection flow of the Engine Combustion Network's Spray G injector were performed using the nine-component PACE-20 fuel with both the unified model and the original HRM. The predicted fuel density in the near-nozzle region matched well with X-ray tomography measurements. The simulation results indicated that, whereas the HRM failed to capture the vaporization due to convective mixing between the fuel and ambient gas, the unified model performed well in predicting the mixing-driven vaporization and the corresponding evaporative cooling. Further comparisons using the nine-component fuel formula and a single-component fuel surrogate demonstrated the unified model's ability to predict preferential vaporization, which affects the predictions of local mixture composition and rate of vaporization. Finally, it is shown that the unified model is capable of representing multiple phase change mechanisms, and the relaxation time factor plays an important role in determining the degree of phase change due to the different mechanisms.

33 ADVANCED PROPULSION SYSTEMS↗

Selection of polynomial chaos bases via Bayesian model uncertainty methods with applications to sparse approximation of PDEs with stochastic inputs

Generalized polynomial chaos (gPC) expansions allow us to represent the solution of a stochastic system using a series of polynomial chaos basis functions. The number of gPC terms increases dramatically as the dimension of the random input variables increases. When the number of the gPC terms is larger than that of the available samples, a scenario that often occurs when the corresponding deterministic solver is computationally expensive, evaluation of the gPC expansion can be inaccurate due to over-fitting. We propose a fully Bayesian approach that allows for global recovery of the stochastic solutions, in both spatial and random domains, by coupling Bayesian model uncertainty and regularization regression methods. It allows the evaluation of the PC coefficients on a grid of spatial points, via (1) the Bayesian model average (BMA) or (2) the median probability model, and their construction as spatial functions on the spatial domain via spline interpolation. The former accounts for the model uncertainty and provides Bayes-optimal predictions; while the latter provides a sparse representation of the stochastic solutions by evaluating the expansion on a subset of dominating gPC bases. Moreover, the proposed methods quantify the importance of the gPC bases in the probabilistic sense through inclusion probabilities. We design a Markov chain Monte Carlo (MCMC) sampler that evaluates all the unknown quantities without the need of ad-hoc techniques. The proposed methods are suitable for, but not restricted to, problems whose stochastic solutions are sparse in the stochastic space with respect to the gPC bases while the deterministic solver involved is expensive. We demonstrate the accuracy and performance of the proposed methods and make comparisons with other approaches on solving elliptic SPDEs with 1-, 14- and 40-random dimensions.

97 MATHEMATICS AND COMPUTING↗

The Single Column Atmosphere Model Version 6 (SCAM6): Not a Scam but a Tool for Model Evaluation and Development

The Single Column Atmosphere Model (SCAM) is a single column model version of the Community Atmosphere Model (CAM). Here we describe the functionality and features of SCAM6, available as part of CAM6 in the Community Earth System Model, version 2 (CESM2). SCAM6 features a wide selection of standard cases, as well as the ability to easily configure a case specified by the user based on a particular point in a CAM 3-D simulation. This work illustrates how SCAM6 reproduces CAM6 results for physical parameterizations, mostly of moisture and clouds. We demonstrate how SCAM6 can be used for model development through different physics selections, as well as with parameter sweep experiments to highlight the sensitivity of cloud properties to the specification of the vapor deposition process in the cloud microphysics. Furthermore, we use SCAM6 to illustrate the sensitivity of CAM6 cloud radiative properties and precipitation to variable drop number (cloud microphysics properties). Finally, we illustrate how SCAM6 can be used to explore critical emergent processes such as cloud feedbacks and show that CAM6 cloud responses to surface warming in stratus and stratocumulus regimes are similar to those in CAM5. CAM6 has a larger response in the shallow cumulus regime than CAM5. CAM6 cloud feedbacks in the shallow cumulus regime are sensitive to turbulence parameters. SCAM6 is thus a valuable tool for model development, evaluation, and scientific analy sis and an important part of the model hierarchy in Community Earth System Model, version 2.

58 GEOSCIENCES↗

A general nonlinear fluid model for reacting plasma-neutral mixtures

A generalized, computationally tractable fluid model for capturing the effects of neutral particles in plasmas is derived. The model derivation begins with Boltzmann equations for singly charged ions, electrons, and a single neutral species. Electron-impact ionization, radiative recombination, and resonant charge exchange reactions are included. Moments of the reaction collision terms are detailed. Moments of the Boltzmann equations for electron, ion, and neutral species are combined to yield a two-component plasma-neutral fluid model. Separate density, momentum, and energy equations, each including reaction transfer terms, are produced for the plasma and neutral equations. The required closures for the plasma-neutral model are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The island coalescence problem: Scaling of reconnection in extended fluid models including higher-order moments

As modeling of collisionless magnetic reconnection in most space plasmas with realistic parameters is beyond the capability of today's simulations, due to the separation between global and kinetic length scales, it is important to establish scaling relations in model problems so as to extrapolate to realistic scales. Recently, large scale particle-in-cell simulations of island coalescence have shown that the time averaged reconnection rate decreases with system size, while fluid systems at such large scales in the Hall regime have not been studied. Here, we perform the complementary resistive magnetohydrodynamic (MHD), Hall MHD, and two fluid simulations using a ten-moment model with the same geometry. In contrast to the standard Harris sheet reconnection problem, Hall MHD is insufficient to capture the physics of the reconnection region. Additionally, motivated by the results of a recent set of hybrid simulations which show the importance of ion kinetics in this geometry, we evaluate the efficacy of the ten-moment model in reproducing such results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modelling of pathologies of the nervous system by the example of computational and electronic models of elementary nervous systems

The paper puts forward principles of action of devices operating similarly to the nervous system and the brain of biological systems. We propose an alternative method of studying diseases of the nervous system, which may significantly influence prevention, medical treatment, or at least retardation of development of these diseases. This alternative is to use computational and electronic models of the nervous system. Within this approach, we represent the brain in the form of a huge electrical circuit composed of active units, namely, neuron-like units and connections between them. As a result, we created computational and electronic models of elementary nervous systems, which are based on the principles of functioning of biological nervous systems that we have put forward. Our models demonstrate reactions to external stimuli and their change similarly to the behavior of simplest biological organisms. The models possess the ability of self-training and retraining in real time without human intervention and switching operation/training modes. In our models, training and memorization take place constantly under the influence of stimuli on the organism. Training is without any interruption and switching operation modes. Training and formation of new reflexes occur by means of formation of new connections between excited neurons, between which formation of connections is physically possible. Connections are formed without external influence. They are formed under the influence of local causes. Connections are formed between outputs and inputs of two neurons, when the difference between output and input potentials of excited neurons exceeds a value sufficient to form a new connection. On these grounds, we suggest that the proposed principles truly reflect mechanisms of functioning of biological nervous systems and the brain. In order to confirm the correspondence of the proposed principles to biological nature, we carry out experiments for the study of processes of formation of connections between neurons in simplest biological objects. Based on the correspondence of function of the created models to function of biological nervous systems we suggest the use of computational and electronic models of the brain for the study of its function under normal and pathological conditions, because operating principles of the models are built on principles imitating the function of biological nervous systems and the brain.

60 APPLIED LIFE SCIENCES↗

Integrated modeling of cryogenic layered highfoot experiments at the NIF

Integrated radiation hydrodynamic modeling in two dimensions, including the hohlraum and capsule, of layered cryogenic HighFoot Deuterium-Tritium (DT) implosions on the NIF successfully predicts important data trends. The model consists of a semi-empirical fit to low mode asymmetries and radiation drive multipliers to match shock trajectories, one dimensional inflight radiography, and time of peak neutron production. Application of the model across the HighFoot shot series, over a range of powers, laser energies, laser wavelengths, and target thicknesses predicts the neutron yield to within a factor of two for most shots. The Deuterium-Deuterium ion temperatures and the DT down scattered ratios, ratio of (10–12)/(13–15) MeV neutrons, roughly agree with data at peak fuel velocities <340 km/s and deviate at higher peak velocities, potentially due to flows and neutron scattering differences stemming from 3D or capsule support tent effects. These calculations show a significant amount alpha heating, 1–2.5× for shots where the experimental yield is within a factor of two, which has been achieved by increasing the fuel kinetic energy. This level of alpha heating is consistent with a dynamic hot spot model that is matched to experimental data and as determined from scaling of the yield with peak fuel velocity. Lastly, these calculations also show that low mode asymmetries become more important as the fuel velocity is increased, and that improving these low mode asymmetries can result in an increase in the yield by a factor of several.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗