The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the penelope physics subroutines to C++, designed as an extension of the Geant4 toolkit, is presented. The Fortran code system penelope performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. Penelope implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of Geant4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original penelope programs. The combined code, named PenG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of Geant4.

Advanced Adiabatic Compressed Air Energy Storage design and modelling accounting for turbomachinery performance

Energy storage plants are going to become a strategic asset in electric grids. This statement is confirmed looking at the increasing shares of renewables composing the energy portfolio of several nations. Therefore the power demand and production mismatches, caused by the intermittent nature of renewables, must be reconciled. Many energy storage solutions are available but Advanced Adiabatic Compressed Air Energy Storage (AA-CAES) plants have potentials similar to pumped hydro systems (PHS). A physical model was developed in Matlab-Simscape to simulate the dynamics of AA-CAES plants, implementing temperature-dependent air properties, efficiency maps for turbomachinery and realistic power ramps. Furthermore, start-up and shut-down phases and energy consumption during idle periods were accounted for. The model embeds a 1D Fortran code to model the detailed behaviour of a packed-bed TES. The grid-to-grid performance of an AA-CAES plant was determined and the assumptions implemented to take into account real turbomachinery behaviour are presented.

Using Software Applications in Teaching of Slope-deflection Method in Subject Static Analysis of Constructions

Static analysis of building structures is an important engineering discipline. The presented paper deals with the teaching process of the slope-deflection method in the subject Static Analysis of Constructions at the Faculty of Civil Engineering of the Technical University in Košice. The slope-deflection method is the method for the analysis of statically indeterminate structures and serves to obtain their inner forces and deformations. For purpose of teaching this subject, teachers prepared programs in the Fortran code and with the help of the output of decisive values obtained from the programs. Teachers know if students have correct solutions of their tasks and can guide students to work to correct mistakes in their work.

SAM++: Porting the E3SM-MMF cloud resolving model using a C++ portability library

In this work, we demonstrate the process for porting the cloud resolving model (CRM) used in the Energy Exascale Earth System Model Multi-Scale Modeling Framework (E3SM-MMF) from its original Fortran code base to C++ code using a portability library. This porting process is performed using the Yet Another Kernel Library (YAKL), a simplified C++ portability library that specializes in Fortran porting. In particular, we detail our step-by-step approach for porting the System for Atmospheric Modeling (SAM), the CRM used in E3SM-MMF, using a hybrid Fortran/C++ framework that allows for systematic reproduction and correctness testing of gradually ported YAKL C++ code. Additionally, analysis is done on the performance of the ported code using OLCF’s Summit supercomputer.

Using Software Applications in Teaching of Slope-Deflection Method in Subject Static Analysis of Constructions

—Static analysis of building structures is an important engineering discipline. The presented paper deals with the teaching process of the slope-deflection method in the subject Static Analysis of Constructions at the Faculty of Civil Engineering of the Technical University in Košice. The slope-deflection method is the method for the analysis of statically indeterminate structures and serves to obtain their inner forces and deformations. For purpose of teaching this subject, teachers prepared programs in the Fortran code and with the help of the output of decisive values obtained from the programs. Teachers know if students have correct solutions of their tasks and can guide students to work to correct mistakes in their work.

Two-loop $$ \mathcal{O} $$((αt + αλ + ακ)2) corrections to the Higgs boson masses in the CP-violating NMSSM

We present our computation of the $$ \mathcal{O} $$ O ((αt + αλ + ακ)2) two-loop corrections to the Higgs boson masses of the CP-violating Next-to-Minimal Supersymmetric Standard Model (NMSSM) using the Feynman-diagrammatic approach in the gaugeless limit at vanishing external momentum. We choose a mixed $$ \overline{\mathrm{DR}} $$ DR ¯ -on-shell (OS) renormalisation scheme for the Higgs sector and apply both $$ \overline{\mathrm{DR}} $$ DR ¯ and OS renormalisation in the top/stop sector. For the treatment of the infrared divergences we apply and compare three different regularisation methods: the introduction of a regulator mass, the application of a small momentum expansion, and the inclusion of the full momentum dependence. Our new corrections have been implemented in the Fortran code NMSSMCALC that computes the Higgs mass spectrum of the CP-conserving and CP-violating NMSSM as well as the Higgs boson decays including the state-of-the-art higher-order corrections. Our numerical analysis shows that the newly computed corrections increase with rising λ and κ, remaining overall below about 3% compared to our previously computed $$ \mathcal{O} $$ O (αt(αt + αs)) corrections, in the region compatible with perturbativity below the GUT scale. The renormalisation scheme and scale dependence is of typical two-loop order. The impact of the CP-violating phases in the new corrections is small. We furthermore show that the Goldstone Boson Catastrophe due to the infrared divergences can be treated in a numerically efficient way by introducing a regulator mass that approximates the momentum-dependent results best for squared mass values in the permille range of the squared renormalisation scale. Our results mark another step forward in the program of increasing the precision in the NMSSM Higgs boson observables.

Mathematical and Numerical Calculation of the Interlayer Slip of a Two-Layer Glued Beam

This study analyses the spatial variation of the deflection, rotation, and slip at the interface of two-layer bamboo scrimber-concrete composite beam simply supported under uniform transverse load on its entire length. The Timoshenko model is considered, and the equations are solved using the analytical methods and direct numerical simulation based on finite differences schemes and a Fortran code. The results obtained show that it is a good agreement between the results from the mathematical calculation and those obtained from the numerical simulation.

Thermosolutal natural convection in a partly porous cavity with sinusoidal wall heating and cooling

Purpose This paper aims to consider numerical analysis of laminar double-diffusive natural convection inside a non-homogeneous closed medium composed of a saturated porous matrix and a clear binary fluid under spatial sinusoidal heating/cooling on one side wall and uniform salting. Design/methodology/approach The domain of interest is a partially square porous enclosure with sinusoidal wall heating and cooling. The fluid flow, heat and mass transfer dimensionless governing equations associated with the corresponding boundary conditions are discretized using the finite volume method. The resulting algebraic equations are solved by an in-house FORTRAN code and the SIMPLE algorithm to handle the non-linear character of conservation equations. The validity of the in-house FORTRAN code is checked by comparing the current results with previously published experimental and numerical works. The effect of the porous layer thickness, the spatial frequency of heating and cooling, the Darcy number, the Rayleigh number and the porous to fluid thermal conductivity ratio is analyzed. Findings The results demonstrate that for high values of the spatial frequency of heating and cooling (f = 7), temperature contours show periodic variations with positive and negative values providing higher temperature gradient near the thermally active wall. In this case, the temperature variation is mainly in the porous layer, while the temperature of the clear fluid region is practically the same as that imposed on the left vertical wall. This aspect can have a beneficial impact on thermal insulation. Besides, the porous to fluid thermal conductivity ratio, Rk, has practically no effect on Shhot wall, contrary to Nuinterface where a strong increase is observed as Rk is increased from 0.1 to 100, and much heat transfer from the hot wall to the clear fluid via the porous media is obtained. Practical implications The findings are useful for devices working on double-diffusive natural convection inside non-homogenous cavities. Originality/value The authors believe that the presented results are original and have not been published elsewhere.

Making legacy Fortran code type safe through automated program transformation

Fortran is still widely used in scientific computing, and a very large corpus of legacy as well as new code is written in FORTRAN 77. In general this code is not type safe, so that incorrect programs can compile without errors. In this paper, we present a formal approach to ensure type safety of legacy Fortran code through automated program transformation. The objective of this work is to reduce programming errors by guaranteeing type safety. We present the first rigorous analysis of the type safety of FORTRAN 77 and the novel program transformation and type checking algorithms required to convert FORTRAN 77 subroutines and functions into pure, side-effect free subroutines and functions in Fortran 90. We have implemented these algorithms in a source-to-source compiler which type checks and automatically transforms the legacy code. We show that the resulting code is type safe and that the pure, side-effect free and referentially transparent subroutines can readily be offloaded to accelerators.

Simulation of Dam Breaks on Dry Bed Using Finite Volume Roe-TVD Method

Dams are one of the most important hydraulic structures. In view of unrecoverable damages occurring after a dam failure, analyzing a dams’ break is necessary. In this study, a dam located in Iran is considered. According to adjacent tourist and entertainment zones, the breaking of the dam could lead to severe problems for the area and bridges downstream of the river. To investigate the issue, a numerical FORTRAN code based on the 2D finite volume Roe-TVD method on a fixed bed is provided to assess the effects of the dam break. Turbulence terms and dry bed conditions were considered in the code. A numerical wave tank (NWT) with a triangular barrier in the bed was numerically modeled and compared with analytical models to verify the capability of the code. Comparing numerical, experimental and analytical results showed that estimated water level and mass conservation in the numerical model is in good agreement with the experimental data and analytical solutions. The 2D approach used has reduced the cost of computing compared to a 3D approach while obtaining accurate results. The code is finally applied to a full-scale dam-break flood. Six KM of the natural river downstream of the dam, including two bridges, B1 and B2, is considered. Flood flow hydrographs and water level variations at bridges B1 and B2 are presented. The results denoted that bridges B1 and B2 will be flooded after 12 and 21 min, respectively, and are at risk of the potential break. Thus, it is necessary to announce and possibly evacuate the resort area alongside the dam in order to decrease losses.