Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeling and Simulation”

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 973 records · Page 54

Evaluating WRF-Chem aerosol indirect effects in Southeast Pacific marine stratocumulus during VOCALS-REx

We evaluate a regional-scale simulation with the WRF-Chem model for the VAMOS (Variability of the American Monsoon Systems) Ocean-Cloud-Atmosphere-Land Study Regional Experiment (VOCALS-REx), which sampled the Southeast Pacific's persistent stratocumulus deck. Evaluation of VOCALS-REx ship-based and aircraft observations focuses on analyzing how aerosol loading affects marine boundary layer (MBL) dynamics and cloud microphysics. We compare local time series and campaign averaged longitudinal gradients, and highlight differences in model simulations with (W) and without wet (NW) deposition processes. The higher aerosol loadings in the NW case produce considerable changes in MBL dynamics and cloud microphysics, in accordance with the established conceptual model of aerosol indirect effects. These include increase in cloud albedo, increase in MBL and cloud heights, drizzle suppression, increase in liquid water content, and increase in cloud lifetime. Moreover, better statistical representation of aerosol mass and number concentration improves model fidelity in reproducing observed spatial and temporal variability in cloud properties, including top and base height, droplet concentration, water content, rain rate, optical depth (COD) and liquid water path (LWP). Together, these help to quantify confidence in WRF-Chem's modeled aerosol-cloud interactions, while identifying structural and parametric uncertainties including: irreversibility in rain wet removal; overestimation of marine DMS and sea salt emissions and accelerated aqueous sulfate conversion. Our findings suggest that WRF-Chem simulates marine cloud-aerosol interactions at a level sufficient for applications in forecasting weather and air quality and studying aerosol climate forcing, including the reliability required for policy analysis and geo-engineering applications.

54 ENVIRONMENTAL SCIENCES↗

Hardware In the Loop for Demand Flexibility (HIL4DF) v1.0

The software package in question is a collection of simulation models in the Modelica language, representing a variety of mechanical system designs and envelope conditions related to LBL's FLEXLAB facility. The collection of models also features multiple controls sequences that can be simulated with the FLEXLAB model to simulate different demand flexibility scenarios. Additionally, this package will feature datasets from 3 experimental tests, used for calibration, validation and comparison against the Modelica models, this includes weather data that can be used to replicate different scenarios in simulation across the same weather conditions experienced in real experiments. Given FLEXLAB high level of instrumentation and available data, the models are calibrated across multiple measurement points, and thus results from the extension of this model to other climate zones or control sequences, would provide high level of confidence.

Huang, Weiping↗

Effect of disturbances on MED-TVC plant characteristics: Dynamic modeling and simulation

This study describes a detailed steady-state analysis and transient operation of a multi-effect desalination plant with a thermal vapor compressor unit (MED-TVC). The model developed is based on solving the basic conservation equations of mass, energy, and salt. It can predict steady-state operation and transient behavior when the plant is subject to abrupt changes in the main input parameters including motive steam flow rate (m ms ), cooling seawater flow rate (m cw ) and cooling seawater temperature (T cw ). These sudden changes simulate the conditions a MED-TVC plant may be exposed to. Both the steady-state and dynamic aspects of the model are validated against published experimental data. Monitoring the variation of four state variables (vapor temperature, brine level, brine and vapor flow rates) in each effect shows that disturbances in the motive steam flow rate and the cooling seawater temperature have a significant effect on the plant performance. On the other hand, the disturbance in the cooling seawater mass flow rate has only a moderate effect. For all applied disturbances, the change in the brine level is the slowest compared to the changes in vapor temperature, and brine and vapor flow rates. Furthermore, the effect of changing seawater salinity shows a slight effect on the total distillate production and the specific heat transfer area.

42 ENGINEERING↗

Prolongation of SMAP to Spatiotemporally Seamless Coverage of Continental U.S. Using a Deep Learning Neural Network

The Soil Moisture Active Passive (SMAP) mission has delivered valuable sensing of surface soil moisture since 2015. However, it has a short time span and irregular revisit schedules. Utilizing a state-of-the-art time series deep learning neural network, Long Short-Term Memory (LSTM), we created a system that predicts SMAP level-3 moisture product with atmospheric forcings, model-simulated moisture, and static physiographic attributes as inputs. The system removes most of the bias with model simulations and improves predicted moisture climatology, achieving small test root-mean-square errors (<0.035) and high-correlation coefficients >0.87 for over 75% of Continental United States, including the forested southeast. As the first application of LSTM in hydrology, we show the proposed network avoids overfitting and is robust for both temporal and spatial extrapolation tests. LSTM generalizes well across regions with distinct climates and environmental settings. With high fidelity to SMAP, LSTM shows great potential for hindcasting, data assimilation, and weather forecasting.

54 ENVIRONMENTAL SCIENCES↗

Precipitation Characteristics and Future Changes Over the Southern Slope of Tibetan Plateau Simulated by a High-Resolution Global Nonhydrostatic Model

Current climate model resolution cannot accurately describe the complex topography over the Tibetan Plateau, which limits our understanding of past and future precipitation over this region. This study investigates the daily precipitation characteristics and its future changes over the Tibetan Plateau, especially for the southern slope of Tibetan Plateau (SSTP), by 14-km Nonhydrostatic ICosahedral Atmospheric Model (NICAM) with explicitly calculated convection for historical and future 30-years period. By comparing with the satellite Global Precipitation Measurement (GPM), NICAM well reproduces the historical precipitation spatial pattern, seasonal cycle, and the extreme precipitation belt over SSTP, but overestimates precipitation amount by~35%. It is found that heavy precipitation probability decreases as elevation becomes higher while the light precipitation generally shows the opposite. For the precipitation changes during June to September from 1979–2008 to 2075–2104, NICAM predicts that mean precipitation will decrease over low-level SSTP but increase over high-level SSTP. Nonprecipitation and heavy precipitation probability will increase while light precipitation probability will decrease in the future over SSTP. The extreme precipitation probability and intensity will increase ~50%/°C and ~8%/°C over SSTP, and this increase is more obvious as elevation becomes higher. The robust increase of extreme precipitation along the SSTP topography is unique and has not been identified by the climate model simulations before. The strong meridional gradient of specific humidity over SSTP is found to be further enhanced under global warming, and this gradient enhancement is suggested to be responsible for the increase in the extreme precipitation over SSTP.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measurement and Monte Carlo simulation of 3D neutrons distribution in a radiotherapy bunker - 25411

Different methods exist to verify bunkers design for radiation therapy medical facilities; analytical methods based on simplified equations and Monte Carlo methods. One of the main purposes of this work is to present the advantages of using Monte Carlo simulation to verify radiotherapy bunker shielding design. This methodology is more accurate and characterizes not only the fluence and dose three-dimensionally, but also the energy spectrum of particles generated by the Linac. The other main purpose of the work is to analyze the importance of the generated photoneutrons when the Linac emits photon beams of high energy (above 8 MeV), since numerous studies have found that they can expose the patient and clinical staff to non-negligible dose. The main novelty introduced by this work, is the creation of more realistic simulation models to represent the radiotherapy facility by using CAD and meshes technologies. Results obtained using these bunker simulation models have been validated experimentally at the Hospital Universitari i Politecnic La Fe de Valencia facilities using two different neutron detectors; the neutron meter LB 6411 (designed between Berthold and Karlsruhe Research Center) and neutron detector model 42 - 41L (Ludlum, Prescila). Results obtained with Monte Carlo and those measured experimentally fit correctly, validating this analysis methodology. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Modeling and simulation of Ultrasonic Guided Waves propagation in the fluid-structure domain by a monolithic approach

The present study introduces the coupled multiphysics model as part of the structural health monitoring (SHM) system. In particular, the Ultrasonic Guided Waves (UGWs) propagation is tracked in order to identify the damage to the structure. For this purpose, a multiphysics mathematical model is proposed. The model constitutes a monolithically coupled system of acoustic and elastic wave equations (WpFSI problem) where the wave signal displacement measurements are analyzed as the UGWs propagates in the solid, fluid and their interface. The ultimate goal of this paper is to explore and develop the efficient numerical solution of the WpFSI problem using the finite element method. A detailed description of the modeling framework and conditions that facilitate the coupling is provided. To couple the displacement-based acoustic wave equations for the isothermal fluid with elastic wave equations for the Saint Venant-Kirchhoff material model, the present study uses a monolithic solution algorithm. In particular, at each time step, wave equations are transformed to a fixed reference configuration via the arbitrary Lagrangian-Eulerian (ALE) mapping and automatically adopt the boundary conditions from the previous time step. The implementation is accomplished via the finite element library deal.II (Goll et al., 2017) based software toolbox DOpElib (Bangerth et al., 2007) which provides modularized higher level algorithms. Beyond SHM systems, the model is relevant for problems from biomechanics and biomedicine, vibromechanics, poroelasticity as well as subsurface and porous media flow. Here, the model developed here serves as a first step towards the on-line SHM system modeling, where structural dynamics are accounted for.

42 ENGINEERING↗

Appending High-Resolution Elevation Data to GPS Speed Traces for Vehicle Energy Modeling and Simulation

Accurate and reliable global positioning system (GPS)-based vehicle use data are highly valuable for many transportation, analysis, and automotive considerations. Model-based design, real-world fuel economy analysis, and the growing field of autonomous and connected technologies (including predictive powertrain control and self-driving cars) all have a vested interest in high-fidelity estimation of powertrain loads and vehicle usage profiles. Unfortunately, road grade can be a difficult property to extract from GPS data with consistency. In this report, we present a methodology for appending high-resolution elevation data to GPS speed traces via a static digital elevation model. Anomalous data points in the digital elevation model are addressed during a filtration/smoothing routine, resulting in an elevation profile that can be used to calculate road grade. This process is evaluated against a large, commercially available height/slope dataset from the Navteq/Nokia/HERE Advanced Driver Assistance Systems product. Results will show good agreement with the Advanced Driver Assistance Systems data in the ability to estimate road grade between any two consecutive points in the contiguous United States.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Shear-induced aggregation or disaggregation in edible oils: Models, computer simulation, and USAXS measurements

The effects of shear upon the aggregation of solid objects formed from solid triacylglycerols (TAGs) immersed in liquid TAG oils were modeled using Dissipative Particle Dynamics (DPD) and the predictions compared to experimental data using Ultra-Small Angle X-ray Scattering (USAXS). The solid components were represented by spheres interacting via attractive van der Waals forces and short range repulsive forces. A velocity was applied to the liquid particles nearest to the boundary, and Lees-Edwards boundary conditions were used to transmit this motion to non-boundary layers via dissipative interactions. The shear was created through the dissipative forces acting between liquid particles. Translational diffusion was simulated, and the Stokes-Einstein equation was used to relate DPD length and time scales to SI units for comparison with USAXS results. The SI values depended on how large the spherical particles were (250 nm vs. 25 nm). Aggregation was studied by (a) computing the Structure Function and (b) quantifying the number of pairs of solid spheres formed. Solid aggregation was found to be enhanced by low shear rates. As the shear rate was increased, a transition shear region was manifested in which aggregation was inhibited and shear banding was observed. Aggregation was inhibited, and eventually eliminated, by further increases in the shear rate. The magnitude of the transition region shear, $\dot{γ}_t$, depended on the size of the solid particles, which was confirmed experimentally.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Heat transfer simulation and thermal measurements of microfabricated x-ray transparent heater stages

A microfabricated amorphous silicon nitride membrane-based nanocalorimeter is proposed to be suitable for an x-ray transparent sample platform with low power heating and built-in temperature sensing. In this work, thermal characterization in both air and vacuum are analyzed experimentally and via simulation. Infrared microscopy and thermoreflectance microscopy are used for thermal imaging of the sample area in air. While a reasonably large isothermal area is found on the sample area, the temperature homogeneity of the entire sample area is low, limiting use of the device as a heater stage in air or other gases. A simulation model that includes conduction, as well as radiation and convection heat loss, is presented with radiation and convection parameters determined experimentally. Simulated temperature distributions show that the homogeneity can be improved by using a thicker thermal conduction layer or reducing the pressure of the gas in the environment but neither are good solutions for the proposed use. A new simple design that has improved temperature homogeneity and a larger isothermal area while maintaining a thin thermal conduction layer is proposed and fabricated. This new design enables applications in transmission x-ray microscopes and spectroscopy setups at atmospheric pressure.

47 OTHER INSTRUMENTATION↗

Review of strongly-coupled composite dark matter models and lattice simulations

We review models of new physics in which dark matter arises as a composite bound state from a confining strongly-coupled non-Abelian gauge theory. We discuss several qualitatively distinct classes of composite candidates, including dark mesons, dark baryons, and dark glueballs. Here, we highlight some of the promising strategies for direct detection, especially through dark moments, using the symmetries and properties of the composite description to identify the operators that dominate the interactions of dark matter with matter, as well as dark matter self-interactions. We briefly discuss the implications of these theories at colliders, especially the (potentially novel) phenomenology of dark mesons in various regimes of the models. Throughout the review, we highlight the use of lattice calculations in the study of these strongly-coupled theories, to obtain precise quantitative predictions and new insights into the dynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Molecular modeling and simulation studies of recombinant laccase from Yersinia enterocolitica suggests significant role in the biotransformation of non-steroidal anti-inflammatory drugs

The YacK gene from Yersinia enterocolitica strain 7, cloned in pET28a vector and expressed in Escherichia coli BL21 (DE3), showed laccase activity when oxidized with 2,2′-azino-bis(3-ethylbenzthiazoline-6-sulfonic acid) (ABTS) and guaiacol. The recombinant laccase protein was purified and characterized biochemically with a molecular mass of ≈58 KDa on SDS-PAGE and showed positive zymogram with ABTS. The protein was highly robust with optimum pH 9.0 and stable at 70 °C upto 12 h with residual activity of 70%. Kinetic constants, K{sub m} values, for ABTS and guaiacol were 675 μM and 2070 μM, respectively, with corresponding Vmax values of 0.125 μmol/ml/min and 6500 μmol/ml/min. It also possess antioxidative property against BSA and Cu{sup 2+}/H{sub 2}O{sub 2} model system. Constant pH MD simulation studies at different protonation states of the system showed ABTS to be most stable at acidic pH, whereas, diclofenac at neutral pH. Interestingly, aspirin drifted out of the binding pocket at acidic and neutral pH, but showed stable binding at alkaline pH. The biotransformation of diclofenac and aspirin by laccase also corroborated the in silico results. This is the first report on biotransformation of non-steroidal anti-inflammatory drugs (NSAIDs) using recombinant laccase from gut bacteria, supported by in silico simulation studies. - Highlights: • Laccase from Yersinia enterocolitica strain 7 was expressed in Escherichia coli BL21 (DE3). • Recombinant laccase was found to be thermostable and alkali tolerant. • The in silico and experimental studied proves the biotransformation of NSAIDs. • Laccase binds to ligands differentially under different protonation state. • Laccase also possesses free radical scavenging property.

60 APPLIED LIFE SCIENCES↗

Deciphering the R Ψ m correlator in search for the chiral magnetic effect in relativistic heavy ion collisions

The chiral magnetic effect (CME) is extensively studied in heavy-ion collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). An azimuthal correlator called R Ψ m was proposed to measure the CME. By observing the same R Ψ 2 and R Ψ 3 (convex) distributions from A Multi-Phase Transport (ampt) model, by contrasting data and model as well as large and small systems and by event shape engineering (ESE), a recent preprint (arXiv:2006.04251v1) from STAR suggests that the R Ψ m observable is sensitive to the CME signal and relatively insensitive to backgrounds, and their Au + Au data are inconsistent with known background contributions. Our ampt results contradict those from a recent preprint by STAR. Our toy model simulations demonstrate that R Ψ 2 is sensitive to both the CME signal and physics backgrounds. Toy model simulations of isobar collisions show similar centrality dependence and magnitudes for the relative R Ψ 2 strengths as well as the relative Δ γ strengths. We conclude that R Ψ 2 and the inclusive Δ γ are essentially the same.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling and Simulation of High-Frequency Solar Irradiance

As the use of solar power as a source of electricity is increasing, so is the interest in modeling radiation at high temporal resolutions. High-dimensional remote sensing data products depend on cloud cover variability, atmosphere aerosol levels, and other atmospheric parameters. Because of weather fronts and aerosols, it is difficult to quantify solar power variability based on distributed solar networks. The global horizontal irradiance (GHI) component in the National Solar Radiation Database (NSRDB) is available at a 30-min time resolution. The algorithm proposed in this paper produces 1-min-resolution GHI samples at locations where NSRDB data are available. The synthetic irradiance datasets are produced using log-additive non-Gaussian mixture models and resampling techniques. The model is trained over historical data, and predicted values are compared with in situ data. This approach allows for estimating the solar irradiance at subhourly temporal resolutions, while featuring variability for locations where measurements are otherwise not available.

14 SOLAR ENERGY↗

A reduced-model (nSOLT) simulation of neutral recycling effects on plasma turbulence in the divertor region of MAST-U

The 2D scrape-off-layer turbulence code (nSOLT), which includes 1D kinetic neutral–plasma interactions, is applied to study effects of neutral recycling on plasma turbulence for parameters illustrative of the MAST-U divertor region. Neutral recycling is modeled by injecting a fraction of the parallel plasma flux to the divertor back into the simulation domain as a source of Franck–Condon-distributed neutrals. Stationary sources, concentrated at the magnetic separatrix, model plasma streaming into the divertor region from the upstream scrape-off-layer and sustain plasma turbulence absent neutral recycling. Starting from one such no-neutrals equilibrium, we initiate recycling in a numerical experiment designed to diagnose and identify the effects of various neutral–plasma interactions on the divertor plasma, divertor turbulence, and plasma exhaust. The onset of recycling triggers an initial burst of enhanced cross field plasma transport that is quelled by ionization cooling and charge–exchange (CX) friction, with growing neutral pressure, leading to a quiescent, turbulence-free state. Diagnosis of this transient burst reveals that (1) the sudden increase in plasma density due to ionization dominates the onset of the burst, (2) electron cooling due to ionization increases collisionality and disconnects blob filaments from the sheath, and (3) CX friction drives tripole polarization of a blob that can dominate the curvature-driven dipole polarization, leading to the stagnation of blob propagation and reduced radial turbulent transport. It is shown that CX friction is negligible compared to sheath physics in determining equilibrium mean flow shearing rates, for parameters considered herein (specifically a short connection length to the divertor target), while it can significantly reduce interchange-instability growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluating Printability of Buried Native EUV Mask Phase Defects through a Modeling and Simulation Approach

The availability of defect-free masks is considered to be a critical issue for enabling extreme ultraviolet lithography (EUVL) as the next generation technology. Since completely defect-free masks will be hard to achieve, it is essential to have a good understanding of the printability of the native EUV mask defects. In this work, we performed a systematic study of native mask defects to understand the defect printability caused by them. The multilayer growth over native substrate mask blank defects was correlated to the multilayer growth over regular-shaped defects having similar profiles in terms of their width and height. To model the multilayer growth over the defects, a novel level-set multilayer growth model was used that took into account the tool deposition conditions of the Veeco Nexus ion beam deposition tool. The same tool was used for performing the actual deposition of the multilayer stack over the characterized native defects, thus ensuring a fair comparison between the actual multilayer growth over native defects, and modeled multilayer growth over regular-shaped defects. Further, the printability of the characterized native defects was studied with the SEMATECH-Berkeley Actinic Inspection Tool (AIT), an EUV mask-imaging microscope at Lawrence Berkeley National Laboratory (LBNL). Printability of the modeled regular-shaped defects, which were propagated up the multilayer stack using level-set growth model was studied using defect printability simulations implementing the waveguide algorithm. Good comparison was observed between AIT and the simulation results, thus demonstrating that multilayer growth over a defect is primarily a function of a defect’s width and height, irrespective of its shape. This would allow us to predict printability of the arbitrarily-shaped native EUV mask defects in a systematic and robust manner.

36 MATERIALS SCIENCE↗

DEPLOYING FAST CHARGING INFRASTRUCTURE FOR ELECTRIC VEHICLES IN URBAN NETWORKS: AN ACTIVITY-BASED APPROACH

This paper explores an important problem under the domain of network modeling, the optimal configuration of charging infrastructure for electric vehicles (EVs) in urban networks considering EV users' daily activities and charging behavior. This study proposes a charging behavior simulation model considering different initial state of charge (SOC), travel distance, availability of home chargers, and the daily schedule of trips for each traveler. The proposed charging behavior simulation model examines the complete chain of trips for EV users as well as the interdependency of trips traveled by each driver. The problem of finding the optimum charging configuration is then formulated as a mixed-integer nonlinear programming problem that considers the dynamics of travel time and travel distance, the interdependency of trips made by each driver, limited range of EVs, remaining battery capacity for recharging, waiting time in queue, and detour to access a charging station. This problem is solved using a metaheuristic approach for a large-scale case network. A series of examples are presented to demonstrate the model efficacy and explore the impact of energy consumption on the final SOC and the optimum charging infrastructure.

Chain of Trips↗

Predicting Deformation of Parts from HP Printers due to Sintering with Modeling and Simulation (CRADA NFE-18-07363)

This work proposes to leverage Oak Ridge National Laboratory’s (ORNL) sintering simulation software developed at the Manufacturing Demonstration Facility to 1) investigate the deformation of sintered parts and 2) to develop fast surrogate-based optimization methods for HP Inc. This includes the development of ORNL’s software on HP systems, integration and testing with HP part geometries and material properties, the development of surrogate models based on the results of the simulation, and modifications to the software based on feedback from HP. The expected impacts of this work are a drastic reduction in print costs, an increase in printer throughput, and the identification of novel, possibly exotic materials for 3D printing. This work will also decrease the time to market for HP’s print technology by decreasing the amount of resources and type dedicated to physical testing.

42 ENGINEERING↗