Biohydrogen is a versatile energy carrier for the generation of electric energy from renewable sources. Hydrogenases can be used in enzymatic fuel cells to oxidize dihydrogen. The rate of electron transfer (ET) at the anodic side between the [NiFe]-hydrogenase enzyme distal iron-sulfur cluster and the electrode surface can be described by the Marcus equation. All parameters for the Marcus equation are accessible from DFT calculations. The distal cubane FeS-cluster has a three cysteine and one histidine coordination [Fe₄S₄](His)(Cys)₃ first ligation sphere. The reorganization energy (inner- and outer-sphere) is almost unchanged upon a histidine-to-cysteine substitution. Differences in rates of electron transfer between the wild-type enzyme and the all-cysteine mutant can be rationalized by a diminished electronic coupling between the donor and acceptor molecules in the [Fe₄S₄](Cys)₄ case. The fast and efficient electron transfer from the distal iron-sulfur cluster is realized by a fine-tuned protein environment which facilitates the flow of electrons. This study enables the design and control of electron transfer rates and pathways by protein engineering.

Biomass gasification acquired great interest over the past decades as an effective and trustable technology to produce energy and fuels with net-zero carbon emissions. Moreover, using biomass waste as feedstock enables to recycle organic wastes and to fit the circular economy goals thus reducing the environmental impacts of waste management. Even if many studies have been already carried on, this kind of process must still be investigated and optimized with the final aim to develop industrial plants for different applications, from the hydrogen production to net-negative emission strategies. Modeling and developing of process simulations became an important tool to investigate the chemical and physical behavior of the plant, allowing to make a first raw optimization of the process and to define heat and material balances of the plant as well as to define optimal geometrical parameters with cost and time effective approaches. The present review paper focuses on the main literature models developed until now to describe biomass gasification process, and in particular on kinetic models, thermodynamic models, and computational fluid dynamic models. The aim of this study is to point out the strengths and the weakness of those models, comparing them and indicating in which situation is better to use an approach instead of another. Moreover, theoretical shortcut models and software simulations, not explicitly addressed by prior reviews, are taken into account. For researchers and designers this review provides a detailed methodology characterization as a guide to develop innovative study or project.

Experiments to measure a single molecule/particle, i.e. an individual molecule/particle, at room temperature or under physiological conditions without immobilization, for example on a surface, or without significant hydrodynamic flow have so far failed. This failure has given impetus to the underlying theory of Brownian molecular motion towards its stochastic due to diffusion. Quantifying the thermodynamic jitter of molecules/particles inspires many and forms the theoretical basis of single-molecule/single-particle biophysics and biochemistry. For the first time, our simulation results for a live cell (cytoplasm) show that the tracks of individual single molecules are localized in Brownian motion, while there is fanning out in fractal diffusion (anomalous diffusion).

Abstract: Although experimental attempts to measure a single molecule/particle at room temperature or under physiological conditions without immobilization, for example on a surface, or without hydrodynamic flow have so far failed, this failure has given impetus to the underlying theory of Brownian motion towards its stochastic due to diffusion. Quantifying the thermodynamic jitter of molecules/particles inspires many and forms the theoretical basis of single-molecule/single-particle biophysics and biochemistry. For the first time, our simulation results for a live cell (cytoplasm) show that the tracks of individual single molecules are localized in Brownian motion, while there is fanning out in fractal diffusion (anomalous diffusion).

An extensive thermodynamic study of N-methylformamide (CAS RN: 123-39-7) and N,N-dimethylformamide (CAS RN: 68-12-2), is presented in this work. The liquid heat capacities of N-methylformamide were measured by Tian–Calvet calorimetry in the temperature interval (250 – 300) K. The vapor pressures for N-methylformamide and N,N-dimethylformamide were measured using static method in the temperature range 238 K to 308 K. The ideal-gas thermo-dynamic properties were calculated using a combination of the density functional theory (DFT) and statistical thermodynamics. A consistent thermodynamic description was developed using the method of simultaneous correlation, where the experimental and selected titerature data for vapor pressures, vaporization enthalpies, and liquid phase heat capacities and calculated ide-al-gas heat capacities are treated together to ensure overall thermodynamic consistency of the results. Resulting vapor pressure equation is valid from the triple point to the normal boiling point temperature.

The attitude dynamics and orbit dynamics of large rigid space structures are coupled with each other under gravity gradient, which may affect the stability of the large rigid space structures. In this paper, the gravity gradient stability of a large rigid space structure is studied under both small and large disturbances. Based on the rigid body dynamics and orbit dynamics, an accurate dynamic model without any linearization of the large rigid space structure is established via the natural coordinate formulation (NCF), which is able to describe the large overall motions of the structures. By using the generalized-α algorithm, the gravity gradient stability of the large rigid space structure is simulated and analyzed via various examples, including the influence of large disturbance angles, the positions at the stabilization and unstabilization regions. Finally, the relationship between spinning stability and gravity gradient stability is also investigated via a large spinning space structure.

A U.S. government health agency conducted epidemiological studies to evaluate whether exposures to drinking water contaminated with volatile organic compounds at U.S. Marine Corps Base Camp Lejeune, North Carolina, were associated with increased health risks to children and adults. These health studies required knowledge of contaminant concentrations in drinking water—at monthly intervals—delivered to family housing, barracks, and other facilities within the study area. Because concentration data were limited or unavailable during much of the period of contamination (1950s–1985), the historical reconstruction process was used to quantify estimates of monthly mean contaminant-specific concentrations. This paper integrates many efforts, reports, and papers into a synthesis of the overall approach to, and results from, a drinking-water historical reconstruction study. Results show that at the Tarawa Terrace water treatment plant (WTP) reconstructed (simulated) tetrachloroethylene (PCE) concentrations reached a maximum monthly average value of 183 micrograms per liter (ug/L) compared to a one-time maximum measured value of 215 ug/L and exceeded the U.S. Environmental Protection Agency’s current maximum contaminant level (MCL) of 5 ug/L during the period November 1957–February 1987. At the Hadnot Point WTP, reconstructed trichloroethylene (TCE) concentrations reached a maximum monthly average value of 783 ug/L compared to a one-time maximum measured value of 1,400 ug/L during the period August 1953–December 1984. The Hadnot Point WTP also provided contaminated drinking water to the Holcomb Boulevard housing area continuously prior to June 1972, when the Holcomb Boulevard WTP came on line (maximum reconstructed TCE concentration of 32 ug/L) and intermittently during the period June 1972–February 1985 (maximum reconstructed TCE concentration of 66 ug/L). Applying the historical reconstruction process to quantify contaminant-specific monthly drinking-water concentrations is advantageous for epidemiological studies when compared to using the classical exposed versus unexposed approach.

Based on DFT (density functional theory) using B3LYP/6-31G method, theoretical investigations applied to demonstrate the structural, electronic properties and stability of (N,Nˊ-(1,4-phenylene)bis(1,8-naphthalimide)) is more stable than the compound (N-(4-amino phenyl)1,8-naphthalimide) by (-1.2762 eV or -29.4299 Kcal.mol-1) depending on the values of HOMO, synthesis reaction of imide is spontaneous and endothermic at temperature 275˚C according to the values of ∆rS, ∆rG and ∆rH.

The discrete element method (DEM) coupled with computational dynamics (CFD) has been considered one of the most sensitive ways of studying fluidized beds. This paper presents gas-solid fluidized bed simulations by use of CFD-DEM. The two-phase interaction is calculated in a way of decoupling. Within the effective neighborhood of the particles, the kernel approximation method is used to calculate the local porosity, thereby calculating the drag force on the particles. At the grid scale, the momentum exchange coefficient is calculated using a two fluid model based method, thereby calculating the source term of the phase interaction in the fluid control equation. The simulated form, size and motion process of the big bubble are all consistent with experimental observations. The simulated solid volume fraction and relative pressure averaged in a horizontal plane at 45 mm above the bottom are in reasonable agreement with the experimental data. The fluctuation of the simulated bed height over time is similar to the experimental measurement results. The reproduction of these qualitative and quantitative results indicates that the proposed method effectively improves simulation performance and accuracy.

The poly lactic acid (PLA) is one of the most used bio-derived thermoplastic polymers in 3D and 4D printing applications. The determination of PLA surface properties is of capital importance in 3D/4D printing technology. (1) Background: The surface thermodynamic properties of PLA polymer were determined by using the inverse gas chromatography (IGC) technique at infinite dilution. The determination of the retention volume of polar and non-polar molecules adsorbed on the PLA particles filling the column allowed us to obtain the dispersive, polar, and Lewis’s acid-base surface properties at different temperatures from 40 °C to 100 °C. (2) Methods: the applied surface method was based on our recent model that used the London dispersion equation and the new chromatographic parameter function of the deformation polarizability and the harmonic mean of the ionization energies of the PLA polymer and organic molecules. (3) Results: The application of this new method led to the determination of the dispersive and polar free surface energy of adsorption of molecules on the polymeric material, as well as the glass transition and the Lewis acid-base constants. (4) Conclusions: four interval temperatures were distinguished showing four zones of the variations of the surface properties of PLA as a function of the temperature before and after the glass transition.

This paper proposes that intelligent processes can be completely explained by thermodynamic principles. They can equally be described by information-theoretic principles that, from the standpoint of the required optimizations, are functionally equivalent. The underlying theory arises from two axioms regarding distinguishability and causality. Their consequence is a theory of computation that applies to the only two kinds of physical processes possible—those that reconstruct the past and those that control the future. Dissipative physical processes fall into the first class, whereas intelligent ones comprise the second. The first kind of process is exothermic and the latter is endothermic. Similarly, the first process dumps entropy and energy to its environment, whereas the second reduces entropy while requiring energy to operate. It is shown that high intelligence efficiency and high energy efficiency are synonymous. The theory suggests the usefulness of developing a new computing paradigm called Thermodynamic Computing to engineer intelligent processes. The described engineering formalism for the design of thermodynamic computers is a hybrid combination of information theory and thermodynamics. Elements of the engineering formalism are introduced in the reverse-engineer of a cortical neuron. The cortical neuron provides perhaps the simplest and most insightful example of a thermodynamic computer possible. It can be seen as a basic building block for constructing more intelligent thermodynamic circuits.

Combining modern thermodynamics theory branches, including finite time thermodynamics or entropy generation minimization, constructal theory and entransy theory, with metallurgical process engineering, this paper provides a new exploration on generalized thermodynamic optimization theory for iron and steel production processes. The theoretical core is to thermodynamically optimize performances of elemental packages, working procedure modules, functional subsystems, and whole process of iron and steel production processes with real finite-resource and/or finite-size constraints with various irreversibilities toward saving energy, decreasing consumption, reducing emission and increasing yield, and to achieve the comprehensive coordination among the material flow, energy flow and environment of the hierarchical process systems. A series of application cases of the theory are reviewed. It can provide a new angle of view for the iron and steel production processes from thermodynamics, and can also provide some guidelines for other process industries.

The present work reports an experimental thermodynamic study of two nitrogen heterocyclic organic compounds, fenclorim and clopyralid, that have been used as herbicides. The sublimation vapor pressures of fenclorim (4,6-dichloro-2-phenylpyrimidine) and of clopyralid (3,6-dichloro-2-pyridinecarboxylic acid) were measured, at different temperatures, using a Knudsen mass-loss effusion technique. The vapor pressures of both crystalline and liquid (including supercooled liquid) phases of fenclorim were also determined using a static method based on capacitance diaphragm manometers. The experimental results enabled accurate determination of the standard molar enthalpies, entropies and Gibbs energies of sublimation for both compounds and of vaporization for fenclorim, allowing a phase diagram representation of the (p,T) results, in the neighborhood of the triple point of this compound. The temperatures and molar enthalpies of fusion of the two compounds studied were determined using differential scanning calorimetry. The standard isobaric molar heat capacities of the two crystalline compounds were determined at 298.15 K, using drop calorimetry. The gas phase thermodynamic properties of the two compounds were estimated through ab initio calculations, at the G3(MP2)//B3LYP level, and their thermodynamic stability was evaluated in the gaseous and crystalline phases, considering the calculated values of the standard Gibbs energies of formation, at 298.15 K.

For paramagnetic and elastic materials in the presence of quasi-static magnetic field, there is correspondingly quasi-static magnetic pressure, in contrast to the radiational electromagnetic pressure associated with electromagnetic momentum density. A thermodynamic pressure demagnetization is demonstrated. For the elastic dielectrics in the presence of quasi-static electrical field, its dual effect can be easily obtained by variable transformation.

Minimal rotating thermodynamic systems are addressed. Particle m placed into the rotating symmetrical double-well potential (bowl), providing binary logical system is considered. The condition providing the transfer of the particle from one frictionless half-well to another, and in this way possibility to record 1 bit of information is derived. The procedure of recording turns out to be irreversible; it is impossible to return the particle to its initial state under rotation about the same axis. The same rotating double-well system exerted to the thermal noise is considered. Minimal rotating thermal engine built of the rotating chamber, movable partition and the particle confined within the chamber is treated. Rotation of the system displaces the partition; thus, enabling erasing of one bit information. Erasing of 1 bit of information is due to the inertia (centrifugal force) acting on the partition. Isothermal expansion of the “minimal gas” expectedly gives rise to the Landauer bound. Compression of the “gas” with the rotation around the same axis is impossible and demands the additional axis of rotation. The interrelation between the possibility of recording/erasing information and the symmetry of the system is considered.

In Newtonian mechanics, time as well as space are perceived as absolute entities. In 2 Einstein’s special relativity, time is frame dependent. Time is also affected by gravitational field 3 and as the field varies in space, time also varies throughout space. In the present article, a 4 thermodynamic-based time is investigated. The entity is called "irreversibility", which is generated 5 when availability (also known as exergy) is destroyed. Since each thermodynamic system may 6 generate different amount of irreversibility, this quantity is system dependent. The time’s arrow is 7 automatically satisfied, since irreversibility generation always proceeds in one direction (toward 8 future). We have demonstrated that, like common time, irreversibility is frame dependent, and 9 affected by gravity in the similar manner as the common time. For this reason, we propose to 10 assign the entity irreversibility of the system as thermodynamic time. A possible application of the 11 thermodynamic time is an interpretation and managing of the aging of biological systems. The 12 metabolic efficiency is related to the irreversibility of the chemical processes and affect the aging of 13 the system. Our sensation of time-flow may be attributed to the flow of availability and destruction 14 of it through the living system. It is shown by other authors that entropy generation (equivalent to 15 irreversibility) is a parameter for the human life span. Since the thermodynamic time is based on a 16 concept of thermodynamics that is universal, further applications to other subjects, such as biological 17 clock, telomere, and cosmology are possible.

In the century and a half since Maxwell first conjured his “finite being” which Lord Kelvin subsequently dubbed a “daemon”, researchers have explored the connections between non-equilibrium thermodynamics, entropy, and information theory. In recent years various Thermodynamic Uncertainty Relations (TURs) have been derived to inform upon the relationship between the entropy production and the precision possible in thermodynamic machines and processes. In this paper the recently derived TURs are applied to the hypothetical thermodynamic economy as described by Parker [1]. The TURs define the lower bound on the total entropy production of the economy. Changes in this entropy production rate has consequences on the stability of the economic system and plays a central role in the business cycle. This new perspective has important implications for policy makers, researchers, and other economic actors.

In the previous work in Ref.[1], we presented a new interaction scheme that describes the interaction of thermodynamic systems. Using this interaction scheme, we showed that the interactions of dark matter and dark energy are chaotic. We also predicted that the results could be generalized to all particle and thermodynamic systems. In addition, we gave physical definitions of the concepts of chaos and self-organization. More importantly, we proposed a new law of physics based on interacting systems. In this study, we will consider two thermodynamic systems interacting under different conditions to prove the universality of the physical law we propose. We will show that these systems also have chaotic interaction dynamics.

Key features of glycogen metabolism in excitable tissues are not well-explained by current concepts. Glycogen stores in brain and skeletal muscle are generally considered to function as local glucose reserves, to be utilized during transient mismatches between glucose supply and demand; however, quantitative measures show that blood glucose supply is likely never rate-limiting for energy metabolism in either brain or muscle under physiological conditions. These tissues nevertheless do normally utilize glycogen during intervals of increased energy demand, despite the availability of free glucose, and despite the ATP cost of cycling glucose through glycogen polymer. This seemingly wasteful shunt can be explained by considering the effect of glycogenolysis on the amount of energy derived from ATP (ΔG’ATP). ΔG’ATP is diminished by elevations in Pi, such as occur at sites of rapid ATP hydrolysis and net phosphocreatine consumption. Glycogen utilization counters this effect by sequestering Pi in glycolytic metabolites (glycogen_n + Pi → glycogen_n-1 + glucose-1-phosphate → phosphorylated glycolytic intermediates), and thereby maintains the amount of energy obtained from ATP at sites of rapid ATP consumption. This thermodynamic effect may be particularly important in the narrow, spatially constricted astrocyte processes that ensheath neuronal synapses. This effect can also explain the co-localization of glycogen and cytosolic phosphocreatine in brain astrocytes, glycolytic super-compensation in brain when glycogen is not available, and aspects of exercise physiology in muscle glycogen phosphorylase deficiency (McArdle’s disease).

The Thermodynamic Formalism provides a rigorous mathematical framework to study quantitative and qualitative aspects of dynamical systems. At its core there is a variational principle and corresponding, in its simplest form, to the Maximum Entropy principle, used as a statistical inference procedure to represent, by specific probability measures (Gibbs measures), the collective behaviour of complex systems. This framework has found applications in different domains of scienThe Thermodynamic Formalism provides a rigorous mathematical framework to study quantitative and qualitative aspects of dynamical systems. At its core there is a variational principle and corresponding, in its simplest form, to the Maximum Entropy principle, used as a statistical inference procedure to represent, by specific probability measures (Gibbs measures), the collective behaviour of complex systems. This framework has found applications in different domains of science, in particular, has been fruitful and influential in neurosciences. In this article, we review how the Thermodynamic Formalism can be exploited in the field of theoretical neuroscience, as a conceptual and operational tool, to link the dynamics of interacting neurons and the statistics of action potentials from either experimental data or mathematical models. We comment on perspectives and open problems in theoretical neuroscience that could be addressed within this formalism.ce, in particular, has been fruitful and influential in neurosciences. In this article, we review how the Thermodynamic Formalism can be exploited in the field of theoretical neuroscience, as a conceptual and operational tool, to link the dynamics of interacting neurons and the statistics of action potentials from either experimental data or mathematical models. We comment on perspectives and open problems in theoretical neuroscience that could be addressed within this formalism.

We present a thermodynamic approach in modeling the evolution of the universe based on a theory that space consists of energy quanta, the spaceons. From wave-particle duality, they can be treated as an ideal gas. The model is similar to the Big Bang but without Inflation. It provides an insight into the nature of dark energy and dark matter, and an explanation for the accelerated expansion of the universe. The universe started from an atomic size volume of spaceons at very high temperature and pressure. Upon expansion and cooling, phase transitions occurred resulting in the formation of fundamental particles, and matter. These nucleate and grow into stars, galaxies, and clusters due to the action of gravity. From the cooling curve of the universe we constructed a thermodynamic phase diagram of cosmic composition, from which we obtained the correlation between dark energy and the energy of space. Using Friedmann’s equations, our model fits well the WMAP data on cosmic composition with an equation of state parameter, w= -0.7. The dominance of dark energy started at 7.25 x 10⁹ years, in good agreement with BOSS measurements. The expansion of space is attributed to a scalar quantum space field. Dark Matter is identified as a plasma form of matter similar to that which existed during the photon epoch, prior to recombination. The thermodynamics of expansion of the universe was adiabatic and decelerating during the first 7 billion years after the Big Bang; it accelerated thereafter. A negative pressure for Dark Energy is required to sustain the latter. This is consistent with the theory of General Relativity and the law of conservation of energy. We propose a mechanism for the acceleration as due to consolidation of matter forming Dark Energy Stars (DES) and other compact objects. The resulting reduction in gravitational potential energy feeds back energy for the expansion. Space will continue to expand and dark energy will undergo phase transition to a Bose-Einstein condensate, a superfluid form of matter. Self-gravitation can cause a bounce and start a new Big Bang. We show how the interplay of gravitational and space forces leads to a cyclic, maybe eternal, universe.

The historical idea of entropy as a property of a body has been reviewed and shown to arise from Clausius’s view of heat as motion. This view of heat, being intermediate between the now defunct idea of heat as substance and the modern view that heat represents an exchange of energy, implied that a body contains a definite quantity of heat. Heat, and therefore entropy, was thus considered a property of a body. These ideas led Clausius to develop his famous inequality and the idea that entropy always increases in irreversible, non-cyclic processes. In this paper, the notion of the entropy as a property of a body is examined in detail. A physical meaning is attached to the thermodynamic entropy by showing that a change in the function Q/T can be understood as representing a change in the number of ways that energy can be distributed among the degrees of freedom active in the system at any given temperature. The idea is illustrated with reference to solids and simple liquids. It is shown that the total thermodynamic entropy is, in general, less than the Boltzmann entropy, except for the case of a monatomic classical ideal gas in which the number of degrees of freedom is independent of temperature. Finally, entropy as a state function is discussed. It is argued that this is entirely mathematical in nature and that the entropy of a state represented in p-V space is not equivalent to a physical property of a physical system in the same thermodynamic state.

The article introduces the reader to the new thermodynamic potential of classical thermodynamics, called the J-potential. For the case of a homogeneous system, it vanishes, but it perfectly serves multiphase (heterogeneous) systems. The J-potential is especially important for the thermodynamics of capillary systems in colloidal science. The J-potential simplifies the derivation of the Neumann and Young equations, and has found application in the thermodynamics of thin films and the theory of strength of materials. In this work, for example, it was used to study the dependence of the chemical potential of a substance of a small bubble on the bubble size.

One open question remains regarding the theory of the generalized variational principle, that is, why the stress-strain relation still be derived from the generalized variational principle while the Lagrangian multiplier method is applied in vain? This study shows that the generalized variational principle can only be understood and implemented correctly within the framework of thermodynamics. As long as the functional has one of the combination $A(\epsilon_{ij})-\sigma_{ij}\epsilon_{ij}$ or $B(\sigma_{ij})-\sigma_{ij}\epsilon_{ij}$, its corresponding variational principle will produce the stress-strain relation without the need to introduce extra constraints by the Lagrangian multiplier method. It is proved herein that the Hu-Washizu functional $\Pi_{HW}[u_i,\epsilon_{ij},\sigma_{ij}]$ and Hu-Washizu variational principle comprise a real three-field functional.

Pathogen-host interaction (virus-host, bacteria-host and fungi-host) is analyzed from the perspective of biothermodynamics. A mechanistic model of pathogen-host interactions, which was earlier applied to virus-host interactions, was applied for the first time to bacteria-host and fungi-host interactions. Through a mechanistic model suggested in this research, it is possible to explain the phenomenon of tropism of microorganisms to develop infections in certain tissues. Gibbs energy of biosynthesis represents the driving force for growth of bacteria inside host organisms. A growing bacterial colony represents a thermodynamic system, while the host organism is its surroundings. The permissiveness coefficient determines the possibility of growth of bacteria in an appropriate environment (tissue).

Systems biology has been established for more than a decade in the post-genomic era. With the help of the computational and mathematical tools, systems biology reconstitutes the entire scenario of the cell, tissue and even organism from the pieces data generated in the past decades. However, the modern biology is mainly focusing on the structure and function of the biomolecule, cell, tissue or organ, which are far from the essence of the life because of missing thermodynamic information. It is doubtable that the current systems biology-based omics is no-how to fully understand the dynamic courses of the structure, function and information in life. For this reason, we promote a novel concept of aquamoleculomics, in which the biological structure and function as well as thermodynamic characteristics and bioinformation of the aquamolecule complexes are included in this theoretical model of systems biology. Water is mother of life, matter and matrix of organism. Indeed, the fundamental roles of H₂O molecules in biological processes might be dramatically underestimated. Extremely speaking, H₂O networks in the living system might be engaged in all the biological processes including building all the biological structures, the residential places of the motherhood molecules as the honeycombs of honeybees.

Because only two variables are needed to characterize a simple thermodynamic system in equilibrium, any such system is constrained on a 2D manifold. Of particular interest are the exact 1-forms on the cotangent space of that manifold, since the integral of exact 1-forms is path-independent, a crucial property satisfied by state variables such as internal energy dE and entropy dS. Our prior work [1] shows that given an appropriate language of vector calculus, a machine can re-discover the Maxwell equations and the incompressible Navier-Stokes equations from data. In this paper, We enhance this language by including differential forms and show that machines can re-discover the equation for entropy dS given data. Since entropy appears in various fields of science in different guises, a potential extension of this work is to use the machinery developed in this paper to let machines discover the expressions for entropy from data in fields other than classical thermodynamics.

The Source Region of the Yellow River (SRYR), renowned as the "Water Tower of the Yellow River", serves as an important water conservation domain in the upper reaches of the Yellow River, significantly influencing water resources within the basin. Based on the Weather Research and Forecasting (WRF) Model Hydrological modeling system (WRF-Hydro) model, the key variables of the atmosphere-land-hydrology coupling processes over the SRYR during the 2013 rainy season are analyzed. The investigation involves a comparative analysis between the coupled WRF-Hydro and the standalone WRF simulations, focusing on the hydrological response to the atmosphere. The results reveal that the WRF-Hydro model's proficiency in depicting streamflow variations over the SRYR, yielding a Nash Efficiency Coefficient (NSE) of 0.44 and of  0.61 during the calibration and validation periods, respectively. Compared to the standalone WRF simulations, the coupled WRF-Hydro model demonstrates enhanced performance in soil heat fluxes simulations, reducing the mean Root Mean Square Error (RMSE) of surface soil temperature by 0.96 K and soil moisture by 0.01 m\textsuperscript{\textsuperscript{3}}/m\textsuperscript{\textsuperscript{3}}. Furthermore, the coupled model adeptly captures the streamflow variation characteristics with the NSE of 0.33. This underscores the significant potential of the coupled WRF-Hydro model for delineating coupled atmosphere-land-hydrology processes in regions characterized by cold climates and intricate topography.

This study investigates the stochastic dynamics of hepatitis B virus (HBV) infection using a newly proposed stochastic model. Contrary to deterministic models that fail to encapsulate the inherent randomness and fluctuations in biological processes, our stochastic model provides a more realistic representation of HBV infection dynamics. It incorporates random variability, thereby acknowledging the changes in viral and cellular populations and uncertainties in parameters such as infection rates and immune responses. We examine the solution's existence, uniqueness, and positivity for the proposed model, followed by a comprehensive stability analysis. We provide the necessary and sufficient conditions for local and global stability, offering deep insight into the infection dynamics. Furthermore, we utilize numerical simulations to corroborate our theoretical results. This research provides a robust tool for understanding the complex behavior of HBV dynamics, contributing significantly to the ongoing quest for more effective HBV control and prevention strategies.

Into the frame of the French TANDEM project (Tsunamis in the Atlantic and the English ChaNnel: Definition of the Effects through numerical Modelling) Principia has been working on the development and qualification of two in-house CFD software’s: the 2D EOLE-SV (Saint Venant) model for simulation of large scale tsunami propagation from the source up to coastal scale and the 3D EOLE-NS (Navier-Stokes) model dedicated to tsunami coastal impact modelling. This paper presents a large range of test cases carried out into the frame of the project and dedicated to the validation of numerical codes in various tsunami wave conditions. The main aspects of phenomena such as wave generation, propagation and coastal impact are investigated on academic situations. A real case simulation is concerned as well, the devastating 2011 Tohoku event which is compared with in-situ data.

Efficiency improvement is the new challenge in all fields of design. In this scenario the reduction of power losses is becoming more and more a main concern also in the design of power transmissions. Appropriate models to predict power losses are therefore from the earliest stages of the design phase. The aim of the project is to carry on lubrication simulations of several variants of a cylindrical-roller-bearing to understand the lubricant distribution and the related churning power losses. Several strategies to reduce the computational effort have been used. Among them the sectorial symmetry and three innovative meshing strategies (purely analytical with and without interfaces and analytical/subtractive) that have been implemented in the OpenFOAM® environment. The results of the different approaches were compared among them and with experimental observations showing good agreement and reasonable savings in terms of computational effort.

In this paper we investigate the geometric structure and control of exponential families depending on additional parameters, called external parameters. These generalized exponential families emerges naturally when one applies the maximum entropy formalism to derive the equilibrium statistical mechanics framework. We study the associated statistical model, compute the Fisher metric and introduce a natural fibration of the parameter space over the external parameter space. The Fisher Riemannian metric allows to endow this fibration with an Ehresmann connection and to study the geometry and control of these statistical models. As an example, we show that horizontal lift of paths in the external parameter space corresponds to an isentropic evolution of the system. We apply the theory to the example of a ideal gas in a rotating rigid container. Most of the results are expressed in local coordinates; in the appendices we hint at possible global extensions of the theory.

A recently developed thermodynamic theory for the determination of the driving force of crystallization and the crystal–melt surface tension is applied to the ice–water system employing the new Thermodynamic Equation of Seawater TEOS-10. The deviations of approximative formulations of the driving force and the surface tension from the exact reference properties are quantified, showing that the proposed simplifications are applicable for low to moderate undercooling and pressure differences to the respective equilibrium state of water. The TEOS-10 based predictions of the ice crystallization rate revealed pressure-induced deceleration of ice nucleation with an increasing pressure, and acceleration of ice nucleation by pressure decrease. This result is in, at least, qualitative agreement with laboratory experiments and computer simulations. Both the temperature and pressure dependencies of the ice–water surface tension were found to be in line with the le Chatelier–Braun principle, in that the surface tension decreases upon increasing degree of metastability of water (by decreasing temperature and pressure), which favors nucleation to move the system back to a stable state. The reason for this behavior is discussed. Finally, the Kauzmann temperature of the ice–water system was found to amount TK=116K, which is far below the temperature of homogeneous freezing. The Kauzmann pressure was found to amount pK=-212MPa, suggesting favor of homogeneous freezing upon exerting a negative pressure on the liquid. In terms of thermodynamic properties entering the theory, the reason for the negative Kauzmann pressure is the higher mass density of water in comparison to ice at the melting point.

We present a precise and universal method for anticipating the thermodynamic characteristics of polyatomic molecules. Unlike conventional computational approaches that rely on numerous experimental spectroscopy or calorimetry data, our model only utilizes a small number of molecular constants. We accurately predicted the values of entropy, enthalpy, specific heat, and Gibbs free energy for boron trifluoride, which were in excellent accordance with the available experimental data. This investigation introduces a new approach to deal with the anharmonic vibrations of polyatomic molecules

The spatial resolution measurement limitation of the position-momentum uncertainty principle is shown to mathematically originate from the Bekenstein entropy bound and the associated second law of thermodynamics, as a special case in which a statistical thermodynamic distribution of energies is specialized to a fixed, definite probe energy equal to the average energy of that distribution. This is used in combination with the Wein displacement law to predict an ultraviolet cutoff for Planck blackbody radiation at about ⅛ of the Wein peak. A new UV photon counting experiment is proposed to test for this. A physical understanding of these results may be provided by a UV-complete, intelligible theory of general relativistic quantum mechanics in which the observation of a blackbody spectrum is simply a remote observation of Hawking radiation emitted from black hole fluctuations in the gravitational vacuum.

Intermediate-luminosity-optical-transients (ILOTs) are stellar outbursts with luminosity between those of classical novae and supernovae. They are divided into a number of sub-groups depending on the erupting progenitor and the properties of the eruption. Many of the ILOTs sit on the slanted Optical Transient Stripe (OTS) in the Energy-Time Diagram (ETD) that shows their total energy vs. duration of their eruption. We describe the different kinds of ILOTs that populate the OTS and other parts of the ETD. We also stand on similarities between Planetary Nebulae (PN) to ILOTs, and suggest that some PNe were formed in an ILOT event. The high energy part of the OTS is reserved to the supernova impostors -- giant eruption of very massive stars. We show results of 3D hydrodynamical simulations of supernova impostors that expose the mechanism behind these giant eruptions, and present new models for recent ILOTs. We stand on the connection between different kinds of ILOTs, and suggest that they are powered by a similar source of energy -- gravitational energy released by mass transfer.

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can provide predictions, reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

In this paper we show how to efficiently implement parallel discrete simulations on Multi-Core and GPU architectures through a real example of application: a cellular automata model of laser dynamics. We describe the techniques employed to build and optimize the implementations using OpenMP and CUDA frameworks. We have evaluated the performance on two different hardware platforms that represent different target market segments: high-end platforms for scientific computing, using an Intel Xeon Platinum 8259CL server with48cores and also an NVIDIA Tesla6V100 GPU, both running on Amazon Web Server (AWS) Cloud, and on a consumer-oriented platform, using an Intel Core i9 9900k CPU and an NVIDIA GeForce GTX 1050 TI GPU. Performance results are compared and analysed in detail. We show that excellent performance and scalability can be obtained in both platforms, and we extract some important issues that imply a performance degradation for them. We also found that current Multi-Core CPUs with large core numbers can bring a performance very near to that of GPUs, even similar in some cases.

Polymers' controlled pyrolysis is an economical and environmentally friendly solution to prepare activated carbon. However, due to the experimental difficulty in measuring the dependence between tissues and pyrolysis parameters at high temperatures, the unknown pyrolysis mechanism hinders access to the target products with desirable morphologies and performances. In this study, we investigate the pyrolysis process of polystyrene under different heating rates and temperatures employing reactive molecular dynamics (ReaxFF-MD) simulations. A clear profile of the generation of pyrolysis products determined by the temperature and heating rate is constructed. It is found that the heating rate affects the type and amount of pyrolysis intermediates and their timing, and low-rate heating helps yield more diverse pyrolysis intermediates. While the temperature affects the pyrolytic tissues of the final equilibrium products, either too low or too high a target temperature is detrimental to generating large areas of graphitized tissue. The established theoretical evolution process matches experiments well, thus contributing to preparing target activated carbons by referring to the regulatory mechanism of pyrolytic tissues.

Eukaryotic cells contain membranes with various curvatures, from the near-plane plasma membrane to the highly curved membranes of organelles, vesicles, and membrane protrusions. These curvatures are generated and sustained by curvature-inducing proteins, peptides, and lipids, and describing these mechanisms is an important scientific challenge. In addition to that, some molecules can sense membrane curvature and a thereby be trafficked to specific locations. The description of curvature-sensing is another fundamental challenge. Curved lipid membranes and their interplay with membrane-associated proteins can be investigated with molecular dynamics (MD) simulations, given the right tools. Various methods for simulating curved membranes with MD are discussed here, including tools for setting up simulation of vesicles, and methods for sustaining membrane curvature. The latter are divided into methods that exploit scaffolding virtual beads, methods that use curvature-inducing molecules, and methods applying virtual forces. The variety of simulation tools allow the researcher to closely match the conditions of experimental studies of membrane curvatures.

There is little doubt that life's origin followed from the known physical and chemical laws of Nature. The most general scientific framework incorporating the laws of Nature and applicable to most known processes to good approximation, is that of thermodynamics and its extensions to treat out-of-equilibrium phenomena. The event of the origin of life should therefore also be amenable to such an analysis. In this paper, I describe the non-equilibrium thermodynamic foundations of the origin of life for the non-expert. This ``Thermodynamic Dissipation Theory for the Origin of Life'' is founded on Classical Irreversible Thermodynamic theory developed by Lars Onsager, Ilya Prigogine, and coworkers.

To embrace the special theory of relativity with classical thermodynamics correct Lorentz transformations of thermodynamic state functions are formulated. All of these Lorentz transformations of state functions equally support the time dilation phenomenon which is experimentally verified fact. Relativistic Boltzmann constant proves entropy to be the Lorentz variant. Twin paradox thought experiments for heat engines prove heat flux to be Lorentz variant. Lorentz transformation of entropy and heat flux proves temperature Lorentz invariant. All thermodynamic state functions are proved to be Lorentz variant which is shown in accord with Lorentz transformation of thermodynamic work.

It is proposed that the primary living cells (“probionts”) cannot emerge of organic substance simply by continuous chemical complication of prebiotic macromolecules and microsystems. The complication must be accompanied by the radical thermodynamic transformation (“jump”) of prebiotic microsystems that resulted in the acquired ability to extract free energy from the environment and export entropy. This transformation is called “the thermodynamic inversion” The inversion may occur by means of the efficient (intensified) response of the microsystems on the oscillations of physic-chemical parameters in hydrothermal environment. In this case the surplus available free energy within a microsystem, when combined with the informational modality, facilitates its conversion into a new microsystem—a living probiont. It is shown the schematic representation of an oscillating prebiotic microsystem that is transforming into a living probiont. A new kind of laboratory and computational experiments on prebiotic chemistry under oscillating conditions is offered to verify the inversion concept.

Purification of biologically-derived therapeutics is a major cost contributor to the production of this rapidly growing class of pharmaceuticals. Monoclonal antibodies comprise a large percentage of these products therefore new antibody purification tools are needed. Small peptides, as opposed to traditional antibody affinity ligands such as Protein A, may have advantages in stability and production costs. Multiple heptapeptides that demonstrate Fc binding behavior that have been identified from a combinatorial peptide library using M13 Phage Display are presented herein. Seven unique peptide sequences of diverse hydrophobicity and charge were identified. All seven peptides showed strong binding to the four major human IgG isotypes, human IgM, as well as binding to canine, rat, and mouse IgG. These seven peptides were also shown to bind human IgG4 from DMEM cell culture media with 5% FCS and 5 g/L ovalbumin present. These peptides may be useful as surface ligands for antibody detection and purification purposes. Molecular docking and classical molecular dynamics (MD) simulations were conducted to elucidate the mechanisms and energetics for the binding of these peptides to the Fc region. The binding site was found to be located between the two glycan chains inside the Fc fragment. Both hydrogen bonding and hydrophobic interactions were found to be crucial for the binding interactions. Excellent agreement for the binding strength was obtained between experimental results and simulations.

Numerical models widely used for hydrocarbon phase behavior and compositional flow simulations are based on assumption of thermodynamic equilibrium. However, it is not uncommon for oil and gas-condensate reservoirs to exhibit essentially non-equilibrium phase behavior, e.g., in the processes of secondary recovery after pressure depletion below saturation pressure, or during gas injection, or for condensate evaporation at low pressures. In many cases the ability to match field data with equilibrium model depends on simulation scale. The only method to account for non-equilibrium phase behavior adopted by the majority of flow simulators is the option of limited rate of gas dissolution (condensate evaporation) in black oil models. For compositional simulations no practical yet thermodynamically consistent method has been presented so far except for some upscaling techniques in gas injection problems. Previously reported academic non-equilibrium formulations have a common drawback of doubling the number of flow equations and unknowns compared to the equilibrium formulation. In the paper a unified thermodynamically-consistent formulation for compositional flow simulations with non-equilibrium phase behavior model is presented. Same formulation and a special scale-up technique can be used for upscaling of an equilibrium or non-equilibrium model to a coarse-scale non-equilibrium model. A number of test cases for real oil and gas-condensate mixtures are given. Model implementation specifics in a flow simulator are discussed and illustrated with test simulations. A non-equilibrium constant volume depletion algorithm is presented to simulate condensate recovery at low pressures in gas-condensate reservoirs. Results of satisfactory model matching to field data are reported and discussed.

With the increased use of laterite nickel ore, the impact of high Al2O3 slag on blast furnace smelting has gradually increased. In this paper, the effects of slag basicity and Al2O3 content on slag viscosity and enthalpy change under constant temperature conditions was investigated. The changes in slag structure were analyzed by activation energy and Fourier Transform Infrared (FT-IR) spectroscopy. The relationship between slag components and slag temperature and viscosity when slag heat is reduced was investigated. The results showed that the viscosity first slightly decreased and then significantly increased with increasing basicity at constant temperature. With the addition of Al2O3 content, the viscosity of the slag increases. The activation energy increases with increasing slag basicity and Al2O3. With increasing basicity, the [SiO4]4- tetrahedral unit trough depth becomes shallow, the [AlO4]5- asymmetric stretching band migrates to lower wave numbers, and the slag structure depolymerizes. With the increase of Al2O3 content, the trough of [SiO4]4- tetrahedra deepens and the center of the symmetric stretching band moves to a higher wave number. The [AlO4]5- asymmetric stretching band becomes obvious, indicating the complexity of the slag structure. When the heat decreases, the slag temperature increases as the basicity increases, and the slag thermal stability is better at the basicity of 0.95-1.05. As the Al2O3 content increases, the thermal stability of the slag becomes worse.

The objective of this study was to evaluate the drying kinetics of the mass of jambu leaves and mass of jambu leaves in foam mat at different temperatures (50, 60 and 70 °C and thickness of 1.0 cm). The physicochemical characterization of the materials was performed before and after drying. Twelve different drying models were fitted to the drying kinetics data and the thermodynamic properties were calculated. The physicochemical parameters for the mass of leaves and foam after drying were moisture content from 2 to 7%, ash from 13 to 17%, protein content from 22 to 30%, lipids from 0.6 to 4% and total titratable acidity from 0.20 to 0.28% of tartaric acid. Wang & Singh and Midilli models were recommended for describing the drying kinetics. The use of foam mat promoted higher values of effective diffusion coefficient and activation energy and lower values of enthalpy and entropy, thus enhancing the reduction of drying time.

Although different assessments and evaluations of Dual-Miller cycle performed, specified output power and thermal performance associated with engine determined. Besides, multi objective optimization of thermal efficiency, Ecological Coefficient of performance ( ) and Ecological function ( ) by the mean of NSGA-II technique and thermodynamic analysis performed. The Pareto optimal frontier obtaining the best optimum solution is chosen by fuzzy Bellman-Zadeh, LINMAP, and TOPSIS decision-making techniques. Based on the results, performances of dual-Miller cycles and their optimization are improved.

Computational modeling (CM) is a versatile scientific methodology used to examine the properties and behavior of complex systems, such as polymeric materials for biomedical bioengineering. CM has emerged as a primary tool for predicting, setting up, and interpreting experimental results. Integrating in silico and in vitro experiments accelerates scientific advancements, yielding quicker results at a reduced cost. While CM is a mature discipline, its recent prominence in Bio-medical Engineering for biopolymer materials has recently gained prominence. In biopolymer biomedical engineering, CM focuses on three key research areas: A) Computer-aided design (CAD/CAM), B) Finite Element Analysis, and C) Molecular Dynamic Simulations. This review centers on Molecular Dynamics (MD) simulations in biopolymers, analyzing structural, functional, and evolutionary aspects of biomolecular systems over time. MD simulations solve Newton's equations of motion for all-atom systems, generating spatial trajectories for each atom, providing insights into properties like water absorption on biopolymer surfaces and interactions with solid surfaces, crucial for biomaterial assessment. This review offers readers a comprehensive overview of the diverse applications of Molecular Dynamics (MD) simulations on biopolymers. Furthermore, it emphasizes the flexibility, robustness, and synergistic relationship between in silico and experimental techniques.

X-rays production through betatron radiation emission from electron bunches is a valuable resource for several research fields. The EuAPS (EuPRAXIA Advanced Photon Sources) project, within the framework of the EuPRAXIA project, aims to provide 1-10 keV photons (soft X-rays) developing a compact plasma-based system designed to exploit self-injection processes that occur in highly nonlinear laser-plasma interaction (LWFA) to drive electron betatron oscillations. Since the emitted radiation spectrum, intensity, angular divergence, and possible coherence strongly depend on the properties of the self-injected beam, accurate preliminary simulations of the process are necessary to evaluate the optimal diagnostic device specifications. Dedicated tools for these tasks are under development: electron trajectories from particle-in-cell (PIC) simulations are currently undergoing numerical analysis through the calculation of retarded fields and spectra for various plasma and laser parameter combinations. At the same time, an analytical approach is being examined to obtain an expression for the emitted spectrum, assuming linear or quadratic electron energy variation. Code structure and some trajectories analysis results are presented.

Modeling the flow of polymer solutions requires knowledge at various length and time scales. The macroscopic behavior is described by the overall velocity, pressure and stress. The polymeric contribution to the stress requires knowledge of the evolution of polymer chains. In this work, we use a microscopic model, the finitely extensible nonlinear elastic (FENE) model, to capture the polymer behavior. The benefit of microscopic models is they remain faithful to the polymer dynamics without information loss via averaging. The downside is the computational cost needed to solve the thousands to millions of differential equations describing the microstructure. We thus describe a multiscale flow solver that utilizes GPUs for massively parallel, efficient simulations. We compare and contrast the microscopic model with its macroscopic counterpart under various flow conditions. In particular, significant differences are observed under nonlinear flow conditions, where the polymers become highly stretched and oriented.

Tsunami modeling and simulation has changed in the past few years more than it had in decades, especially so with respect to coastal inundation. Among other things, this change is supported by the approaching era of exa-scale computing, whether via GPU or more likely forms of hybrid computing whose presence is growing across the geosciences. For reasons identified across this review, exa-scale computing efforts will impact the on-shore, highly turbulent r\'egime to higher degree than the 2D shallow water equations used to model tsunami propagation in the open ocean. This short review describes the different approaches to tsunami modeling from generation to impact and underlines the limits of each model based on the flow r\'egime. Moreover, from the perspective of a future comprehensive multi-scale modeling infrastructure to simulate a full tsunami, we underline the current challenges associated with this approach and review the few efforts that are currently underway to achieve this goal. A table of existing tsunami software packages is provided along with an open Github repository to allow developers and model users to update the table with additional models as they are published and help with model discoverability.

The Meissner effect is a special property of superconductors. In this paper, based on the characteristics of the Meissner effect, a superconducting phase transition model of a superconductor under a changing magnetic field is proposed. A thermodynamic analysis was performed on the model to calculate both the changes in the thermodynamic parameters of the superconductor during the phase transition and the changes in the overall work and energy of the model. For Model 1 in the present study, the thermodynamic energy of the superconductor was not conserved due to the phase change of the superconductor under a changing magnetic field.

Friction, one of the oldest technologies harnessed by humanity, has thus far eluded a theoretical model capable of predicting kinetic friction forces across a broad range of sliding velocities. This study employs the first principle of thermodynamics to delve into sliding friction. For the first time, a general expression for the friction force has been derived, revealing that the damping of sliding friction is inversely proportional to the relative sliding velocity between the two interfaces in motion, contrary to the assumption in traditional models that it is directly proportional to the relative velocity. To simplify the general result, the concept of a friction boundary layer is introduced. An approximate expression is found to reveal the mechanism for reducing friction, which involves increasing the sliding velocity while decreasing the vibration velocity of the object perpendicular to the sliding interface.

Considering the significance of thermodynamic cycles in the global energy system, it is necessary to develop new general classes of thermodynamic cycles to relieve current energy and environmental problems. Inspired by the relationship between power cycles and refrigeration cycles, we realize that general classes of thermodynamic cycles should occur in pairs with opposite functions. Here we reverse class 1 heating cycles to obtain another new general class of thermodynamic cycles named class 2 heating cycles (HC-2s). HC-2s have two basic forms, and each contains six thermodynamic processes. HC-2s present the simplest and most general approach to utilizing the temperature difference between a medium-temperature heat source and a low-temperature heat sink to achieve efficient high-temperature heating. HC-2s fill the gaps that have existed since the origin of thermal science, and they will play significant roles in the global sustainable energy system.

Thermodynamic cycles are not only the core concepts of thermal science, but also key approaches to energy conversion and utilization. So far, power cycles and refrigeration cycles have been the only two general classes of thermodynamic cycles. While diverse types of systems have been developed to perform thermodynamic cycles, no new general classes of thermodynamic cycles have been proposed. Based on the basic principles of thermodynamics, here we propose and analyze a new general class of thermodynamic cycles named class 1 heating cycles (HC-1s). Two basic forms of HC-1s are obtained by connecting six essential thermodynamic processes in the proper order and forming a thermodynamic cycle. HC-1s present the simplest and most general approach to utilizing the temperature difference between a high-temperature heat source and a medium-temperature heat sink to achieve efficient medium-temperature heating and/or low-temperature cooling. HC-1s fill the gaps that have existed since the origin of thermal science, and they will play significant roles in energy conservation and emission reduction.

In this study, a new structure is proposed based on the body-centered cubic (BCC) lattice structure: adding a cubic truss in the center of the BCC structure, which is denoted by TLC (Truss-lattice-cube). The different dimensions of the central cube can affect the mechanical properties of the lattice structure notably. With a fixed length (15mm) of a unit cell, the optimal size for the central cube is determined to be 5 mm. Quasi-static compressive tests are performed on specimens made of Polylactic Acid (PLA) using additive manufacturing technology. The deformation characteristics of the new structure is analyzed in detail by experiments and numerical simulations. Compared to the BCC structure, the mechanical properties of the TLC structure are significantly improved. The initial flow stress of the TLC has been increased by 122% at a strain of 0.1; the specific strength has been enhanced by 293% at a strain of 0.5; and the specific energy absorption has been improved by 312% at the strain of 0.6. Printing defects in the lattice structure could deteriorate the mechanical property remarkably. In this work, incorporation of microcracks into the finite element model makes the simulation capture the influence of printing defects and improve the predictive accuracy of simulation significantly.

Turbulent dissipation is a central issue in the star and galaxy formation process. Its fundamental property of space-time intermittency, well characterised in incompressible laboratory experiments, remains elusive in cosmic turbulence. Progress requires the combination of state-of-the-art modelling, numerical simulations and observations. The power of such a combination is illustrated here where the statistical method intended to locate extrema of velocity shears in a turbulent field is applied to numerical simulations of compressible magneto-hydrodynamical (MHD) turbulence dedicated to dissipation scales and to a nearby turbulent diffuse molecular cloud. We demonstrate that short-spacing increments of observables can detect strongly dissipative structures. In our simulations, we compute structure functions of various synthetic observables and show that they verify Extended Self-Similarity. This allows to compute their intermittency exponents and we show how they could help constraining some properties of the underlying three-dimensional turbulence. In observations of a turbulent cloud close to the Sun in our Galaxy, a remarkable coherent structure of velocity shear extremum is disclosed . At the location of the largest velocity shear, this coherent structure is found to foster 10 × more carbon monoxide molecules than standard diffuse molecular gas, an enrichment supported by models of non-equilibrium warm chemistry triggered by turbulent dissipation. The power of the combination between modelling and observations is also illustrated by observations of the CH+ cation that provide unique quantitative informations on the kinetic energy trail in the massive, multi-phase and turbulent circum-galactic medium of a galaxy group at redshift z = 2.8.

Shale gas and coalbed methane are energy sources that partly or even mainly consist of methane stored in an adsorbed state in the pores of the organic-rich rock and coal seams. In this study, the graphene nanoslit model is employed to model the nanometer slit pores in shale and coal. The Grand Canonical Monte Carlo (GCMC) and Molecular Dynamics (MD) modelling methods are used to investigate the mechanisms of adsorption and displacement of methane in the slit pore. It was found that the CVFF forcefield resulted in the largest adsorption amount, while the PCFF forcefield resulted in the least. The COMPASS and COMPASS II force fields led to similar results. As the width of the slit pore increases, the adsorption amount of gas molecules increases, and the number density profile of adsorbed methane molecules alters from a single adsorption layer to multi-adsorption layers. The minimum slit pore width at which methane molecules can penetrate the slit pore was found to be 0.7 nm. Moreover, it is demonstrated that by lowering the temperature, the adsorption rate of the methane increases since the adsorption is an exothermic process. Enhancing methane recovery was investigated by the injection of gases such as CO2 and N2 to displace the adsorbed methane. The comparison of adsorption isotherms of gas molecules provided the following order in terms of the amount of adsorption: CO2>CH4>N2.

This study explores wind energy resources in different locations through the Gulf of Oman and also their future variability due climate change impacts. In this regard, EC-EARTH near-surface wind outputs obtained from CORDEX-MENA simulations are used for historical and future projection of the energy. The ERA5 wind data are employed to assess the suitability of the climate model. Moreover, the ERA5 wave data over the study area are applied to compute sea surface roughness as an important variable for converting near-surface wind speeds to those of wind speed at turbine hub height. Considering the power distribution, bathymetry and distance from the coats, some spots as tentative energy hotspots to provide a detailed assessment of directional and temporal variability and also to investigate climate change impact studies. RCP8.5 is a common climatic scenario is used to project and extract future variation of the energy in the selected sites. The results of this study demonstrate that the selected locations have a suitable potential for wind power turbine plans and constructions.

This study has developed a Matlab application for simulating statistical models project (SMp) probabilistic distributions that are similar to binomial and Poisson, which were created by mathematical procedures. The simulated distributions are graphically compared with these popular distributions. The application allows to obtain many probabilistic distributions, and shows the trend (τ ) for n trials with success probability p, i.e. the maximum probability as τ=np. While the Poisson distribution PD(x;µ) is a unique probabilistic distribution, where PD=0 in x=+∞, the application simulates many SMp(x;µ,Xmax) distributions, where µ is the Poisson parameter and value of x with generally the maximum probability, and Xmax is the upper limit of x with SMp(x;µ,Xmax) ≥ 0 and limit of the stochastic region of a random discrete variable. It is shown that by simulation via, one can get many and better probabilistic distributions than by mathematical one.

This study has developed a Matlab application for simulating statistical models project (SMp) probabilistic distributions that are similar to binomial and Poisson, which were created by mathematical procedures. The simulated distributions are graphically compared with these legendary distributions. The application allows to obtain many probabilistic distributions, and shows the trend (τ ) for n trials with success probability p, i.e. the maximum probability as τ=np. While the Poisson distribution PD(x;µ) is a unique probabilistic distribution, where PD=0 in x=+∞, the application simulates many SMp(x;µ,Xmax) distributions, where µ is the Poisson parameter and value of x with generally the maximum probability, and Xmax is upper limit of x with SMp(x;µ,Xmax) ≥ 0 and limit of the stochastic region of the random discrete variable X. It is shown that by simulation via, one can get many and better probabilistic distributions than by mathematical one.

Transient analysis of gas flow in pipeline was studied. Finite Element Method based on ANSYS FLOTRAN was used to account for changes in pressure, temperature and flow rate. Compressibility factor function of temperature and pressure was considered. For non-isothermal transient results, the pressure and the flowrate gave wave propagation as a result of slow transients created by demand condition at the outlet end of custody transfer. Results obtained were in agreement with the demand restrictions at the outlet end of custody transfer indicating that the predictions are accurate and reliable. The results demonstrated that the Finite Element Method gave accurate prediction of pressure, temperature and flowrate in transient gas studies. For steady state non-isothermal model, results showed that the magnitude of the average pressure drop was higher when pressure was predicted with constant compressibility factor, but the same average pressure drop was reduced when the pressure was predicted with variable compressibility factor, z. Since compressibility factor is a function of temperature and pressure, the above findings signifies that in the case when gas temperature does not stabilize, the prediction of pressure with isothermal model and constant compressibility factor will lead to significant errors.

It is always a difficult task to assign the peaks recorded from vibrational spectrum. Herein we explored a new pathway of density functional theory (DFT) simulation to present three kinds of spectra of ice XIV and can be referenced as inelastic neutron scattering (INS), IR, and Raman experimental spectrum. The INS spectrum is proportional to phonon density of states (PDOS) while the photon scattering signals reflect the normal vibration frequencies near the Brillouin zone (BZ) center. Based on good agreements with experimental data, we identified the relative frequency and made scientific assignments by normal vibration modes analysis. The debating two hydrogen bond (H-bond) peaks among ice phases from INS have been discussed and the dynamic process of H-bond vibrations were found to be classified into two basic modes. we deduce that two H-bond modes are a general rule among ice family and more works are ongoing to investigate this subject.

The efficiencies of conventional triple organic Rankine cycle (TORC) configurations are restricted by the low critical temperatures of common working fluids (＜320 °C). This paper proposes a high-temperature TORC system. A near-azeotropic mixture biphenyl/diphenyl oxide (BDO), which has stellar track record in the high-temperature ORC applications, is innovatively adopted as the top and middle ORC fluid simultaneously. The top cycle evaporation temperature is up to 400 °C in view of BDO’s high stability. Five conventional organic fluids are chosen for the bottom ORC. Fundamentals of the TORC system are illustrated, and mathematical models are built. Thermo-economic performance is investigated. The results hint that the optimal thermodynamic property is achieved on the use of benzene as the bottom fluid. The maximum thermal and exergy efficiencies are respectively 40.86% and 74.14%. The largest irreversible loss occurs inside the heat exchanger coupling the middle and bottom ORCs, accounting for above 30% of the total entropy generation. The levelized energy cost (LEC) is 0.0284 $/kWh. Given the same heat source condition, the TORC system can boost efficiency by 1.02% and drive down LEC by 0.0024 $/kWh compared with a BDO mixture-based cascade ORC.

Ni-Al alloys have good corrosion-resistant properties making them preferable as a material of choice compared to ordinary metals. A good understanding of the phase diagram and thermodynamic properties of the Ni-Al alloy binary system is essential in determining relevant industrial applications of the system. The Ni-Al phase diagram is being reassessed since those found in the literatures are based on very old experimental studies. The phase diagram and thermodynamic properties of the Ni-Al alloy is determined using the Thermo-Calc 2022b software over a temperature range of 200 K to 2000 K and 1 atm of pressure. A total of nineteen fields (including elemental Al and Ni) comprising of five single solid phases, one liquid phase, and twelve double regions was obtained from the equilibrium calculation. The thermodynamic calculation has six invariant reactions comprising of two eutectic reactions (914.83 ± 0.025 K and 1641.80 ± 0.025 K), three peritectic reactions (1123.55 ± 0.025 K, 1400.72 ± 0.025 K and 1642.65 ± 0.025 K), and one peritectoid reaction (913.96 ± 0.025 K). The solubility of Al in Ni is also determined. The Thermo-Calc temperature – composition of the Ni-Al alloy is compared with the existing experimental studies to identify the physical properties of the different thermodynamic properties of the Ni-Al alloy binary system.

Background and motivation: Combining Deep Reinforcement Learning (Deep RL) and Health Systems Simulations has significant potential, for both research into improving Deep RL performance and safety, and in operational practice. While individual toolkits exist for Deep RL and Health Systems Simulations, no framework to integrate the two has been established. Aim: Provide a framework for integrating Deep RL Networks with Health System Simulations, and to ensure this framework is compatible with Deep RL agents that have been developed and tested using OpenAI Gym. Methods: We developed our framework based on the OpenAI Gym framework, and demonstrate its use on a simple hospital bed capacity model. We built the Deep RL agents using PyTorch, and the Hospital Simulation using SimPy. Results: We demonstrate example models using a Double Deep Q Network or a Duelling Double Deep Q Network as the Deep RL agent. Conclusion: SimPy may be used to create Health System Simulations that are compatible with agents developed and tested on OpenAI Gym environments. GitHub repository of code: https://github.com/MichaelAllen1966/learninghospital

Blazars are highly variable and display complex characteristics. A key characteristic is the flux distribution or flux PDF whose shape depends upon the form of the underlying physical process driving variability. The BL Lacertae, Mrk 421 is one of the brightest and most variable blazars across the electromagnetic spectrum. It has been reported to show hints of lognormality across the spectrum from observed histograms of radio to gamma-ray fluxes. This would imply that the underlying mechanisms may not conform to the "standard" additive, multi-zone picture, but could potentially have multiplicative processes. This is investigated by testing the observed lightcurves at different wavelengths with time-series simulations. We find that the simulations reveal a more complex scenario, than a single lognormal distribution explaining the multiwavelength lightcurves of Mrk 421.

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. There is evidence that some ligands that bind to the GPCRs activate different downstream signaling pathways. G protein activation or -arrestin biased signaling involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. Molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators' interactions between receptors and ligands. Here, we highlight recent simulation studies and methods used to study biased G protein-coupled receptor signaling and their conformational dynamics. We also highlight applications of MD simulations to drug discovery.

Laboratory experiments and investigations into natural phenomena in this research series have unveiled the presence of a gravitational repulsion force that permeates our environment on both microscopic and macroscopic scales, contingent upon the thermal energy present. This paper presents a study exploring novel properties of gravitational forces among gas molecules, employing principles of thermodynamics. A model has been devised based on the interactions between pairs of gas molecules. While traditional models treat gravitational interaction as a singular force, our experimental approach validates it as the composite of two distinct forces: gravitational repulsion and attraction. By utilizing established experimental data on gas thermodynamics, our model demonstrates robust performance both analytically and experimentally. It verifies the coexistence of gravitational repulsion and attraction forces among gas molecules, showcases their Inverse-Cube relationship with distance, and elucidates the direct proportionality of the repulsion force to absolute temperature. This bridges a crucial gap between energy and fundamental forces. The order of magnitudes of gravitational repulsion and attraction forces are found to be considerably large, contrasting with the low values predicted by classical theory, which results from their amalgamation. Recognizing these forces as substantial in magnitude promises unprecedented outcomes and advancements.

Symmetry-breaking is a phenomenon that is observed in various contexts, from the early universe to complex organisms, and is considered a key puzzle in understanding the emergence of life. The presence of enantiomeric amino acids and proteins highlights its critical role. However, the origin of symmetry-breaking has yet to be comprehensively explained, particularly from an energetic standpoint. This article explores a novel approach by considering energy dissipation, specifically lost free energy, as a crucial factor in elucidating symmetry-breaking. Through a thorough thermodynamic analysis applicable to enantiomeric systems at various scales, from elementary particles to aggregated matter like crystals, we present experimental findings that establish a direct correlation between non-equilibrium free energy and energy dissipation during formation processes. Results emphasize the pivotal role of energy dissipation not only as an outcome but as the trigger for symmetry-breaking. This insight suggests that understanding the origins of complex systems, from cells to living beings and the universe itself, requires a lens focused on non-equilibrium processes, ultimately facilitating the formation of non-racemic systems.

The entropy of a system gives a powerful insight into its microscopic degrees of freedom, however standard experimental ways of measuring entropy through heat capacity are hard to apply to nanoscale systems, as they require the measurement of increasingly small amounts of heat. Two alternative entropy measurement methods have been recently proposed for nanodevices: through charge balance measurements and transport properties. We describe a self-consistent thermodynamic framework for treating few-electron nanodevices which incorporates both existing entropy measurement methods, whilst highlighting several ongoing misconceptions. We show that both methods can be described as special cases of a more general relation and prove its applicability in systems with complex microscopic dynamics – those with many excited states of various degeneracies.

In computational neuroscience, spiking neurons are often analyzed as computing devices that register bits of information, with each action potential carrying at most one bit of Shannon entropy. Here, I question this interpretation by using Landauer's principle to estimate an upper limit for the quantity of thermodynamic information that can be dissipated by a single action potential in a typical mammalian neuron. I show that an action potential in a typical mammalian cortical pyramidal cell can carry up to approximately 3.4e11 natural units of thermodynamic information, or about 4.9e11 bits of Shannon entropy. This result suggests that an action potential can process much more information than a single bit of Shannon entropy.

The Boltzmann-Gibbs (BG) entropy has been used in a wide variety of problems for more than a century. It is well known that BG entropy is extensive, but for certain systems such as those dictated by long-range interactions, the entropy must be non-extensive. Tsallis entropy possesses non-extensive characteristics, which is parametrized by a variable q (q = 1 being the classic BG limit), but unless q is determined from microscopic dynamics, the model remains but a phenomenological tool. To this date very few examples have emerged in which q can be computed from first principles. This paper shows that the space plasma environment, which is governed by long-range collective electromagnetic interaction, represents a perfect example for which the q parameter can be computed from micro-physics. By taking the electron velocity distribution function measured in the heliospheric environment into account, and considering them to be in quasi equilibrium state with electrostatic turbulence known as the quasi-thermal noise, it is shown that the value corresponding to q = 9/13 = 0.6923 may be deduced. This prediction is verified against observation made by spacecraft, and it is shown to be in excellent agreement.

Computational simulations of three-dimensional flow around a NACA 0018 wing with an aspect ratio of AR=5 are carried out using the Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations with the Shear-Stress Transport turbulence model closure. Simulations were performed to capture the aerodynamic stall hysteresis using a developed Pseudo-Transient Continuation (PTC) method based on a dual time step approach in the CFD code OpenFOAM. The flow was characterized by an incompressible Mach number of M=0.12 and a moderate Reynolds number Re=0.67×10^6. The results obtained indicate the presence of a noticeable aerodynamic hysteresis in the static dependencies of the force and moment coefficients, as well as the manifestation of bi-stable flow separation patterns, accompanied by the development of asymmetry in the stall zone. The URANS simulation results are in good agreement with the experimental data obtained for the NACA 0018 finite aspect ratio wing from the low-speed wind tunnel under the same test conditions. A new phenomenological bifurcation model of aerodynamic stall hysteresis under static and dynamic conditions is formulated, which is proven to be able to closely match the experimental data.

Fires are complex multi-physics problems that span wide spatial scale ranges. Capturing this complexity in computationally affordable numerical simulations for process studies and “outer-loop” techniques (e.g., optimization and uncertainty quantification) is a fundamental challenge in reacting flow research. Further complications arise for propagating fires where a priori knowledge of the fire spread rate and direction is typically not available. In such cases, static mesh refinement at all possible fire locations is a computationally inefficient approach to bridging the wide range of spatial scales relevant to fire behavior. In the present study, we address this challenge by incorporating adaptive mesh refinement (AMR) in fireFoam, an OpenFOAM solver for simulations of complex fire phenomena involving pyrolyzing solid surfaces. The AMR functionality in the extended solver, called fireDyMFoam, is load balanced, models gas, solid, and liquid phases, and allows us to dynamically track regions of interest, thus avoiding inefficient over-resolution of areas far from a propagating flame. We demonstrate the AMR capability and computational efficiency for fire spread on vertical panels, showing that the AMR solver reproduces results obtained using much larger statically refined meshes, but at a substantially reduced computational cost. We then leverage the computational efficiency of the AMR solver to demonstrate an optimization framework for fire suppression based on the open-source Dakota toolkit.

Water pollution is regarded as a dangerous problem that needs to be resolved right away. This is largely due to the positive correlation between the increase in global population and waste production, especially food waste. Hydrogel beads based on sodium alginate (Alg) and pomegranate fruit peels (PP) were developed for the adsorption of Safranin O dye (SO) in aqueous solutions. The obtained Alg-PP beads were widely characterized. The effects of the contact time (0-180 min) and initial concentration (10-300 mg/l), the initial pH (2-10), the adsorbent dosage (1-40 g/l) and the temperature (293-333 K) were investigated through batch tests. The data proved that the adsorption kinetics of SO reached equilibrium within 30 min and up to 180 min. The adsorption data of SO onto synthesized beads follow the pseudo second-order model. The experimental data fitted very well to both Freundlich and Langmuir models. The maximum adsorption capacity of Alg-PP could reach 8.99 mg/g at 293 K. Calculation of Gibbs free energy and enthalpy indicated that adsorption of SO onto Alg-PP is spontaneous and endothermic. Finally, analysis of diffusion and mass transport phenomena were presented. The reported results revealed that the Alg-PP beads could be used as a novel natural adsorbent for the removal of cationic dye pollutants from aqueous environments.

The interactions of metal ions with fulvic acids resulted in changes in the particle size and charge. Particle size distribution was trimodal for fulvic acids, bimodal for fulvic complexes with calcium and magnesium, and copper-fulvic complexes had only one size fraction. The compensation of negative charge of carboxylic and phenolic functional groups by positive charged metal ions resulted in the increase in zeta potential which was getting closer to zero in the case of copper-fulvic complexes. However, all metal humic complexes behaved as colloidally unstable which resulted in visually observable sedimentation. Calorimetric measurements provided positive values of changes in enthalpy which indicated endothermic processes. In contrast quantum chemical calculations as well as experiments with model compounds provided negative values indicated exothermic processes. Changes in Gibbs energy were determined as negative and changes in entropy as positive which meant spontaneous processes.

The removal of metal ions by biosorption on inexpensive materials is still a challenge for environmental engineering research. In this study, marine green algae biomass (Ulva lactuca sp.) and the biochars obtained from this biomass, at 320 C (BC-320) and 550 C (BC-550), were used as biosorbents for the removal of Cu(II) ions from aqueous solution. In addition to comparing the biosorption capacities, the determination of the thermodynamic parameters allows the choice of the most suitable material for the biosorption processes. The experimental results, obtained for Cu(II) ions biosorption on each biosorbent (algae biomass (AB), BC-320 and BC-550), at three different temperatures (10, 30 and 50 C) were analyzed using Langmuir and Freundlich isotherm models, while pseudo-first order, pseudo-second order and intra-particle diffusions models were used to model the kinetic data. The biosorption of Cu(II) ions is best described by the Langmuir model and the pseudo-second kinetic model, regardless of the type of biosorbent. The maximum biosorption capacity (qmax, mg/g) depends on the temperature, but also on the type of biosorbent, and follow the order: BC-320 < AB < BC-550. Using the experimental isotherms, the thermodynamic parameters (G0, H0 and S0) for the biosorption of Cu(II) ions on each biosorbent were calculated. The analysis of the obtained values constitutes the main arguments in choosing BC-550 as the most effective biosorbent for the removal of Cu(II) ions from aqueous media.

Li3FeN2 material was synthesized by two-step solid-state method from Li3N (adiabatic camera) and FeN2 (tube furnace) powders. Phase investigation of Li3N, FeN2 and Li3FeN2 were carried out. Discharge capacity of Li3FeN2 is 343 mAh g-1, that is about 44.7% of theoretic capacity. The molar heat capacity of Li3FeN2 at constant pressure in the temperature range 298-900 K should be calculated as Cp,m = 77,831 + 0,130 × T – 6,289 × T-2, where T is absolute temperature, . Thermodynamic characteristics of Li3FeN2 were determined as next: entropy S0298 = 116.2 J mol-1 K-1, molar enthalpy of dissolution ΔdHLFN = ˗ 206,537 ± 2,8 kJ mol−1, the standard enthalpy of formation ΔfH0 = ˗ 291.331 ± 5.7 kJ mol−1, entropy S0298 = 113.2 J mol-1 K-1 (Neumann-Kopp rule) and 116.2 J mol-1 K-1 (W.Herz rule), the standard Gibbs free energy of formation ∆f G0298 = ˗276,7 kJ mol-1.

A triangular lattice model for pattern formation by core-shell particles at fluid interfaces is introduced and studied for the particle to core diameter ratio equal to 3. Repulsion for overlapping shells and attraction at larger distances due to capillary forces are assumed. Ground states and thermodynamic properties are determined analytically and by Monte Carlo simulations for soft outer- and stiffer inner shells, with different decay rates of the interparticle repulsion. We find that thermodynamic properties are qualitatively the same for slow and for fast decay of the repulsive potential, but the ordered phases are stable for temperature range depending strongly on the shape of the repulsive potential. More importantly, there are two types of patterns formed for fixed chemical potential - one for a slow and another one for a fast decay of the repulsion at small distances. In the first case two different patterns - for example clusters or stripes - occur with the same probability for some range of the chemical potential. For fixed concentration, an interface is formed between two ordered phases with the closest concentration, and the surface tension takes the same value for all stable interfaces. In the case of degeneracy, a stable interface cannot be formed for one out of four combinations of the coexisting phases, because of a larger surface tension. Our results show that by tuning the architecture of a thick polymeric shell, many different patterns can be obtained for sufficiently low temperature.

As known all physical properties of solids are described well by the system of quantum linear harmonic oscillators. It is shown in the present paper that the system consisting of classical linear harmonic oscillators having temperature dependent masses or (and) frequencies has the same partition function as the system consisting of quantum linear harmonic oscillators having temperature independent masses and frequencies while the means of the square displacements of the positions of the oscillators from their mean positions for the system consisting of classical linear harmonic oscillators having: the temperature dependent masses; temperature dependent frequencies; and temperature dependent masses and frequencies differ from each other and from that of the system consisting of quantum linear harmonic oscillators, and hence, the system consisting of classical linear harmonic oscillators describes well the thermodynamic properties of the system consisting of quantum linear harmonic oscillators and solids.

Drilling operation is the major cost factor for the industry. Appropriately designed operations are essential for successful drilling. Optimized drilling operations also allow enhanced drilling performance and reduce the drilling cost. This is achieved by increasing the bit life (minimizing premature bit wear), drilling faster that allows reducing drilling time, and also reducing tripping operations. The paper presents two parts. The first part compares the parametric physic-based swab and surge simulation results obtained from Bingham Plastic, Power Law, and Robertson and Stiff models. The aim is to show how the model predictions deviate from each other. Two 80-20 Oil Water ratios (OWR) and two 90-10 OWR oil-based drilling fluids with 1.96 sg and 2.0 sg were considered in vertical- and deviated wells. Simulation results analysis revealed that the deviations depend on the drilling fluid’s physical and rheological parameters as well as the well trajectory. Moreover, the model’s predictions were inconsistent. Data-driven machine learning (ML) modeling makes up the second section. Data-driven modeling was done using both software-generated datasets and field datasets. Results show that the Random Forest Regressor (RF), Artificial Neural Network (ANN), Long Short-Term Memory (LSTM), LightGBM, XGboost, and Multivariate regression models predict the training and test datasets with higher R-squared and minimum root means square values. Deploying the ML model in real-time applications and the planning phase would have the potential application of artificial intelligence for well planning and optimization processes.

Large office buildings are responsible for a substantial portion of energy consumption in urban districts. However, thorough assessments regarding the Nordic countries are still lacking. In this paper we analyse the largest dataset to date for a Nordic office building, by considering a case study located in Stockholm, Sweden, that is occupied by nearly a thousand employees.Distinguishing the lighting and occupants’ appliances energy use from heating and cooling, we can estimate the impact of occupancy without any schedule data. A standard frequentist analysis is compared with Bayesian inference, and the according regression formulas are listed in tables that are easy to implement into building performance simulations (BPS). Monthly as well as seasonal correlations are addressed, showing the critical importance of occupancy.A simple method, grounded on the power drain measurements aimed at generating boundary conditions for the BPS, is also introduced; it shows how, for this type of data and number of occupants, no more complexities are needed in order to obtain reliable predictions. For an average year, we overestimate the measured cumulative consumption by only 4.7%. The model can be easily generalised to a variety of datasets.

NMR Spectroscopy is used, in the temperature range 180-350K, to study local order and transport in liquid water (pure and confined) and its solutions with glicerol and methanol at different molar fractions. Being the liquid water thermodynamic dominated by polymorphism (two coexisting liquid phases: high- and low-density HDL and LDL) - with the LDL due to the hydrophilic HB interactions, originating in the supercooled regime the tetrahedral networking and the liquid-liquid transition – we focused our interest to hydrophobic effects (HE) on these. Nowadays, if compared to the hydrophilicity, little is known about hydrophobicity so that the main purpose of this study is to gain new information on it. We measured the relaxation times (T1 and T2) and the self-diffusion (DS). From the times we took advantage of the NMR property to follow the behaviors of each molecular component (the hydrophilic and hydrophobic groups) separately; they are studied directly and DS in terms of the Adam-Gibbs model: obtaining the configurational entropy (Sconf ) and the specific heat contributions (CP,conf ). Due to the HE all the studied quantities, behave differently. For water-glycerol the HB interaction is dominant for all the conditions, whereas for water-methanol are observable two different T-regions above and below 265 K, dominated respectively by the hydrophilicity and hydrophobicity. A situation linked to the water polymorphism. Below this temperature, where the LDL phase and the HB networking develops and grows, the times and CP,conf change behaviors leading to maxima and minima. Above it, where the HB becomes weak and less stable and the HDL dominates, the hydrophobicity determines the solution properties.

In the case of vascular plants the process of water loss by leafs and water absorption by the root is well known. There is agreement on the passive nature of long-distance moisture movement in the dead cells of the xylem; however, controversy exists focusing on the long-distance water transport principle. Hales (1726) founded a view of bulk flow based on water suction after experiments with cut twigs. The previous doctrine of long-distance water transport within vessel elements and tracheid of the xylem of intact plants – the relevant cohesion theory in text books – was developed mainly by Boehm (1893), Renner (1911) and Dixon (1914) with plant artefacts. Water movement according to this theory is based on an assumed hydrodynamic bulk fluid flow in xylem in continuous water columns (free of water vapour space), under tension, according to the law of Poiseuille (see e.g. Dixon 1914). Physically hydrodynamics is part of fluid mechanics, as a result Poiseuille’s law is usually valid only for hydrodynamic bulk flow in ideal capillaries (Sutera & Skalak, 1993). Besides the basic requirement for transport, according to cohesion theory, the existence of ideal capillaries is not compatible with either: “Because vessel elements and tracheid do not stand as ideal capillaries. …” (Bresinsky et al. 2008, translated from German). Unlike ideal capillaries, the walls of vessel elements and tracheid interact with the transported water. These walls are able to function as a source or as a sink for the transported water because of interaction with the cell walls. With the interaction, vessel elements and tracheid, part of the xylem, can shrink and swell, unlike ideal capillaries. Because the xylem (in woody plants part of the wood) is inconsistent with the basic law of fluid flow, the equation of mass balance (Zimmermann et al. 2004) and cohesion theory are not strictly followed.Many plant physiologists view the cohesion theory as appropriate, however, this theory remains controver­sial, i.e. by Eisenhut (1988), Laschimke (1990) and Hahn (1993). Nultsch (1996) gives doubts referring to the present doctrine of plant water transport. Zimmermann et al. (2004) reject the cohesion theory and conclude: “... that the arguments of the proponents of the Cohesion Theory are completely misleading” (Zimmermann et al. 2004). Hence cohesion theory can be treated as inapplicable and the question arises: how does water transport in fact function? In the following, it is gone into in more detail. A sorption hypothesis of actual water transport, based on empirical fact, shall be addressed in this paper.

A new technique, the sessile drop accelerometry (SDACC) for the study and measurement of the interfacial energies of solid-liquid-gas systems is tested and discussed in this study. The instrument –a combination of a goniometer with high speed camera and a laboratory drop tower- and the evaluation method, were designed to evaluate the interfacial energies as a function of the geometrical changes of a sessile drop shape due to the effect of “switching off” gravity during the experiment. The method bases on the thermodynamic equilibrium of the system interfaces and not on the balance of bi-dimensional tensors on the solid-liquid-gas contour line. A comparison of the model with the widely accepted Young`s equation is discussed in this study.

A pivotal strategy in decarbonizing the transport sector is by transitioning from internal combustion engine vehicles (ICEVs) to electric vehicles (EVs), especially if they are charged with renewable energy. To provide an incentive to the common populace to switch to EVs, cities should set up - or enable or even encourage the private sector to do so - electric vehicle charging stations (EVCS) at multiple and suitably accessible locations. This paper aims to present a digitaltwin-enabled methodology study for selecting the most suitable possible locations to establish EVCS with consideration of the low voltage (LV) network operation in terms of power losses, active and reactive power in EV charging demand and the effect of the EVCS on the transformers, based on examples in Newcastle upon Tyne (England, UK). It also evaluates the most suitable type of charger that could be installed at these locations by studying the estimated traffic, cost and impact of the EVCS on the power grid, both at peak and off hours. This paper demonstrates how digitaltwin technologies, such as the web-based tool AutoDesign by Northern Powergrid and EA Technology and open source Pandapower can be used for identifying the most viable and cost-effective locations to connect EVCSs to a Distribution Network Operators (DNO)’s LV network. The results of the analysis indicate that with the installation of the most suitable type of EVCS in a particular location, the local power grid would not be significantly affected at any hour of the day.

3D bioprinting, like 3D printing, is a process that builds structures by depositing material layer by layer. In the specific case of bioprinting, these layers are composed of a biocompatible material together with living cells, allowing the creation of three-dimensional structures that recreate functional tissues and organs. Bioprinting cell-laden structures is complicated by the high rate of cell damage and stress that occurs during the process, caused by the high pressures and stresses to which they are subjected. To minimise damage during the process, it is important to study and optimise certain bioprinting conditions beforehand and to analyse how they affect the cells by means of computational simulations. In this review, both quantitative and qualitative data are collected to improve the nozzle geometry through computational simulation studies. Optimal ranges for nozzle diameter (0.2-1mm) and length (8-10mm; 300-900μm) have been defined, and recommendations for improving nozzle performance during bioprinting, such as the provision of an inner angle of 20-30°, an inner EDTA coating and a shear stress of less than 10kPa, have been gathered. Finally, based on the data collected, a design of experiments is proposed to obtain an optimal bioprinting configuration for a particular bioink.

This study investigates the magnetization of Cobalt (Co) nanoparticles with two different crystal structures: a centered cubic (fcc) and an hexagonal close packed (hcp), under the magnetic field of the Earth. For that reason, the hysteresis loops of the two material samples are exported and compared with the corresponding value by experimental procedure in the laboratory under the Earth’s magnetization. For both samples, the nanoparticle number variation and the nearest neighbor distance effect, in the Co hysteresis loop are investigated, structural information that is useful for the sample preparation. For the simulations the Finite Element Method (FEM) was used, while experiments are conducted using a measurement set-up based on fluxgate magnetometers. From the simulation, it is proven that, for both samples, the nanoparticle number effect does not change the hysteresis loop. In contrast, when the distance between the nanoparticles increases, the magnetization obtains smaller value. This study aims to build a novel, simple and economical measurement arrangement enabling in situ magnetic measurements for ferromagnetic nano-material characterization. As the need for continuous knowledge of the condition and the phase of ferromagnetic nanoparticles that is crucial in many applications, like catalytic processes, this work aims to introduce an alternative approach that will add another dimension to conventional chemical reaction engineering this of the magnetic state of the material resolved in time and space and to contribute in the field of the monitoring and the early fault diagnosis of nanoparticle materials.

Air pollution is responsible for any adverse effects on human beings. Thermal discomfort, on the other hand, is able to overload the human body and eventually provoke health implications due to the heat imbalance. Methods: The aim of the present work is to study the behavior of two bio-climatic indexes and statistical characteristics of the air quality index for Sofia city - the capital of Bulgaria for the period 2008 - 2014. The study is based on WRF-CMAQ model system simula-tions with a spatial resolution of 1 kilometer. The air quality is estimated by the air quality index, taking into account the influence of different pollutants and the thermal conditions by two indi-ces, respectively, for hot and cold weather. Results: It was found that half of the heat and cold index categories are present in the simulations. Their distribution has some spatial features. All air quality categories are present in the domain, with dominance only of the O3 and PM10. Conclu-sions: It was found that Sofia is not so hot and air polluted place, but in some situations, people have to have some concerns when intend to be outdoors for a prolonged time.

The availability of high quality surface observations of precipitation and volume observations by polarimetric operational radars make it possible to constrain, evaluate and validate numerical models with a wide variety of microphysical schemes. In this article, a novel particle-based Monte-Carlo microphysical model (called “McSnow”) is used to simulate the outer rain bands of Hurricane Dorian which traversed the densely instrumented precipitation research facility operated by NASA at Wallops Island, Virginia. The rain bands showed steady stratiform vertical profiles with radar signature of dendritic growth layers near −15 °C and peak reflectivity in the bright band of 55 dBZ along with polarimetric signatures of wet snow with sizes inferred to exceed 15 mm. A 2D-video disdrometer measured frequent occurrences of large drops >5 mm and combined with an optical array probe the drop size distribution was well-documented in spite of uncertainty for drops <0.5 mm due to high wind gusts and turbulence. The 1D McSnow control run and four numerical “experiments” were conducted and compared with observations. One of the main findings is that even at the moderate rain rate of 10 mm/h collisional breakup is essential for the shape of the drop size distribution.

In this manuscript, an automated framework dedicated to design space exploration and design optimization studies is presented. The framework integrates a set of numerical simulation, computer-aided design, numerical optimization, and data analytics tools using scripting capabilities. The tools used are open-source and freeware, and can be deployed on any platform. The main feature of the proposed methodology is the use of a cloud-based parametrical computer-aided design application, which allows the user to change any parametric variable defined in the solid model. We demonstrate the capabilities and flexibility of the framework using computational fluid dynamics applications; however, the same workflow can be used with any numerical simulation tool (e.g., a structural solver or a spread-sheet) that is able to interact via a command line interface or using scripting languages. We conduct design space exploration and design optimization studies using quantitative and qualitative metrics, and to reduce the high computing times and computational resources intrinsic to these kinds of studies, concurrent simulations and surrogate-based optimization are used.

The present studies deliver the computational investigations of a 10 MW turbine with a diameter of 205.8 m developed within the framework of the AVATAR (Advanced Aerodynamic Tools for Large Rotors) project. The simulations were carried out using two methods with different fidelity levels, namely the computational fluid dynamics (CFD) and blade element and momentum (BEM) approaches. For this purpose, a new BEM code namely B-GO was developed employing several correction terms and three different polar and spatial interpolation options. Several flow conditions were considered in the simulations, ranging from the design condition to the off-design condition where massive flow separation takes place, challenging the validity of the BEM approach. An excellent agreement is obtained between the BEM computations and the 3D CFD results for all blade regions, even when massive flow separation occurs on the blade inboard area. The results demonstrate that the selection of the polar data can influence the accuracy of the BEM results significantly, where the 3D polar datasets extracted from the CFD simulations are considered the best. The BEM prediction depends on the interpolation order and the blade segment discretization.

The total enthalpies of the 16 different spin configurations that can be realized in the unit cell of the archetype spin crossover complex [Fe(phen)₂(NCS)₂] (phen=1,2-phenanthroline) were calculated applying periodic density functional theory combined with the Hubbard model and the Grimme-D2 dispersion correction (DFT+U/D2). The obtained enthalpy differences between the individual spin configurations were used to determine spin couplings of an Ising-like model and subsequent Monte Carlo simulations for this model allowed to estimate the phenomenological interaction parameter Г of the Slichter-Drickamer model, which is commonly used to describe the cooperativity of the spin transition. The calculational procedure described here, which led to an estimate of about 3 kJ/mol⁻¹ for Г, in good agreement with experiment, may be used to predict from first principles how modifications of spin crossover complexes can change the character of the spin transition from gradual to abrupt and vice versa.

We consider the quantum thermodynamic origin of the gravitational force of matter by applying the spacetime entanglement entropy and the Unruh effect originating from vacuum quantum fluctuations. By analyzing both the local thermal equilibrium and quasi-static processes of a system, we may get both the magnitude and direction of Newton's gravitational force in our theoretical model. Our work shows the possibility that the elusive Unruh effect has already shown its manifestation through gravitational force.

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this mini-review, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions and, on other hand, highly simplified structure-based Elastic Network Models of protein flexibility. In contrast to classical Molecular Dynamics the modeling strategies discussed here allow quite accurate modeling of much larger systems and longer time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.

Microplastics and radionuclides pose significant challenges to the sustainable management of water systems. The interaction of uranium-232 and americium-241 with polyurethane (PU) and polylactic acid (PLA) microplastics has been investigated in aqueous laboratory and environ-mental solutions (e.g., seawater and wastewater) as a function of temperature in various pH (4, 7, 9). The temperature increase affects positively the binding of uranium-232 and americium-241. The highest adsorption efficiency for uranium and americium is observed at the neutral and al-kaline pH region, respectively. In environmental water samples (pH ~8) the adsorption efficiency decreases significantly due to the competitive adsorption of other metals present in natural wa-ters (e.g., Ca2+) as well as the stabilization of the actinides (particularly uranium) in solution (e.g., UO2(CO3)34-). The solution composition which governs both the actinide speciation, and the type of surface-active sites is strongly associated with the surface adsorption thermodynamics and de-termines the values of the associated parameters (ΔΗo and ΔSo). Generally, the values of ΔΗo and ΔSo are positive indicating an entropy-driven reaction. However, in the case of the U(VI) adsorp-tion by PLA in seawater samples both ΔΗo and ΔSo values become negative suggesting an enthal-py-driven binding mechanism associated with a decline in randomness at the surface upon ad-sorption.

The use of biodiesel fuels in compression ignition engines leads to decrease CO, PM, HC and smoke opacity. However, NOx emissions increase importantly. Various methods to reduce NOx are used namely the EGR, the water injection method and the steam injection method. In this study, the steam injection method has been used instead of the other methods because of its benefits. This study examines the effects of steam injection on combustion products of diesel and biodiesel fuels by considering chemical equilibrium in order to determine the equilibrium combustion products involving 10 combustion products. A developed simulation code determing the equilibrium mole fractions and thermodynamic properties of combustion products is used for diesel and biodiesel fuels. It can be used for any fuel consisting of carbon, hydrogen or any oxygenated fuel. The results show that the mole fraction of CO2 and CO decrease with the steam injection ratios. NO mole fractions decrease with steam injection for lean mixtures but they increase slightly in the case of rich mixtures. The specific heat of combustion products increase with the steam injection ratios. Thus, the engine performance can be improved using the method. The model has been validated by comparing model results with the ones of NASA CEA software and GASEQ software using the methane as fuel. The relative errors of equilibrium mole fractions and thermodynamic properties of combustion products are less than 0.98 %.

Point defects induced by doping rare earth elements RE (Nd, Er) into the magnesium oxide host were investigated by classical atomistic simulations utilising the General Utility Lattice Program (GULP). Formation and association energies were calculated for the potential defect structures. Both isolated defects and defect complexes were considered. The most energetically favourable structures of defect complexes were found for rare earth doped and Li co-doped systems. The correlation between the association energy and the structure of the defect complex was investigated. The influences of Li were revealed with respect to energy and structure. The simulation results contribute to the understanding of the point defects of doped MgO and how Li influences the doping of rare earth elements in the MgO host.