Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled process modeling”

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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DiscoFluxM version 1.x

DiscoFluxM is software for modeling the deformation behavior of crystalline materials within the third-party finite-element based software, MOOSE. The software enables full-field crystal plasticity finite element simulations of deformation processes of metals and other crystalline materials. DiscoFluxM is an implementation of the DiscoFlux theoretical and algorithmic framework for coupling single-crystal plasticity with dislocation transport. The continuum mesoscale modeling framework is founded upon a nonlocal dislocation-density based crystal plasticity theory. The nonlocal theory couples continuum dislocation transport with an otherwise spatially-local single-crystal model employing nonlinear thermoelasticity and crystallographic plasticity. Dislocation transport is modeled by enforcing dislocation conservation at a slip-system level through the solution of advection-diffusion equations. The configuration of geometrically necessary dislocation density gives rise to a back-stress that inhibits or accentuates the flow of dislocations.

Luscher, Darby↗
NASA NTRS2005

Applications and Improvement of a Coupled, Global and Cloud-Resolving Modeling System

Recently Grabowski (2001) and Khairoutdinov and Randall (2001) have proposed the use of 2D CFWs as a "super parameterization" [or multi-scale modeling framework (MMF)] to represent cloud processes within atmospheric general circulation models (GCMs). In the MMF, a fine-resolution 2D CRM takes the place of the single-column parameterization used in conventional GCMs. A prototype Goddard MMF based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM) is now being developed. The prototype includes the fvGCM run at 2.50 x 20 horizontal resolution with 32 vertical layers from the surface to 1 mb and the 2D (x-z) GCE using 64 horizontal and 32 vertical grid points with 4 km horizontal resolution and a cyclic lateral boundary. The time step for the 2D GCE would be 15 seconds, and the fvGCM-GCE coupling frequency would be 30 minutes (i.e. the fvGCM physical time step). We have successfully developed an fvGCM-GCE coupler for this prototype. Because the vertical coordinate of the fvGCM (a terrain-following floating Lagrangian coordinate) is different from that of the GCE (a z coordinate), vertical interpolations between the two coordinates are needed in the coupler. In interpolating fields from the GCE to fvGCM, we use an existing fvGCM finite- volume piecewise parabolic mapping (PPM) algorithm, which conserves the mass, momentum, and total energy. A new finite-volume PPM algorithm, which conserves the mass, momentum and moist static energy in the z coordinate, is being developed for interpolating fields from the fvGCM to the GCE. In the meeting, we will discuss the major differences between the two MMFs (i.e., the CSU MMF and the Goddard MMF). We will also present performance and critical issues related to the MMFs. In addition, we will present multi-dimensional cloud datasets (i.e., a cloud data library) generated by the Goddard MMF that will be provided to the global modeling community to help improve the representation and performance of moist processes in climate models and to improve our understanding of cloud processes globally (the software tools needed to produce cloud statistics and to identify various types of clouds and cloud systems from both high-resolution satellite and model data will be also presented).

Tao, W.-K.↗
DOE OSTI2021

Incorporation of Thermal Hydraulic Models for Thermal Power Dispatch into a PWR Power Plant Simulator

This report describes the development, modeling, and results of a generic pressurized water reactor power plant simulator that incorporates coupled electrical and thermal power dispatch to an industrial process located approximately one kilometer from the nuclear power plant. The simulator is a commercial PWR simulator that has been modified to include thermal power dispatch as described in past milestone reports [ , ]. The commercial PWR simulator is a generic simulator available from GSE SYSTEMS® (Sykesville, MD, USA) that is built using RELAP5-HDTM Real-Time Solution and in-house software developed by GSE Systems. This generic PWR (GPWR) simulator performs real-time simulation of the complete power plant from the reactor neutronics to the electricity generation and distribution. All primary, secondary, and auxiliary systems are modeled including all control logic in order to provide the most accurate representation of actual nuclear power plant (NPP) operation, and the simulator results have been rigorously verified by an actual NPP operating at approximately 1 GWe. This report is a continuation of worked performed in previous years, and supplemental information from previous reports is included in the appendix for reference.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

Evaluation of anticipatory signal to steam generator pressure control program for 700 MWe Indian pressurized heavy water reactor

700 MWe Indian Pressurized Heavy Water Reactor (IPHWR) is horizontal channel type reactor with partial boiling at channel outlet. Due to boiling, it has a large volume of vapor present in the primary loops. It has two primary loops connected with the help of pressurizer surge line. The pressurizer has a large capacity and is partly filled by liquid and partly by vapor. Large vapor volume improves compressibility of the system. During turbine trip or load rejection, pressure builds up in Steam Generator (SG). This leads to pressurization of Primary Heat Transport System (PHTS). To control pressurization of SG and PHTS, around 70% of the steam generated in SG is dumped into the condenser by opening Condenser Steam Dump Valves (CSDVs) and rest of the steam is released to the atmosphere by opening Atmospheric Steam Discharge Valves (ASDVs) immediately after sensing the event. This is accomplished by adding anticipatory signal to the output of SG pressure controller. Anticipatory signal is proportional to the thermal power of reactor and the proportionality constant is set so that SG pressure controller's output jacks up to ASDV opening range when operating at 100% FP. To simulate this behavior for 700 MWe IPHWR, Primary and secondary heat transport system is modeled. SG pressure control and other process control program have also been modeled to capture overall plant dynamics. Analysis has been carried out with 3-D neutron kinetics coupled thermal hydraulic computer code ATMIKA.T to evaluate the effect of the anticipatory signal on PHT pressure and over all plant dynamics during turbine trip in 700 MWe IPHWR. This paper brings out the results of the analysis with and without considering anticipatory signal in SG pressure control program during turbine trip. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Hydrological analysis in watersheds with a variable-resolution global climate model (VR-CESM)

Traditionally, watershed-scale hydrology is simulated by distributed hydrological models with offline meteorological forcing data, or by regional regional climate models that link atmospheric and land hydrology interactions. Global climate models (GCMs) are rarely used to study watershed-scale hydrology due to the relatively coarse grid resolution, computationally expensive downscaling, and simplified physical processes. Recently, however, watershed-scale hydrology analysis has become possible in GCMs due to the development of variable-resolution GCMs that dynamically couple the hydrological processes between atmosphere and land systems at fine resolutions in selected regions and coarse resolution elsewhere. In this study, we used the variable-resolution Community Earth System Model (VR-CESM) with refined-resolution (14 km) in the western U.S. and eastern China to evaluate watershed-scale hydrology. Furthermore, we compared the historical VR-CESM outputs with gauge measurements and other hydrological models (e.g., National Water Model in the U.S.) and calibrated the subsurface runoff capacities in four mountainous watersheds. An RCP8.5 projection from 2007 to 2050 is used to estimate the impact of changing precipitation and snow climatology on watershed hydrology. We also analyzed the long-term runoff variability and the possibility of extreme runoff events as simulated by the VR-CESM. Although calibration is not possible in larger-scale watersheds, VR-CESM simulates the long-term annual variability of watersheds and provides insights on climate change impacts on hydrology. We conclude that refined-resolution VR-CESM can be used for watershed-scale hydrology analysis to understand water resources and water balance, in addition to traditional watershed-scale hydrological models. It enables hydrological analysis at multiple watersheds in one simulation and can help understand the two-way dynamics between land surface hydrology and atmospheric processes, and is especially practical for projecting climate change impacts. However, it is challenging to apply VR-CESM for hydrologic analysis in regulated watersheds as human factors (e.g., pumping, irrigation, water diversion) have not been fully addressed in VR-CESM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Numerical parameter optimizations of the Ignition and Growth model for a HMX plastic bonded explosive

We present here a novel method for optimizing parameters of the Ignition and Growth (I&G) reactive flow model for high explosives. The I&G model predicts the shock initiation response of explosives subjected to dynamic loading. However, calibrating the model is a time-consuming process because it has several parameters, typically requires multiple data sets, and there are no direct correlations between parameters and data. In this study, we couple a scalable optimization algorithm to simulations of shock initiation experiments in the multi-physics code ALE3D. We develop four I&G model parameter sets for the HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane)-based explosive LX-07 [90 wt. % HMX, 10 wt. % Viton A] based on minimizing the difference between calculations and measurements at the embedded pressure gauges in 1D gun shock experiments. Our study shows that a cost function based on both shock time of arrival and pressure pulse shape demonstrates the best agreement with experimental data.

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

Particle-in-cell Monte Carlo-collision modeling of non-ideal effects in wave-heated dense microplasmas

A computational model for non-ideal plasma effects during the time evolution of a second-stage laser-heated discharge at high pressures is presented. The model extends a classical one-dimensional particle-in-cell Monte Carlo-collision (PIC-MCC) approach coupled with Maxwell's equations for the laser-heating process of a xenon plasma at 300 K temperature and 10–100 bar pressure. Plasma non-ideality resulting from Coulomb coupling at high plasma densities is manifested as a depression in the effective ionization potential of atoms and enhanced collision cross sections. These non-ideal effects are represented using the Ecker–Kröll model in the context of the PIC-MCC approach. We find that full ionization of the plasma is obtained on the picosecond timescale, starting from the skin layer and quickly expanding throughout the domain through an anomalous extension of the skin depth. More critically, we show that the inclusion of the non-ideal plasma effects results in more rapid ionization compared to an ideal plasma, especially at higher pressures. The ionization delay reduction is of the order of a fraction of a picosecond, corresponding to a 16% decrease at 100 bar. As the amplitude of the wave field is lowered, the ionization rate is lowered, making the plasma non-ideality effects more prominent.

Solmaz, Evrim (ORCID:0000000211830881)↗
DOE OSTI2012

Measurement of the parity-violating asymmetry in inclusive electroproduction of π - near the Delta 0 resonance

The parity-violating (PV) asymmetry of inclusive π - production in electron scattering from a liquid deuterium target was measured at backward angles. The measurement was conducted as a part of the G0 experiment, at a beam energy of 360 MeV. The physics process dominating pion production for these kinematics is quasi-free photoproduction off the neutron via the Δ 0 resonance. In the context of heavy-baryon chiral perturbation theory (HBχPT), this asymmetry is related to a low energy constant d Δ - that characterizes the parity-violating γNΔ coupling. Zhu et al. calculated d Δ - in a model benchmarked by the large asymmetries seen in hyperon weak radiative decays, and predicted potentially large asymmetries for this process, ranging from A γ - = -5.2 to +5.2 ppm. The measurement performed in this work leads to A γ - = -0.36 ± 1.06 ± 0.37 ± 0.03 ppm (where sources of statistical, systematic and theoretical uncertainties are included), which would disfavor enchancements considered by Zhu et al. proportional to V ud /V us . The measurement is part of a program of inelastic scattering measurements that were conducted by the G0 experiment, seeking to determine the N-Δ axial transition form-factors using PV electron scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Dynamic Data Reconciliation and Validation of a Dynamic Model for Solvent-Based CO 2 Capture Using Pilot-Plant Data

This paper is focused on the development and validation of a dynamic model for an MEA-based CO 2 capture unit. Starting with a rigorous steady-state process model, the dynamic model is developed in Aspen Plus Dynamics ® using a modified Murphree-efficiency approach and is validated using dynamic data collected from the National Carbon Capture Center (NCCC) in Wilsonville, Alabama. The dynamic test runs were designed with due consideration of process excitation and process nonlinearities constrained by the available time to conduct the test runs and implement the experimental designs in the plant control system. The experimental data include solvent composition and loading. Since the experimental data were found to violate mass and energy balances, filtering algorithms coupled with dynamic data reconciliation techniques were used. Without adjusting any model parameters, the dynamic model satisfactorily predicted dynamic response of CO 2 capture due to step changes in the solvent flowrate, flue gas flowrate, and steam flowrates. Transient studies show that the process gain and time constants can considerably change depending on the direction of the disturbance or manipulated variables, indicating process nonlinearities. Due to a storage tank between the absorber and stripper, the time constant of the full plant can be considerably longer than the time constants of the absorber or stripper individually. Furthermore, impact of the control system on the process efficiency while following an optimal CO 2 capture schedule was also studied.

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

Charge Delocalization and Vibronic Couplings in Quadrupolar Squaraine Dyes

Squaraines are prototypical quadrupolar charge-transfer chromophores that have recently attracted much attention as building blocks for solution-processed photovoltaics, fluorescent probes with large two-photon absorption cross sections, and aggregates with large circular dichroism. Their optical properties are often rationalized in terms of phenomenological essential state models, considering the coupling of two zwitterionic excited states to a neutral ground state. As a result, optical transitions to the lowest S 1 excited state are one-photon allowed, whereas the next higher S 2 state can only be accessed by two-photon transitions. A further implication of these models is a substantial reduction of vibronic coupling to the ubiquitous high-frequency vinyl-stretching modes of organic materials. Here, in this study, we combine time-resolved vibrational spectroscopy, two-dimensional electronic spectroscopy, and quantum-chemical simulations to test and rationalize these predictions for nonaggregated molecules. We find small Huang–Rhys factors below 0.01 for the high-frequency, 1500 cm –1 modes in particular, as well as a noticeable reduction for those of lower frequency modes in general for the electronic S 0 → S 1 transition. The two-photon allowed state S 2 is well separated energetically from S 1 and has weak vibronic signatures as well. Thus, the resulting pronounced concentration of the oscillator strength in a narrow region relevant to the lowest electronic transition makes squaraines and their aggregates exceptionally interesting for strong and ultrastrong coupling of excitons to localized light modes in external resonators with chiral properties that can largely be controlled by the molecular architecture.

36 MATERIALS SCIENCE↗
NASA NTRS1974

Collisional coupling in counterstreaming laser-produced plasmas

The collisional processes which transfer momentum between counterstreaming plasmas are reviewed and applied to the example of a laser-produced plasma expanding into a partially ionized background. Experimental measurements of the dependence of the ion flow field on collisional momentum transfer demonstrate the validity of the simplified treatment of collision processes which have been adopted. A numerical model which simulates the laser-plasma interaction with the background confirms the importance of collisions in previous experimental studies of momentum coupling, and provides some insight into the distinction between collisional and collisionless flow regimes.

Koopman, D. W.↗
DOE OSTI2011

On the Use of Lambert W Function for Modeling of Quantum Resonances

The Lambert W function is used to represent data (bound state energies) for the process of analytical continuation in the coupling constant. It is shown how the presence of the singularity of the W function restricts the accuracy of determination of the branch point. The accuracy of the determination of resonance parameters is also discussed. It is shown that even for moderately accurate input data the process of analytical continuation yields resonance parameters of accuracy surpassing modern experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

Day and night models of the Venus thermosphere

A model atmosphere of Venus for altitudes between 100 and 178 km is presented for the dayside and nightside. Densities of CO2, CO, O, N2, He, and O2 on the dayside, for 0800 and 1600 hours local time, are obtained by simultaneous solution of continuity equations. These equations couple ionospheric and neutral chemistry and the transport processes of molecular and eddy diffusion. Photodissociation and photoionization J coefficients are presented to facilitate the incorporation of chemistry into circulation models of the Venus atmosphere. Midnight densities of CO2 CO, O, N2, He, and N are derived from integration of the continuity equations, subject to specified fluxes. The nightside densities and fluxes are consistent with the observed airglow of NO and O2(1 Delta). The homopause of Venus is located near 133 km on both the dayside and nightside.

Massie, S. T.↗
NASA NTRS1990

Effects of source vs. interactive actuator assumptions on modelling of a smart structure

Embedded or internal actuators in a structure are being utilized to perform vibration attenuation, pointing maneuvers, damage control, and other adaptive tasks. The controls-structures-interaction (CSI) design process depends strongly on the coupled dynamic response of the structure and the actuators. A bond graph approach is used to guide this analytical study to uncover the effects of certain generic actuator types on the closed-loop, initial condition control of a single-mode structural model and an eighteen-state model of an active truss. So-called 'flow' and 'effort' source actuator types were studied and were shown to have significant differences in closed-loop response for ranges of initial conditions and system open-loop natural frequencies.

Robertshaw, Harry H.↗
DOE OSTI2019

Computational Evidence for Kinetically Controlled Radical Coupling during Lignification

Lignin is an alkyl-aromatic biopolymer that, despite its abundance, is underutilized as a renewable feedstock because of its highly complex structure. A method to overcome this challenge that has gained prominence in recent years leverages the plasticity and malleability of lignin biosynthesis to tune lignin structure in planta through genetic approaches. An improved understanding of lignin biosynthesis can thus provide fundamental insights critical for the development of effective tailoring and valorization strategies. Although it is widely accepted that lignin monomers and growing chains are oxidized enzymatically into radicals that then undergo kinetically controlled coupling in planta, direct experimental evidence has been scarce because of the difficulty of exactly replicating in planta lignification conditions. Here, we computationally investigate a set of radical reactions representative of lignin biosynthesis. We demonstrate that, contrary to the notion that radical coupling reactions are usually barrierless and dynamically controlled, the computed activation energies can be qualitatively consistent with key structural observations made empirically for native lignin in a variety of biomass types. We also rationalize the origins of regioselectivity in coupling reactions through structural and activation strain analyses. Our findings lay the groundwork for first-principles lignin structural models and more detailed multiscale simulations of the lignification process.

09 BIOMASS FUELS↗
NASA NTRS2022

Effects of Increasing the Category Resolution of the Sea Ice Thickness Distribution in a Coupled Climate Model on Arctic and Antarctic Sea Ice Mean State

Many modern sea ice models used in global climate models represent the subgrid-scale heterogeneity in sea ice thickness with an ice thickness distribution (ITD), which improves model realism by representing the significant impact of the high spatial heterogeneity of sea ice thickness on thermodynamic and dynamic processes. Most models default to five thickness categories. However, little has been done to explore the effects of the resolution of this distribution (number of categories) on sea-ice feedbacks in a coupled model framework and resulting representation of the sea ice mean state. Here, we explore this using sensitivity experiments in CESM2 with the standard 5 ice thickness categories and 15 ice thickness categories. Increasing the resolution of the ITD in a run with preindustrial climate forcing results in substantially thicker Arctic sea ice year-round. Analyses show that this is a result of the ITD influence on ice strength. With 15 ITD categories, weaker ice occurs for the same average thickness, resulting in a higher fraction of ridged sea ice. In contrast, the higher resolution of thin ice categories results in enhanced heat conduction and bottom growth and leads to only somewhat increased winter Antarctic sea ice volume. The spatial resolution of the ICESat-2 satellite mission provides a new opportunity to compare model outputs with observations of seasonal evolution of the ITD in the Arctic (ICESat-2; 2018–2021). Comparisons highlight significant differences from the ITD modeled with both runs over this period, likely pointing to underlying issues contributing to the representation of average thickness.

Madison M. Smith↗
DOE OSTI2023

Influence of strong Coulomb coupling on diffusion in atmospheric pressure plasmas

Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion–ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions (x i ≲ 10 —6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions (10 —6 ≲ x i ≲ 10 —2 ) there is a transition from ion-neutral to ion–ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion–ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion–ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions (x i ≳ 10 —2 ), strongly correlated ion–ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. Importantly, the results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.

fast neutral gas heating↗
DOE OSTI2025

Coupling to rotational manifolds to improve gas-phase pump–probe spectroscopic models

The physical picture of gas-phase optical transitions is normally presented as an isolated two-level system balanced by upward and downward processes. Isolated models assume a phenomenological treatment of collisional dephasing but do not strictly account for collisional population exchange with the rotational baths. While this assumption is valid under low-intensity conditions, where excitation is rate-limiting, isolated models can deviate from Beer’s Law at sufficient pressures and monochromatic intensities when both collisional broadening and power broadening are comparable to (or greater than) lifetime broadening, which are not uncommon conditions for cavity enhanced spectroscopies in the mid-IR spectral range. Although this problem has been addressed by rate-equation models for linear absorption measurements, a general treatment for multi-level quantum mechanical models suitable for non-linear absorption measurements (two-photon/two-color/pump–probe) is lacking. Isolated models require physical parameter inputs that disagree with expected values by at least an order of magnitude. These non-physical models undermine the ability to predict non-linear signal strengths under untested conditions and thereby limit the potential to optimize the sensitivity of non-linear spectroscopies and to expand their analytical applications (e.g., new analytes and/or buffer gases, changes in cavity free-spectral-range, changes in intracavity powers or wavelengths, and accurate investigation of physical phenomena). In this study, we derive bath-coupled models for gaseous pump–probe spectroscopy by application of the quantum Lindblad equation and detailed balance. Bath-coupled models are shown to fit data consistently across variations in intensity and agree with all physically expected values.

Cavity ring-down spectroscopy↗