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 667 records · Page 37

On the spectral stability of time integration algorithms for a class of constrained dynamics problems

Incomplete field formulations have recently been the subject of intense research because of their potential in coupled analysis of independently modeled substructures, adaptive refinement, domain decomposition, and parallel processing. This paper discusses the design and analysis of time-integration algorithms for these formulations and emphasizes the treatment of their inter-subdomain constraint equations. These constraints are shown to introduce a destabilizing effect in the dynamic system that can be analyzed by investigating the behavior of the time-integration algorithm at infinite and zero frequencies. Three different approaches for constructing penalty-free unconditionally stable second-order accurate solution procedures for this class of hybrid formulations are presented, discussed and illustrated with numerical examples. The theoretical results presented in this paper also apply to a large family of nonlinear multibody dynamics formulations. Some of the algorithms outlined herein are important alternatives to the popular technique consisting of transforming differential/algebraic equations into ordinary differential equations via the introduction of a stabilization term that depends on arbitrary constants and that influences the computed so1ution.

Farhat, Charbel↗

Powder Bed Fusion Laser Beam Metals Additive Manufacturing: Process Monitoring Approaches for Qualification and Certification

The use of in-situ process monitoring is of interest to lower the cost of inspection for the qualification of powder bed fusion laser beam metal (PBF-LB/M) additively manufactured (AM) parts. Precise monitoring of the PBF-LB/M AM build process constitutes a multi-scale and multi-discipline task. There are several significant challenges to the in-situ approach: the synchronization of sensor signals to process steps; the physical interpretation and classification of sensor signals; managing very large datasets; and comparing the inputs with the observed monitoring signals. At NASA Langley Research Center, a configurable architecture additive testbed has been developed to monitor the build process with synchronized sensors. The philosophy and method adopted for the synchronization of the cameras with laser power and position throughout a complex PBF-LB/M AM build will be described. The synchronized in-situ monitoring signals are compared with ex-situ nondestructive inspection, x-ray computed tomography (XCT). Such comparisons permit a better understanding of how the sequential process actions of LPBF-AM can affect build quality. The multi-scale and complex process of printing additively manufactured (AM) parts can have unexpected, but predictable, build conditions that result in material microstructure variability. This presentation will describe an additive manufacturing model-based process metric (AM-PM) computational method that is a fully parallel reduced order modeling approach developed to evaluate the evolution of AM processes. This method couples the known sequence of the AM process with a physically informed nearest neighbors’ calculation to map the conditions of a part-scale build. The result is a map of the build that is derived directly from build files or in-situ process monitoring sensors. The methodology of the approach will be described and mapped to the porosity observed from XCT for a complex PBF-LB/M build. Such comparative results develop understanding of how the sequential process actions can affect the PBF-LB/M AM build quality and microstructure variability.

Laser Powder Bed Fusion↗

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↗

Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site

Numerical simulations are essential tools for understanding the complex hydrologic response of Arctic regions to a warming climate. However, strong coupling among thermal and hydrological processes on the surface and in the subsurface and the significant role that subtle variations in surface topography have in regulating flow direction and surface storage lead to significant uncertainties. Careful model evaluation against field observations is thus important to build confidence.We evaluate the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations from polygonal tundra at the Barrow Environmental Observatory. ATS couples a multiphase, 3D representation of subsurface thermal hydrology with representations of overland nonisothermal flows, snow processes, and surface energy balance. We simulated thermal hydrology of a 3D ice-wedge polygon with geometry that is abstracted but broadly consistent with the surface microtopography at our study site. The simulations were forced by meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of parameters appearing in the soil evaporation model, the 3-year simulations agreed reasonably well with snow depth, summer water table elevations in the polygon center, and high-frequency soil temperature measurements at several depths in the trough, rim, and center of the polygon. Upscaled evaporation is in good agreement with flux tower observations. The simulations were found to be sensitive to parameters in the bare soil evaporation model, snowpack, and the lateral saturated hydraulic conductivity.Timing of fall freeze-up was found to be sensitive to initial snow density, illustrating the importance of including snow aging effects. The study provides new support for an emerging class of integrated surface/subsurface permafrost simulators.

54 ENVIRONMENTAL SCIENCES↗

Modeling and Observational Framework for Diagnosing Local Land-Atmosphere Coupling on Diurnal Time Scales

Land-atmosphere interactions play a critical role in determining the diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture states. The degree of coupling between the land surface and PBL in numerical weather prediction and climate models remains largely unexplored and undiagnosed due to the complex interactions and feedbacks present across a range of scales. Further, uncoupled systems or experiments (e.g., the Project for Intercomparison of Land Parameterization Schemes, PILPS) may lead to inaccurate water and energy cycle process understanding by neglecting feedback processes such as PBL-top entrainment. In this study, a framework for diagnosing local land-atmosphere coupling is presented using a coupled mesoscale model with a suite of PBL and land surface model (LSM) options along with observations during field experiments in the U. S. Southern Great Plains. Specifically, the Weather Research and Forecasting (WRF) model has been coupled to the Land Information System (LIS), which provides a flexible and high-resolution representation and initialization of land surface physics and states. Within this framework, the coupling established by each pairing of the available PBL schemes in WRF with the LSMs in LIS is evaluated in terms of the diurnal temperature and humidity evolution in the mixed layer. The co-evolution of these variables and the convective PBL is sensitive to and, in fact, integrative of the dominant processes that govern the PBL budget, which are synthesized through the use of mixing diagrams. Results show how the sensitivity of land-atmosphere interactions to the specific choice of PBL scheme and LSM varies across surface moisture regimes and can be quantified and evaluated against observations. As such, this methodology provides a potential pathway to study factors controlling local land-atmosphere coupling (LoCo) using the LIS-WRF system, which will serve as a testbed for future experiments to evaluate coupling diagnostics within the community.

Santanello, Joseph A., Jr.↗

Theory of Plasmonic Hot-Carrier Generation and Relaxation

Hot-carrier (HC) generation from (localized) surface plasmon decay has recently attracted much attention due to its promising applications in physical, chemical, materials, and energy science. However, the detailed mechanisms of plasmonic HC generation, relaxation, and trapping are less studied. In this work, we developed and applied a quantum-mechanical model and coupled master equation method to study the generation of HCs from plasmon decay and their following relaxation processes with different mechanisms treated on equal footing. First, a quantum-mechanical model for HC generation is developed. Its connection to existing semiclassical models and time-dependent density functional theory (TDDFT) is discussed. Second, the relaxation and lifetimes of HCs are investigated in the presence of electron–electron and electron–phonon interactions. A GW-like approximation is introduced to account for the electron–electron scattering. The numerical simulations on the Jellium nanoparticles with a size up to 1.6 nm demonstrate the electron–electron scattering and electron–phonon scattering dominate different time scale in the relaxation dynamics. We also generalize the model to study the extraction of HCs to attached molecules. The quantum yield of extracting HCs for other applications is found to be size-dependent. In general, the smaller size of NP improves the quantum yield, which is in agreement with recent experimental measurements. Even though we demonstrate this newly developed theoretical formalism with Jellium model, the theory applies to any other atomistic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of numerical optimization techniques to control system design for nonlinear dynamic models of aircraft

Control system design for general nonlinear flight dynamic models is considered through numerical simulation. The design is accomplished through a numerical optimizer coupled with analysis of flight dynamic equations. The general flight dynamic equations are numerically integrated and dynamic characteristics are then identified from the dynamic response. The design variables are determined iteratively by the optimizer to optimize a prescribed objective function which is related to desired dynamic characteristics. Generality of the method allows nonlinear effects to aerodynamics and dynamic coupling to be considered in the design process. To demonstrate the method, nonlinear simulation models for an F-5A and an F-16 configurations are used to design dampers to satisfy specifications on flying qualities and control systems to prevent departure. The results indicate that the present method is simple in formulation and effective in satisfying the design objectives.

Lan, C. Edward↗

Soil Evaporation Stress Determines Soil Moisture‐Evapotranspiration Coupling Strength in Land Surface Modeling

Model‐based estimates of soil moisture (SM)‐evapotranspiration (ET) coupling strength(ρ) vary widely and are prone to bias. Here we apply numerical modeling and remote sensing to identify the process‐level source of modeledρbias with the goal of improving the fidelity of current Earth system models. Results illustrate that modeledρis most strongly determined by soil evaporation (E) stress, and (generally positive)ρmodeling bias is attributable to the oversimplification of soil texture impacts on E stress. Based on new remotely sensed estimates ofρ, we demonstrate that removingρbias via a single optimized E stress parameter leads to improved ET accuracy and resolves a well‐known modeling bias in the partitioning of ET into E and T. As such, we highlight the importance of the stress function relating E and SM and its central role in regulating land‐atmosphere coupling processes impacting local climate.

Jianzhi Dong↗

Fuel Matrix Degradation Model Development Update: Alloy Corrosion Rates and Hydrogen Generation

This work is being performed as part of the DOE NE Spent Fuel and Waste Science and Technology Campaign, Argillite and Crystalline Rock R&D work packages: SF-19AN01030101 and SF-19AN01030201. This document meets the milestone M4SF-19AN010301013 for Argillite R&D and the milestone M4SF-19AN010302013 for Crystalline R&D. The fuel matrix degradation (FMD) model calculates the degradation rate of spent UO2 fuel based on fundamental electrochemical theory and thermodynamics. The model provides the quantitative basis for radionuclide source term predictions in the generic disposal system analysis (GDSA) performance assessment (PA) model. The FMD model has been implemented in a manner that facilitates its integration with the GDSA-PA code and a preliminary, integrated FMD-GDSA model has been successfully tested (Jerden et al., 2017). The specific focus of on-going work is to accurately and quantitatively represent the generation of H2 in a breached waste package and its effect on the spent fuel degradation rate so that this key process can be accurately represented in PA models. It has been shown experimentally that millimolar concentrations of dissolved H2 in contact with spent fuel will inhibit oxidative dissolution and decrease the fuel degradation rates by as much as three orders of magnitude (e.g., Röllin et al., 2001, Ollila, 2008). This will lead to a significant decrease in the radionuclide source term values used in repository PA models. Sensitivity studies using the FMD model have shown that the dissolved H2 concentration is the dominant environmental variable affecting the UO2 spent fuel dissolution rate as long as a source of H2 persists in the breached waste package (Jerden et al., 2015). Therefore, on-going experimental and modeling efforts focus on H2 production, accumulation and interfacial reactions within a breached waste package. The anoxic corrosion of metallic engineering materials within a breached waste package will be the main source of H2 in crystalline and argillite rock repository systems. The most important H2-producing alloys present within a waste package include stainless steel, carbon steel, and aluminum alloys and the Zircaloy fuel cladding. A steel corrosion module was added to the FMD model to account for these H2 sources and couple the H2 generation rate with fuel degradation processes (Jerden et al., 2017). As discussed below, there is a need for experimental data from electrochemical corrosion tests to quantify the H2 generation rates due to the corrosion of these waste package alloys as a function of key environmental variables (Eh, pH, groundwater chemistry). This report presents results from series of electrochemical corrosion tests with materials present in waste packages that provide information needed to parameterize and validate the FMD model. The alloys tested include 316L stainless steel (316SS), AISI 4320 carbon steel (C-steel) and Zircaloy-4 fuel cladding (unirradiated). Potentiodynamic polarization and potentiostatic polarization tests were performed in buffer solutions at pH 4, 7 and 10 with NaCl concentrations of 0, 4 and 100 millimolar (mM). The potentiostatic tests were performed at voltages varying from 0.5 volts down to -0.66 volts, vs the standard hydrogen electrode (SHE), and all tests were performed at laboratory ambient temperatures (~22 °C). The corrosion rates measured during the electrochemical tests were used as direct input to a new FMD/in-package chemistry model that is being developed to provide accurate long-term spent fuel degradation rates as chemical conditions (Eh, pH, speciation) evolve within a breached waste package. This new model was developed by coupling the FMD model with the reactive transport code X1t, which is a module within the Geochemist’s Workbench (GWB) reaction path modeling code (Jerden et al., 2019). The reactive transport model consists of a 1D domain discretized with 21 reaction diffusion cells and includes a single waste packag

Jerden, James↗

Phase-field model of precipitation processes with coherency loss

A phase-field model is proposed to simulate coherency loss coupled with microstructure evolution. A special field variable is employed to describe the degree of coherency loss of each particle and its evolution is governed by a Ginzburg-Landau type kinetic equation. For the sake of computational efficiency, a flood-fill algorithm is introduced that can drastically reduce the required number of field variables, which allows the model to efficiently simulate a large number of particles sufficient for characterizing their statistical features during Ostwald ripening. The model can incorporate size dependence of coherency loss, metastability of coherent particles, and effectively incorporate the underlying mechanisms of coherency loss by introducing a so-called differential energy criterion. The model is applied to simulate coarsening of Al 3 Sc precipitates in aluminum alloy and comprehensively compared with experiments. Our results clearly show how the particle size distribution is changed during coherency loss and affects the coarsening rate.

36 MATERIALS SCIENCE↗

Environment-assisted non-Markovian speedup dynamics control

For the damped Jaynes–Cummings model in the weak system–environment coupling regime, a scheme to manipulate the non-Markovian dynamics has been previously proposed in the controllable environment. We here mainly focus on the two-Lorentzian model and the band-gap model, which are composed by adding a new Lorentzian spectrum to the original single-Lorentzian spectrum. By manipulating the parameters of the new added Lorentzian spectrum, i.e., changing the weights of the two Lorentzians, and tuning the spectral width of the new second Lorentzian, the dynamical crossover from Markovian to non-Markovian is achieved in the weak system–environment coupling regime. In addition, according to the controllable non-Markovian dynamics behavior, the roles of the parameters of the environment on the speedup evolution of the system are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Pore Resolved Simulations of Joule Heating in Fibrous Media using an Embedded Boundary Method

Joule heating has been regarded as an energy-efficient and sustainable method for heating materials and gases at large scales. The modeling of local temperature effects at pore-resolved scales for such systems, however, has been difficult to achieve due to challenges in coupling thermo-chemical processes in complex porous media and in large representative volume elements (RVEs). To this end, we developed an electro-thermal model at the pore scale to study Joule heating effects in large heterogeneous systems with different microstructures. This was achieved using the level set method to implicitly delineate distinct regions within the domain, and an embedded boundary method to facilitate heat exchange across the fluid-solid interface. Moreover, we applied this method to investigate unsteady non-linear electro-thermal effects in non-woven fibrous graphite conductors for RVEs with characteristic lengths of 2 mm, with different fiber orientations, porosity (80% – 90%) and fiber diameters (10 – 20µm). The coupled equations were solved numerically and they produced peak temperatures greater than 2000 K resulting in heating rates as high as 80,000 K/s. Moreover, the results depended strongly on the microstructure of the fiber skeleton and current density. Geometries with large fibers (∼ 20µm) had the highest average and peak temperatures with the mean temperature increasing by 3.9 % while the peak temperature increased by 9.9 %. Anisotropic domains on the other hand had the lowest mean and peak temperatures with peak and mean temperatures of 2293 K and 1437.7K respectively representing a corresponding 12.1% and 5.1% drop in the temperatures. An increase in porosity from 80% to 90%, however, led to an increase in the peak temperature by 5.1%.

Joule heating↗

Simulation of the Onset of the Southeast Asian Monsoon During 1997 and 1998: The Impact of Surface Processes

The onset of the southeast Asian monsoon during 1997 and 1998 was simulated with a coupled mesoscale atmospheric model (MM5) and a detailed land surface model. The rainfall results from the simulations were compared with observed satellite data fiom the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The simulation with the land surface model captured basic signatures of the monsoon onset processes and associated rainfall statistics. The sensitivity tests indicated that land surface processes had a greater impact on the simulated rainfall results than that of a small sea surface temperature change during the onset period. In both the 1997 and 1998 cases, the simulations were significantly improved by including the land surface processes. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation; the southwest low-level flow over the Indo- China peninsula and the northern cold front intrusion from southern China. The surface sensible and latent heat exchange between the land and atmosphere modified the lowlevel temperature distribution and gradient, and therefore the low-level. The more realistic forcing of the sensible and latent heat from the detailed land surface model improved the monsoon rainfall and associated wind simulation.

Wang, Yansen↗

Simulation of the Onset of the Southeast Asian Monsoon During 1997 and 1998: The Impact of Surface Processes

The onset of the southeast Asian monsoon during 1997 and 1998 was simulated with a coupled mesoscale atmospheric model (MM5) and a detailed land surface model. The rainfall results from the simulations were compared with observed satellite data from the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The simulation with the land surface model captured basic signatures of the monsoon onset processes and associated rainfall statistics. The sensitivity tests indicated that land surface processes had a greater impact on the simulated rainfall results than that of a small sea surface temperature change during the onset period. In both the 1997 and 1998 cases, the simulations were significantly improved by including the land surface processes. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation; the southwest low-level flow over the Indo-China peninsula and the northern cold front intrusion from southern China. The surface sensible and latent heat exchange between the land and atmosphere modified the low-level temperature distribution and gradient, and therefore the low-level. The more realistic forcing of the sensible and latent heat from the detailed land surface model improved the monsoon rainfall and associated wind simulation. The model results will be compared to the simulation of the 6-7 May 2000 Missouri flash flood event. In addition, the impact of model initialization and land surface treatment on timing, intensity, and location of extreme precipitation will be examined.

Tao, W.-K.↗

Simulation of the Onset of the Southeast Asian Monsoon during 1997 and 1998: The Impact of Surface Processes

The onset of the southeast Asian monsoon during 1997 and 1998 was simulated with a coupled mesoscale atmospheric model (MM5) and a detailed land surface model. The rainfall results from the simulations were compared with observed satellite data from the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The simulation with the land surface model captured basic signatures of the monsoon onset processes and associated rainfall statistics. The sensitivity tests indicated that land surface processes had a greater impact on the simulated rainfall results than that of a small sea surface temperature change during the onset period. In both the 1997 and 1998 cases, the simulations were significantly improved by including the land surface processes. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation; the southwest low-level flow over the Indo-China peninsula and the northern cold front intrusion from southern China. The surface sensible and latent heat exchange between the land and atmosphere modified the low-level temperature distribution and gradient, and therefore the low-level. The more realistic forcing of the sensible and latent heat from the detailed land surface model improved the monsoon rainfall and associated wind simulation. The model results will be compared to the simulation of the 6-7 May 2000 Missouri flash flood event. In addition, the impact of model initialization and land surface treatment on timing, intensity, and location of extreme precipitation will be examined.

Tao, W.-K.↗

Self-Consistent Magnetosphere-Ionosphere Coupling and Associated Plasma Energization Processes

Magnetosphere-Ionosphere (MI) coupling and associated with this process electron and ion energization processes have interested scientists for decades and, in spite of experimental and theoretical research efforts, are still ones of the least well known dynamic processes in space plasma physics. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is large scale ring current (RC) electrodynamic coupling that includes calculation of the magnetospheric electric field that is consistent with the ring current (RC) distribution. A general scheme for numerical simulation of such large-scale magnetosphere-ionosphere coupling processes has been presented earlier in many works. The mathematical formulation of these models are based on "modified frozen-in flux theorem" for an ensemble of adiabatically drifting particles in the magnetosphere. By tracking the flow of particles through the inner magnetosphere, the bounce-averaged phase space density of the hot ions and electrons can be reconstructed and the magnetospheric electric field can be calculated such that it is consistent with the particle distribution in the magnetosphere. The new a self-consistent ring current model has been developed that couples electron and ion magnetospheric dynamics with calculation of electric field. Two new features were taken into account in addition to the RC ions, we solve an electron kinetic equation in our model, self-consistently including these results in the solution. Second, using different analytical relationships, we calculate the height integrated ionospheric conductances as the function of precipitated high energy magnetospheric electrons and ions as produced by our model. This results in fundamental changes to the electric potential pattern in the inner magnetosphere, with a smaller Alfven boundary than previous potential formulations would predict but one consistent with recent satellite observations. This leads to deeper penetration of the plasma sheet ions and electrons into the inner magnetosphere and more effective ring current ions and electron energization.

Khazanov, G. V.↗

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A comprehensive modeling method for proton exchange membrane electrolyzer development

Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.

08 HYDROGEN↗