Search NASA⌕ Search

SEARCH · Search NASA

Results for “stability”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Insights into dual-functional modification for water stability enhancement of mesoporous zirconium metal–organic frameworks

The stability of metal–organic frameworks (MOFs) in water affects their ability to function as chemical catalysts, their capacity as adsorbents for separations in water vapor presence, and their usefulness as recyclable water harvesters. Here, we have examined water stability of four node-modified variants of the mesoporous MOF, NU-1000, namely formate-, Acac-, TFacac-, and Facac-NU-1000, comparing these with node-accessible NU-1000. These NU-1000 variants present ligands grafted to NU-1000's hexa-Zr(IV)-oxy nodes by displacing terminal aqua and hydroxo ligands. Facac-NU-1000, containing the most hydrophobic ligands, showed the greatest water stability, being able to undergo at least 20 water adsorption/desorption cycles without loss of water uptake capacity. Computational studies revealed dual salutary functions of installed Facac ligands: (1) enhancement of framework mechanical stability due to electrostatic interactions; and (2) transformation and shielding of the otherwise highly hydrophilic nodes from H-bonding interactions with free water, presumably leading to weaker channel-stressing capillary forces during water evacuation – consistent with trends in free energies of dehydration across the NU-1000 variants. Water harvesting and hydrolysis of chemical warfare agent simulants were examined to gauge the functional consequences of modification and mechanical stabilization of NU-1000 by Facac ligands. The studies revealed a harvesting capacity of ~1.1 L of water vapor per gram of Facac-NU-1000 per sorption cycle. They also revealed retention of catalytic MOF activity following 20 water uptake and release cycles. This study provides insights into the basis for node-ligand-engendered stabilization of wide-channel MOFs against collapse during water removal.

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

Critical role of oxide anchoring groups in organic electron transport layers for perovskite solar cell stability

Charge transport layers (CTLs) play a critical role in the performance and long-term stability of perovskite solar cells (PSCs) by facilitating efficient charge extraction and providing interfacial stabilization. For organic CTLs, different functional groups have been designed to modulate the interactions with the substrate. Anchoring groups are used to establish strong adhesion to transparent conductive oxides (TCOs) but are not routinely incorporated in organic CTL design. In this study, we investigate the influence of anchoring groups in organic electron transport layers (ETLs) on PSC performance and stability by synthesizing and evaluating two naphthalene diimide (NDI)-based ETLs: NDI-(PhPA) 2 , functionalized with polar phenyl phosphonic acid (PhPA) groups, and NDI-(PhBr) 2 , modified with nonpolar bromophenyl (PhBr) groups. X-ray photoelectron spectroscopy confirms that PhPA functionalization leads to NDIs anchored onto fluorine-doped tin oxide (FTO) following chemical bath deposition, while NDI-(PhBr) 2 shows no detectable surface presence after deposition. We show that NDI-(PhPA) 2 deposits as a thin film on FTO, exhibiting strong adhesion, high solvent resistance, and enhanced light and thermal stability, achieving a maximum device efficiency of 14.3%. Stability assessments conducted under continuous illumination at 25 °C for 200 hours and at 65 °C for 100 hours demonstrate performance on par with TiO 2 -based ETLs. In contrast, NDI-(PhBr) 2 , which lacks strong anchoring interactions, fails to adhere to the FTO surface, resulting in a rapid decrease in efficiency over a few hours, comparable to devices without an ETL. These findings highlight the importance of anchoring groups in organic CTLs for achieving long-term PSC stability under operational conditions.

Kim, Sanggyun [Georgia Institute of Technology, At↗
DOE OSTI2026

Stability Constants of Lanthanide-nitrate Complexes in Aqueous Solutions: A Theoretical Study

Calculating stability constants for lanthanide--nitrate complexes in aqueous solution is challenging due to the complex free-energy landscapes of the participating species. In this work, we compare cluster-continuum solvation and condensed-phase approaches using universal machine learning interatomic potentials (MLIPs) for determining the stability constants of lanthanide--nitrate complexes in aqueous solutions. Within the cluster--continuum solvation framework at the B3LYP level of theory, reactions involving lanthanide coordination numbers of both 8 and 9 are found to be relevant. After an empirical linear free-energy correction, the cluster--continuum results fall on the same order-of-magnitude scale as the experimental stability constants. By contrast, MACE-MP0 MLIP underestimates lanthanide hydration numbers and gives PMF-derived stability constants with large deviations from experiment, whereas MACE-MATPES-R2SCAN improves both hydration structure and the stability-constant scale but still does not quantitatively reproduce the detailed lanthanide trend. Across both approaches, nitrate binding is best viewed as a labile coordination motif rather than a fixed mono- or bidentate structure, with hydration-shell structure influencing which configurations are favored. Overall, the cluster--continuum calculations provide a practical semi-quantitative baseline for the experimental stability-constant scale, while the explicit-solvent MLIP benchmarks show clear progress from MACE-MP0 to MACE-MATPES-r2SCAN but also highlight the need for lanthanide-targeted training, fine-tuning, or improved long-range and polarization treatments to obtain predictive thermodynamics for complex aqueous lanthanide chemistry.

Dinpajooh, Mohammadhasan↗
DOE OSTI2020

Predicting metal–metal interactions. I. The influence of strain on nanoparticle and metal adlayer stabilities

Strain-engineering of bimetallic nanomaterials is an important design strategy for developing new catalysts. In this paper, we introduce an approach for including strain effects into a recently introduced, density functional theory (DFT)-based alloy stability model. The model predicts adsorption site stabilities in nanoparticles and connects these site stabilities with catalytic reactivity and selectivity. Strain-based dependencies will increase the model’s accuracy for nanoparticles affected by finite-size effects. In addition to the stability of small nanoparticles, strain also influences the heat of adsorption of epitaxially grown metal-on-metal adlayers. In this respect, we successfully benchmark the strain-including alloy stability model with previous experimentally determined trends in the heats of adsorption of Au and Cu adlayers on Pt (111). For these systems, our model predicts stronger bimetallic interactions in the first monolayer than monometallic interactions in the second monolayer. We explicitly quantify the interplay between destabilizing strain effects and the energy gained by forming new metal–metal bonds. While tensile strain in the first Cu monolayer significantly destabilizes the adsorption strength, compressive strain in the first Au monolayer has a minimal impact on the heat of adsorption. Hence, this study introduces and, by comparison with previous experiments, validates an efficient DFT-based approach for strain-engineering the stability, and, in turn, the catalytic performance, of active sites in bimetallic alloys with atomic level resolution.

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

The theory of kinetic effects on resistive wall mode stability in tokamaks

Tokamak fusion plasmas benefit from high pressures but are then susceptible to modes of instability. These magnetohydrodynamic (MHD) modes are macroscopic distortions of the plasma, but certain collective motions of individual particles can provide stabilizing effects opposing them. The presence of a resistive wall slows the mode growth, converting a kink to a resistive wall mode (RWM). A kinetic MHD model includes Maxwell's equations, ideal MHD constraints, and kinetic effects included through the pressure tensor, calculated with the perturbed drift-kinetic distribution function of the particles. The kinetic stabilizing effects on the RWM arise through resonances between the plasma rotation and particle drift motions: precession, bounce, and transit. A match between particle motions and the mode allows efficient transfer of energy that would otherwise drive the growth of the mode, thus damping the growth. The first approach to calculating RWM stability is to write a set of equations for the complex mode frequency in terms of known quantities and then to solve the system. The “energy principle” approach, which has the advantage of clarity in distinguishing the various stabilizing and destabilizing effects, is to change the force balance equation into an equation in terms of changes of kinetic and potential energies, and then to write a dispersion relation for the mode frequency in terms of those quantities. These methods have been used in various benchmarked codes to calculate kinetic effects on RWM stability. Importantly, the theory has illuminated the important roles of plasma rotation, energetic particles, and collisions in RWM stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effect of anisotropic thermal transport on tearing mode stability in negative versus positive triangularity plasmas

The combined effects of anisotropic thermal transport and the plasma shaping, including negative triangularity, on the n = 1 (n is the toroidal mode number) tearing mode (TM) stability are numerically investigated utilizing the MARS-F code [Liu et al., Phys. Plasmas 7, 3681–3690 (2000)]. While varying the plasma boundary triangularity, the TM stability is found to be dictated by the competing effects of the Shafranov shift induced stabilization and the bad-curvature induced destabilization. The negative triangularity shape increases the Shafranov shift (stabilizing) in the plasma core but also enlarges bad-curvature regions (destabilizing) near the plasma edge, with the net effect being largely destabilizing for the TM as compared to the positive triangularity counter-part. Large negative triangularity however can also lead to more stabilization for the plasma core-localized TM. Anisotropic thermal transport reduces the stabilizing effect on the TM associated with the favorable averaged curvature, resulting in more unstable core-localized TMs in both negative and positive triangularity plasmas. But the opposite effect can also take place for the edge-localized TM in finite-pressure plasmas with negative triangularity.

Physics↗
DOE OSTI2020

Stabilizing effect of enhanced resistivity on peeling-ballooning instabilities on EAST

Previous stability analysis of NSTX equilibrium with lithium-conditioning demonstrates that enhanced resistivity due to the increased effective charge number Z eff (i.e. increased impurity level) can provide a stabilizing effect on low-n edge localized modes. This paper extends the resistivity stabilizing effect to the intermediate-n peeling-ballooning (PB) instabilities with the linear stability analysis of EAST high-confinement mode equilibria in NIMROD two-fluid calculations. However, in this work, the resistivity stabilizing effect on PB instabilities in the EAST tokamak appears weaker than that found in NSTX. This work may provide better insight into the physical mechanism behind the beneficial effects of impurities on pedestal stability.

EAST tokamak↗
DOE OSTI2021

A comparative study of internal kink stability in EU DEMO designs with negative and positive triangularity

Abstract Internal kink (IK) instability is investigated for European demonstration fusion reactor (EU DEMO) plasmas in both negative triangularity (NT) and positive triangularity (PT) configurations. For NT plasmas, the IK becomes more unstable as an ideal conformal wall moves away from the plasma boundary, with the mode growth rate saturating at the wall radial location of about b / a = 1.5 , where a is the plasma minor radius and b the wall radial location. The plasma resistivity destabilizes the IK mode. The effect of sub-sonic toroidal plasma flow is sufficiently weak and can thus be ignored for these EU DEMO equilibria. These results are consistent with those for PT plasmas, albeit with larger mode growth rate in the NT configuration. Both perturbative and self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid calculations predict (partial) stabilization of the IK modes in both NT and PT configurations, with inclusion of various kinetic contributions. Precessional drift motion of trapped fusion-born alphas in EU DEMO produces weak stabilization to the IK mode. Stronger stabilization occurs with the toroidal precession of trapped thermal particles (ions and electrons) and the bounce-transit motion of thermal ions. The stabilization is similar between the NT and PT configurations, due to the similarity of the mode eigenfunction (occupying a nearly circular region in the plasma core) despite the sign difference in the triangularity. The non-perturbative MHD-kinetic hybrid model predicts much less stabilization of the mode than the perturbative model, primarily due to the self-consistent determination of the mode eigenvalue in the former. Generally, no significant difference in the IK mode stability is found between the NT and PT plasmas in EU DEMO.

Physics↗
DOE OSTI2024

Small-Signal Stability Impacts of Load and Network Dynamics on Grid-Forming Inverters

This paper presents several stability analyses for grid-forming inverters and synchronous generators considering the dynamics of transmission lines and different load models. Load models are usually of secondary importance compared to generation source models, but as the results show, they play a crucial role in stability studies with the introduction of inverter-based resources. Given inverter control time scales, the implications of considering or neglecting electromagnetic transients of the network are very relevant in the stability assessments. In this paper, we perform eigenvalue analyses for inverter-based resources and synchronous machines connected to a load and explore the effects of multiples models under different network representations. We explore maximum loadability of inverter-based resources and synchronous machines, while analyzing the effects of load and network dynamic models on small-signal stability. The results show that the network representation plays a fundamental role in the stability of the system of different load models. The resulting stability regions are significantly different depending on the source and load model considered.

grid forming inverters↗
DOE OSTI2025

Compatibilization Strategy and Mechanism for Co-stabilizing Commingled Plastics and Pyrolyzed Rubber in Asphalt

Hot mix asphalt mixture is considered the ideal approach to reuse waste plastics in high-value applications because of its very high amount of usage in highway construction. However, the differences in polarity and density between polymers and asphalt lead to polymer coalescence and therefore the poor storage stability of modified asphalt. These challenges are exalted when recycling commingled plastics. This study introduced an innovative compatibilization strategy and mechanism for co-stabilizing commingled plastics and pyrolyzed rubber in asphalt. Commingled plastics were first grafted with maleic anhydride for surface activation, followed by reactive kneading with pyrolyzed rubber and crosslinking agent to form an integrated thermoplastic elastomer (ITPE) for asphalt modification. The mechanical, thermal, and interfacial behaviors of the ITPE were evaluated through tensile testing, thermogravimetric analysis, and scanning electron microscopy. The storage stability and rheological properties of the modified binder blends were evaluated through the cigar tube test and dynamic shear rheometer testing. Results demonstrated a successful formation of imide bonds in the ITPE, which can improve the strength, ductility, and thermal stability of rubber–plastic composites. Appropriate utilization of crosslinking agents can improve both rutting and fatigue resistance of ITPE-modified asphalt with good storage stability because of the co-existence of rigid plastic and soft rubbery regimes and the formation of a crosslink network. Furthermore, excessive content of crosslinker led to severe phase separation and reduced storage stability of modified binder blends. Extra crosslinker tended to float in asphalt because of its low density and caused an excessive formation of the crosslink network in the top section of the asphalt.

asphalt binder modifiers↗
DOE OSTI2023

Enhancing Oxygen Stability in Low-Cobalt Layered Oxide Cathode Materials by Three-Dimensional Targeted Doping

In this project, we propose to develop a new concept and a generic platform that can lead to the greatly enhanced stabilization of all high-energy cathode materials, and in particular high-nickel (Ni) and low-cobalt (Co) oxides. The new concept is a 3D doping technology that hierarchically combines surface and bulk doping. We will use surface doping to stabilize the surface of primary particles and also introduce dopants in the bulk to further enhance oxygen stability, conductivity, and structural stability in low-Co oxides under high voltage and deep discharging operating conditions. This new concept not only will deliver a low-cost, high-energy cathode but also will provide a generic method that can stabilize all high-energy cathodes. The proposed novel 3D doping approach is poised to resolve some longstanding challenges in fundamental doping effects on battery materials as well as to reduce Li-ion batteries’ cost and improve their safety, energy density, and lifetime. To tackle this problem, we have formed a highly complementary multi-university/national labs/industry team to enable a doping-central and systematic investigation of low-Co materials and create a knowledge base for many electrode materials to be used in advanced electric vehicles. The successful execution of the proposed project relies on five components that can be carried out by the complementary team members: (1) a theoretical investigation of the surface and bulk stabilizing dopants (Persson), (2) a precise synthesis of materials with targeted doping (Lin and Xin), (3) development of electrolytes for high-Ni low-Co oxides (Xu), (4) multi-scale characterization of the structures and their interfaces by scanning transmission electron microscopy (SEM) and synchrotron X-ray imaging and spectroscopy tools (Xin and Lin), and (5) pouch cell-level integration (Fan). The UCI-led project will enable a doping-central and systematic investigation of low-Co materials and create a knowledge base for many electrode materials to be used in advanced electric vehicles.

25 ENERGY STORAGE↗
DOE OSTI2021

In-Home Mobility Frequency and Stability in Older Adults Living Alone With or Without MCI: Introduction of New Metrics

Background: Older adults spend a considerable amount of time inside their residences; however, most research investigates out-of-home mobility and its health correlates. We measured indoor mobility using room-to-room transitions, tested their psychometric properties, and correlated indoor mobility with cognitive and functional status. Materials and Methods: Community-dwelling older adults living alone ( n = 139; age = 78.1 ± 8.6 years) from the Oregon Center for Aging & Technology (ORCATECH) and Minority Aging Research Study (MARS) were included in the study. Two indoor mobility features were developed using non-parametric parameters (frequency; stability): Indoor mobility frequency (room-to-room transitions/day) was detected using passive infrared (PIR) motion sensors fixed on the walls in four geographic locations (bathroom; bedroom; kitchen; living room) and using door contact sensors attached to the egress door in the entrance. Indoor mobility stability was estimated by variances of number of room-to-room transitions over a week. Test-retest reliability (Intra-class coefficient, ICC) and the minimal clinically important difference (MCID) defined as the standard error of measurement (SEM) were generated. Generalized estimating equations models related mobility features with mild cognitive impairment (MCI) and functional status (gait speed). Results: An average of 206 days (±127) of sensor data were analyzed per individual. Indoor mobility frequency and stability showed good to excellent test-retest reliability (ICCs = 0.91[0.88–0.94]; 0.59[0.48–0.70]). The MCIDs of mobility frequency and mobility stability were 18 and 0.09, respectively. On average, a higher indoor mobility frequency was associated with faster gait speed (β = 0.53, p = 0.04), suggesting an increase of 5.3 room-to-room transitions per day was associated with an increase of 10 cm/s gait speed. A decrease in mobility stability was associated with MCI (β = −0.04, p = 0.03). Discussion: Mobility frequency and stability in the home are clinically meaningful and reliable features. Pervasive-sensing systems deployed in homes can objectively reveal cognitive and functional status in older adults who live alone.

Wu, Chao-Yi↗
DOE OSTI2022

High Temperature, High Pressure Stability of Aqueous Foams for Potential Application in Enhanced Geothermal System (EGS)

Enhanced Geothermal Systems (EGS) are typically constructed by injecting high-pressure water into deep hot dry rocks (HDR) under carefully controlled conditions to create new or re-open existing fractures, which usually uses an immense quantity of water. In lieu, a more sustainable technique is to utilize aqueous foams as fracturing fluid to reduce water consumption and waste-water treatments of conventional hydraulic fracturing. Although foam-based fracturing has shown promising results in oil and gas industries, its feasibility is not demonstrated in EGS conditions that usually involve high temperature and high pressures. One potential barrier of utilizing foams in EGS applications is that foams are thermodynamically unstable and will become more unstable with increasing temperature due to liquid drainage, bubble coarsening, and coalescence. This work focuses on evaluation of the stability of selected aqueous foams under high-temperature and high-pressure conditions. Specifically, foams generated with surfactant alfa olefin sulfonate (AOS) were studied at temperature up to 200oC, while the pressure ranged between 100 psi and 1000 psi. The effect of additional stabilizing agents was also examined, including guar gum, bentonite clay, borate salt crosslinker, silicon dioxide nanoparticles (SiO2), and graphene oxide (GO) dispersion. Results showed that the stabilizing agents can enhance the foam thermal stability. Foams made with AOS and the borate salt crosslinker exhibited the longest half-life of 20 min at 200°C when the pressure was at 1000 psi. Data fitting showed that foam stability decreased exponentially as temperature increased. On the other hand, pressure increased the foam half-life which followed a power model. This study indicates that it may be possible to obtain highly stable foams at high-temperature and high-pressure conditions with appropriate stabilizing agents.

Thakore, Viren↗
DOE OSTI2024

Influence of Joule Heating on the Stability of High Temperature Irradiation-resistant Thermocouples

Development of in-core instrumentation is driven by the pursuit of safer, more economic energy production from the perspective of both existing nuclear reactors and Generation IV reactor designs. Idaho National Laboratory has developed high-temperature irradiation-resistant thermocouples (HTIR-TCs) for temperature sensing inside Generation IV nuclear reactors. These thermocouples are composed of phosphorus-doped niobium (Nb-P) and lanthana-doped molybdenum (Mo-LaO) thermoelements, an alumina (Al2O3) insulation, and a niobium sheath. HTIR-TCs require an initial heat treatment exceeding the maximum service temperature to stabilize the generated electromotive force (EMF) signal. The mechanism behind the stabilization of the HTIR-TCs through traditional heat treatment methods is understood; however, the traditional heat treatment method is expensive, time consuming, and results in a heterogeneous microstructure. Therefore, we investigated a rapid method for stabilization and microstructure homogeneity, through Joule heating. This work evaluates the impact of Joule heating on the thermoelements’ microstructures, chemical stability, and mechanical properties so as to determine the mechanisms by which stabilization of the EMF signal occurs. Accordingly, during the Joule heat treatment, a secondary Nb3P phase coarsened along the length of the Nb-P thermoelement, along with the formation of an interaction region at the Al2O3/niobium interface. The interaction between the alumina insulation and the Nb-P thermoelement was also observed within the Mo-LaO thermoelement. Joule heating induced stability within the generated HTIR-TC EMF signal via the formation of secondary phases within the Nb-P and the interaction between the alumina insulation and the thermoelements.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2023

Deep Sigma Point Processes-Assisted Chance-Constrained Power System Transient Stability Preventive Control

Here this paper proposes a deep sigma point processes (DSPP)-assisted chance-constrained power system transient stability preventive control method to deal with uncertain renewable energy and loads-induced stability risk. The traditional transient stability-constrained preventive control is reformulated as a chance-constrained optimization problem. To deal with the computational bottleneck of the time-domain simulation-based probabilistic transient stability assessment, the DSPP is developed. DSPP is a parametric Bayesian approach that allows us to predict system transient stability with high computational efficiency while accurately quantifying the confidence intervals of the predictions that can be used to inform system instability risk. To this end, with a given preset confidence probability, we embed DSPP into the primal dual interior point method to help solve the chance-constrained preventive control problem, where the corresponding Jacobian and Hessian matrices are derived. Comparison results with other existing methods show that the proposed method can significantly speed up preventive control while maintaining high accuracy and convergence

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Reversed Impedance-Based Stability Criterion for IBR Grids: Preprint

The existing impedance-based stability criterion is effective for analyzing local control interactions; however, it is difficult to scale the existing criterion to analyze wide-area control interactions among numerous IBRs through a complex power system network. The scaled version of the existing criterion requires the impedance response of each IBR in the system as well as of the network looking from all the IBRs. It is quite challenging to obtain all these impedance responses because of the computational effort and the requirement of separately scanning the impedance of the network and the IBRs. We propose a reversed criterion for the impedance-based stability analysis to address these problems. In contrast to the existing criterion, the reversed criterion analyzes the stability of a power system when an IBR is disconnected from the system. The reversed criterion estimates the impact of an IBR on the frequency and damping of power system oscillation modes using the impedance scans of only the IBR and the grid at its terminal. It can be sequentially applied at different IBRs to evaluate their impact on the power system stability. In addition to scalability, the reversed criterion gives flexibility to focus only on a few selected IBRs, depending on their rating, the magnitude of oscillations observed at their terminals, and the vendor support available for implementing stabilizing control system updates. The reversed criterion is demonstrated on a 14-bus power system with 100% IBRs.

control interactions↗
DOE OSTI2023

Improving Stability of Triple-Cation Perovskite Solar Cells under High-Temperature Operation

Metal halide perovskite photovoltaic performance required for commercial technology encompasses both efficiency and stability. Advances in both these parameters have recently been reported; however, these strategies are often difficult to directly compare due to differences in perovskite composition, device architecture, fabrication methods, and accelerated stressors applied in stability tests. In particular, it is found that there is a distinct lack of elevated temperature, operational (light and bias) stability data. Furthermore, significant testing is required to understand the interactions when combinations are used (e.g., additives used with posttreatments). Herein, individual and combined additive, posttreatment, and contact layer strategies from recent literature reports under standardized operational stability tests of p-i-n CsMAFA perovskites at 70 degrees C are evaluated. Through analysis of over 1000 devices, it is concluded that the hole-transport layer (HTL) is the most significant component impacting elevated temperature operational stability. This analysis motivates future development of high-performance HTLs.

14 SOLAR ENERGY↗
DOE OSTI2020

Stability Analysis of Interface Conditions for Ocean-Atmosphere Coupling

In this paper, we analyze the stability of different coupling strategies for multidomain PDEs that arise in general circulation models used in climate simulations. We focus on fully coupled ocean–atmosphere models that are needed to represent and understand the complicated interactions of these two systems, becoming increasingly important in climate change assessment in recent years. Numerical stability issues typically arise because of different time-stepping strategies applied to the coupled PDE system. In particular, the contributing factors include using large time steps, lack of accurate interface flux, and single-iteration coupling. We investigate the stability of the coupled ocean–atmosphere models for various interface conditions such as the Dirichlet–Neumann condition and the bulk interface condition, which is unique to climate modeling. By analyzing a simplified model, we demonstrate here how the parameterization of the bulk condition and other numerical and physical parameters affect the coupling stability and establish stability conditions for different coupling strategies.

97 MATHEMATICS AND COMPUTING↗