The assessment of the deep geothermal potential is an essential task during the early phases of any geothermal project. The well-known “Heat-In-Place” volumetric method is the most widely used technique to estimate the available stored heat and the recoverable heat fraction of deep geothermal reservoirs at the regional scale. Different commercial and open-source software packages have been used to date to estimate these parameters. However, these tools are either not freely available, can only consider the entire reservoir volume or a specific part as a single-voxel model, or are restricted to certain geographical areas. The 3DHIP-Calculator tool presented in this contribution is an open-source software designed for the assessment of the deep geothermal potential at the regional scale using the volumetric method based on a stochastic approach. The tool estimates the Heat-In-Place and recoverable thermal energy using 3D geological and 3D thermal voxel models as input data. The 3DHIP-Cal... [more]

Engineering optimization is the subject of interest for many scientific research teams on a global scale; it is a part of today’s mathematical modelling and control of processes and systems. The attention in this article is focused on optimization modelling of technological processes of surface treatment. To date, a multitude of articles are devoted to the applications of mathematical optimization methods to control technological processes, but the situation is different for surface treatment processes, especially for anodizing. We perceive their lack more, so this state has stimulated our interest, and the article contributes to filling the gap in scientific research in this area. The article deals with the application of non-linear programming (NLP) methods to optimise the process of anodic oxidation of aluminium using MATLAB toolboxes. The implementation of optimization methods is illustrated by solving a specific problem from engineering practice. The novelty of this article lies i... [more]

According to the IPCC, a global carbon budget of 400 GtCO2 is required to limit the temperature rise to 1.5 °C with a 67% likelihood by 2050. The finance industry is increasingly committed to ambitious climate targets. In this article, we describe the detailed methodology and energy model architecture of a MATLAB-based integrated energy assessment model for industry-specific 1.5 °C pathways, with a high technical resolution of target parameters as key performance indicators (KPIs). The additionality of OECM 2.0 is the high technical resolution in terms of the level of detail of industry-specific energy demand and supply parameters that can be modeled—a prerequisite to define industry-specific KPIs. We found that a database of industry-sector-specific energy demands and energy intensities, with a consistent methodology, is required to improve the accuracy of calculations in future research. We supplement the technical documentation with the results for a transport scenario.

The development and utilization of marine renewable energy is an important measure for achieving energy conservation, emissions reduction and carbon neutrality. Ocean thermal energy is the most stable energy among all the types of marine renewable energy. This paper built a simulation model of an ocean thermal energy conversion system based on actual device specifications by Aspen and MATLAB and put forward a corresponding control strategy. The opening control signal of the control valve at the turbine inlet was the condenser inlet pressure in this paper, and the frequency control of the working fluid pump depended on the evaporating pressure and flow rate of the working fluid. This paper analyzed the key operating parameter changes of the system under different working conditions. According to the analysis results, the turbogenerator in this system was able to generate 50 kW power for about 8 months per year. The highest net output power of the Organic Rankine Cycle was 47.3 kW; the h... [more]

In this paper, we propose a sophisticated battery model for vanadium redox flow batteries (VRFBs), which are a promising energy storage technology due to their design flexibility, low manufacturing costs on a large scale, indefinite lifetime, and recyclable electrolytes. Primarily, fluid distribution is analysed using computational fluid dynamics (CFD) considering only half-cells. Based on the analysis results, a novel model is developed in the MATLAB Simulink environment which is capable of identifying both the steady-state and dynamic characteristics of VRFBs. Unlike the majority of published studies, the inherent characteristics of the flow battery, such as shunt current, ion diffusion, and pumping energy consumption, are considered. Furthermore, simplified charge transfer resistance (CTR) is taken into account based on electrochemical impedance spectroscopy (EIS) measurement results. The accuracy of the model was determined by comparing the simulation results generated by the equiv... [more]

This paper presents research on the model developed in the Matlab environment for simulating effects of overvoltage in an unmanned aerial vehicle (UAV) upon lightning discharge. They are based on transmittance obtained from voltage surge impulse measured in drone circuits. Overvoltage waveforms were measured at the input and output of different parts of the machine. It was then possible to calculate the transmittance of those (chosen) circuits. The motors, supply system, communication lines, and sensors were primarily tested. Both positive and negative polarization of the surge pulse were used and compared. The shape of pulse, is standardized by international norms for avionics tests (RTCA DO-160). The special surge generators were used to prepare the same repetition of each pulse (for all measurements). The simplified model of surge pulses propagation in drone circuits was prepared in Matlab. The differential between direct and step-by-step paths of pulse propagation in some connected... [more]

The growing production of plastic waste and improper dumping after use has become a worldwide challenge. This waste is a substantial source of petroleum and can be effectively converted into pyrolytic oil and other useful products. A statistical prediction of the rate constants is essential for optimizing pyrolysis process parameters, such as activation energy (Ea), frequency factor (Ao), temperature (T), and kinetic rate constants (k). In this research, we utilized Box−Behnken using RSM with Design Expert software to predict statistical rate constants at 500 °C and 550 °C. The efficiency of the predicted rate constants was investigated and compared to the findings of experimental rate constants extracted from the literature. At 500 °C, the estimated rate constants did not reveal a significant rise in the oil output since these constants promoted high gas yield. Compared to the experimental rate constants, statistically predicted rate constants at 550 °C demonstrated substantially high... [more]

Bubble size and its distribution are the important parameters which have a direct impact on mass transfer in bubble column reactors. For this, a new robust image processing technique was presented for investigating hydrodynamic aspects and bubble behavior in real chemical or biochemical processes. The experiments were performed in a small-scale bubble column. The study was conducted for the wide range of clear liquid heights and superficial gas velocities. However, a major challenge in image analysis techniques is identification of overlapping or cluster bubbles. This problem can be overcome with the help of the proposed algorithm. In this respect, large numbers of videos were recorded using a high-speed camera. Based on detailed experiments, the gas−liquid dispersion area was divided into different zones. A foam region width was found as inversely proportional to the clear liquid height. An entry region width was found as directly proportional to the clear liquid height. Hydrodynamic... [more]

This is the source code for the STORE (STochasic Omentum REsponse) model. The model is used to simulate how naive T-cells in the omentum will prime and multiply during the expansion period (8 days) after a T-cell vaccination in a mouse.

This is the matlab source code used in the following paper:

Christian DA, Adams TA II, Smith TA, Shallberg LA, Theisen DJ, Phan AT, Abraha M, Perry J, Ruthel G, Clark JT, Murphy KM, Kedl, RM, Hunter CA. cDC1 coordinate innate and adaptive responses in the omentum required for T cell priming and memory. Science Immunology 7, eabq7432 (2022).

This is fixed legacy code used for the paper for scientific auditing and reproduction purposes. See paper for documentation.

The paper presents the concept of a hybrid power system with additional energy storage to support electric vehicles (EVs) charging stations. The aim is to verify the possibilities of mutual cooperation of individual elements of the system from the point of view of energy balances and to show possibilities of utilization of accumulation for these purposes using mathematical modeling. The description of the technical solution of the concept is described by a mathematical model in the Matlab Simulink programming environment. Individual elements of the assembled model are described in detail, together with the algorithm of the control logic of charging the supporting storage system. The resulting model was validated via an actual small-scale hybrid system (HS). Within the outputs of the mathematical model, two simulation scenarios are presented, with the aid of which the benefits of the concept presented were verified.

Ion exchange resin is used to remove potentially corrosive impurities from coolant in the first circuit of a nuclear power plant. After one operational cycle, the used and unused resin in the mixed bed is discarded as solid waste. The aim of this work is to create a mathematical model to predict the operational cycle time of the mixed bed resin for reducing unused resin discharge. A partial differential equation (PDE) was set up with the conservation of matter. A finite difference method was used to solve the PDE. Matlab was the programming and calculating tool used in this work. The data from solution were obtained at different time and space nodes. The model was then verified experimentally using different ions on exchange columns. Concentrations of K+, Mn2+, and Cl- were calculated to verify the validation of the model by comparing it with experimental data. The calculated values showed good consistency with the experimental value.

While topics in the chemical engineering curriculum generally involve the macroscale, e.g., focusing on process design and optimizations, many of the grand challenges facing chemical engineers today involve molecular level chemistry. Having an appreciation for molecular-level phenomena will better prepare chemical engineers to handle problems in energy, water, food, catalysis, advanced materials, and other fields. In this presentation, we discuss homework/project problems that we have developed for use in the Material and Energy Balances, Separations, Kinetics, and Thermodynamics courses that explore molecular level phenomena and molecular-level simulation approaches. Specifically, these problems focus on intermolecular interactions, connections between molecular configurations and macroscopic observables, how molecular level properties control material performance, and simulation strategies that enable exploration of the molecular level. An aim of the assignments is to demonstrate sim... [more]

Source code for the case study presented in the paper "Approximating Nonlinear Relationships for Optimal Operation of Natural Gas Transport Networks". The case study involves solving the compressor fuel cost minimization problem (FCMP) on three simple natural gas networks. For each gas network three different formulations of the FCMP are tested: a common simplified FCMP model (FCMP_S), the novel approximation FCMP model (FCMP_N) that is developed in the paper, and a partially rigorous FCMP model (FCMP_PR) that models components of the model using their most rigorous calculations where feasible. The FCMP for each of these tests was optimized using GAMS, for which the code is provided. The accuracy of each of the three models was then assessed by comparing them to a rigorous simulation. The rigorous simulation was coded in Matlab and is provided, where separate files are used to calculate the rigorous gas pressure drop along a pipeline, and the energy input required for gas compression... [more]