Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled 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 109 records · Page 6

Coupled Climate Simulations With E3SM‐MMF

Simulations of the recent historical period from 1950 to 2014 are conducted with E3SM‐MMF, which uses an embedded 2D cloud resolving model that runs efficiently on GPUs in place of traditional parameterizations for cloud and turbulence. Analysis of the climate and variability reveal several aspects where E3SM‐MMF produces smaller biases compared to E3SMv2, including better agreement with the observed evolution of global mean surface temperature, although the representation of ENSO is too weak and fast. Three idealized abrupt CO 2 experiments were also conducted to assess climate sensitivity and feedbacks. These yield three estimates of effective climate sensitivity (4.38, 5.21, and 6.06 K), with a corresponding spread in the shortwave cloud feedbacks. These estimates are on the higher end of sensitivity estimates from CMIP ensembles, and the spread indicates substantial state‐dependent feedbacks. These results demonstrate how multiscale modeling framework (MMF) models can be used for climate relevant experiments and projections by leveraging modern GPU enabled computational platforms. The unique qualities of E3SM‐MMF shown in previous literature are largely still present, but various instances of reduced biases suggest that MMF models have utility in improving future projections.

E3SM↗

Describing the Mechanism of Instability Suppression Using a Central Pilot Flame With Coupled Experiments and Simulations

Pilot flames are commonly used to extend combustor operability limits and suppress combustion oscillations in low-emissions gas turbines. Combustion oscillations, a coupling between heat release rate oscillations and combustor acoustics, can arise at the operability limits of low-emissions combustors where the flame is more susceptible to perturbations. While the use of pilot flames is common in land-based gas turbine combustors, the mechanism by which they suppress instability is still unclear. In this study, we consider the impact of a central jet pilot on the stability of a swirl-stabilized flame in a variable-length, single-nozzle combustor. Previously, the pilot flame was found to suppress the instability for a range of equivalence ratios and combustor lengths. We hypothesize that combustion oscillation suppression by the pilot occurs because the pilot provides hot gases to the vortex breakdown region of the flow that recirculate and improve the static, and hence dynamic, stability of the main flame. This hypothesis is based on a series of experimental results that show that pilot efficacy is a strong function of pilot equivalence ratio but not pilot flow rate, which would indicate that the temperature of the pilot products as well as the combustion intensity of the pilot flame play more of a role in oscillation stabilization than the length of the pilot flame relative to the main flame. Further, the pilot-flame efficacy increases with pilot-flame equivalence ratio until it matches the main-flame equivalence ratio; at pilot equivalence ratios greater than the main equivalence ratio, the pilot-flame efficacy does not change significantly with pilot equivalence ratio. To understand these results, we use large-eddy simulation (LES) to provide a detailed analysis of the flow in the region of the pilot flame and the transport of radical species in the region between the main flame and pilot flame. Furthermore, the simulation, using a flamelet/progress variable-based chemistry tabulation approach and standard eddy viscosity/diffusivity turbulence closure models, provides detailed information that is inaccessible through experimental measurements.

42 ENGINEERING↗

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗

Machine Learning Calibration of Groundwater Table Depth in ELM: Impact on Land Surface Hydrology and Land‐Atmosphere Fluxes

Accurate representation of groundwater table depth (GWTD) is crucial for simulating hydrological cycling in Earth system models (ESM). Nevertheless, there is a notable gap in the literature regarding the validation of GWTD simulations in ESMs and their subsequent impact on downstream hydrological components. This study explores the calibration of parameterization of global GWTD using machine learning within the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Despite achieving significant gains in simulating GWTD through calibration, offline ELM simulations unexpectedly show that these improvements do not translate to substantial enhancements in model performance for other key hydrological variables, including soil moisture (SM), runoff, groundwater contribution to runoff or base flow index (BFI), and evapotranspiration and its partitioning. The performance in SM and runoff was even degraded in some regions, while BFI was mostly overestimated. Although there is significant improvement in GWTD within the critical range of 1–5 m, where groundwater traditionally influences land surface energy fluxes, these improvements occurred mostly in humid areas where the impact of GWTD on surface processes is minimal. Although the impacts of model calibration are generally small in offline ELM simulations, coupled land-atmosphere simulations exhibit much stronger responses to GWTD calibration, highlighting the role of land-atmosphere feedbacks in Earth system modeling. These findings underscore the need for integrated calibration strategies that simultaneously optimize multiple hydrological variables. However, if a single-variable approach is necessary, it is crucial to establish clear priorities for calibration, identifying the most critical variables that have the greatest impact on overall model performance.

Fang, Yilin [Pacific Northwest National Laboratory↗

A Three-Dimensional Coupled Internal/External Simulation of a Film-Cooled Turbine Vane

A three-dimensional Navier-Stokes simulation has been performed for a realistic film-cooled turbine vane using the LeRC-HT code. The simulation includes the flow regions inside the coolant plena and film cooling holes in addition to the external flow. The vane is the subject of an upcoming NASA Glenn Research Center experiment and has both circular cross-section and shaped film cooling holes. This complex geometry is modeled using a multi-block grid which accurately discretizes the actual vane geometry including shaped holes. The simulation matches operating conditions for the planned experiment and assumes periodicity in the spanwise direction on the scale of one pitch of the film cooling hole pattern. Two computations were performed for different isothermal wall temperatures, allowing independent determination of heat transfer coefficients and film effectiveness values. The results indicate separate localized regions of high heat transfer coefficient values, while the shaped holes provide a reduction in heat flux through both parameters. Hole exit data indicate rather simple skewed profiles for the round holes, but complex profiles for the shaped holes with mass fluxes skewed strongly toward their leading edges.

Heidmann, James D.↗

Microscopic Foundation and Simulation of Coupled Carrier-Temperature Diffusions in Semiconductor Lasers

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.

Li, J.↗

Computational Study of Modeling Fully-coupled Combustor-Turbine Interactions by the Open National Combustion Code (OpenNCC)

We investigate the flow fields inside the combustor and the high-pressure turbine including the stator and moving rotors, from the Energy Efficient Engine (E3) program. We then compare the fully-coupled simulations results with the prior results from the sequential approach, where key performance metrics for the combustor and HPT are compared. The results will provide one example of the impact of performing fully-coupled combustor/HPT simulations relative to sequential (one-way coupled) combustor and HPT simulations.

CFD↗

Computational Study of Modeling Fully-coupled Combustor-Turbine Interactions by the Open National Combustion Code (OpenNCC)

We investigate the flow fields inside the combustor and the high-pressure turbine including the stator and moving rotors, from the Energy Efficient Engine (E3) program. We then compare the fully-coupled simulations results with the prior results from the sequential approach, where key performance metrics for the combustor and HPT are compared. The results will provide one example of the impact of performing fully-coupled combustor/HPT simulations relative to sequential (one-way coupled) combustor and HPT simulations.

CFD↗

Numerical Study of Combustor-Turbine Interactions Using Open-National Combustion Code (OpenNCC)

We investigate the flow fields inside the combustor and the high-pressure turbine including the stator and moving rotors, from the Energy Efficient Engine (E3) program. We then compare the fully-coupled simulations results with the prior results from the sequential approach, where key performance metrics for the combustor and HPT are compared. The results will provide one example of the impact of performing fully-coupled combustor/HPT simulations relative to sequential (one-way coupled) combustor and HPT simulations.

CFD↗

Carbon–concentration and carbon–climate feedbacks in CMIP6 models and their comparison to CMIP5 models

Results from the fully and biogeochemically coupled simulations in which CO 2 increases at a rate of 1 % yr –1 (1pctCO2) from its preindustrial value are analyzed to quantify the magnitude of carbon–concentration and carbon–climate feedback parameters which measure the response of ocean and terrestrial carbon pools to changes in atmospheric CO 2 concentration and the resulting change in global climate, respectively. The results are based on 11 comprehensive Earth system models from the most recent (sixth) Coupled Model Intercomparison Project (CMIP6) and compared with eight models from the fifth CMIP (CMIP5). The strength of the carbon–concentration feedback is of comparable magnitudes over land (mean ± standard deviation = 0.97 ± 0.40 PgC ppm –1 ) and ocean (0.79 ± 0.07 PgC ppm –1 ), while the carbon–climate feedback over land (–45.1 ± 50.6 PgC °C –1 ) is about 3 times larger than over ocean (–17.2 ± 5.0 PgC °C –1 ). The strength of both feedbacks is an order of magnitude more uncertain over land than over ocean as has been seen in existing studies. These values and their spread from 11 CMIP6 models have not changed significantly compared to CMIP5 models. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The transient climate response to cumulative emissions (TCRE) from the 11 CMIP6 models considered here is 1.77 ± 0.37 °C EgC –1 and is similar to that found in CMIP5 models (1.63 ± 0.48 °C EgC –1 ) but with somewhat reduced model spread. The expressions for feedback parameters based on the fully and biogeochemically coupled configurations of the 1pctCO2 simulation are simplified when the small temperature change in the biogeochemically coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters is used to gain insight into the reasons for differing responses among ocean and land carbon cycle models.

54 ENVIRONMENTAL SCIENCES↗