Dynamic Numerical Model Research Articles

Research on autocannon firing dispersion based on bond space method

This paper is devoted to constructing a dynamic numerical model for the firing dispersion caused by the autocannon dynamic characteristics. The buffer dynamics and projectile-barrel coupling are considered. First, the simulation of autocannon firing dispersion using commercial software usually leads to expensive computational costs. Second, autocannon dynamic design includes multiple subsystem models with nonlinearities. The conventional design method makes it difficult to describe the dynamic response of such complex systems with a unified model. To end these, a dynamic junction bond space method is proposed for analyzing muzzle vibration and firing dispersion under continuous firing loads, where the gene expression programming (GEP) method is adopted to construct the surrogate model for the buffer flow field. For the coupling analysis of the projectile and barrel, the projectile load is applied at a moving junction, which coincides with the flexible node of the barrel. By this, the dynamic numerical model for autocannon firing dispersion is established, and then the system state equation is obtained for each time step. Moreover, an autocannon standing target shooting example is presented to demonstrate the validity of the proposed method; the results show that the firing dispersion from the bond space model is consistent with the test.

Numerical stabilization methods for level-set-based ice front migration

Abstract. Numerical modeling of ice sheet dynamics is a critical tool for projecting future sea level rise. Among all the processes responsible for the loss of mass of the ice sheets, enhanced ice discharge triggered by the retreat of marine-terminating glaciers is one of the key drivers. Numerical models of ice sheet flow are therefore required to include ice front migration in order to reproduce today's mass loss and to be able to predict their future. However, the discontinuous nature of calving poses a significant numerical challenge for accurately capturing the motion of the ice front. In this study, we explore different stabilization techniques combined with varying reinitialization strategies to enhance the numerical stability and accuracy of solving the level-set function, which tracks the position of the ice front. Through rigorous testing on an idealized domain with a semicircular and a straight-line ice front, including scenarios with diverse front velocities, we assess the performance of these techniques. The findings contribute to advancing our ability to model ice sheet dynamics, specifically calving processes, and provide valuable insights into the most effective strategies for simulating and tracking the motion of the ice front.

Multi-field and multi-scale dynamic numerical modeling of supercritical CO2 and fly ash mineralization

Fly ash-CO2 mineralization technology can effectively control gas–solid waste emissions and improve the climate environment. However, rationally describing the evolution mechanism of the fly ash-supercritical CO2 mineralization reaction and accurately predicting the efficiency of the mineralization reaction are of great significance for improving and evaluating the supercritical CO2 mineralization sequestration capacity. Fly ash micro- and nanopore multi-scale effects are the key to influencing the supercritical CO2 seepage-mineralization reaction mechanism. None of the current theoretical models consider the pore multi-scale effects, thus the validity and accuracy of predicting and describing the evolutionary mechanism of supercritical CO2 mineralization rates and efficiencies are questionable. To address this scientific issue, a multi-scale continuous pore structure theoretical model was established based on the structural characteristics of fly ash micro- and nanopores. Additionally, a new model of supercritical CO2 multi-field and multi-scale dynamic mineralization theory was developed by considering parameters such as temperature, pressure, and humidity. The accuracy of the model was verified through the mineralization experiments. The new model effectively describes the dynamic evolution mechanism of supercritical CO2 mineralized fly ash. The effects of parameters such as temperature, pressure, water vapor diffusion coefficient, supercritical CO2 permeability, and pore multi-scale effect on the efficiency and mineralization rate were investigated through numerical simulation. The results showed the following: The mineralization efficiency exhibited a parabolic increase with temperature, with 72℃ being the turning point for the growth rate of mineralization efficiency, and 82℃ being the optimal effective temperature to promote the mineralization efficiency of supercritical CO2. A 3.5-fold increase in supercritical CO2 pressure changes the mineralization efficiency by a maximum of only 0.0493 %. Increasing the water vapor diffusion coefficient by 5 orders of magnitude changes the mineralization efficiency by only 8.3 %. The pore attenuation coefficient increases by 3 orders of magnitude, and the supercritical CO2 mineralization efficiency changes by a factor of 2.27. A model that does not account for pore multi-scale effects overestimates the fly ash- supercritical CO2 mineralization efficiency by a factor of 2.47 at a permeability of 3 × 10-10 m2. Therefore, neglecting the pore continuum multi-scale effect can lead to a serious overestimation of the mineralization efficiency of supercritical CO2.

Numerical modeling and validation of dike-induced water flow dynamics using OpenFOAM

ABSTRACT The basic purpose of the present research is to numerically elucidate the flow around a single dike in four different geometric configurations (CN: no cylinders, C5: 5 cylinders, C10: 10 cylinders, C15: 15 cylinders), under conditions of a constant flow rate and subcritical flow. This involved substituting the impermeable dike with varying numbers of piles to validate the observed experimental flow patterns, particularly the conditions leading to reverse flow formation. OpenFOAM, an open-source algorithm, was utilized for the simulation, incorporating the volume of fluid (VOF) method. Accurately depicting the reverse flow formation with minimal numerical diffusion, a strategy involving multizone meshing and the application of mesh refinement zones was utilized. These refinements were particularly focused on the section between 1.8 m and 2.2 m within the numerical domain. The water surface profile in front section of the dike and the averaged velocity at five different locations, which were then compared with the numerical results. The numerically predicted streamwise velocity showed a percentage error ranging between 0.5% and 4.5%. The numerical results from this study are pivotal in dike design, aiming to mitigate accelerated wear during flood events and to protect both the dike head and the adjacent bank from high energy flow failures.

Numerical Modeling of Plasma Dynamics and Neutron Generation in Z-pinch at the ANGARA-5-1 Facility

Numerical Modeling of Plasma Dynamics and Neutron Generation in Z-pinch at the ANGARA-5-1 Facility

An Eigenvalue‐Based Framework for Constraining Anisotropic Eddy Viscosity

AbstractEddy viscosity is employed throughout the majority of numerical fluid dynamical models, and has been the subject of a vigorous body of research spanning a variety of disciplines. It has long been recognized that the proper description of eddy viscosity uses tensor mathematics, but in practice it is almost always employed as a scalar due to uncertainty about how to constrain the extra degrees of freedom and physical properties of its tensorial form. This manuscript borrows techniques from outside the realm of geophysical fluid dynamics to consider the eddy viscosity tensor using its eigenvalues and eigenvectors, establishing a new framework by which tensorial eddy viscosity can be tested. This is made possible by a careful analysis of an operation called tensor unrolling, which casts the eigenvalue problem for a fourth‐order tensor into a more familiar matrix‐vector form, whereby it becomes far easier to understand and manipulate. New constraints are established for the eddy viscosity coefficients that are guaranteed to result in energy dissipation, backscatter, or a combination of both. Finally, a testing protocol is developed by which tensorial eddy viscosity can be systematically evaluated across a wide range of fluid regimes.

Characterization and optimization of continuous ohmic thermal sterilization based on the development of a predictive computational toolbox

Continuous thermal processing (CTP) is a common method for sterilizing food. However, it can result in an uneven temperature distribution, which can lead to a varying degree of processing intensity. Ohmic heating (OH) can be advantageous in this regard, as it enables volumetric heating for more homogenous treatments. However, evaluating the processing intensity distribution inside the equipment for OH is challenging due to the complex interaction between electrical, mechanical and thermal phenomena. Furthermore, the comparison of OH and conventional heating treatments often lack a profound basis of comparable treatment intensity considerations. To gain a deeper mechanistic understanding of the technology, a numerical computational fluid dynamics model for the OH sterilization of a clear carrot juice from the heating region to the cooling process was developed. The model was validated with thermal and electrical measurements and showed an error rate below 2.5% in its prediction capacities. Moreover, the model was implanted for the validation of the products sterilization and compared to a conventional validation approach, reviling a 33.3% underestimation of the thermal load by conventional manners, which can lead to faulty sterilization of the food product. Additionally, the model was expanded to also be able to predict the microbial inactivation ratio of the system with an average error of 1.10±0.74%. In addition, results indicate that the numerical calculation of the F0 values and their validation with the microbial inactivation ratio have a notable potential for localization and evaluation of hotspots in OH simulations. Therefore, it can be seen as a promising step for establishing a foundation for computer-assisted optimization of CTP and targeted processing.

Experimental and numerical studies on moisture adsorption/desorption characteristics across the circular fin tube desiccant coated heat exchanger

To improve the energy transfer abilities, enhance the quality of air, and curtail the growth of microorganisms, the desiccant coated heat exchangers are a capable substitute compared to conventional heat exchangers due to their capability of decoupling the sensible and latent loads and utilizing low-grade thermal energy. The concept of desiccant coated heat exchangers (DCHE) can be employed in various heat exchange/transport systems like heat pumps, atmospheric water harvesters, adsorption chillers, and air conditioners. Thus, in this research work, a circular fin tube desiccant coated heat exchanger (DCCHE) has been proposed and designed and a dynamic numerical model is established for simulating the combined heat and mass transport occurring across the adsorption/desorption process of a DCCHE employing finite element approach. The moisture and temperature distribution and flow of fluid are considered in the transient numerical model covering several domains such as water flow, air flow, and the coated desiccant. The desiccant material utilized in the current analysis is silica gel. The validation study showed good agreement, and the maximum error between the experimental result and simulated model for the humidity ratio of air and temperature of outlet air were observed to be ± 8 % and ± 7 %, respectively. A parametric analysis is carried out to assess the effect of varying the different input parameters on DCCHE’s performance. The field emission scanning electron microscopy results concluded that a uniform coating is formed on the circular fin surface. The experimental examination on the sorption kinetics of the silica gel showed that the moisture uptake capability was 0.34 g/g. Moreover, the proposed circular fin tube desiccant coated heat exchanger has been compared with the conventional plate type fin tube DCHE to assess its performance in terms of thermal coefficient of performance and specific dehumidification energy.

A thermal fluid dynamic model for the melt region during the laser powder bed fusion of polyamide 11 (PA11)

Studying the physical characteristics of the melt region in polymer-based laser powder bed fusion (PBF-LB/P) offers good comprehension that is essential for quality control, process optimization, understanding of material behavior, and preserving process consistency. Among the few numerical modelling investigations of PBF-LB/P, this study focuses on 3D multiphysics simulation of melt pool morphologies, beginning with powder bed preparation incorporating the actual particle size distribution using the discrete element method (DEM). A thermal fluid dynamic numerical model is then developed based on the finite volume method (FVM) using Flow-3D Weld to simulate PBF-LB while using polyamide 11 (PA11) as the primary material. The cooling after sintering plays a crucial role for the formation of the melt region due to the high viscosity of the polymer material effectively resulting in creeping flow. Here, we propose a model for both sintering as well as the subsequent cooling step, a two-step model capable of simulating the finished structure accurately with relatively low computational resources. The model's accuracy is verified through a mesh dependency analysis and its predictive power via validation against internal single-track experiments providing dimension data on the melt region width and depth. The simulation results allow for a detailed exploration of the effects of laser settings on melt region morphologies and the occurrence of defects. The model's scope is further extended to multi-track simulation, revealing the connection between surface roughness and laser settings. In general, this study provides significant contributions to the comprehension of PBF-LB/P methodologies which in turn can be used for further process improvement.

Effects of windbreak types on aerodynamics of high-speed trains traversing from flat ground to semi-cutting and semi-embankment under crosswinds

Under the operation of strong crosswinds, the aerodynamic performance of high-speed trains (HSTs) will be seriously deteriorated when the transition section of flat ground and semi-cutting and semi-embankment (FGSCSE) is traversed, and the setting of windbreaks will help to slow down the impact of strong crosswinds on the trains. In this study, a three-dimensional coupled computational fluid dynamics numerical model to assess the aerodynamic performance of train–windbreak–FGSCSE–air system is developed. A comparative assessment is carried out to identify the variations in aerodynamic performance on the train carriage: no windbreak (NW), 50% ventilation windbreak (VW), and solid windbreak (SW), and the reasons for these variations are elucidated by examining the flow field structure's evolution. Furthermore, the operational safety of the train is discussed based on the indicator of wheel unloading ratio (fΔQ). Across the three distinct conditions, significant abrupt changes in aerodynamic load coefficients (ALCs) and the shedding of vortex structures are experienced by HSTs traversing the FGSCSE transition sections. Compared to the VW condition, the NW and SW conditions result in a greater number of shedding vortices on the leeward side and the tail of the train, and the VW condition results in the smallest magnitude of ALCs fluctuation. The power spectral density peak values of the aerodynamic loads follow the order: SW > NW > VW. Upon the train fully enters the subsequent operational environment, the VW condition has the smallest standard deviation of these coefficients. The standard deviations of CFy, CFz, CMx, CMy, and CMz for the head train in the VW condition are only 57.17% (46.81%), 55.85% (54.15%), 72.74% (34.62%), 57.99% (51.92%), and 44.60% (43.82%) of the corresponding values in the NW (SW) condition, respectively. In the NW, VW, and SW conditions, the fΔQ exceeds 0.9 when the wind speeds reach 30, 40, and 35 m/s, respectively. The windbreak with a ventilation rate of 30% performs the best, providing the most effective safety and stability for train operation.

Aerodynamic interference between road vehicle and train during meeting on a single-level rail-cum-road bridge under crosswind

A single-level rail-cum-road bridge is a relatively new type of bridge structure that accommodates both roadway and railway traffic on the same level. The aerodynamic interactions between road vehicles, trains, and the bridge deck under strong crosswind are significant, potentially posing safety risks to both vehicles and trains. This study aims to elucidate the aerodynamic interference between a road vehicle and a train during meeting on a single-level rail-cum-road bridge under crosswind. To achieve this, a three-dimensional, incompressible, unsteady Reynolds averaged Navier–Stokes method is utilized to simulate the meeting process between a van and a train on a prototype single-level rail-cum-road bridge under crosswind. Using an established computational fluid dynamics numerical model, the flow structure and aerodynamic loads of the van–train–bridge system under crosswind are studied. The results show that the auxiliary facilities on the bridge deck (such as pedestrian guardrails and anti-glare barriers) significantly alter the flow field around the van and train, affecting their aerodynamic characteristics. Additionally, the aerodynamic interference between the van and the train during their meeting under crosswind is substantial. The aerodynamic coefficients of the van and train, in terms of both magnitude and fluctuation, increase significantly with their driving speeds during the meeting process under crosswind. Furthermore, this aerodynamic interference intensifies as the lateral distance between the van and the train decreases during the meeting under crosswind.

The effect of fastener clip fatigue for high-speed railway on vehicle-track dynamic interaction: Numerical analysis and probabilistic evaluation

To investigate the vehicle-induced fatigue damage evolution of fastener clips and their impact on the vehicle-track (VT) system, a numerical dynamic model is established to achieve probabilistic evaluation. The “rain flow” counting method is adopted to describe the stress cycles of the hazard area of the clips, and the S-N curve is used to predict the fatigue life of the clips. The solution of the dynamic equation of motion is obtained by the time-domain Green's function method, which is an effective method to improve computational efficiency. A versatile iterative algorithm is adopted to solve the nonlinear system considering the damaged clips. In the probabilistic evaluation framework, a novel random field method is employed to depict the spatial variability of clips’ strength. The applicability of the model established in this work in aspects of expressing random damage of fasteners and evaluating the dynamic behavior of the system is demonstrated by different numerical examples. The results show that the dynamic responses of the train and track system significantly exacerbate as the number of train cycles increases, which is closely related to the variability and non-uniform damage of the fasteners. The probability of cumulative damage to the fasteners induced by the train follows a three-stage rule, and the fatigue damage of the non-uniform fastener system caused by the train is greater than that of the uniform situation. In addition, the short wave components of track irregularity will significantly exacerbate the damage to fasteners and therefore degrade the status of the railway line. As the driving speed increases, the damage to the fasteners also increases, but the impact of vehicle speed on the fastener damage is smaller than that of track irregularity, especially in the range of fatigue damage greater than 0.9. Therefore, maintaining the geometric smoothness of the track, controlling the uniformity of the mechanical properties of the fasteners, and reducing the driving speed are feasible measures to extend the service life of the fastener system, thereby ensuring the operation of high-speed trains.

Hydrodynamic mechanism for dynamical structure formation of a system of rotating particles

Based on the hydrodynamic mechanism, which takes into account the interaction of all particles, a numerical simulation of the formation of a dynamical structure in a viscous fluid was carried out. This structure is a result of the collective dynamics of rotating particles in the fluid. It is supposed that the particles have a magnetic moment and are driven into rotation by an external variable uniform magnetic field. The results of numerical modeling of collective dynamics are presented for three initial structures that can be formed by interacting dipole particles in the absence of an external magnetic field. Such equilibrium structures are a straight chain, a closed chain, and a periodic structure in the form of a flat system of particle chains. The rotation of particles sets the surrounding fluid in motion, whose flow creates hydrodynamic forces and moments that move the particles. The collective dynamics of a system of rotating particles leads to the formation of a new dynamical structure from the original one, and this new structure has its own characteristic features for each case considered. A qualitative comparison of the results of the dynamics for a particles’ system set in motion due to the action of an external moment or an external force is carried out. The proposed hydrodynamic mechanism for the formation of a dynamical structure as a result of the collective dynamics of a rotating particles’ system can be used to control structure formation in a liquid-particle system.

Study on fault-slip process and seismic mechanism under dynamic loading of hard roof fracture disturbance

The dynamic stress disturbance induced by mining is a significant factor triggering fault-slip. To better understand the impact of mining-induced hard roof fracture on fault-slip, combined with the geological conditions of a longwall face in contact with the fault, a dynamic and static numerical model was established using FLAC3D. The influence of dynamic loading from hard roof fracture on fault stability was investigated. The characteristics of changes in fault shear stress, shear displacement, slip area, and seismic moment under dynamic loading were analyzed. Simultaneously, numerical results were compared and analyzed against microseismic (MS) monitoring records in the mining, validating the accuracy of the numerical results. The results indicate that dynamic loading is a critical factor in inducing fault-slip. As a geological weak surface, the fault impedes the transmission of shock waves and causes their attenuation. Under dynamic loading, when the fluctuation of fault shear stress and friction coefficient exceeds the threshold, the fault-slip. During the fault-slip process, accompanied by a sharp increase in shear displacement, the seismic moment increases. Finally, the influence of changes in seismic parameters (e.g. longwall–fault distance, seismic energy, and seismic location) on the fault-slip mechanism was analyzed. This study contributes to a better understanding of the disturbance mechanism of mining-induced dynamic loading on fault-slip and provides a theoretical basis for fault-slip rock bursts.

Effect of microbial growth and electron competition on nitrous oxide source and sink function of hyporheic zones

The hyporheic zones (HZs) are key sites of the production of nitrous oxide (N2O), a potent ozone-depleting greenhouse gas. Denitrification is the primary process of N2O production in HZs, including four reduction steps (NO3−→NO2−→NO→N2O→N2). Electron competition occurs between the four reduction steps and can significantly impact the production of N2O. However, denitrification was typically considered as simplified two step reactions for investigating the release of N2O, neglecting the electron competition among the four-step reactions. Moreover, the N2O production/consumption patterns are regulated by both hydraulic and biogeochemical conditions in HZs. Dynamic microbial growth cannot only mediate the biogeochemical reactions, but also change hydraulic properties spatiotemporally by bioclogging. But microbial growth is rarely considered for investigating N2O dynamics of HZs. To assess these effects on hyporheic N2O dynamics and source-sink function, we establish a novel numerical model of N2O dynamics of HZs, coupling porous flow, reactive transport, electron competition, microbial growth and bioclogging. The results show that the weak electron competitiveness of N2O reductase results in a less allocation of electrons to the N2O reduction process, particularly in situations with limited carbon sources, thus increasing the release of N2O into the rives. Microbial growth significantly influences N2O release from HZs into rivers, increasing by more than two orders of magnitude on average compared to the model neglecting microbial dynamics. In contrast to the classical knowledges that HZs in coarse sediments tending to short residence time cannot act as sources of N2O, dynamic microbial growth obviously increases the potential for N2O release from HZs in coarse sediments to the rivers. The global Monte Carlo regional sensitivity analyses indicate that microbial biomass is the most critical factor determined the hyporheic source-sink function for N2O, followed by carbon oxidation rate and residence time. These are significantly different from previous knowledge that the residence time and oxygen/nitrogen uptake rate are the most sensitive parameters, which may lead to misunderstanding of the key controlling factors of N2O release from HZs. In addition, we propose a new Damköhler number (DaO2∗) of dissolved oxygen by multiplying the classical DaO2 with a dimensionless microbial modification factor for identifying N2O source-sink function of HZs, with DaO2∗ < 1 for N2O sink, while DaO2∗ > 1 for N2O source.

Dynamic Numerical Simulation and Transfer Learning-Based Rapid Rock Identification during Measurement While Drilling (MWD)

In constructing rapid rock identification models for measurement while drilling (MWD) via neural network methods, collecting actual drilling data to train the model is extremely time-consuming and labor-intensive. This requires extensive drilling experiments in various rock types, resulting in limited neural network training data for rock identification that covers a limited range of rock types. To suitably address this issue, a dynamic numerical simulation model for rock drilling is established that generates extensive drilling data. The input parameters for the simulations include torque, drill bit rotation speed, and drilling speed. A neural network model is then developed for rock classification using large datasets from dynamic numerical simulations, specifically those of granite, limestone, and sandstone. Building upon this model, transfer learning is appropriately applied to store the knowledge obtained in the rock identification based on the neural network model. Further training through transfer learning is conducted with smaller datasets obtained during actual drilling, making the model suitable for practical rock identification and prediction in the drilling processes. The neural network rock classification model, incorporating dynamic numerical simulation and transfer learning, achieves a prediction accuracy of 99.36% for granite, 99.53% for sandstone, and 99.82% for limestone. This reveals an enhancement in prediction accuracy of up to 22.94% compared to the models without transfer learning.

Dynamic subsidence in the Colorado basin, offshore Argentina, South Atlantic

The Colorado basin in the SW Atlantic is one of the key exploratory frontier basins along the Argentine shelf. The stratigraphy shows a typical rift basin, divided into two major seismic successions, a Jurassic-Lower Cretaceous synrift followed by an Upper Cretaceous to Cenozoic postrift. They represent, respectively, the breakup of Gondwana and passive margin evolution. While the synrift basin formation was associated with the development of structural grabens, the preservation of thick (thousands of meters) postrift sequences are still debated because they cannot be accounted for only by thermal subsidence. In this work, we conducted a multimethodological subsidence analysis across the Colorado basin, including part of the continental slope. We performed a 1D backstripping analysis on 10 hydrocarbon exploration wells and 32 pseudo-wells, covering the whole basin area and a uniform stretching model, which were supported using dynamic topography data from the study region. From the backstripping analyses, considering the six major postrift sedimentary successions described in the literature, we could recognize three subsidence stages. An initial rapid subsidence stage during the post-rift (125–83 Ma), which is restricted to the main depocenters, followed by a stage of generalized and low subsidence focused on the main axis of the basin (83–33 Ma), with an increase in the last 33 my (33–0 Ma). The stretching models exhibit a typical asymptotic curve, depicting significant subsidence at the onset of the postrift stage, gradually diminishing to minimal values within a few tens of millions of years. Upon comparing the backstripping curves with the stretching subsidence models, particularly during the final flat segments, a residual negative value emerges over time (locally exceeding −500 m), indicating considerable residual subsidence during the Cenozoic postrift evolution along the Colorado Basin. This could be interpreted as dynamic subsidence. This residual sinking in the Argentina shelf matches with dynamic topography numerical models in the study region as well as the occurrence of accumulated slabs in the mantle detected by seismic tomography studies. We propose that the accumulation of slabs along this segment of the South American subduction zone might be the drivers of dynamic subsidence, which assisted to the preservation of a large thickness of post-rift strata along the SW Atlantic margin.

Ammonia as a hydrogen carrier: Simulation study of the ammonia cracker on thermo-fluid stability and key factors for scale-up

A typical technical problem in a hydrogen society is the difficulty in securing the economic efficiency of hydrogen storage and long-distance transportation. Hydrogen carriers that can convert hydrogen into other compounds are required to overcome this issue. Ammonia has many practical advantages as a hydrogen carrier and is drawing attention as a suitable substance. Therefore, a cracking technology that converts ammonia into high-purity hydrogen is necessary. In this study, a numerical computational fluid dynamics model of an ammonia cracker was developed. The flow and heat transfer characteristics were analyzed using a numerical model. The results showed that a nonuniform temperature distribution caused local cracking of ammonia and adversely affected the hydrogen production process. This results of this study can be used to help design large-scale ammonia crackers and efficient hydrogen production processes and contribute to building infrastructure for hydrogen storage/transport/utilization systems.

State of art of hydrogen utilization for building sector and set-up with preliminary experimental results of 1 kWel solid oxide fuel cell installed in a nearly zero energy house

The diffusion on large scale of the hydrogen-fueled technologies could be a strategy for boosting the expected decarbonisation of the civil sector. However, there are still a lot of unanswered questions about the real feasibility and the problems of the integration in a pre-existing building-plants system. Despite some conflicting opinions, the research about the application of the hydrogen system in the building sector must continue to make the technology more attractive and easily available at low cost. Therefore, first of all, the paper analyses advantages and issues through a broad overview. The main limits of the available studies regard both the approach and the technology integration. Indeed, there are few data about in-field performance analysis since the simulative/numerical investigations are usually based on in-lab prototypes. The most diffused systems referred to integrated energy solutions with proton exchange membrane fuel cell technology and in any case the integration in very low energy buildings such as the nearly zero energy (nZEB) ones is not discussed. Then basing on the discussed research gap, the paper focuses on the set-up of a new installation of micro-combined cooling, heating and power unit based on solid oxide fuel cell technology fueled by green hydrogen. The complexity in the design of the integration with other innovative plant configurations is discussed considering as real case study a nZEB in Mediterranean climate. Preliminarily experimental results show that the combination of fuel cell and photovoltaic system could bring to positive energy with surplus energy, during autumn season, of about 23.8 kWh. Moreover, it is also proposed the validation of a dynamic numerical model developed by means of TRNSYS (Transient System Simulation Tool) that will be used for further sensitivity analysis.

Dynamic numerical modeling and performance optimization of solar and wind assisted combined heat and power system coupled with battery storage and sophisticated control framework

Incorporating distributed renewable energy sources into heating, power, and cooling systems is facilitating the drive toward intelligent building energy solutions. However, the inherent uncertainties associated with controllable renewable resources pose challenges for the thermo-economic scheduling of smart buildings. Nonetheless, the flexibility inherent in smart buildings can be leveraged through various market mechanisms. Hence, this research introduces a sustainable energy-building system driven by an autonomous solar/dish Stirling engine (SDSE) and wind turbine combined with a sophisticated control strategy and battery storage, which is designed to provide heat and power requirements. The proposed control strategy is also devised to optimize energy generation, ensuring prudent conservation and minimal waste, thereby amplifying overall efficiency and financial viability. The meticulous design of the system is accomplished through advanced MATLAB/Simulink® modeling. A thorough technical sensitivity analysis meticulously hones design parameters, revealing optimal operational thresholds. Simulation outcomes unveil a consistent Stirling engine (SE) efficiency, achieving a pinnacle of 35 %, while the SDSE attains over 28 %, respectively. The horizontal axis wind turbine, encompassing 100 kW and 500 kW modules, demonstrates power coefficients spanning from 0.18 to 0.09, with corresponding land area requirements ranging from 4.22 m2 to 21.10 m2, emphasizing the pivotal roles of module power and land area optimization. This paper also casts a spotlight on the environmental repercussions of the system, illustrating its potential to avert up to 30,000 kg of CO2 emissions per kWh/year. Moreover, the levelized cost of electricity ranges from 0.13 to 0.16 $/kWh, accompanied by an hourly cost that fluctuates between 3 $/h and 40 $/h, respectively. Conclusively, the developed modeling can be regarded as a significant stride in the realm of hybrid renewable energy systems, replacing the conventional photovoltaic/wind models with cutting-edge solar/wind configurations managed by a sophisticated control strategy and battery storage system.