Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated current”

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 397 records · Page 22

Data driven methods to recognize patterns in EIC weak-strong simulation

Beam-Beam simulations are currently being studied in preparation for future EIC experiments to study beam-beam effects and, in turn, maximize luminosity. Weak-strong methods are studied for single-particle dynamics during collision. 1 million macro-particles for 1 million turns are typically tracked, corresponding to only 10 seconds in the EIC. The goal of this study is to predict beam properties over the scale of hours. A potential solution focuses on using data-driven methods such as machine learning methods to analyze and extend the insights of the beam properties such as long-term nonlinear effects. This would aid in long-term predictions where results would be more efficiently acquired than a typical tracking simulation. Some limitations such as inaccurate predictions and spatial complexity are also discussed. These methods can then be applied to strong-strong simulations in the future studies.

Accelerator Physics↗

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems Using MOOSE-based Multiphysics Simulations

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of an Efficient Binaural Simulation for the Analysis of Structural Acoustic Data

Applying binaural simulation techniques to structural acoustic data can be very computationally intensive as the number of discrete noise sources can be very large. Typically, Head Related Transfer Functions (HRTFs) are used to individually filter the signals from each of the sources in the acoustic field. Therefore, creating a binaural simulation implies the use of potentially hundreds of real time filters. This paper details two methods of reducing the number of real-time computations required by: (i) using the singular value decomposition (SVD) to reduce the complexity of the HRTFs by breaking them into dominant singular values and vectors and (ii) by using equivalent source reduction (ESR) to reduce the number of sources to be analyzed in real-time by replacing sources on the scale of a structural wavelength with sources on the scale of an acoustic wavelength. The ESR and SVD reduction methods can be combined to provide an estimated computation time reduction of 99.4% for the structural acoustic data tested. In addition, preliminary tests have shown that there is a 97% correlation between the results of the combined reduction methods and the results found with the current binaural simulation techniques

Johnson, Marty E.↗

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗

Nominal and adversarial synthetic PMU data for standard IEEE test systems

GridSTAGE (Spatio-Temporal Adversarial scenario GEneration) is a framework for the simulation of adversarial scenarios and the generation of multivariate spatio-temporal data in cyber-physical systems. GridSTAGE is developed based on Matlab and leverages Power System Toolbox (PST) where the evolution of the power network is governed by nonlinear differential equations. Using GridSTAGE, one can create several event scenarios that correspond to several operating states of the power network by enabling or disabling any of the following: faults, AGC control, PSS control, exciter control, load changes, generation changes, and different types of cyber-attacks. Standard IEEE bus system data is used to define the power system environment. GridSTAGE emulates the data from PMU and SCADA sensors. The rate of frequency and location of the sensors can be adjusted as well. Detailed instructions on generating data scenarios with different system topologies, attack characteristics, load characteristics, sensor configuration, control parameters are available in the Github repository - https://github.com/pnnl/GridSTAGE. There is no existing adversarial data-generation framework that can incorporate several attack characteristics and yield adversarial PMU data. The GridSTAGE framework currently supports simulation of False Data Injection attacks (such as a ramp, step, random, trapezoidal, multiplicative, replay, freezing) and Denial of Service attacks (such as time-delay, packet-loss) on PMU data. Furthermore, it supports generating spatio-temporal time-series data corresponding to several random load changes across the network or corresponding to several generation changes. A Koopman mode decomposition (KMD) based algorithm to detect and identify the false data attacks in real-time is proposed in https://ieeexplore.ieee.org/document/9303022. Machine learning-based predictive models are developed to capture the dynamics of the underlying power system with a high level of accuracy under various operating conditions for IEEE 68 bus system. The corresponding machine learning models are available at https://github.com/pnnl/grid_prediction.

99 GENERAL AND MISCELLANEOUS↗

Influence of the Biosphere on Precipitation: July 1995 Studies with the ARM-CART Data

Ensemble sets of simulation experiments were conducted with a single column model (SCM) using the Goddard GEOS II GCM physics containing a recent version of the Cumulus Scheme (McRAS) and a biosphere based land-fluxes scheme (SSiB). The study used the 18 July to 5 August 1995 ARM-CART (Atmospheric Radiation Measurement-Cloud Atmospheric Radiation Test-bed) data, which was collected at the ARM-CART site in the mid-western United States and analyzed for single column modeling (SCM) studies. The new findings affirm the earlier findings that the vegetation, which increases the solar energy absorption at the surface together with soil and soil-moisture dependent processes, which modulate the surface, fluxes (particularly evapotranspiration) together help to increase the local rainfall. In addition, the results also show that for the particular study period roughly 50% of the increased evaporation over the ARM-CART site would be converted into rainfall with the Column, while the remainder would be advected out to the large-scale. Notwithstanding the limitations of only one-way interaction (i.e., the large-scale influencing the regional physics and not vice versa), the current SCM simulations show a very robust relationship. The evaporation-precipitation relationship turns out to be independent of the soil types, and soil moisture; however, it is weakly dependent on the vegetation cover because of its surface-albedo effect. Clearly, these inferences are prone to weaknesses of the SCM physics, the assumptions of the large-scale being unaffected by gridscale (SCM-scale) changes in moist processes, and other limitations of the evaluation procedures.

Sud, Y. C.↗

Projections of Declining Surface-Water Availability for the Southwestern United States

16 of the CMIP5 models had all the data needed for this work for at least one simulation that was continuous from 1950 to 2040. Details of the models analyzed here are provided in Table S1. The model data analyzed here are available at http://strega.ldeo.columbia.edu:81/expert/home/.naomi/.AR5/.v2/.historical:rcp85/.mmm16/ a. Assessing the climatology of the models Despite increases in horizontal resolution of many models compared to their CMIP3 counterparts none of these models can adequately resolve the topography of the south west United States, such as the Sierra Nevada and Rocky Mountains and the associated orographic precipitation. This requires that caution be used when interpreting the results presented here. To assess the ability of the models to simulate the current hydroclimate, in Figure S1 we show the observed (from the Global Precipitation Climatology Centre gridded rain gauge data, (1)) monthly climatology of precipitation and the same for all the models and the multimodel mean for the California-Nevada, Colorado headwaters and Texas regions. The GPCC data uses rain gauges only and interpolates to regular grids of which we used the 1◦ by 1◦ one. Details of the data set can be found in (2). While the models apparently overestimate precipitation in California and Nevada the seasonal cycle with wet winters and dry summers is very well represented. It is also possible that the rain gauge observations are biased low by inadequately sampling the higher mountain regions. How ever the models might also be expected to underestimate orographic precipitation due to inadequate horizontal resolution. The 25 models are also too wet in the Colorado headwaters region but correctly represent the quite even distribution though the year. The bimodal distribution of precipitation in Texas, with peaks in May and September, and the absolute amounts, are well modeled but with the September peak too weak. The positive precipitation bias translates into a positive runoff bias for the Colorado headwaters as also shown in Figure S1. Here the observed runoff values are taken from simulations of the Variable Infiltration Capacity (VIC) land surface-hydrology model (3) forced by observed meteorology (5) that were conducted as part of the North American Land Data Assimilation System project phase 2 ( (NLDAS-2), http://www.emc.ncep.noaa.gov/mmb/nldas/. Runoff for California-Nevada is better simulated but there is a positive bias over Texas despite no strong precipitation bias. To check whether regional climate models better simulate P and runoff in these regions we analyzed the historical simulation with the Regional Climate Model version 3 driven by the National Centers for Environmental Prediction-Department of Energy Reanalysis 2 available from the North American Regional Climate Change Assessment Program (http://www.narccap.ucar.edu). This model configuration retained these biases in P and runoff although they were reduced in amplitude. Given these varying biases we plot P and P − E changes in actual values but apply the simplest bias correction possible to the runoff and soil moisture values and show the modeled changes in terms of percentages of the 20th Century model climatologies. A thorough assessment of the simulation of North American climate in CMIP5 models is conducted in Sheffield at al. (North American Climate in CMIP5 Experiments. Part I: Evaluation of 20th Century Continental and Regional Climatology, manuscript submit ted to J. Climate, available at http://www.climate.noaa.gov/index.jsp?pg=./cpo pa/ mapp/cmip5 publications.html). Sheffield et al. analyze the climatology of precipitation, surface air temperature, low level winds, moisture fluxes, runoff etc. and conclude that the main features of the hydrological cycle, including characteristics of the atmospheric moisture balance and its seasonality, are captured in the CMP5 models subject to biases in total precipitation amounts. We chose to use all available models instead of selecting some and rejecting others based on an assessment of model realism. This is in accord with the suggestions of Mote et al. for CMIP3 (4) but future work needs to revisit this matter for the case of the CMIP5 ensemble.

Seager, Richard↗

Development of a Detailed Surface Chemistry Framework in DSMC

Many of the current direct simulation Monte Carlo (DSMC) codes still employ only simple surface catalysis models. These include only basic mechanisms such as dissociation, recombination, and exchange reactions, without any provision for adsorption and finite rate kinetics. Incorporating finite rate chemistry at the surface is increasingly becoming a necessity for various applications such as high speed re-entry flows over thermal protection systems (TPS), micro-electro-mechanical systems (MEMS), surface catalysis, etc. In the recent years, relatively few works have examined finite-rate surface reaction modeling using the DSMC method.In this work, a generalized finite-rate surface chemistry framework incorporating a comprehensive list of reaction mechanisms is developed and implemented into the DSMC solver SPARTA. The various mechanisms include adsorption, desorption, Langmuir-Hinshelwood (LH), Eley-Rideal (ER), Collision Induced (CI), condensation, sublimation, etc. The approach is to stochastically model the various competing reactions occurring on a set of active sites. Both gas-surface (e.g., ER, CI) and pure-surface (e.g., LH, desorption) reaction mechanisms are incorporated. The reaction mechanisms could also be catalytic or surface altering based on the participation of the bulk-phase species (e.g., bulk carbon atoms). Marschall and MacLean developed a general formulation in which multiple phases and surface sites are used and we adopt a similar convention in the current work. Microscopic parameters of reaction probabilities (for gas-surface reactions) and frequencies (for pure-surface reactions) that are require for DSMC are computed from the surface properties and macroscopic parameters such as rate constants, sticking coefficients, etc. The energy and angular distributions of the products are decided based on the reaction type and input parameters. Thus, the user has the capability to model various surface reactions via user-specified reaction rate constants, surface properties and parameters.

Gas↗

High Level Architecture (HLA) as a Research Tool for Advanced Air Mobility (AAM)

In this document, a NASA investigation into the High-Level Architecture (HLA) standard for distributed simulation as a means of furthering Advanced Air Mobility (AAM) research is illustrated along with a description of the evolution of distributed simulation leading to the HLA standard. HLA enables efficient definition and insertion of components and participants in simulation sessions. Novel concepts can be added or removed as the work evolves. A general history of the HLA standard is described to show past needs for improved methods of distributed simulation, particularly among participating laboratories in separate locations, as well as previous applications of HLA simulation to airspace research. These examples from the past illustrate how new and emerging concepts can be simulated by layering notional ideas over simulated current-day systems and capabilities in a controlled environment for examination. In the 1990s, the Federal Aviation Administration (FAA) and the MITRE Corporation performed joint research into then-emergent capabilities in the National Airspace System (NAS). In the early 2000s, the MITRE Corporation developed an open, extendable set of objects and protocols designed for global and collaborative aviation research called AviationSimNet, which has HLA as its foundation. This led to NASA and the FAA using AviationSimNet to connect laboratories for AAM research, resulting in the NASA/FAA Laboratory Integrated Test Environment (NFLITE). The underlying capabilities of HLA enabled testing a notional interaction between AAM and legacy ATC systems that resulted in AAM vehicles receiving beacon codes for identification on ATC screens in Class C airspace. The potential next steps in this work include defining and examining notional interactions between AAM and legacy ATC systems within Class B airspace.

UAM↗

On the shear-current effect: toward understanding why theories and simulations have mutually and separately conflicted

ABSTRACT The shear-current effect (SCE) of mean-field dynamo theory refers to the combination of a shear flow and a turbulent coefficient β21 with a favourable negative sign for exponential mean-field growth, rather than positive for diffusion. There have been long-standing disagreements among theoretical calculations and comparisons of theory with numerical experiments as to the sign of kinetic ($\beta ^u_{21}$) and magnetic ($\beta ^b_{21}$) contributions. To resolve these discrepancies, we combine an analytical approach with simulations, and show that unlike $\beta ^b_{21}$, the kinetic SCE $\beta ^u_{21}$ has a strong dependence on the kinetic energy spectral index and can transit from positive to negative values at $\mathcal {O}(10)$ Reynolds numbers if the spectrum is not too steep. Conversely, $\beta ^b_{21}$ is always negative regardless of the spectral index and Reynolds numbers. For very steep energy spectra, the positive $\beta ^u_{21}$ can dominate even at energy equipartition urms ≃ brms, resulting in a positive total β21 even though $\beta ^b_{21}\lt 0$. Our findings bridge the gap between the seemingly contradictory results from the second-order-correlation approximation versus the spectral-τ closure, for which opposite signs for $\beta ^u_{21}$ have been reported, with the same sign for $\beta ^b_{21}\lt 0$. The results also offer an explanation for the simulations that find $\beta ^u_{21}\gt 0$ and an inconclusive overall sign of β21 for $\mathcal {O}(10)$ Reynolds numbers. The transient behaviour of $\beta ^u_{21}$ is demonstrated using the kinematic test-field method. We compute dynamo growth rates for cases with or without rotation, and discuss opportunities for further work.

79 ASTRONOMY AND ASTROPHYSICS↗

Ulysses: A functional description and simulation software system

Current design tools for digital circuits and systems are not well-integrated among the behavioral, gate, and transistor levels of design. Ulysses is a prototype software system that consists of a description language, a description compiler, and a simulator that make no distinction among these levels. The language is uniform over the entire range of logical descriptions, the description is hierarchical with no fundamental restrictions on depth or mixing of levels, and the simulator is fully integrated with the description. The structure of the language, compiler, and simulator are described in terms of their relationships to the abstractions of physical systems that are made in order to create logical descriptions and models of behavior.

Griswold, T. W.↗

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗

Spectral form factor in sparse SYK models

We investigate the spectral form factor of the sparse Sachdev-Ye-Kitaev model. We use numerical methods to establish that at intermediate times the connected part of the spectral form factor is the dominant one. These connected contributions arise from fluctuations around the disconnected geometry, not from a new saddle point. A similar effect was previously conjectured in SYK but required a value of N out of reach of current numerical simulations.

2D Gravity↗

Athermal glass work at the nanoscale: Engineered electron-beam-induced viscoplasticity for mechanical shaping of brittle amorphous silica

Amorphous silica deforms viscoplastically at elevated temperatures, which is common for brittle glasses. The key mechanism of viscoplastic deformation involves interatomic bond switching, which is thermally activated. Here, for this study, we precisely control the mechanical shaping of brittle amorphous silica at the nanoscale via engineered electron–matter interactions without heating. We observe a ductile plastic deformation of amorphous silica under a focused scanning electron beam with low acceleration voltages (few to tens of kilovolts) during in-situ compression studies, with unique dependence on the acceleration voltage and beam current. By simulating the electron–matter interaction, we show that the deformation of amorphous silica depends strongly on the volume where inelastic scattering occurs. The electron–matter interaction via e-beam irradiation alters the Si–O interatomic bonds, enabling the high-temperature deformation behavior of amorphous silica to occur athermally. Finally, by systematically controlling the electron–matter interaction volume, it is possible to mechanically shape the brittle amorphous silica on a small scale at room temperature to a level comparable to glass shaping at high temperatures. The findings can be extended to develop new fabrication processes for nano- and microscale brittle glasses.

36 MATERIALS SCIENCE↗

A machine learning approach to predict thermal expansion of complex oxides

Although it is of scientific and practical importance, the state-of-the-art of predicting the thermal expansion of oxides over broad temperature and composition ranges by physics-based atomistic simulations is currently limited to qualitative agreements. We present an emerging machine learning (ML) approach to accurately predict the thermal expansion of cubic oxides with a dataset consisting of experimentally measured lattice parameters while using the metal cation polyhedron and temperature as descriptors. High-fidelity ML models that can accurately predict temperature- and composition-dependent lattice parameters of cubic oxides with isotropic thermal expansions have been successfully trained. The ML-predicted thermal expansions of oxides not included in the training dataset have shown good agreement with available experiments. The limitations of the current approach and challenges to go beyond cubic oxides with isotropic thermal expansion are also briefly discussed.

36 MATERIALS SCIENCE↗

Thickness threshold study of polyaniline-based lightning strike protection coating for carbon/glass fiber reinforced polymer composites

Different polyaniline (PANI) based lightning strike protection (LSP) coating thicknesses were tested against simulated lightning strike tests to evaluate the threshold thickness of the coating. A series of lightning strike tests were performed on carbon/glass fiber reinforced polymer (CFRP/GFRP) composites subjected to simulated lightning currents ranging from 40 to 100 kA. Five LSP-coated specimens were considered herein: two CFRP and three GFRP composites. A minimum of ~ 345 μm thick PANI-based LSP coating on CFRP and GFRP composites under a 40 kA peak current retained 83 % and 100 % residual strength, respectively. Additionally, GFRP specimens tested against 60 and 100 kA peak currents exhibited 100 % residual strength. Finally, this study proposes that a PANI-based LSP coating is suitable for protecting GFRP composites, while the careful design is necessary for CFRP composites to avoid any catastrophic failure.

36 MATERIALS SCIENCE↗

A New Ecosystem Model for Arctic Phytoplankton Phenology From Ice‐Covered to Open‐Water Periods: Implications for Future Sea Ice Retreat Scenarios

Abstract To predict ecosystem change in the Arctic Ocean, understanding Arctic phytoplankton phenology is essential. We develop a marine ecosystem model focusing on phytoplankton dynamics and the competition for multiple resources based on knowledge from in situ data obtained in the Chukchi Sea. The model is designed to include the unique ecological characteristics of the Arctic Ocean. To the best of our knowledge, this is the first model successfully simulating the current paradigm for Arctic phytoplankton phenology, including both under‐ice blooms and succession of phytoplankton groups. Sensitivity experiments show that the dominance of diatoms can be sustained by lateral transport of a high‐silicate water mass (Pacific Winter Water). Experiments for future scenarios show that sea ice retreat drives a decrease in large diatoms and an increase in the other relatively small non‐diatom groups. These results suggest that dominant phytoplankton groups can be shifted from diatoms to non‐diatoms under ongoing Arctic Ocean sea ice retreat.

Choi, Jang‐Geun↗