Conservation Equations Research Articles

Assisting and opposing mixed convective flow over a solid sphere at its bottom stagnation point with a constant surface heat, mass and motile microorganism flux

This study examined the laminar mixed convective flow around the lower point of stagnation for a solid sphere with a steady flux of gyrotactic microbial motility, mass and temperature. Biological convection phenomena have been considered in this study due to their intermixing traits and capacity in order to enhance mass transit in numerous biotechnology and biological systems. The first objective of this study is to examine and establish the mathematical formulation comprising partial differential equations for energy, momentum, mass conservation and mobile microorganism conservation balances. To do numerical calculations, the controlling partial differential equations are initially transformed by similarity transformations into a collection of interconnected non-linear ordinary differential equations. These equations are then solved by using the MATLAB Bvp4c scheme. The results indicate that several controlling parameters, particularly bioconvection parameters such as the Lewis number ( Le ), Bioconvection Lewis parameter ( Lb ) and Bioconvection Peclet number ( Pe ), have significant effects on flow transfer rates. A dual solution is observed in a certain region of mixed convection λ , where the mixed convection parameter is within the range of 0.41–0.65. The findings reveal significant results from dual solution analysis, highlighting the occurrence of unstable and stable phenomena for specific values of bioconvection parameters exclusively in the case of opposing flow. The assisting flow demonstrates only a single solution that is physically stable and unique. Observing lower stagnation point flow is essential for predicting and analysing complex fluid–microorganism interactions in various scientific and engineering applications such as designing spherical-shaped microbial fuel cell. Optimising flow conditions can enhance the efficiency of processes like wastewater treatment, where microorganisms play a crucial role in breaking down pollutants.

Three-point functions of higher-spin supercurrents in 4D N=1 SCFT: General formalism for arbitrary superspins

We analyze the general structure of the three-point functions involving conserved higher-spin “vectorlike” supercurrents Js(z)≔Jα(s)α˙(s)(z) in four-dimensional N=1 superconformal field theory. Using the constraints of superconformal symmetry and superfield conservation equations, we utilize a computational approach to analyze the general structure of the three-point function ⟨Js1(z1)Js2′(z2)Js3′′(z3)⟩ and provide a general classification of the results. We demonstrate that the three-point function is fixed up to 2 min(si)+2 independent conserved structures, which we propose to hold for arbitrary superspins. In addition, we show that the conserved structures can be classified as parity-even or parity-odd in superspace based on their transformation properties under superinversion. Published by the American Physical Society 2024

A new methodology to calculate interface thermal resistance between intumescent coating and steel substrate

Interface thermal resistance (ITR) plays a very important role in heat transfer and thermomechanical coupling response from intumescent coating to steel substrate in fire. However, the dynamic ITR between intumescent coating and steel substrate during burning is seldom reported. In this paper, a new methodology is presented to calculate dynamic ITR between intumescent coating and steel substrate under high incident heat flux based on the interface conservation equation. The new methodology focuses on measuring the temperature distribution and thermal conductivity of intumescent coating and steel substrate, with fewer variables. In this way, the calculation of ITR can be simplified. Firstly, the temperature variation function curves inside the intumescent coating and the steel substrate were fitted separately based on the temperature measurement data. Secondly, the interface temperatures (Tci,Tsi) were determined separately at the interface position by temperature variation function curves inside the intumescent coating or the steel substrate. Finally, the ITR calculation formula was deduced and the key parameters required for ITR calculation were obtained by cone calorimetry et al. An illustrative example was studied and validate the applicability of the new methodology.

Enhancing building energy efficiency: Numerical analysis of a hybrid thermoelectric ventilation system and earth-air heat exchange with photovoltaic-thermoelectric power integration for heating and cooling loads

The lack of studies on hybrid systems combining thermoelectric ventilation (TEV) and earth-air heat exchangers with a photovoltaic-thermoelectric generator (PV-TEG) as the power supply for the hybrid system is addressed in this research. This study aims to investigate the effectiveness and performance of such a combined system in providing the required thermal load of a building. To simulate the performance of TEV and PV-TEG systems, energy conservation equations for various components of the system were formulated for the quasi-state mode. The performance of the earth-air heat exchanger system was evaluated using the effectiveness-number of transfer units method. The resulting set of algebraic equations was then solved using a MATLAB code. Results indicated significant preheating (5.41–15.03 °C) and precooling (0.7–10.22 °C) effects during the daily hours of 15-February and 15-July, respectively. Additionally, the TEG component effectively reduced PV temperatures by 3.54 °C–16.99 °C in 15-February and 3.99 °C–21.54 °C in 15-July. The system also demonstrated a high monthly useful thermal energy output (1215–1584 kWh) in the cold months (December to March), contrasting with higher monthly useful thermal exergy (712–763 kWh) in the summer months. Furthermore, increasing the number of TEGs from 10 to 50 resulted in significant improvements in the yearly average energy and exergy proficiency indexes (PIya,en and PIya,ex) by 164.97 % and 459.82 %, respectively. The supply air mass flow rate and the length and depth of the earth-air duct were identified as the second and third most influential parameters affecting these indexes. Moreover, the system at CPV concentration ratio of 3 provided the highest PIya,en and PIya,ex of 1.512 and 14.65, respectively.

Quark Stars in $f(R,T)$ Gravity: Mass-to-Radius Profiles and Observational Data

Abstract In this paper, we explore the $f(R,T)$ gravity theory, which introduces a coupling between matter and curvature, through the simplest linear functional form $f(R,T)= R + 2\beta T$. We derive the modified Einstein field equations and conservation equations for this theory. We then apply this framework to study the structural properties of quark stars (QSs) composed of interacting quark matter, considering perturbative QCD corrections and color superconductivity. By solving the modified Tolman-Oppenheimer-Volkoff equations, we investigate the mass-radius relation, stability criteria, and energy conditions of QSs. Our results indicate that the $f(R,T)$ gravity significantly influences the properties of QSs, leading to deviations from General Relativity. The analysis shows consistency with recent observational data, suggesting that the modified gravity framework could provide viable models for the study of compact stars.

A space-time meshfree method for heat transfer analysis in porous material

Abstract Porous material batteries are a new type of battery technology that uses porous materials as electrode materials, with advantages such as long lifespan. In the present study, we manily propose a space-time meshfree method for heat transfer problems in porous material energy storage battery. The thermal behavior is demonstrated based on a numerical solution of the energy conservation equation. A simple space-time meshless collocation scheme, which is based on a new type space-time radial basis function, is presented to get the approximate solution of the energy conservation equation. The energy conservation equation was transferred into a system of linear algebraic equations. By introducing prescribed boundary conditions, the heat transfer analysis in porous material energy storage battery can be shown with numerical results. Numerical discussions show that the proposed meshless collocation method is appropriate to simulate thermal behavior in porous material energy storage battery, while the traditional finite element method failed to provide the reasonable results for the tested examples in this paper.

The analytical model of flame characteristics of hydrogen–air through wall and gas interaction analysis

AbstractIn this article, the effect of different boundary conditions and different thermal and physical properties of walls and gas on flame characteristics and stability of hydrogen–air mixture are investigated using an analytical method. This method solves the gas–wall energy equation, and the hydrogen mass conservation equations. The jump conditions are obtained by integrating the energy and mass equation into a small control volume around the flame. For validation of this model, the temperature distribution on the outer surface of the wall is compared with experimental data that show the maximum relative error of 3.5% for Q = 400 mL/min and 4.9% for Q = 200 mL/min. The maximum variation of gas temperature is nearly 6.5 times of wall temperature variation. The wall can be considered one‐dimensional for conventional wall materials with K > 10. For the existence of combustion inside the chamber, when the value of K is greater than 10, the Péclet number should also be considered greater than 10. In a constant equivalence ratio, increasing the medium temperature increases flame stability.

A transient network model for cross-regional layered fire smoke diffusion in high-rise buildings

At present, research on smoke diffusion in high-rise building fires mainly focuses on the single smoke layer propagation within a limited area, with less attention paid to the smoke propagation under the joint action of the fire floor and the vertical shaft area. Moreover, research on fire spread often focuses on the composition of combustible materials and building structures, neglecting natural ventilation. To address the above issues, we propose a transient network model for layered smoke diffusion (CLDTN) based on multi-region integration and complex control parameters. This model considers the smoke stratification mode of the fire floor affected by the fire source and ventilation, as well as the vertical stratification mode of the shaft affected by the stack effect, achieving cross-regional smoke stratification and diffusion modeling through smoke coupling. Then, a time-dependent transient network model is used to analyze furniture materials and open ventilation conditions of doors and windows, and a differential equation for energy conservation in multiple scene layers is constructed to analyze the dynamic behavior of smoke propagation.The comparative experiment of smoke propagation between CLDTN and the Fire Dynamics Simulator (FDS) shows that CLDTN is more efficient in computation time and the model is more accurate.

A unifying bubble-layer-based theoretical model for flow boiling and application to a rod bundle

Current analyses of flow boiling primarily utilize two−/three-dimensional numerical simulations or one-dimensional empirical correlations/models. To develop a quasi-two-dimensional theoretical model, a unifying bubble-layer-based theoretical model is proposed, which can be applicable to the entire flow boiling process, from the single-phase liquid flow to the dryout point. This model considers bubble behaviors in each region and heat and mass transfer at the interfaces of two regions, consisting of conservation equations of mass, momentum, and energy, supplemented by a few empirical correlations. A sensitivity analysis of these correlations is conducted, revealing that the slip ratio correlation and bubble condensation rate correlation have the most significant impacts. Using experimental data on void fraction, the most accurate combination of empirical correlations is identified, with a total absolute error of 27.9 %. The verification ranges are 0.1-8Mpa for pressure, 182-1700 kW/m2 for heat flux, 240-2200 kg/m2s for mass flux, 0.83–1.76 for Prandtl number, and 4.4e3–2.98e5 for Reynold number. Finally, this model is applied to the analysis of the NHR200-II rod bundle, revealing the variations of several two-phase parameters under micro-boiling conditions and the influence of operating conditions on the length proportion of different flow sections, serving as important references for analyzing the impact of boiling on reactor reactivity.

Existence of solutions to k‐Wave models of nonlinear ultrasound propagation in biological tissue

AbstractWe investigate models for nonlinear ultrasound propagation in soft biological tissue based on the one that serves as the core for the software package k‐Wave. The systems are solved for the acoustic particle velocity, mass density, and acoustic pressure and involve a fractional absorption operator. We first consider a system that incorporates additional viscosity in the equation for momentum conservation. By constructing a Galerkin approximation procedure, we prove the local existence of its solutions. In view of inverse problems arising from imaging tasks, the theory allows for the variable background mass density, speed of sound, and the nonlinearity parameter in the systems. Second, under stronger conditions on the data, we take the vanishing viscosity limit of the problem, thereby rigorously establishing the existence of solutions for the limiting system as well.

The Impact of Flow-Thermal Characteristics in Ship-Board Solid Oxide Fuel Cells

Hydrogen is increasingly recognized as a clean and reliable energy vector for decarbonization in the future. In the marine sector, marine solid oxide fuel cells (SOFCs) that employ hydrogen as an energy source have already been developed. In this study, a multi-channel plate-anode-loaded SOFC was taken as the research object. A three-dimensional steady-state computational fluid dynamics (CFD) model for anode-supported SOFC was established, which is based on the mass conservation, energy conservation, momentum conservation, electrochemical reactions, and charge transport equations, including detailed geometric shapes, model boundary condition settings, and the numerical methods employed. The polarization curves calculated from the numerical simulation were compared with experimental results from the literature to verify the model’s accuracy. The curved model was applied by enlarging the flow channels or adding blocks. Numerical calculations were employed to obtain the current density, temperature distribution, and component concentration distribution under the operating conditions of the SOFC. Subsequently, the distribution patterns of various physical parameters during the SOFC operation were analyzed. Compared to the classical model, the temperature of the curved model was reduced by 1.3%, and the velocities of the cathode and anode were increased by 4.9% and 5.0%, respectively, with a 2.42% enhancement in performance. The findings of this study provide robust support for research into and the application of marine SOFCs, and offer they insights into how we may achieving “dual carbon” goals.

Prediction of wellbore gas hydrate generation based on temperature-pressure coupling model in deepwater testing

With the advancement of oil and gas exploration and development technology, the development of the world's oil and gas resources is gradually advancing from land and shallow waters to deep waters. During the development of gas wells in deepwater regions, the temperature of the wellbore decreases from the bottom of the well to the wellhead, and therefore, common well conditions of high pressure, low temperature, and the presence of free water are encountered, which can lead to gas hydrate plugging of the pipeline and affect the transportation of natural gas. Simulation of temperature and pressure along the wellbore and identification of gas hydrate-generating zones in the wellbore can be used to take timely preventive measures and minimize losses. In this study, a wellbore gas hydrate generation and prediction model is developed by combining the mass, momentum, and energy conservation equations, and the model is solved by the fourth-order Runge-Kutta method, taking into account the temperature-pressure coupling in the wellbore. The effects of gas production, specific heat capacity, thermal conductivity, tubing size, test duration, and ground temperature gradient on the location of gas hydrate generation and the amount of inhibitor used were analyzed in a case study well in the South China Sea. Some suggestions for practical engineering operations, such as increasing the gas production in moderation and choosing smaller tubing sizes within the acceptable range, were given. The study's results can provide a theoretical basis for flow assurance for hydrocarbon production in the wellbore of deepwater gas wells.

A Numerical Framework of Simulating Flow-Induced Deformation during Liquid Composite Moulding

Fibre deformation (or shearing of yarns) can develop during the liquid moulding of composites due to injection pressures or polymerisation (cross-linking) reactions (e.g., chemical shrinkage). On that premise, this may also induce potential residual stress–strain, warpage, and design defects in the composite part. In this paper, a developed numerical framework is customised to analyse deformations and the residual stress–strain of fibre (at a micro-scale) and yarns (at a meso-scale) during a liquid composite moulding (LCM) process cycle (fill and cure stages). This is achieved by linking flow simulations (coupled filling–curing simulation) to a transient structural model using ANSYS software. This work develops advanced User-Defined Functions (UDFs) and User-Defined Scalers (UDSs) to enhance the commercial CFD code with extra models for chemorheology, cure kinetics, heat generation, and permeability. Such models will be hooked within the conservation equations in the thermo-chemo-flow model and hence reflected by the structural model. In doing so, the knowledge of permeability, polymerisation, rheology, and mechanical response can be digitally obtained for more coherent and optimised manufacturing processes of advanced composites.

A two-fluid multiphase flow model coupled with the population balance equations for the flow of a multidispersed particle system

SummaryWe present a two-fluid multiphase flow model that can be applied to the flow of polydisperse particle systems. In the model, the flow of polydisperse particle systems with different N particle sizes is considered only as two phases of fluid and particle phases, without considering the individual N phases, and the motion of particles with different diameters is represented by the motion characteristics of fluid and granular phases. The equations for the granular phase are obtained by averaging the Euler equations for the particles of different particle sizes. And the velocity of particles of different particle sizes is calculated by an explicit approximation, similar to determining the slip velocity in a mixed model of multiphase flow without solving the Eulerian momentum equation. Also, the conservation equation of particles of different particle sizes was obtained by adding the source term to the population balance equation, which is caused by the difference in the particle velocity of the particles with different particle sizes. The proposed method was compared with the results of experimental studies presented in the previous study to verify the validity of the method. The method has the advantage of not only visualizing the particle size fraction distribution but also of being stable in calculation and not increasing the computational effort significantly in studying the flow of polydisperse particle systems with particle size distribution.

Mixed volumes of networks with binomial steady-states

The steady-state degree of a chemical reaction network is the number of complex steady-states for generic rate constants and initial conditions. One way to bound the steady-state degree is through the mixed volume of an associated steady-state system. In this work, we show that for partitionable binomial chemical reaction systems, whose resulting steady-state systems are given by a set of binomials and a set of linear (not necessarily binomial) conservation equations, computing the mixed volume is equivalent to finding the volume of a single mixed cell that is the translate of a parallelotope. Additionally, we give a coloring condition on cycle networks to identify reaction systems with binomial steady-state ideals. We highlight both of these theorems using a class of networks referred to as species-overlapping networks and give a formula for the mixed volume of these networks.

Development of a smoothed particle hydrodynamics model for porous media flows with enhanced volume conservation and the revisit of the mass conservation equation

Porous media exist extensively in hydraulic and coastal engineering structures, while the modeling of wave/flow interaction with porous media remains challenging. This work develops a smoothed particle hydrodynamics (SPH) model for accurately simulating wave/flow interaction with porous media. The mass and momentum conservation equations incorporating the mixture theory are adopted. The resistant forces of the solid skeleton of porous media on fluid flows are described by the nonlinear empirical formula. The research contributions of the work lie in two aspects. First, two categories of mass conservation equations for porous media flow are revisited and analyzed to examine the influences of the local time derivative term of fluid volume fraction on simulation results. Second, the Volume Conservation Shifting scheme is, for the first time, introduced into SPH to enhance volume conservation for simulating porous media flows. The developed SPH model is validated by an analytical case of seepage flows in a U-tube with porous media and then applied to study four benchmark examples involving both saturated and unsaturated porous media, i.e., dam-break flow through a crushed stone dam, rapid seepage flow through a rockfill dam, solitary wave propagation over a porous seabed, and solitary wave propagation over a submerged porous breakwater. The morphological features and dynamic pressure heads of the porous media flows have been satisfactorily predicted, demonstrating the good accuracy and enhanced volume conservation of the developed SPH model.

A parallel compositional reservoir simulator for large-scale CO2 geological storage modeling and assessment

Storing CO2 in deep aquifers and depleted gas reservoirs is an effective way to achieve carbon neutrality. However, the numerical simulation of CO2 storage in these formations is challenging due to the complexity of gases-brine systems. The number of gas species included in the gases-brine fluid models of existing simulators cannot meet the rapidly evolving CO2 sequestration scenarios. To address this intricate issue, we developed a three-dimensional fully implicit parallel CO2 geological storage simulator (PRSI-CGCS) on distributed-memory computers based on our in-house parallel platform. This simulator uses a compositional fluid model with a diverse range of gas species, including CO2, C1 ~ C3, N2, H2S, as well as newly added gases H2 and O2, which may be encountered in geological CO2 storages. Besides, we provide more suitable scaling factors for different gases in the stability analysis bypassing (SAB) method to accelerate the gases-brine phase equilibrium calculations. PRSI-CGCS does not incorporate energy conservation equations, and salt precipitation or dissolution is also not considered. Numerical experiments show that our simulator is scalable, robust and validated to simulate large-scale CO2 storage problems with hundreds of millions of grid blocks on a parallel supercomputer cluster. Besides, after our modification on scaling factors, the SAB method can reduce the number of stability analyses by 61.39 % to 88.71 %, thereby reducing simulation time. Furthermore, case studies indicate that injecting O2 and H2 along with CO2 reduces the stability or capacity of CO2 storage and increases the pressure required for injection. However, this impact is not significant when the impurity content is less than 10 %.

Influence of concentration on thermophoresis of spherical aerosol particles within a Brinkman medium

Abstract We are examining the thermophoretic movement of a uniform mixture of spherical aerosol particles, all with the same properties, as they are situated within a porous material. These particles can have various thermal conductivity and surface characteristics. This analysis focuses on situations where the Péclet and Reynolds numbers are small. The influence of particle interactions is carefully considered by using a unit cell model, a well-established method known for its accurate predictions in the context of sedimentation for monodisperse suspensions of spherical particles. The porous medium is represented as a Brinkman fluid characterized by a Darcy permeability, which can be determined directly from experimental observations. This medium is considered to be uniform and isotropic, and the solid matrix is in thermal equilibrium with the fluid flowing through the voids of the medium. The Knudsen number is assumed to be low, enabling the description of fluid flow through the porous medium using a continuum model that includes temperature jump, thermal creep, frictional slip, and thermal stress slip at the aerosol particle’s surface. The conservation equations for energy and momentum are individually tackled within each cell. In this model, each cell represents a spherical particle enclosed by a concentric shell of surrounding fluid. The thermophoretic particle migration velocity is determined across different cases. We derive analytical expressions for this average particle velocity, expressing it in terms of the particle volume fraction. It is observed that different cell models yield somewhat varied results for particle velocity. Generally, with a fixed permeability parameter characterizing the porous medium, an increase in the thermal stress slip coefficient tends to decrease the normalized thermophoretic velocity across the different cell models. The results are in good agreement with the available data as documented in the existing literature. Additionally, a parallel examination of aerosol sphere sedimentation is provided.

Numerical Approach of the First Instability Appearance in Inclined Taylor–Couette System

The present study numerically investigates the three-dimensional forced and mixed convection heat transfer in an inclined Taylor–Couette system (η=0.5 and Γ=20) submitted to a radial temperature gradient. The main objective of this study is to determine the effect of the angle inclination duct on the thermal and dynamic fields. Several cases have been dealt with depending on the inclination angle to detect the critical Reynolds number in each case. The model of the conservation equations with their boundary conditions is numerically solved by the finite volume method with a second-order spatiotemporal discretization. The results show that, in forced convection, the effect of the inclination angle is inexistent on the velocity field. However, with the presence of buoyancy effects, which impact flow stability and the transition to turbulence, the inclination influences both velocity and temperature fields. It also shows that selecting the vertical position of the annulus is preferable to obtain hydrodynamic stability in mixed convection. At the same time, from the thermal point of view, it is preferable to select the horizontal position to get dynamic and thermal stability.

Inundation Processes with Active Sediment Transportation in the Floodplain of West Rapti River, Nepal

The present study discusses flood hazard characteristics in the lower reaches of the West Rapti River based on the results obtained from field surveys and numerical computations using depth averaged 2-D numerical models for flood flow and associated sediment transportation. To evaluate the inundation process with sediment erosion and deposition in the floodplain, a new erosion term was introduced into the mass conservation equations for suspended sediment and bed sediment. The results obtained from numerical computations indicated that field data on the spatial distributions of depths for inundation, and sediment erosion and deposition can be evaluated by the numerical model. Thus, numerical predictions were performed for the inundated areas, and the accumulated volumes of sediment erosion/deposition for floods with return periods of 50, 100, and 200 years, as a step towards damage assessment and risk assessment due to floods with active sediment transportation in the floodplains.