Cycle-to-cycle variations (CCV) in spark-ignited (SI) engines impose performance limitations and in the extreme limit can lead to very strong, potentially damaging cycles. Thus, CCV force sub-optimal engine operating conditions. A deeper understanding of CCV is key to enabling control strategies, improving engine design and reducing the negative impact of CCV on engine operation. This paper presents a new simulation strategy which allows investigation of the impact of individual physical quantities (e.g., flow field or turbulence quantities) on CCV separately. As a first step, multi-cycle unsteady Reynolds-averaged Navier−Stokes (uRANS) computational fluid dynamics (CFD) simulations of a spark-ignited natural gas engine are performed. For each cycle, simulation results just prior to each spark timing are taken. Next, simulation results from different cycles are combined: one quantity, e.g., the flow field, is extracted from a snapshot of one given cycle, and all other quantities are ta... [more]

Geothermal energy piles (GEPs) are an environmentally friendly heat exchange technology that dualizes the role of the structural foundation pile for load support and in meeting the building heating/cooling need. Energy loops made from high-density polyethylene, which allow heat carrier fluid circulation, are fitted into the pile foundation elements to extract or inject and store heat energy in the soil surrounding the pile. This paper reports the results of a numerical study investigating the long-term behaviour of a group of energy piles embedded in unsaturated soils (sand and clay) under continuous cyclic heating and cooling load. Additionally, two scenarios were investigated where: (1) the whole GEPs were heated and cooled collectively; (2) alternate piles were heated and cooled. It was found that the trend of temperature magnitude at all the observed locations decreases with time as a result of the continuous heating and cooling cycles. Furthermore, subjecting alternate GEPs to the... [more]

Battery Cell design and control have been widely explored through modeling and simulation. On the one hand, Doyle’s pseudo-two-dimensional (P2D) model and Single Particle Models are among the most popular electrochemical models capable of predicting battery performance and therefore guiding cell characterization. On the other hand, empirical models obtained, for example, by Machine Learning (ML) methods represent a simpler and computationally more efficient complement to electrochemical models and have been widely used for Battery Management System (BMS) control purposes. This article proposes ML-based ensemble models to be used for the estimation of the performance of an LIB cell across a wide range of input material characteristics and parameters and evaluates 1. Deep Learning ensembles for simulation convergence classification and 2. structured regressors for battery energy and power predictions. The results represent an improvement on state-of-the-art LIB surrogate models and indic... [more]

The estimation of the subsurface acoustic impedance (AI) model is an important step of seismic data processing for oil and gas exploration. The full waveform inversion (FWI) is a powerful way to invert the subsurface parameters with surface acquired seismic data. Nevertheless, the strong nonlinear relationship between the seismic data and the subsurface model will cause nonconvergence and unstable problems in practice. To divide the nonlinear inversion into some more linear steps, a 2D AI inversion imaging method is proposed to estimate the broadband AI model based on a broadband reflectivity. Firstly, a novel scheme based on Gaussian beam migration (GBM) is proposed to produce the point spread function (PSF) and conventional image of the subsurface. Then, the broadband reflectivity can be obtained by implementing deconvolution on the image with respect to the calculated PSF. Assuming that the low-wavenumber part of the AI model can be deduced by the background velocity, we implemented... [more]

To enhance grid reliability, weak points must be monitored. One of the weaknesses is the cable joints, which are prone to failure and can cause great losses from both a technical and economical point of view. Joints failures are usually caused by impurities unintentionally added during installation that cause partial discharges (PDs). In time, these discharges erode the insulation and generate treeing up to a destructive discharge between the conductor and the ground shield. For this reason, a method for the early detection of defects in joint installation and their online monitoring is required. A previously developed sensor was improved by adapting it for this purpose. It is based on the measurement of the induced current on a conductor due to a charge variation. It was experimentally tested on an actual joint in which defects were artificially introduced. Results show that the sensor is able to detect partial discharges. Moreover, a method for PD localization was developed. The firs... [more]

This paper proposes a model of a coupled inductor which takes into account the influence of frequency, temperature, and a constant component, IDC, of currents in the windings on the parameters of the considered element. A description of the model and methods of measuring parameters of the inductor using an impedance analyzer and a chamber for thermal measurements is given. The obtained results of measurements are compared with the results of calculations proving a satisfactory match.

Automatic proportional integral derivative control techniques are applied in a single-stage solar absorption cooling system, showing 3.8 kW (~1 ton) cooling capacity, with a coefficient of performance of 0.6 and −4.1 °C evaporator cooling temperature. It is built with plate heat exchangers as main components, using ammonia−water as the working mixture fluid and solar collectors as the main source of hot water. Control tuning was verified with a dynamical simulation model for a solution regarding mass flow stability and temperature control in the solar absorption cooling system. The controller improved steady thermodynamic state and time response. According to experimental cooling temperatures, the system could work in ranges of refrigeration or air-conditioning end-uses, whose operation makes this control technique an attractive option to be implemented in the solar absorption cooling system.

Battery electric vehicles offer many advantages in terms of performance and zero-emission pollutants, but their limited range for long-distance trips compromises their large-scale market penetration. The problem of range can be solved with a dense network of fast-charging stations and an increase in embedded battery capacity. Simultaneously, improvements in high-power charging point units offer range gains of hundreds of kilometers in a mere 20 min. One risk remains: The travel time depends on the availability of charging stations, which can drop during rush hours, due to long queues, or power grid constraints. These situations could significantly affect the user experience. In this paper, we presented an approach to coordinate EV charging station choices in the case of long-distance trips. This system relies on vehicle-to-infrastructure communications (V2X). The objective is to enhance the use of the infrastructure by improving the distribution of vehicles between the different chargi... [more]

A design method was developed for automated, systematic design of hydrokinetic turbine rotor blades. The method coupled a Computational Fluid Dynamics (CFD) solver to estimate the power output of a given turbine with a surrogate-based constrained optimization method. This allowed the characterization of the design space while minimizing the number of analyzed blade geometries and the associated computational effort. An initial blade geometry developed using a lifting line optimization method was selected as the base geometry to generate a turbine blade family by multiplying a series of geometric parameters with corresponding linear functions. A performance database was constructed for the turbine blade family with the CFD solver and used to build the surrogate function. The linear functions were then incorporated into a constrained nonlinear optimization algorithm to solve for the blade geometry with the highest efficiency. A constraint on the minimum pressure on the blade could be set... [more]

The dynamic two-phase flow in porous media was theoretically developed based on mass, momentum conservation, and fundamental constitutive relationships for simulating immiscible fluid-fluid retention behavior and seepage in the natural geomaterial. The simulation of transient two-phase flow seepage is, therefore, dependent on both the hydraulic boundaries applied and the immiscible fluid-fluid retention behavior experimentally measured. Many previous studies manifested the velocity-dependent capillary pressure−saturation relationship (Pc-S) and relative permeability (Kr-S). However, those works were experimentally conducted on a continuum scale. To discover the dynamic effects from the microscale, the Computational Fluid Dynamic (CFD) is usually adopted as a novel method. Compared to the conventional CFD methods solving Naiver−Stokes (NS) equations incorporated with the fluid phase separation schemes, the two-phase Lattice Boltzmann Method (LBM) can generate the immiscible fluid-fluid... [more]

Rotary heat exchangers have been widely used in paint shops, combustion power plants, and in heating, ventilation, and air conditioning systems in buildings. For these processes, many types of heat exchangers are available in the market: Tube-shell heat exchangers, plate heat exchangers, and rotary heat exchangers, among others. For the rotary heat exchangers, the problem is that there is no net present value method and lifecycle assessment method-based optimization found in the literature. In this work, we address this issue: An optimization is carried out with help of an empirically validated simulation model, a life-cycle assessment model, an economical assessment, and an optimization algorithm. The objective function of the optimization simultaneously considers economic and environmental aspects by using different CO2 pricing. Different CO2 pricing scenarios lead to different optimization results. The ambient air empty tube velocity va, 2.1 optimum was found at 1.2 m/s, which... [more]

Phase change materials (PCMs) have a large number of applications for thermal energy storage (TES) and temperature reduction in buildings due to their thermal characteristics and latent heat storage capabilities. The thermal mass of typical brick walls can be substantially increased using a suitable PCM primarily based on phase change temperature and heat of fusion for different weather conditions in summer and winter. This study proposed a novel dual-layer PCM configuration for brick walls to maintain human comfort for hot and cold climatic conditions in Islamabad, Pakistan. Numerical simulations were performed using Ansys Fluent for dual PCMs layered within a brick wall for June and January with melting temperatures of 29 °C and 13 °C. This study examined and discussed the charging and discharging cycles of PCMs over an extended period (one month) to establish whether the efficacy of PCMs is hindered due to difficulties in discharging. The results show that the combined use of both P... [more]

The problem of freshwater deficit in the last decade has progressed, not only in Africa or Asia, but also in European countries. One of the possible solutions is to obtain freshwater from drifting icebergs. The towing of large icebergs is the topic analyzed in various freshwater supply projects conducted in different zone-specific regions of the world. These projects show general effects of iceberg transport efficiency but do not present a detailed methodology for the calculation of their mass losses. The aim of this article is to develop the methodology to calculate the mass losses of icebergs transported on a selected route. A multi-agent simulation was used, and the numerical model to estimate the melting rate of the iceberg during its energy-efficient towing was developed. Moreover, the effect of towing speed on the iceberg’s mass loss was determined. It was stated that the maximum use of ocean currents, despite longer route and increased transport time, allows for energy-efficient... [more]

This article investigates the design, modeling, and fabrication of small-size (150 × 90 × 1.6 mm) broadband printed biconical antenna. The proposed antenna is intended for use a reference antenna for electromagnetic interference measurement inside the EMC chamber. The reflection coefficient (S11-parameter) is verified by modeling the equivalent circuit of the structure in terms of lumped elements. This structure offers a −10 dB impedance bandwidth (from 0.65 GHz to 2.3 GHz) with the tapered balun feeding method. Therefore, it has a high probability of estimating the electromagnetic waves emitted from several applications such as GSM, LTE, UMTS, 3G, Wi-fi, Bluetooth, ZigBee and more. The simulated standard antenna parameters are compatible with the measured parameters results. Furthermore, azimuth omnidirectional radiation pattern and well-realized gain (3.8 dBi) are achieved, reflecting good values of antenna factor compared to the commercial design.

In order to find the temperature field in the ground with a heat exchanger, it is necessary to determine temperature responses of the ground caused by heat sources and the influence of the environment. To determine the latter, a new model of heat transfer in the ground under natural conditions was developed. The heat flux of the evaporation of moisture from the ground was described by the relationship taking into account the annual amount of rainfall. The analytical solution for the equations of this model is presented. Under the conditions for which the calculations were performed, the following data were obtained: the average ground surface temperature Tsm = 10.67 °C, the ground surface temperature amplitude As = 13.88 K, and the phase angle Ps = 0.202 rad. This method makes it possible to easily determine the undisturbed ground temperature at any depth and at any time. This solution was used to find the temperature field in the ground with an installed slinky-coil heat exchanger tha... [more]

A power transformer is an essential device for stable and reliable power transfer to customers. Therefore, accurate modeling of transformers is required for simulation-based analysis with the model. The paper proposes an efficient and straightforward parameter estimation of power transformers based on sweep frequency response analysis (SFRA) test data. The method first develops a transformer model consisting of repetitive RLC sections and mutual inductances and then aligns the simulated SFRA curve with the measured one by adjusting parameters. Note that this adjustment is based on individual parameter impacts on the SFRA curve. After aligning the two curves, the final transformer model can be obtained. In this paper, actual single-phase, three-winding transformer model parameters were estimated based on field SFRA data, showing that SFRA curves simulated from the estimated model are consistent with the measured data.

This article discusses the optimization of railway transition curves, through the application of polynomials of 9th and 11th degrees. In this work, the authors use a 2-axle rail vehicle model combined with mathematically understood optimization methods. This model is used to simulate rail vehicle movement negotiating both a transition curve and circular arc. Passenger comfort is applied as the criterion to assess which transition is actually is the best one. The 4-axle vehicle was also used to verify the results obtained using the 2-axle vehicle. Our results show that the traditionally used in a railway engineering transition—3rd degree parabola—which is not always the optimum curve. This fact is especially valid for the longest curves, with lengths greater than 150 m. For such cases, the transition curves similar to standard curves of 9th and 11th degrees is the optimum ones. This result is confirmed by the use of the 4-axle vehicle.

The negative impact of the displacement currents on the operation of electromagnetic converters results in additional losses and faster insulation degradation, as well as the self-resonance phenomenon. Effective measurement of the dielectric displacement currents in converters is quite complex; thus, advanced simulation programs should be used. However, currently, they do not enable the analysis of the systems in terms of the displacement currents distribution. In order to elaborate an effective tool for analyzing the distribution of the displacement currents by means of the Finite Element Method, we have decided to supplement the well-known reluctance-conductance network model with an additional capacitance model. In the paper, equations for the linked reluctance-conductance-capacitance network model have been presented and discussed in detail. Moreover, we introduce in the algorithm the Harmonic Balance Finite Element Method (HBFEM) and the Fixed-Point Method. This approach enables u... [more]

Making an accurate estimate of the CO2 storage capacity before the commencement of a carbon capture and storage (CCS) project is crucial to the project design and feasibility investigation. We present herein a numerical modelling study on the CO2 storage capacity in depleted gas reservoirs. First, we show a simple volumetric equation that gives the CO2 storage capacity in a depleted gas reservoir, which considers the same volume of CH4 at reservoir pressure and temperature conditions produced from the reservoir. Next, the validity and the limitations of this equation are investigated using a numerical reservoir simulation with the various reservoir characteristics of reservoir heterogeneity, aquifer water encroachment, and rock compaction and its reversibility. Regardless of the reservoir heterogeneity, if a reservoir is subjected to a weak or moderate aquifer support, the volumetric equation provides an estimate of the CO2 storage capacity as structurally trapped gas within 1% of that... [more]

Airflow occurring in a ventilation duct is characterized by low velocity and hence low Reynolds number. In these conditions, either a laminar, transitional or turbulent flow will occur. Different flow conditions result in different values of the friction coefficient. To achieve the transitional flow in numerical simulation, a modified algebraic model for bypass transition (modified k−ω) was used. Numerical simulation was validated using Particle Tracking Velocimetry (PTV) in the circular channel. The modified algebraic model consists of only two partial differential equations, which leads to much faster calculation than the shear stress transport model. Results of the modified algebraic model are largely consistent with either the measurement and shear stress transport model considering laminar and transitional flow. Consistency slightly decreased in turbulent flow in relation to the model using shear stress transport method.

Lithium zirconate is a candidate material in the design of electrochemical devices and tritium breeding blankets. Here we employ an atomistic simulation based on the classical pair-wise potentials to examine the defect energetics, diffusion of Li-ions, and solution of dopants. The Li-Frenkel is the lowest defect energy process. The Li-Zr anti-site defect cluster energy is slightly higher than the Li-Frenkel. The Li-ion diffuses along the c axis with an activation energy of 0.55 eV agreeing with experimental values. The most favorable isovalent dopants on the Li and Zr sites were Na and Ti respectively. The formation of additional Li in this material can be processed by doping of Ga on the Zr site. Incorporation of Li was studied using density functional theory simulation. Li incorporation is exoergic with respect to isolated gas phase Li. Furthermore, the semiconducting nature of LZO turns metallic upon Li incorporation.

The short-term effects of temperature on deammonification sludge were evaluated in a laboratory-scale sequencing batch reactor (SBR). Mathematical modeling was used for further evaluations of different intermittent aeration strategies for achieving high and stable deammonification performance at decreasing temperatures. As for the biomass cultivated at high temperatures (e.g., 30 °C), a higher temperature dependency (the adjusted Arrhenius coefficient θ for 11−17 °C = 1.71 vs. θ for 17−30 °C = 1.12) on the specific anammox growth rates was found at lower temperatures (11−17 °C) in comparison with higher temperatures (17−30 °C). Further evaluations of recovering the nitrogen removal efficiency at decreasing temperatures with the mathematical model by modifying the intermittent aeration strategies (aeration frequency (F) and the ratio (R) between non-aerated (non-aer) phase and aerated (aer) phase durations) indicated that intermittent aeration with a prolonged non-aerated phase (e.g., R... [more]

The turbulent pulsating flow and heat transfer in straight and 90° curved square pipes are investigated in this study. Both experimental temperature field measurements at the cross-sections of the pipes and conjugate heat transfer (CHT) simulation were performed. The steady turbulent flow was investigated and compared to the pulsating flow under the same time-averaged Reynolds number. The time-averaged Reynolds number of the pulsating flow, as well as the steady flow, was approximately 60,000. The Womersley number of the pulsating flow was 43.1, corresponding to a 30 Hz pulsating frequency. Meanwhile, the Dean number in the curved pipe was approximately 31,000. The results showed that the local heat flux of the pulsating flow was greater than that of the steady flow when the location was closer to the upstream pulsation generator. However, the total heat flux of the pulsating flow was less than that of the steady flow. Moreover, the instantaneous velocity and temperature fields of the... [more]

The purpose of this study is to compare the load calculation results by a model using the air changes per hour (ACH) method and a model using an airflow network (AFN) and to ascertain what causes the difference between the two models. In the basic case study, the difference in the heat transfer distribution of the model in the interior space was investigated. The most significant difference between the two models is the heat transfer that results from infiltration. Parameter analysis was performed to investigate the relationship between the difference and the environmental variables. The result shows that the greater the difference is between the air temperature inside the balcony and the outdoor air temperature, and the greater the air flows from the balcony to the residential area, and the greater the heating and cooling load difference occurs. The analysis using the actual weather files of five domestic cities in South Korea rather than a virtual case shows that the differences are... [more]

Important factors controlling the effective utilization of geothermal energy are favorable reservoir properties of rock formations, which determine both the availability and the transfer opportunities of reservoir fluids. Hence, crucial to the successful utilization of a given reservoir is the preliminary recognition of distribution of reservoir parameters as it enables the researchers to select the prospective areas for localization of future geothermal installations and to decide on their characters. The objectives of this paper are analyses and discussion of the properties of quartz sandstones buried down to a depth interval from about 3000 to under 5000 m below surface. These sandstones belong to Ediacaran−Lowery Cambrian Łeba, Kluki and Żarnowiec formations. The source data from the Słupsk IG-1 provided the basis for 1D reconstruction of burial depth and paleothermal conditions as well as enabled the authors to validate of the results of 2D models. Then, porosity distribution with... [more]