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 829 records · Page 46

A high-resolution large-eddy simulation framework for wildland fire predictions using TensorFlow

Background: Wildfires are becoming more severe, so we need improved tools to predict them over a wide range of conditions and scales. One approach towards this goal entails the use of coupled fire/atmosphere modelling tools. Although significant progress has been made in advancing their physical fidelity, existing tools have not taken full advantage of emerging programming paradigms and computing architectures to enable high-resolution wildfire simulations. Aims: The aim of this study was to present a new framework that enables landscape-scale wildfire simulations with physical representation of combustion at an affordable cost. Methods: We developed a coupled fire/atmosphere simulation framework using TensorFlow, which enables efficient and scalable computations on Tensor Processing Units. Key Results: Simulation results for a prescribed fire were compared with experimental data. Predicted fire behavior and statistical analysis for fire spread rate, scar area, and intermittency showed overall reasonable agreement. Scalability analysis was performed, showing close to linear scaling. Conclusions: While mesh refinement was shown to have less impact on global quantities, such as fire scar area and spread rate, it benefits predictions of intermittent fire behavior, buoyancy-driven dynamics, and small-scale turbulent motion. Implications: This new simulation framework is efficient in capturing both global quantities and unsteady dynamics of wildfires at high spatial resolutions.

54 ENVIRONMENTAL SCIENCES↗

Dimensionality effects on multicomponent ionic transport and surface complexation in porous media

Coulombic interactions between charged species in pore water and at surface/solution interfaces are of pivotal importance for multicomponent reactive transport in porous media. In this study, we investigate the impact of domain dimensionality on electrostatically coupled dispersion and surface-solution reactions during transport of acidic plumes and major ions in porous media. Column and quasi two-dimensional flow-through experiments were performed, with identical silica porous media and under the same advection-dominated conditions. Equal mass fluxes of different electrolyte solutions (i.e., HCl - pH ~ 2.8, NaBr - 100 mM, HCl - pH ~ 2.8 plus NaBr - 100 mM) were continuously injected in the 1-D and 2-D experiments and breakthrough curves of pH and major ions were measured at the outlet of the domains. The presence of pronounced ionic strength gradients in the transverse direction in the 2-D setup caused distinct retardation and transport behaviors of protons and major ions which were not observed in the one-dimensional column experiments. Furthermore, in the cases of salt electrolytes injection, considerably enhanced release of H + (>61%) from the quartz surface was observed in the multidimensional system compared to the one-dimensional setup. Reactive transport modeling was performed to reproduce the experimental outcomes and to analyse the coupling between transport processes, based on the Nernst-Planck formulation of diffusive/dispersive fluxes and on surface complexation reactions at the solid-solution interface. Electrostatic interactions between Na + , B – , and H + , and deprotonation of the quartz surface upon the formation of sodium outer-sphere complexes, are the primary controllers of the spatial and temporal features displayed by the pH and major ions measurements. The reactive transport simulations allowed us to interpret the experimental observations, to visualize the distribution and spatio-temporal evolution of dissolved and solid species, to identify a spatially heterogeneous zonation of Coulombic interactions with distinct behavior at the fringe and core of the injected plumes in the multidimensional setup, and to quantify the different components of the Nerst-Planck fluxes of the charged solutes. This study demonstrates that the domain dimensionality directly affects electrostatic interactions between charged aqueous species in the pore water and surface complexation reactions at the solid-solution interface. The non-trivial effects of dimensionality on multicomponent ionic transport result in a significantly different behavior in 1-D and 2-D systems.

58 GEOSCIENCES↗

Investigation of turbulent processes in magnetospheric boundary layers

A self-consistent non-evolving two dimensional slab model of a viscous low-latitude boundary layer (LLBL) coupled to the ionosphere was developed by Phan, et al., (1989). Numerical results from the model and possible use of observations to determine the model parameters are discussed. The dynamical model developed by Lotko, et al., (1987) was used by Lotko and Shen (1991) to examine dynamical processes relevant to the LLBL with particular application to post-noon auroral shear layers. Initial results from a magnetohydrodynamic study of flank-side mangetopause boundary configuration are described. Effects of compressibility, scalar viscosity, and electrical resistivity are included in the MHD equations.

Lotko, William↗

Use of a Hybrid Edge Node-Centroid Node Approach to Thermal Modeling

A recent proposal submitted for an ESA mission required that models be delivered in ESARAD/ESATAN formats. ThermalDesktop was the preferable analysis code to be used for model development with a conversion done as the final step before delivery. However, due to some differences between the capabilities of the two codes, a unique approach was developed to take advantage of the edge node capability of ThermalDesktop while maintaining the centroid node approach used by ESARAD. In essence, two separate meshes were used: one for conduction and one for radiation. The conduction calculations were eliminated from the radiation surfaces and the capacitance and radiative calculations were eliminated from the conduction surfaces. The resulting conduction surface nodes were coincident with all nodes of the radiation surface and were subsequently merged, while the nodes along the edges remained free. Merging of nodes on the edges of adjacent surfaces provided the conductive links between surfaces. Lastly, all nodes along edges were placed into the subnetwork and the resulting supernetwork included only the nodes associated with radiation surfaces. This approach had both benefits and disadvantages. The use of centroid, surface based radiation reduces the overall size of the radiation network, which is often the most computationally intensive part of the modeling process. Furthermore, using the conduction surfaces and allowing ThermalDesktop to calculate the conduction network can save significant time by not having to manually generate the couplings. Lastly, the resulting GMM/TMM models can be exported to formats which do not support edge nodes. One drawback, however, is the necessity to maintain two sets of surfaces. This requires additional care on the part of the analyst to ensure communication between the conductive and radiative surfaces in the resulting overall network. However, with more frequent use of this technique, the benefits of this approach can far outweigh the additional effort.

Peabody, Hume L.↗

Use of a Hybrid Edge Node-Centroid Node Approach to Thermal Modeling

A recent proposal submitted for an ESA mission required that models be delivered in ESARAD/ESAT AN formats. ThermalDesktop was the preferable analysis code to be used for model development with a conversion done as the final step before delivery. However, due to some differences between the capabilities of the two codes, a unique approach was developed to take advantage of the edge node capability of ThermalDesktop while maintaining the centroid node approach used by ESARAD. In essence, two separate meshes were used: one for conduction and one for radiation. The conduction calculations were eliminated from the radiation surfaces and the capacitance and radiative calculations were eliminated from the conduction surfaces. The resulting conduction surface nodes were coincident with all nodes of the radiation surface and were subsequently merged, while the nodes along the edges remained free. Merging of nodes on the edges of adjacent surfaces provided the conductive links between surfaces. Lastly, all nodes along edges were placed into the subnetwork and the resulting supernetwork included only the nodes associated with radiation surfaces. This approach had both benefits and disadvantages. The use of centroid, surface based radiation reduces the overall size of the radiation network, which is often the most computationally intensive part of the modeling process. Furthermore, using the conduction surfaces and allowing ThermalDesktop to calculate the conduction network can save significant time by not having to manually generate the couplings. Lastly, the resulting GMM/TMM models can be exported to formats which do not support edge nodes. One drawback, however, is the necessity to maintain two sets of surfaces. This requires additional care on the part of the analyst to ensure communication between the conductive and radiative surfaces in the resulting overall network. However, with more frequent use of this technique, the benefits of this approach can far outweigh the additional effort.

Peabody, Hume L.↗

Frequency and Time Domain Nuclear-Electronic Orbital Equation-of-Motion Coupled Cluster Methods: Combination Bands and Electronic-Protonic Double Excitations

The accurate description of excited vibronic states is important for modeling a wide range of photoinduced processes. The nuclear-electronic orbital (NEO) approach, which treats specified protons on the same level as the electrons, can describe excited electronic-protonic states. Herein the multicomponent equation-of-motion coupled cluster with singles and doubles (NEO-EOM-CCSD) method and its time-domain counterpart, TD-NEO-EOM-CCSD, are developed and implemented. The application of these methods to the HCN molecule highlights their capabilities. These methods predict qualitatively reasonable energies and intensities for a combination band corresponding to simultaneous excitation of two vibrational modes, as well as an overtone. These methods also describe states with double excitation character, such as excited electronic-protonic states corresponding to the simultaneous excitation of an electron and a proton. The ability of the NEO-EOM-CCSD method and its time-dependent counterpart to describe combination bands, overtones, and double excitations will enable a wide range of photochemical applications.

Pavosevic, Fabijan↗

Effects of heat exchanger tubes on hydrodynamics and CO 2 capture of a sorbent-based fluidized bed reactor

In virtual design and scale up of pilot-scale carbon capture systems, the coupled reactive multiphase flow problem must be solved to predict the adsorber’s performance and capture efficiency under various operation conditions. This paper focuses on the detailed computational fluid dynamics (CFD) modeling of a pilot-scale fluidized bed adsorber equipped with vertical cooling tubes. Multiphase Flow with Interphase eXchanges (MFiX), an open-source multiphase flow CFD solver, is used for the simulations with custom code to simulate the chemical reactions and filtered models to capture the effect of the unresolved details in the coarser mesh for simulations with reasonable simulations and manageable computational effort. Previously developed two filtered models for horizontal cylinder drag, heat transfer, and reaction kinetics have been modified to derive the 2D filtered models representing vertical cylinders in the coarse-grid CFD simulations. The effects of the heat exchanger configurations (i.e., horizontal or vertical) on the adsorber’s hydrodynamics and CO2 capture performance are then examined. The simulation result subsequently is compared and contrasted with another predicted by a one-dimensional three-region process model.

bubbling bed↗

The evolution from weak to strong geomagnetic activity - An interpretation in terms of deterministic chaos

An analogue of the magnetosphere developed on the basis of Shaw's (1984) dripping faucet model was used to model the mechanisms of the magnetospheric response to energy transfer from the solar wind. It is demonstrated that geomagnetic activity results from nonlinearly coupled physical processes and that the strength and the nature of the coupling changes dramatically as the magnetosphere is driven harder and harder by increasing energy input. Based on initial results obtained from the model, is is suggested that a chaotic transition takes place in the analogue system as the loading rate is increased beyond a critical value. This model is able to explain many of the features in the results of linear prediction filtering techniques.

Baker, D. N.↗

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation, surface-to-surface radiation exchange, and flowfield coupling. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗

Use of satellite data and modeling to asses the influence of stratospheric processes on the troposphere

The research is comprised of the following tasks: use of simple analytical and numerical models of a coupled troposphere-stratosphere system to examine the effects of radiation and ozone on planetary wave dynamics and the tropospheric circulation; use of satellite data obtained from the Nimbus 7 Limb Infrared Monitor of the Stratosphere (LIMS) instrument and Solar Backscattered Ultraviolet (SBUV) experiment, in conjunction with National Meteorological Center (NMC) data, to determine the planetary wave vertical structures, dominant wave spectra, ozone spectra, and time variations in diabatic heating rate; and synthesis of the modeling and observational results to provide a better understanding of the effects that stratospheric processes have on tropospheric dynamics.

Nathan, Terrence R.↗

Development of the GEOS-MIT gcm Atmosphere-Ocean Model for Coupled Data Assimilation

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MIT-GCM (General Circulation Model)) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (approximately 1/8 degrees) configuration of the GEOS atmosphere to a latitude-longitude-cap-1080 (approximately 1/12 degrees) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗

Modeling the processes of soil moisture in regulating microbial and carbon-nitrogen cycling

Soil carbon (C) and nitrogen (N) cycles and their complex responses to hydro-climatic forcing have gained increasing attention. While the temperature effects have been intensively studied, soil moisture response functions (SMRFs) are not well documented for various microbial and enzymatic processes due to the difficulties in directly measuring and differentiating the moisture effects on various processes. In this work we extended our C-only Microbial-ENzyme Decomposition (MEND) model to the C-N coupled MEND model with flexible element stoichiometry. Our model calibration showed good agreement between simulated and observed C:N ratios in soil organic matter and microbial biomass, as well as the ammonium and nitrate concentrations. We show that the selection of SMRFs for specific biogeochemical processes could result in significant differences in model simulated microbial and C-N processes. In particular, it is essential to account for the soil moisture effects on microbial dormancy and resuscitation, as the changes in microbial physiology under favorable or stressful conditions will exert strong controls on soil C and N dynamics. We also advocate the utilization of dynamic (time-variant) data (e.g., litter input, N deposition, soil temperature and moisture), instead of time-invariant data, to drive model simulations and analyses. Dynamic forcing data (particularly dynamic soil moisture) better represent the real-world climate and environmental conditions, which could facilitate more realistic modeling and understanding of soil C and nutrient cycling in a changing world.

54 ENVIRONMENTAL SCIENCES↗

NASA Systems Engineering Research Consortium: Defining the Path to Elegance in Systems

The NASA Systems Engineering Research Consortium was formed at the end of 2010 to study the approaches to producing elegant systems on a consistent basis. This has been a transformative study looking at the engineering and organizational basis of systems engineering. The consortium has engaged in a variety of research topics to determine the path to elegant systems. In the second year of the consortium, a systems engineering framework emerged which structured the approach to systems engineering and guided our research. This led in the third year to set of systems engineering postulates that the consortium is continuing to refine. The consortium has conducted several research projects that have contributed significantly to the understanding of systems engineering. The consortium has surveyed the application of the NASA 17 systems engineering processes, explored the physics and statistics of systems integration, and considered organizational aspects of systems engineering discipline integration. The systems integration methods have included system exergy analysis, Akaike Information Criteria (AIC), State Variable Analysis, Multidisciplinary Coupling Analysis (MCA), Multidisciplinary Design Optimization (MDO), System Cost Modelling, System Robustness, and Value Modelling. Organizational studies have included the variability of processes in change evaluations, margin management within the organization, information theory of board structures, social categorization of unintended consequences, and initial looks at applying cognitive science to systems engineering. Consortium members have also studied the bidirectional influence of policy and law with systems engineering.

Watson, Michael D.↗

NASA Systems Engineering Research Consortium: Defining the Path to Elegance in Systems

The NASA Systems Engineering Research Consortium was formed at the end of 2010 to study the approaches to producing elegant systems on a consistent basis. This has been a transformative study looking at the engineering and organizational basis of systems engineering. The consortium has engaged in a variety of research topics to determine the path to elegant systems. In the second year of the consortium, a systems engineering framework emerged which structured the approach to systems engineering and guided our research. This led in the third year to set of systems engineering postulates that the consortium is continuing to refine. The consortium has conducted several research projects that have contributed significantly to the understanding of systems engineering. The consortium has surveyed the application of the NASA 17 systems engineering processes, explored the physics and statistics of systems integration, and considered organizational aspects of systems engineering discipline integration. The systems integration methods have included system energy analysis, Akaike Information Criteria (AIC), State Variable Analysis, Multidisciplinary Coupling Analysis (MCA), Multidisciplinary Design Optimization (MDO), System Cost Modeling, System Robustness, and Value Modeling. Organizational studies have included the variability of processes in change evaluations, margin management within the organization, information theory of board structures, social categorization of unintended consequences, and initial looks at applying cognitive science to systems engineering. Consortium members have also studied the bidirectional influence of policy and law with systems engineering.

Watson, Michael D.↗

Connecting dark matter annihilation to the vertex functions of Standard Model fermions

We consider here scenarios in which dark matter is a Majorana fermion which couples to Standard Model fermions through the exchange of charged mediating particles. The matrix elements for various dark matter annihilation processes are then related to one-loop corrections to the fermion-photon vertex, where dark matter and the charged mediators run in the loop. In particular, in the limit where Standard Model fermion helicity mixing is suppressed, the cross section for dark matter annihilation to various final states is related to corrections to the Standard Model fermion charge form factor. These corrections can be extracted in a gauge-invariant manner from collider cross sections. While current measurements from colliders are not precise enough to provide useful constraints on dark matter annihilation, improved measurements at future experiments, such as the International Linear Collider, could improve these constraints by several orders of magnitude, allowing them to surpass the limits obtainable by direct observation

79 ASTRONOMY AND ASTROPHYSICS↗

On the Importance of Using Event-Specific Wave Diffusion Rates in Modeling Diffuse Electron Precipitation

A few to tens of keV electron precipitation that carries substantial energy source down to the upper atmosphere to create aurora is manifested as an important magnetosphere-ionosphere coupling process. The precipitation is usually caused by scattering processes associated with plasma waves in the magnetosphere. The scattering process is often quantified by wave diffusion rates that indicate how fast an electron is scattered. Global models commonly use diffusion coefficients that are derived from statistical wave models. However, due to the statistical nature, many localized, transient features could be smeared out. In this study, we investigate electron precipitation using event-specific diffusion coefficients that are obtained based on simultaneous in-situ measured/inferred, rather than statistical, chorus wave dynamics. We find that the application of the event-specific diffusion coefficients associated with a more dynamic and intense chorus wave model leads more electrons, particularly at several to tens of keV in the dawn-to-noon sector at L > 3, to precipitate than using statistical coefficients. Here, the new simulation roughly captures both the intensity and variability of the precipitating flux as detected by the NOAA/POES satellites. Ionospheric electron density in the lower E region (100–120 km) observed by the mid-latitude Millstone Hill radar is also much better reproduced, while the case using statistical diffusion coefficients underestimates the ionization rate. This study implies the importance of using event-specific diffusion rates in simulating the diffuse electron precipitation and understanding the magnetosphere-ionosphere coupling.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear gyrokinetic predictions of SPARC burning plasma profiles enabled by surrogate modeling

Multi-channel, nonlinear predictions of core temperature and density profiles are performed for the SPARC tokamak accounting for both kinetic neoclassical and fully nonlinear gyro-kinetic turbulent fluxes. A series of flux-tube, nonlinear, electromagnetic simulations using the CGYRO code with six gyrokinetic species are coupled to a nonlinear optimizer using Gaussian process regression techniques. The simultaneous evolution of energy sources, including alpha heat, radiation, and energy exchange, coupled with these high fidelity models and techniques, leads to a converged solution in electron temperature, ion temperature and electron density channels with a minimal number of expensive gyrokinetic simulations without compromising accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simultaneous Decomposition and Dewetting of Nanoscale Alloys: A Comparison of Experiment and Theory

In this work, we consider the coupled process of phase separation and dewetting of metal alloys of nanoscale thickness deposited on solid substrates. The experiments involve applying nanosecond laser pulses that melt the Ag40Ni60 alloy films in two setups: either on thin supporting membranes or on bulk substrates. These two setups allow for extracting both temporal and spatial scales on which the considered processes occur. The theoretical model involves a longwave version of the Cahn–Hilliard formulation used to describe spinodal decomposition, coupled with an asymptotically consistent longwave-based description of dewetting that occurs due to destabilizing interactions between the alloy and the substrate, modeled using the disjoining pressure approach. Careful modeling, combined with linear stability analysis and fully nonlinear simulations, leads to results consistent with the experiments. In particular, we find that the two instability mechanisms occur concurrently, with the phase separation occurring faster and on shorter temporal scales. The modeling results show a strong influence of the temperature dependence of relevant material properties, implying that such a dependence is crucial for the understanding of the experimental findings. The agreement between theory and experiment suggests the utility of the proposed theoretical approach in helping to develop further experiments directed toward formation of metallic alloy nanoparticles of desired properties.

36 MATERIALS SCIENCE↗