Flux Calculations Research Articles

Self-forced inspirals with spin-orbit precession

We develop the first model for extreme mass-ratio inspirals with misaligned angular momentum and primary spin, and zero eccentricity—also known as quasispherical inspirals—evolving under the influence of the first-order in mass-ratio gravitational self-force. The forcing terms are provided by an efficient spectral interpolation of the first-order gravitational self-force in the outgoing radiation gauge. In order to speed up the calculation of the inspiral, we apply a near-identity (averaging) transformation to eliminate all dependence of the orbital phases from the equations of motion while maintaining all secular effects of the first-order gravitational self-force at postadiabatic order. The resulting solutions are defined with respect to “Mino time”; thus, we perform a second averaging transformation, so the inspiral is parametrized in terms of Boyer-Lindquist time, which is more convent for LISA data analysis. We also perform a similar analysis using the two-timescale expansion and find that using either approach yields self-forced inspirals that can be evolved to subradian accuracy in less than a second. The dominant contribution to the inspiral phase comes from the adiabatic contributions, so we further refine our self-force model using information from gravitational wave flux calculations. The significant dephasing we observe between the lower and higher accuracy models highlights the importance of accurately capturing adiabatic contributions to the phase evolution. Published by the American Physical Society 2024

Thin-Film Composite Membrane Compaction: Exploring the Interplay among Support Compressive Modulus, Structural Characteristics, and Overall Transport Efficiency.

Water scarcity has driven the demand for water production from unconventional sources and the reuse of industrial wastewater. Pressure-driven membranes, notably thin-film composite (TFC) membranes, stand as energy-efficient alternatives to the water scarcity challenge and various wastewater treatments. While pressure drives solvent movement, it concurrently triggers membrane compaction and flux deterioration. This necessitates a profound comprehension of the intricate interplay among compressive modulus, structural properties, and transport efficacy amid the compaction process. In this study, we present an all-encompassing compaction model for TFC membranes, applying authentic structural and mechanical variables, achieved by coupling viscoelasticity with Monte Carlo flux calculations based on the resistance-in-series model. Through validation against experimental data for multiple commercial membranes, we evaluated the influence of diverse physical parameters. We find that support polymers with a higher compressive modulus (lower compliance), supports with higher densities of "finger-like" pores, and "sponge-like" pores with optimum void fractions will be preferred to mitigate compaction. More importantly, we uncover a trade-off correlation between steady-state permeability and the modulus for identical support polymers displaying varying porosities. This model holds the potential as a valuable guide in shaping the design and optimization for further TFC applications and extending its utility to biological scaffolds and hydrogels with thin-film coatings in tissue engineering.

DUGKS-GPU: An efficient parallel GPU code for 3D turbulent flow simulations using Discrete Unified Gas Kinetic Scheme

This paper presents a parallel implementation of the Discrete Unified Gas Kinetic Scheme (DUGKS) on the GPU system using the CUDA Fortran and CUDA C++ programming languages. Firstly, we conducted an extensive revision of our original CPU-based code, resulting in a threefold decrease in memory usage. This new implementation is also paired with a novel approach to compute cell face flux using trilinear interpolation. It is shown analytically that the interpolation-based approach to flux calculation is more accurate compared to the one used in the original DUGKS. The initial simulation results using this new approach suggest that trilinear interpolation can reduce numerical errors on a coarse mesh. For example, in the case of the decaying Taylor-Green vortex flow at a 1283 mesh resolution, the relative numerical error in the energy dissipation rate at t⁎=2, using the spectral simulation result as the benchmark, is approximately 30% lower than that of the original implementation. The improved GPU DUGKS method is applied to laminar and turbulent flows in periodic and wall-bounded boundary configurations. A performance comparison of the GPU implementation is also presented and compared to the previous CPU implementation. A maximum speedup of 7.64x was achieved on a desktop-level GPU compared to a 32-core CPU. The strong scaling test, conducted on an eight-GPU node, demonstrated the efficient utilization of available multiple GPU resources by the code. Program summaryProgram Title: DUGKS-GPUCPC Library link to program files:https://doi.org/10.17632/yykv5s9g2n.1Developer's repository link:https://github.com/kairzhan/DUGKS-GPUCode Ocean capsule:https://codeocean.com/capsule/3b1e4f74-cdd3-4781-8923-8514d5923dfb/Licensing provisions: GNU General Public License 3Programming language: CUDA C++, CUDA FortranSupplementary material:Nature of problem: DUGKS-GPU is an advanced numerical code designed to accelerate the simulation of turbulent fluid flows through the application of the kinetic method known as the Discrete Unified Gas Kinetic Scheme (DUGKS). This computational tool comprises a CUDA Fortran version, tailored for a single GPU, as well as a CUDA C++ version developed for multi-GPU platforms.Solution method: The code employs a second-order finite-volume discretization of the discrete velocity Boltzmann equation with the BGK collision model. It can be used to simulate small to medium-scale problems on modern multi-GPU platforms with CUDA architecture.

The effects of changes in HITRAN and the water vapor continuum model on infrared radiative transfer calculations and remote sensing applications

This study examines the implications of recent updates to the High Resolution Transmission (HITRAN) database and the MT_CKD water vapor continuum model on infrared radiative transfer calculations and satellite remote sensing applications. Specifically, it assesses the impact of the latest versions of the HITRAN database (i.e. HITRAN2016 and HITRAN2020) and the MT_CKD WV continuum model (i.e. versions 3.2 and 4.1.1) on clear-sky infrared flux calculations and the simulation of satellite infrared channels, specifically those of the Geostationary Operational Environmental Satellite (GOES) Advanced Baseline Imager (ABI). The results indicate that updates to the MT_CKD model tend to reduce atmospheric opacity, while updates to the HITRAN database have the opposite effect. The most significant differences noted are attributed to updates in the MT_CKD, leading to an increase in upward flux of 0.13 Wm−2 for the TRO and a decrease in downward flux of 0.28 Wm−2 for the MLS profile. However, these impacts are not uniform across the spectrum, atmospheric height, or among different atmospheric profiles. Notably, changes in the WV continuum model have significant effects in the far-infrared region. The updates in HITRAN exhibit pronounced differences associated with changes in the spectroscopic parameters of H2O, CO2 and O3. The results of GOES’ ABI channels indicate that changes in HITRAN lead to an average brightness temperature decrease of 0.3 K for channel 12, which exceeds the channel’s sensitivity. This research underscores the importance of continuous updates and evaluations of spectroscopic databases and radiative transfer models to improve the accuracy of atmospheric remote sensing data.

Modelling of sap flux density of oak in a humid region in China

AbstractTranspiration plays a vital role in determining the watershed water cycle. However, we still have little knowledge of the characteristics of tree transpiration in the Hanjiang River Basin, which is the water source for the middle route of South‐to‐North water diversion project. Here, we measured sap flux density of oak trees (Quercus, the dominant species here) at the 10‐min resolution for 2 years and explored its response to the environmental conditions. The incoming short‐wave radiation and vapour pressure deficit well explained the variation of daytime sap flux density, and a statistical model was then proposed to calculate the daytime sap flux density correspondingly; then a nighttime sap flux density module was proposed based on the daytime sap flux density calculation. Sap flux density showed clear counter‐clockwise hysteresis response to incoming short‐wave radiation, and clockwise hysteresis response to vapour pressure deficit. Our sap flux density model can well reproduce the corresponding hysteresis response to incoming short‐wave radiation and vapour pressure deficit. This study unravelled the environmental controls of sap flux density of the oak trees in the Hanjiang River Basin, proposed an efficient model for the sap flux density simulation, provided important knowledge for understanding the corresponding forests' water use, which is of critical significance in determining the water availability for the middle route of the South‐to‐North water diversion project.

Subgrid-Scale Topographic Effects on Radiation for Global Weather Forecast Models

The incoming solar radiation arriving the Earth’s surface is strongly influenced by surface terrain. Conventional global weather forecast models, with grid scales of about 10 km, lack the resolution to accurately capture terrain-induced variations. We devised a new parameterization method to incorporate high-resolution subgrid-scale terrain data into grid-scale radiative flux calculations without averaging or smoothing topographic features. Utilizing a 15″ digital elevation model from the Shuttle Radar Topography Mission, we computed subgrid-scale data and transformed them onto the Korean Integrated Model’s cubed sphere grid using a Voronoi diagram to maintain geographical accuracy. The new scheme was initially evaluated through offline ideal tests and case studies. The results demonstrated that the scheme accurately captured the variations in downward shortwave flux, keeping the mean flux on a global scale nearly constant. The global mean flux difference in all skies was less than 0.01%. Statistical analyses demonstrated improved temperature and geopotential height predictions compared to reanalysis data. The anomaly correlation coefficient for East Asia at 850 hPa increased by 0.036 at 240 forecast hours. Overall, the anomaly correlation coefficient and root mean square error of geopotential height and temperature showed enhancements, particularly in the Northern Hemisphere and tropics. Importantly, the scheme introduces negligible additional memory and CPU requirements, making it suitable for both regional and global models. Only a 0.58% increase in CPU time was observed for the 10-day forecast.

Dry and wet periods determine stem and soil greenhouse gas fluxes in a northern drained peatland forest

Greenhouse gas (GHG) fluxes from peatland soils are relatively well studied, whereas tree stem fluxes have received far less attention. Simultaneous year-long measurements of soil and tree stem GHG fluxes in northern peatland forests are scarce, as previous studies have primarily focused on the growing season. We determined the seasonal dynamics of tree stem and soil CH4, N2O and CO2 fluxes in a hemiboreal drained peatland forest. Gas samples for flux calculations were manually collected from chambers at different heights on Downy Birch (Betula pubescens) and Norway Spruce (Picea abies) trees (November 2020–December 2021) and analysed using gas chromatography. Environmental parameters were measured simultaneously with fluxes and xylem sap flow was recorded during the growing season. Birch stems played a greater role in the annual GHG dynamics than spruce stems. Birch stems were net annual CH4, N2O and CO2 sources, while spruce stems constituted a CH4 and CO2 source but a N2O sink. Soil was a net CO2 and N2O source, but a sink of CH4. Temporal dynamics of stem CH4 and N2O fluxes were driven by isolated emissions' peaks that contributed significantly to net annual fluxes. Stem CO2 efflux followed a seasonal trend coinciding with tree growth phenology. Stem CH4 dynamics were significantly affected by the changes between wetter and drier periods, while N2O was more influenced by short-term changes in soil hydrologic conditions. We showed that CH4 emitted from tree stems during the wetter period can offset nearly half of the soil sink capacity. We presented for the first time the relationship between tree stem GHG fluxes and sap flow in a peatland forest. The net CH4 flux was likely an aggregate of soil-derived and stem-produced CH4. A dominating soil source was more evident for stem N2O fluxes.

Methods and metrics to assess the accuracy of heat flux data from a spark ignition IC engine

The measurement of in-cylinder heat transfer can be a valuable diagnostic tool for quantifying and improving IC engine performance. Increased heat transfer out of the combustion chamber negatively impacts overall efficiency, while heat transfer from the metal to the charge can increase knock propensity. The heat flux at various locations in an engine cylinder can be measured using heat flux probes consisting of two thermocouples – a surface thermocouple exposed to the changing in-cylinder temperatures, and a far field thermocouple which has a constant temperature at steady state operating conditions. The buildup of carbon deposits on the surface thermocouple can affect the heat flux measurements and lead to errors in the data. Heat flux measurements were performed on a proprietary single-cylinder research engine with instrumented heat flux probes in the cylinder head and cylinder liner. The engine was run rich with a high PMI fuel to expedite the buildup of carbon deposits. A metric was developed based on the output of the surface thermocouple to determine the cleanliness of the heat flux probes non-intrusively. Heat flux calculations were done using the FFT method and validated using the Cook-Felderman method. An uncertainty analysis was conducted on the heat flux data to verify that the negative heat flux observed was not due to errors in the distance between the two thermocouples or the measurement of temperature. Finally, to bolster confidence in the results, heat flux measurements were made for motoring conditions with varying coolant and intake air temperatures to gauge whether the trends in the output were in the expected direction.

Modeling of the dissolution phenomena of solid particles into an infinitely large medium from the analytical solutions of reaction–diffusion equations

Mathematical solutions of the concentration of dissolved component in surrounding fluid were obtained for analyses of the dissolution phenomena of spherical, cylindrical, and slab-type particles. Calculation of the surface flux enabled material balance equations to be derived for the prediction of transient particle size. For the spherical particle, the resulting differential equation was solved numerically to compare with the analytical solution. For the cylindrical particle, the Weber transform was applied to predict the radius of the infinitely long cylinder as a function of time. After obtaining an analytical solution for the thickness of the dissolving slab, the result was compared with those of the sphere and cylinder. The spherical particle exhibited the shortest complete dissolution time (CDT), due to having the greatest surface-to-volume ratio. When external film resistance was important for the mass transfer, the Biot number was included in the governing equation, and to predict the particle size during dissolution, the concentration in the fluid phase could be derived by the combination of variables, or the Weber-type transform. When first-order or zero-order reaction occurs in the surrounding fluid, the concentration could also be predicted by the Laplace transform, and the particle size could be calculated for the solid particles. Similar to the simple diffusion process, CDT was predicted to increase in the order sphere < cylinder < slab. The experimental data of the dissolution of polymeric fiber was compared with the modeling result to estimate the diffusivity and solubility factor, which varied with the pH of the surrounding fluid.

Seasonal variability of coastal pH and CO2 using an oceanographic buoy in the Canary Islands

Ocean acidification, caused by the absorption of carbon dioxide (CO2) from the atmosphere into the ocean, ranks among the most critical consequences of climate change for marine ecosystems. Most studies have examined pH and CO2 trends in the open ocean through oceanic time-series research. The analysis in coastal waters, particularly in island environments, remains relatively underexplored. This gap in our understanding is particularly important given the profound implications of these changes for coastal ecosystems and the blue economy. The present study focuses on the ongoing monitoring effort that started in March 2020 along the east coast of Gran Canaria, within the Gando Bay, by the CanOA-1 buoy. This monitoring initiative focuses on the systematic collection of multiple variables within the CO2 system, such as CO2 fugacity (fCO2), pH (in total scale, pHT), total inorganic carbon (CT), and other hydrographic variables including sea surface salinity (SSS), sea surface temperature (SST) and wind intensity and direction. Accordingly, the study allows the computation of the CO2 flux (FCO2) between the surface waters and the atmosphere. During the study period, stational (warm and cold periods) behavior was found for all the variables. The lowest SST values were recorded in March, with a range of 18.8-19.3°C, while the highest SST were observed in September and October, ranging from 24.5-24.8°C. SST exhibited an annual increase with a rate of 0.007°C yr-1. Warmer months increased SSS, while colder periods, influenced by extreme events like tropical storms, led to lower salinity (SSS=34.02). The predominant Trade Winds facilitated the arrival of deeper water, replenishing seawater. The study provided insights into atmospheric CO2. Atmospheric fCO2 averaged 415 ± 4 µatm (2020-2023). Surface water fCO2sw presented variability, with the highest values recorded in September and October, peaking at 437 µatm in September 2021. The lowest values for fCO2sw were found in February 2021 (368 µatm). From 2020 to 2023, surface water fCO2sw values displayed an increasing rate of 1.9 µatm yr-1 in the study area. The assessment of fCO2sw decomposition into thermal and non-thermal processes revealed the importance of SST on the fCO2sw. Nevertheless, in the present study, it is crucial to remark the impact of non-thermal factors on near-shallow coastal regions. Our findings highlight the influence of physical factors such as tides, and wind effect to horizontal mixing in these areas. The CT showed a mean concentration of 2113 ± 8 μmol kg-1 and pH at in-situ temperature (pHT,IS) has a mean value of 8.05 ± 0.02. The mean FCO2 from 2020 to 2023 was 0.34 ± 0.04 mmol m-2 d-1 (126 ± 13 mmol m-2 yr-1) acting as a slight CO2 source. In general, between May and December were the months when the area was a source of CO2. Extrapolating to the entire 6 km2 of Gando Bay, the region sourced 33 ± 4 Tons of CO2 yr-1.

The Relationship between Vertical Kinematic Eddy Heat Flux, Air Temperature and Turbulent Kinetic Energy in Atmospheric Boundary Layer: Baghdad City

Studying the heat flux in the boundary layer and understanding its relationship to the various atmospheric variables reflects positively on understanding the nature of turbulence in this layer and thus understanding the nature of the spread and movement of pollutants and the transmission and distribution of energy in this layer, also, the vertical heat flux is considered a significant influence on the movement of buoyancy and stability in the boundary layer and its effect on the horizontal wind movement, and therefore it is considered one of the important studies indirectly involved in the estimates of wind energy production, pollutant diffusion, and turbulence. This study involved the calculation of the Eddy Heat Flux and turbulent kinetic energy across Baghdad city. Our investigation revealed a correlation between the Eddy Heat Flux, temperature, and Turbulent Kinetic Energy (TKE). The observations have been made for wind speed with three components (u, v, w) and temperature by using a fast-response anemometer. As for the atmospheric pressure data, it was obtained from the automatic weather station located in the department of Atmospheric Sciences at Mustansiriyah University. The range of observations extended through the duration of thirty days for 24 hours from 1st July 2016 to 30th July 2016 every second. The maximum Eddy Heat Flux value was 0.092 J/(m2.s) at 10:00 on 23th July 2016, while the minimum Eddy Heat Flux value was -0.013 J/(m2.s) at 21:00 hour on 25th July 2016, the negative sign refer to a change in the direction of heat transfer,. It was also found that there is a positive relationship between the eddy heat flux and temperatures with a correlation coefficient of 0.93, as well as between the eddy heat flux and the turbulent kinetic energy with a correlation coefficient of 0.8.

Heterogeneous CPU-GPU parallelization for modeling supersonic reacting flows with detailed chemical kinetics

Accurate simulations of complex combustion phenomena associated with supersonic reacting flows require the use of detailed chemical kinetic mechanisms. However, detailed chemical mechanisms may consist of a large number of species and reactions resulting in extremely high computation cost. In order to accelerate the simulations of supersonic reacting flows with detailed chemical mechanisms, a heterogeneous CPU-GPU parallel algorithm and its portable software implementation are presented. The parallel algorithm was broken down into two parts: the CPU handled the computation of fluid dynamics while the GPU evaluated the chemical source terms and gas physical properties. The use of overlapping computations of chemical source terms on GPU and calculations of viscous flux on CPU is also presented. A study of performance tests was conducted. The performance results show that evaluating chemical source terms and gas physical properties on 2 GPUs are about 157.8× and 78.5× faster than running on the 16-core CPU when using the most complex mechanism on a grid of 3.3 million cells, respectively, resulting in an excellent speedup of the whole iteration up to 47. The significant performance improvement provided by the parallel algorithm can provide a significant perspective for designing heterogeneous CPU-GPU algorithms for applications in simulating supersonic reacting flows with detailed chemical kinetics. Program summaryProgram Title: OpenHurricaneCPC Library link to program files:https://doi.org/10.17632/m9pphg9cjj.1Licensing provisions: GPLv3Programming language: C++, CUDANature of problem: The computation cost in evaluating chemical source terms and theirs diagonal Jacobian matrices, and in calculating gas physical properties dominates the simulations of supersonic reacting flows with detailed chemical kinetic models. And the parallelization of these computations can yield great benefits from GPU acceleration.Solution method: The parallel implementation of evaluating the chemical source terms and theirs diagonal Jacobian matrices, and calculating gas physical properties were conducted on GPU to speed up the computation.

A Comprehensive Numerical Overview of the Performance of Godunov Solutions Using Roe and Rusanov Schemes Applied to Dam-Break Flow

As open channel simulations are of great economic and human significance, many numerical approaches have been developed, with the Godunov schemes showing particular promise. To evaluate, confirm, and extend the simulation results of others, a variety of first- and second-order FVMs are available, with Rusanov and Roe schemes being used here to simulate the demanding case of 1D and 2D flows following a dam break. The virtual boundary cells approach is shown to achieve a monotonic solution for both interior and boundary cells, and while flux computation is employed at boundary cells, a refinement is only rarely used in existing models. A number of variations are explored, including the TVD MUSCL-Hancock (monotone upwind scheme for conservation laws) numerical scheme with several slope limiters in a quest to avoid spurious oscillations. The sensitivity of the results to both channel length and the ratio of downstream to initial upstream water depth is explored using 1D and 2D models. The Roe scheme with a Van Leer limiter as a slope limiter is shown to be both fast and slightly more accurate than other slope limiters for this problem, but the Rusanov scheme with different slope limiters works well for 1D simulations. Significantly, the selection of an appropriate slope limiter is shown to be best based on the ratio of the downstream to upstream water depth. However, this study focuses on the special case where the ratio of the initial depth downstream to upstream of the dam is equal to or less than 0.5, and these outcomes are compared to theoretical results. The 2D dam-break problem is used to further explore first- and second-order methods using different slope limiters, and the results show that the Superbee limiter can be problematic due to an observed large dispersion in depth contours. However, the most promising approaches from previous studies are confirmed to deserve the high regard given to them by many researchers.

Motor magnetic field analysis using the edge-based smooth finite element method (ES-FEM)

This paper proposes a smooth finite element method (S-FEM) for efficient and accurate analysis of motor magnetic fields. The edge-based smooth finite element (ES-FEM) formulations are derived for two-dimensional triangular element meshes suitable for multi-material motor structures, and then a magnetic flux density calculation method appropriate for S-FEM to calculate nonlinear electromagnetic fields is introduced. In the calculation of motor magnetic field, ES-FEM has a proper softening effect on the stiffness matrix of the finite elements, which makes the calculation more accurate and more efficient than FEM. The magnetic flux intensity calculation method introduced in this paper is not only appropriate for ES-FEM to calculate the nonlinear electromagnetic field of motor, but also can improve the calculation accuracy of magnetic flux density. In the calculation of nonlinear magnetic field of motor, ES-FEM can also improve the calculation accuracy of magnetic field analysis compared with FEM. Several numerical examples demonstrate the improving effect of the ES-FEM discussed above on the motor's magnetic field calculation.

Non‐turbulent motion identified from properties of transport and its influence on the calculation of turbulent flux

AbstractTurbulent fluxes play a critical role in atmospheric science and are usually calculated from the eddy covariance system. However, the presence of motions of larger scale and the bias from observational instruments often affects the procurement of turbulent fluxes. Based on the Hilbert–Huang transform, the properties of transport can be defined and used to distinguish non‐turbulent motions from the observational data, and a new way of decomposing the variables is thereby put forward to reconstruct the turbulent flux series. To quantify the influence of non‐turbulent motions on the calculation of turbulent flux, the non‐turbulent motions extracted from the five‐level turbulence data of the Tianjin 255‐m meteorological tower from July 1 to August 31, 2017 are examined. The results reveal that the presence of non‐turbulent motion can lead to a universal overestimation of turbulent flux, and the degree of overestimation increases with the complexity and the intensity of non‐turbulent motion. According to the consequences above, an empirical relationship, together with the corresponding coefficients, is given to provide a guidance on the correction of turbulent fluxes in practical use, such as the simulation of atmospheric turbulence and the parameterization in meteorological and climate models.

A Conservative Solid-Gas Coupling Scheme with Overlapping Meshes for Combustion Simulation of Propellants with Voids

ABSTRACT A conservative coupling scheme with overlapping meshes for the combustion simulation of propellants with voids is established. The two main ingredients of this scheme are the conservative calculation of fluxes of mass, momentum, energy, and species and the robust treatment of voids. Specifically, for the discrete flux calculation, firstly the solid-gas interfaces are captured by the level-set method and then sampling points are generated between the interfaces at consecutive time steps, and these points are used to ensure the conservation of these quantities up to each grid of the gas domain. For the second ingredient, the void before exposure is treated as “solid” with properties of the air. Once an exposure condition is met, the topological change of the void is realized by reinitializing the level set function inside and on the boundary of the void. Overall, the scheme is conservative and robust, which is validated with experimental results.

A centered limited finite volume approximation of the momentum convection operator for low‐order nonconforming face‐centered discretizations

AbstractWe propose in this article a discretization of the momentum convection operator for fluid flow simulations on quadrangular or generalized hexahedral meshes. The space discretization is performed by the low‐order nonconforming Rannacher–Turek finite element: the scalar unknowns are associated with the cells of the mesh while the velocities unknowns are associated with the edges or faces. The momentum convection operator is of finite volume type, and its expression is derived, as in MUSCL schemes, by a two‐step technique: computation of a tentative flux, here, with a centered approximation of the velocity, and limitation of this flux using monotonicity arguments. The limitation procedure is of algebraic type, in the sense that its does not invoke any slope reconstruction, and is independent from the geometry of the cells. The derived discrete convection operator applies both to constant or variable density flows and may thus be implemented in a scheme for incompressible or compressible flows. To achieve this goal, we derive a discrete analogue of the computation (with the velocity, one of its component, the density, and assuming that the mass balance holds) and discuss two applications of this result: first, we obtain stability results for a semi‐implicit in time scheme for incompressible and barotropic compressible flows; second, we build a consistent, semi‐implicit in time scheme that is based on the discretization of the internal energy balance rather than the total energy. The performance of the proposed discrete convection operator is assessed by numerical tests on the incompressible Navier–Stokes equations, the barotropic and the full compressible Navier–Stokes equations and the compressible Euler equations.

Lessons Learned from the Updated GEWEX Cloud Assessment Database

AbstractSince the first Global Energy and Water Exchanges cloud assessment a decade ago, existing cloud property retrievals have been revised and new retrievals have been developed. The new global long-term cloud datasets show, in general, similar results to those of the previous assessment. A notable exception is the reduced cloud amount provided by the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) Science Team, resulting from an improved aerosol–cloud distinction. Height, opacity and thermodynamic phase determine the radiative effect of clouds. Their distributions as well as relative occurrences of cloud types distinguished by height and optical depth are discussed. The similar results of the two assessments indicate that further improvement, in particular on vertical cloud layering, can only be achieved by combining complementary information. We suggest such combination methods to estimate the amount of all clouds within the atmospheric column, including those hidden by clouds aloft. The results compare well with those from CloudSat-CALIPSO radar–lidar geometrical profiles as well as with results from the International Satellite Cloud Climatology Project (ISCCP) corrected by the cloud vertical layer model, which is used for the computation of the ISCCP-derived radiative fluxes. Furthermore, we highlight studies on cloud monitoring using the information from the histograms of the database and give guidelines for: (1) the use of satellite-retrieved cloud properties in climate studies and climate model evaluation and (2) improved retrieval strategies.

Flux coupling approach on an exchange grid for the IOW Earth System Model (version 1.04.00) of the Baltic Sea region

Abstract. In this article the development of a high-resolution Earth System Model (ESM) for the Baltic Sea region is described. In contrast to conventional coupling approaches, the presented model features an additional (technical) component, the flux calculator, which calculates fluxes between the model components on a common exchange grid. This approach naturally ensures conservation of exchanged quantities, a locally consistent treatment of the fluxes, and facilitates interchanging model components in a straightforward manner. The main purpose of this model is to downscale global reanalysis or climate model data to the Baltic Sea region as typically, global model grids are too coarse to resolve the region of interest sufficiently. The regional ESM consists of the Modular Ocean Model 5 (MOM5) for the ocean and the COSMO model in CLimate Mode (CCLM, version 5.0_clm3) for the atmosphere. The bi-directional ocean–atmosphere coupling allows for a realistic air–sea feedback that outperforms the traditional approach of using uncoupled standalone models, as typically pursued with the EURO-CORDEX protocol. In order to address marine environmental problems (e.g., eutrophication and oxygen depletion), the ocean model is internally coupled with the marine biogeochemistry model, ERGOM, set up for the Baltic Sea's hydrographic conditions. The regional ESM can be used for various scientific questions such as climate sensitivity experiments, reconstruction of ocean dynamics, study of past climates, and natural variability, as well as investigation of ocean–atmosphere interactions. Therefore, it can serve for a better understanding of natural processes via attribution experiments that relate observed changes to mechanistic causes.

CPU-GPU-coupled acceleration method for point flux calculation in Monte Carlo particle transport.

In Monte Carlo particle transport simulations, point flux tallying is a variance reduction technique that performs well with small detectors and finds broad application in source-detector problems and local point dose calculations. However, its use in large-scale point flux tallying computation adds substantial computational time. To address this issue, we propose a CPU-GPU-coupled acceleration method, which separates the complex logic and computationally intensive parts of particle transport calculation and assigns them to the CPU and GPU, respectively. This proposed method greatly enhances the efficiency of large-scale point flux tallies, providing significant convenience for subsequent dose calculations and other related steps. We validated our method by comparing the performance of a pure CPU program with our CPU-GPU accelerated program using the NUREG/CR-6115 PWR benchmark problem. The results indicate identical outcomes for photon point flux estimation, with the accelerated program being ~50 times faster.