Coupling the Lorentz Integral Transform (LIT) and the Coupled Cluster (CC) Methods: A way towards continuum spectra of "not-so-few-body" systems
Explore the source record for details and available documents.
Orlandini, Giuseppina↗
SEARCH · Engineering Papers
Results for “COUPLINGS”
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.
Spin-orbit coupling and magnetic long-range dipolar interaction for coupled spin-lattice dynamics.
Abstract not provided.
Tranchida, Julien↗
Air‐Ice‐Ocean Coupling during a Strong Mid‐Winter Cyclone: Observing Coupled Dynamic Interactions across Scales
Explore the source record for details and available documents.
54 ENVIRONMENTAL SCIENCES↗
A Study on the Integration between SIMT and Scalar Cores: Loosely Coupled to Tightly Coupled
Explore the source record for details and available documents.
Ramshanker, Abinands↗
Validation of a Hybrid Domain Overlapping Coupling Between SAM and CFD Against the TALL-3D Transients
The System Thermal Hydraulics (STH) code SAM has been coupled to the Computational Fluid Dynamics (CFD) code Simcenter STAR-CCM+ utilizing a hybrid domain overlapping method with an explicit coupling in time. The coupling aims to extend the STH code’s applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional (3D) mixing. The coupling method’s numerical stability was verified in the past against two closed-loop configurations, and it was validated against a double T-junction experiment with 3D scalar mixing. In the present work, the coupling method is validated against the TALL-3D STH/CFD coupling benchmark facility. TALL-3D is a three-legged, liquid-metal facility with a large, pool-type enclosure (test section) that exhibits 3D flow effects to be modeled by a CFD code. The rest of the system exhibits approximately 1D behavior well-predicted by an STH code. First, the present STAR-CCM+ CFD model of the 3D test section is validated against experimental data. Then, the SAM-STARCCM+ coupled model is validated against six different TALL-3D steady states, including SAM standalone model results for comparison. Lastly, the SAM-STARCCM+ coupled model is validated against two TALL-3D transients, one exhibiting flow reversal in the test section and one exhibiting nonlinear, Limit Cycle Oscillations (LCO). For the first transient, the SAM-STARCCM+ coupled model properly predicts an increase in the test section’s inlet temperature during flow reversal, and this is not predicted by the SAM standalone model. Following flow reversal, the SAM-STARCCM+ coupled model better-predicts the initial flow recovery and following oscillations as the system approaches a final natural circulation state. For the second transient, no true final steady state is observed due to LCO. Neither the SAM-STARCCM+ coupled model nor the SAM standalone model can perfectly capture the experiment’s changing oscillation frequency during the transient. However, the SAM-STARCCM+ coupled model does reproduce the oscillatory feedback observed in the system. This is a significant achievement as the SAM-STARCCM+ coupled model only uses an explicit coupling in time, as opposed to a semi-implicit coupling. In comparison, previous STH/CFD coupling efforts of the TALL-3D facility required semi-implicit coupling to obtain similar results.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Building energy simulation coupled with CFD for indoor environment: A critical review and recent applications
This paper presents a comprehensive review of the open literature on motivations, methods and applications of linking stratified airflow simulation to building energy simulation (BES). First, we reviewed the motivations for coupling prediction models for building energy and indoor environment. This review classified various exchanged data in different applications as interface data and state data, and found that choosing different data sets may lead to varying performance of stability, convergence, and speed for the co-simulation. Second, our review shows that an external coupling scheme is substantially more popular in implementations of co-simulation than an internal coupling scheme. The external coupling is shown to be generally faster in computational speed, as well as easier to implement, maintain and expand than the internal coupling. Third, the external coupling can be carried out in different data synchronization schemes, including static coupling and dynamic coupling. In comparison, the static coupling that performs data exchange only once is computationally faster and more stable than the dynamic coupling. However, concerning accuracy, the dynamic coupling that requires multiple times of data exchange is more accurate than the static coupling. Furthermore, the review identified that the implementation of the external coupling can be achieved through customized interfaces, middleware, and standard interfaces. The customized interface is straightforward but may be limited to a specific coupling application. The middleware is versatile and user-friendly but usually limited in data synchronization schemes. The standard interface is versatile and promising, but may be difficult to implement. Current applications of the co-simulation are mainly energy performance evaluation and control studies. Finally, we discussed the limitations of the current research and provided an overview for future research.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Effects on calculated half-widths and shifts from the line coupling for asymmetric-top molecules
The refinement of the Robert-Bonamy formalism by considering the line coupling for linear molecules developed in our previous studies [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013); 140, 104304 (2014)] have been extended to asymmetric-top molecules. For H{sub 2}O immersed in N{sub 2} bath, the line coupling selection rules applicable for the pure rotational band to determine whether two specified lines are coupled or not are established. Meanwhile, because the coupling strengths are determined by relative importance of off-diagonal matrix elements versus diagonal elements of the operator −iS{sub 1} − S{sub 2}, quantitative tools are developed with which one is able to remove weakly coupled lines from consideration. By applying these tools, we have found that within reasonable tolerances, most of the H{sub 2}O lines in the pure rotational band are not coupled. This reflects the fact that differences of energy levels of the H{sub 2}O states are pretty large. But, there are several dozen strongly coupled lines and they can be categorized into different groups such that the line couplings occur only within the same groups. In practice, to identify those strongly coupled lines and to confine them into sub-linespaces are crucial steps in considering the line coupling. We have calculated half-widths and shifts for some groups, including the line coupling. Based on these calculations, one can conclude that for most of the H{sub 2}O lines, it is unnecessary to consider the line coupling. However, for several dozens of lines, effects on the calculated half-widths from the line coupling are small, but remain noticeable and reductions of calculated half-widths due to including the line coupling could reach to 5%. Meanwhile, effects on the calculated shifts are very significant and variations of calculated shifts could be as large as 25%.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Analysis of the Temporal Response of Coupled Asymmetrical Zero-Power Subcritical Bare Metal Reactor Systems
The behavior of symmetrical coupled-core systems has been extensively studied, yet there is a dearth of research on asymmetrical systems due to the increased complexity of the analysis of such systems. In this research, the multipoint kinetics method is applied to asymmetrical zeropower, subcritical, bare metal reactor systems. Existing research on asymmetrical reactor systems assumes symmetry in the neutronic coupling; however, it will be shown that this cannot always be assumed. Deep subcriticality adds another layer of complexity and requires modification of the multipoint kinetics equations to account for the effect of the external neutron source. A modified set of multipoint kinetics equations is derived with this in mind. Subsequently, the Rossi-alpha equations are derived for a two-region asymmetrical reactor system. The predictive capabilities of the radiation transport code MCNP6 for neutron noise experiments are shown in a comparison to the results of a series of Rossi-alpha measurements performed by J. Mihalczo utilizing a coupled set of symmetrical bare highly-enriched uranium (HEU) cylinders. The ptrac option within MCNP6 can generate time-tagged counts in a cell (list-mode data). The list-mode data can then be processed similarly to measured data to obtain values for system parameters such as the dual prompt neutron decay constants observable in a coupled system. The results from the ptrac simulations agree well with the historical measured values. A series of case studies are conducted to study the effects of geometrical asymmetry in the coupling between two bare metal HEU cylinders. While the coupling behavior of symmetrical systems has been reported on extensively, that of asymmetrical systems remains sparse. In particular, it appears that there has been no previous research in obtaining the coupling time constants for asymmetrically-coupled systems. The difficulty in observing such systems is due in part to the inability to determine the individual coupling coefficients from measurement: unlike the symmetrical cases, only the product of the values can be obtained. A method is proposed utilizing MCNP6 tally ratios to separate the coupling coefficients for such systems. This work provides insight into the behavior of asymmetrically-coupled systems as the separation distance between the two cores is changed and also as the asymmetry is increased. As the asymmetry increases, both the slower and the faster observable prompt neutron decay constants increase in magnitude. The coupling time constants are determined from the measured decay constants. As the separation distance increases, both coupling coefficients decrease as expected. Based on these findings, an effective computational method utilizing MCNP6 and the Rossialpha technique can be applied to the prediction of asymmetrical coupled system measurements.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Effects on Calculated Half-Widths and Shifts from the Line Coupling for Asymmetric-Top Molecules
The refinement of the Robert-Bonamy formalism by considering the line coupling for linear molecules developed in our previous studies [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013); 140, 104304 (2014)] have been extended to asymmetric-top molecules. For H2O immersed in N2 bath, the line coupling selection rules applicable for the pure rotational band to determine whether two specified lines are coupled or not are established. Meanwhile, because the coupling strengths are determined by relative importance of off-diagonal matrix elements versus diagonal elements of the operator -iS1 -S2, quantitative tools are developed with which one is able to remove weakly coupled lines from consideration. By applying these tools, we have found that within reasonable tolerances, most of the H2O lines in the pure rotational band are not coupled. This reflects the fact that differences of energy levels of the H2O states are pretty large. But, there are several dozen strongly coupled lines and they can be categorized into different groups such that the line couplings occur only within the same groups. In practice, to identify those strongly coupled lines and to confine them into sub-linespaces are crucial steps in considering the line coupling. We have calculated half-widths and shifts for some groups, including the line coupling. Based on these calculations, one can conclude that for most of the H2O lines, it is unnecessary to consider the line coupling. However, for several dozens of lines, effects on the calculated half-widths from the line coupling are small, but remain noticeable and reductions of calculated half-widths due to including the line coupling could reach to 5%. Meanwhile, effects on the calculated shifts are very significant and variations of calculated shifts could be as large as 25%.
formalism↗
Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations
The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.
CDA↗
Exploiting the Liquid Phase to Enhance the Cross-Coupling of Alcohols over Nanoporous Gold Catalysts
Controlling selectivity in the cross-coupling of two similar compounds is a classic challenge in heterogeneous catalysis. Here it is shown that the phase in which the catalysis is performed has a dramatic impact on the selectivity of the oxidative coupling of alcohols to esters over unsupported nanoporous gold catalysts, affording nearly optimal cross-coupling to a single ester at equimolar concentrations in the liquid. Operation in the liquid vs the gas phase affects 1) the relative C-H activation rates of methoxy and 1-propoxy, 2) the fraction of C-H activation which leads to esters vs aldehydes, and 3) the fraction of esters which result from cross-coupling vs self-coupling. While activation of the critical adsorbed reactive alkoxy intermediate is faster for 1-propoxy than methoxy in both phases, the liquid phase is more effective in coupling the resulting aldehyde with adsorbed methoxy or 1-propoxy to yield an ester. Additionally, operation in the liquid phase promotes cross-coupling to methyl propionate, whereas in the gas-phase self-coupling of 1-propanol to propyl propionate is favored. The promotion of self-coupling in the gas phase results from the stabilization of larger alkoxides on the surface by Au-alkyl van der Waals forces. However, such forces do not appear dominant in the liquid phase, as evidenced by similar cross-coupling selectivities of methanol with ethanol, 1-propanol, and 1-butanol. Introducing steric hindrance into the higher alcohol (i.e., 2-methyl-1-propanol) further promotes cross-coupling. Here, this promotion is attributed to a kinetic preference for an aldehyde to couple with less-hindered alkoxides. Altogether, these findings demonstrate that alcohol cross-coupling selectivities are strongly impacted by the phase in which the catalysis is conducted, thus altering the phase provides opportunities for selective and efficient chemical syntheses.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The first-principle coupled calculations using TMCC and CFX for the pin-wise simulation of LWR
The coupling of neutronics and thermal-hydraulics plays an important role in the reactor safety, core design and operation of nuclear power facilities. This paper introduces the research on the coupling of Monte Carlo method and CFD method, specifically using TMCC and CFX. The methods of the coupling including the coupling approach, data transfer, mesh mapping and transient coupling scheme are studied firstly. The coupling of TMCC and CFX for the steady state calculations is studied and described for the single rod model and the 3 x 3 Rod Bundle model. The calculation results prove that the coupling method is feasible and the coupled calculation can be used for steady state calculations. However, the oscillation which occurs during the coupled calculation indicates that this method still needs to be improved for the accuracy. Then the coupling for the transient calculations is also studied and tested by two cases of the steady state and the lost of heat sink. The preliminary results of the transient coupled calculations indicates that the transient coupling with TMCC and CFX is able to simulate the transients but instabilities are occurring. It is also concluded that the transient coupling of TMCC and CFX needs to be improved due to the limitation of computational resource and the difference of time scales. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Present and Future Changes in Land‐Atmosphere Coupling of Water and Energy Over Extratropical Forest Regions
Abstract Couplings between land and the near surface atmosphere are modulated by interactions between soil conditions, vegetation dynamics, turbulent fluxes, and atmospheric properties. How the land‐atmosphere coupling responds to warming and elevated CO 2 are important for understanding the land surface carbon, energy, and water cycles. In this work, we documented this coupled land‐atmosphere network based on observations and the Energy Exascale Earth System Model (E3SM) simulations over extratropical forest ecosystems. We employed a transfer entropy approach and novel network metrics to reveal patterns and strength of the land‐atmosphere coupling under historical conditions and a future high emission scenario (SSP585). We found that, in observations, the present‐day extratropical forest coupling network has high network connectivity (72%–88% of the targeted processes are coupled). E3SM reasonably captured the extratropical forest coupling network (modeled network connectivity was 81%–96%) and predicted that the coupling strength would significantly increase by 28% (±3%) under warming and elevated CO 2 conditions. Furthermore, E3SM factorial coupled experiments suggested that warming enhanced soil nitrogen mineralization favoring plant nitrogen uptake and vegetation growth were responsible for the strengthening future land‐atmosphere coupling. This work provides new metrics to analyze and document complex couplings for coupled earth system processes and highlights the important roles soil nutrient availability and biogeochemistry have on land‐atmosphere coupling.
54 ENVIRONMENTAL SCIENCES↗
Assessment of the sea surface temperature diurnal cycle in CNRM-CM6-1 based on its 1D coupled configuration
A single-column version of the CNRM-CM6-1 global climate model has been developed to ease development and validation of the boundary layer physics and air–sea coupling in a simplified environment. This framework is then used to assess the ability of the coupled model to represent the sea surface temperature (SST) diurnal cycle. To this aim, the atmospheric–ocean single-column model (AOSCM), called CNRM-CM6-1D, is implemented in a case study derived from the CINDY2011/DYNAMO campaign over the Indian Ocean, where large diurnal SST variabilities have been well documented. Comparing the AOSCM and its uncoupled components (atmospheric SCM and oceanic SCM, called OSCM) highlights the fact that the impact of coupling in the atmosphere results from both the possibility to take into account the diurnal variability of SST, which is not usually available in forcing products, and the change in mean state SST as simulated by the OSCM, with the ocean mean state not being heavily impacted by the coupling. This suggests that coupling feedbacks in the 3D model do not arise from the coupling of ocean and atmosphere vertical column physics but are more due to the large-scale dynamics resolved by the 3D model. Additionally, a sub-daily coupling frequency is needed to represent the SST diurnal variability, but the choice of the coupling time step between 15 min and 3 h does not impact the diurnal temperature range simulated much. The main drawback of a 3 h coupling is delaying the SST diurnal cycle by 5 h in asynchronous coupled models. Overall, the diurnal SST variability is reasonably well represented in CNRM-CM6-1 with a 1 h coupling time step and the upper-ocean model resolution of 1 m. This framework is shown to be a very valuable tool to develop and validate the boundary layer physics and the coupling interface. It highlights the interest to develop other atmosphere–ocean coupling case studies.
54 ENVIRONMENTAL SCIENCES↗
Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well
Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.
deformation flow↗
Synthetic and Mechanistic Implications of Chlorine Photoelimination in Nickel/Photoredox C(sp 3 )–H Cross-Coupling
In recent years, the development of light-driven reactions has contributed numerous advances in synthetic organic chemistry. A particularly active research area combines photoredox catalysis with nickel catalysis to accomplish otherwise inaccessible cross-coupling reactions. In these reactions, the photoredox catalyst absorbs light to generate an electronically excited charge-transfer state that can engage in electron or energy transfer with a substrate and the nickel catalyst. Our group questioned whether photoinduced activation of the nickel catalyst itself could also contribute new approaches to cross-coupling. Over the past 5 years, we have sought to advance this hypothesis for the development of a suite of mild and site-selective C(sp 3 )–H cross-coupling reactions with chloride-containing coupling partners via photoelimination of a Ni–Cl bond. On the basis of a report from the Nocera laboratory, we reasoned that photolysis of a Ni(III) aryl chloride species, generated by single-electron oxidation of a typical Ni(II) intermediate in cross-coupling, might allow for the catalytic generation of chlorine atoms. Combining this with the ability of Ni(II) to accept alkyl radicals, we hypothesized that photocatalytically generated chlorine atoms could mediate hydrogen atom transfer (HAT) with C(sp 3 )–H bonds to generate a substrate-derived alkyl radical that is captured by the Ni center in cross-coupling. A photoredox catalyst was envisioned to promote the necessary single-electron oxidation and reduction of the Ni catalyst to facilitate an overall redox-neutral process. Overall, this strategy would offer a visible-light-driven mechanism for chlorine radical formation enabled by the sequential capture of two photons. As an initial demonstration, we developed a Ni/photoredox-catalyzed α-oxy C(sp 3 )–H arylation of cyclic and acyclic ethers. This method was extended to a mild formylation of abundant and complex aryl chlorides through selective 2-functionalization of 1,3-dioxolane. Seeking to develop a suite of reactions that introduce carbon at all different oxidation states, we explored C(sp 3 )–H cross-coupling with trimethyl orthoformate, a common laboratory solvent. We found that trimethyl orthoformate serves as a source of methyl radical for a methylation reaction via β-scission from a tertiary radical generated upon chlorine-mediated HAT. Since chlorine radical is capable of abstracting unactivated C(sp 3 )–H bonds, our efforts have also been directed at cross-coupling with a range of feedstock chemicals, such as alkanes and toluenes, along with late-stage intermediates, using chloroformates as coupling partners. Overall, this platform enables access to valuable synthetic transformations with (hetero)aryl chlorides, which despite being the most ubiquitous and inexpensive aryl halide coupling partners, are rarely reactive in Ni/photoredox catalysis. Little is known about the photophysics and photochemistry of organometallic Ni complexes relevant to cross-coupling. We have conducted mechanistic investigations, including computational, spectroscopic, emission quenching, and stoichiometric oxidation studies, of Ni(II) aryl halide complexes common to Ni/photoredox reactions. These studies indicate that chlorine radical generation from excited Ni(III) is operative in the described C(sp 3 )–H functionalization methods. More generally, the studies illustrate that the photochemistry of cross-coupling catalysts cannot be ignored in metallaphotoredox reactions. We anticipate that further mechanistic understanding should facilitate catalyst design and lead to the development of new synthetic methods.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Conservative tightly-coupled simulations of stochastic multiscale systems
Multiphysics problems often involve components whose macroscopic dynamics is driven by microscopic random fluctuations. The fidelity of simulations of such systems depends on their ability to propagate these random fluctuations throughout a computational domain, including subdomains represented by deterministic solvers. When the constituent processes take place in nonoverlapping subdomains, system behavior can be modeled via a domain-decomposition approach that couples separate components at the interfaces between these subdomains. Its coupling algorithm has to maintain a stable and efficient numerical time integration even at high noise strength. We propose a conservative domain-decomposition algorithm in which tight coupling is achieved by employing either Picard's or Newton's iterative method. Coupled diffusion equations, one of which has a Gaussian white-noise source term, provide a computational testbed for analysis of these two coupling strategies. Fully-converged (“implicit”) coupling with Newton's method typically outperforms its Picard counterpart, especially at high noise levels. This is because the number of Newton iterations scales linearly with the amplitude of the Gaussian noise, while the number of Picard iterations can scale superlinearly. At large time intervals between two subsequent inter-solver communications, the solution error for single-iteration (“explicit”) Picard's coupling can be several orders of magnitude higher than that for implicit coupling. Increasing the explicit coupling's communication frequency reduces this difference, but the resulting increase in computational cost can make it less efficient than implicit coupling at similar levels of solution error, depending on the communication frequency of the latter and the noise strength. This trend carries over into higher dimensions, although at high noise strength explicit coupling may be the only computationally viable option.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Results from Tight and Loose Coupled Multiphysics in Nuclear Fuels Performance Simulations using BISON
The behavior of nuclear fuel in the reactor environment is affected by multiple physics, most notably heat conduction and solid mechanics, which can have a strong influence on each other. To provide credible solutions, a fuel performance simulation code must have the ability to obtain solutions for each of the physics, including coupling between them. Solution strategies for solving systems of coupled equations can be categorized as loosely-coupled, where the individual physics are solved separately, keeping the solutions for the other physics fixed at each iteration, or tightly coupled, where the nonlinear solver simultaneously drives down the residual for each physics, taking into account the coupling between the physics in each nonlinear iteration. In this paper, we compare the performance of loosely and tightly coupled solution algorithms for thermomechanical problems involving coupled thermal and mechanical contact, which is a primary source of interdependence between thermal and mechanical solutions in fuel performance models. The results indicate that loosely-coupled simulations require significantly more nonlinear iterations, and may lead to convergence trouble when the thermal conductivity of the gap is too small. We also apply the tightly coupled solution strategy to a nuclear fuel simulation of an experiment in a test reactor. Studying the results from these simulations indicates that perhaps convergence for either approach may be problem dependent, i.e., there may be problems for which a loose coupled approach converges, where tightly coupled won’t converge and vice versa.
97 MATHEMATICS AND COMPUTING↗