We are witnessing a progressive divestment of some institutions with strong traditions and skills in physical modelling and their consequent impoverishment, to the detriment of numerical modelling. For many reasons, the economic imperatives and the exponential growth of computational means and numerical methods should certainly not be excluded. In this work, we aimed to highlight the new requirements of the recent sophisticated developments in physical modelling, precisely due to the new needs imposed on them by mathematical and numerical modelling and the growing risks in civil construction works. In this context, reflections are reported, justified by scientific and real-world examples, on the need for maintenance and reinforcement of investments in physical modelling, both to support the scientific community and to design buildings of significant economic, social and environmental impact.

Since the first case of 2019 novel coronavirus disease (COVID-19) detected on Jan 30, 2020, in India, the number of cases rapidly increased to 3819 cases including 106 deaths as of 5 April 2020. Taking this into account, in the present work, we are studying a Bats-Hosts-Reservoir-People transmission fractional-order COVID-19 model for simulating the potential transmission with the thought of individual social response and control measures by the government. The real data available about infectious cases from $14^{th}$ March to $26^{th}$ March 2020 is analysed and accordingly various parameters of the model are estimated or fitted. The Picard successive approximation technique and Banach's fixed point theory have been used for verification of the existence and stability criteria of the model. Numerical computations are done utilizing the iterative Laplace transform method. In the end, we illustrate the obtained results graphically. The purpose of this study is to estimate the effectiveness of preventive measures, predicting future outbreaks and potential control strategies using the mathematical model.

Automotive aerodynamics comprises of the study of aerodynamics of road vehicles. Its main goals are reducing drag, minimizing noise emission, improving fuel economy, preventing undesired lift forces and minimizing other causes of aerodynamic instability at high speeds. The Ahmed body has the form of a highly simplified car, consisting of a blunt nose with rounded edges fixed onto a box-like middle section and a rear end that has an upper slanted surface, the angle of which can be varied. It retains vital features of real vehicles in order to study the flow fields around it and the related turbulence models which characterizes the actual flow at elevated Reynolds number. In the present study, the aerodynamic behavior of this body is investigated numerically by the aid of commercial CFD tool: Ansys Fluent. The results of the simulation are validated with available experimental data and results of the simulations from other literatures. The numerical data were obtained for a fixed free stream velocity of 25 m/s at the inlet. The simulations were performed at a fixed slant angle of 25 degree and zero yaw angle. The present study focuses on how local refinement of mesh inside the concerned body and the outside, helps affect the results and for which grid dependency test is the primary objective of this paper. The present study also helps demonstrate how the drag of the body behaves, which is mainly the effect of pressure drag force generated at the rear portion of the body. The study also focuses on important properties like the velocity magnitude at different locations for different meshing cases, and to capture the flow pattern in the front or near the wake region. The study can be further helpful to future researchers in determining resistance, fuel efficiency etc. helping designers to optimize in specialized areas for better efficiency.

While the knowledge about the development of understanding positive integers is rapidly growing, the development of understanding zero is not well-known. Here we tested several components of preschoolers’ understanding zero: whether they can use empty sets in numerical tasks, whether they can use empty sets soon after they understand the cardinality principle, whether they know what the word “zero” refers to, and whether they categorize zero as a number. The results show that preschoolers can handle empty sets in numerical tasks as soon as they understand the cardinality principle or even earlier, and some of them know that these sets are labeled as “zero.” However, they are unsure whether zero is a number. These results identify three components of knowledge about zero: operational knowledge, linguistic knowledge, and meta knowledge. To account for these results we propose that preschoolers might understand numbers as the properties of items or objects in a set. In this view, zero cannot be a number, because an empty set does not include any items, and the missing items cannot have any property, excluding also the number property. This model might explain why zero is handled correctly in numerical tasks, while it is not regarded to be a number.

Small and micro hydropower represents an attractive solution for electricity generation, with low cost and low environmental impact. The pump-as-turbine (PAT) approach has promise in this application owing to its low purchase and maintenance costs. In this paper, a new method to predict the inverse characteristic of industrial centrifugal pumps is presented. This method is based on results of simulations performed with commercial three-dimensional CFD software. Model results have been first validated in pumping mode using data supplied by pump manufacturers. Then, results have been compared to experimental data for a pump running in reverse condition. Experimentation has been performed on a dedicated test bench installed in the Department of Civil Construction and Environmental Engineering of the University of Naples Federico II. Three different pumps, with different specific speeds, have been analyzed. Using the model results, the inverse characteristic and the best efficiency point have been evaluated. Finally, results of this methodology have been compared to prediction methods available in the literature.

A solution for the Prandtl-Blasius equation is essential to all kinds of boundary layer problems. This paper revisits this classic problem and presents a general Maple code as its numerical solution. The solutions were obtained from the Maple code, using the Runge-Kutta method. The study also considers convergence radius expanding and an approximate analytic solution is proposed by curve fitting. Similarly, the study resolves some boundary layer related problems and provide relevant Maple codes for these.

Rapidly spreading disease, COVID-19 is classified as the human-to-human transmission-able disease and currently is a pandemic in the Globe. In this paper, we propose conceptual mathematical models for COVID-19 outbreak and it's control measurement; quarantine, hospitalization and the effect of panic and anxiety. In this situation, mathematical models are a important tool to employ an effective strategy in order to fight against this pandemic. We establish the positivity and boundedness of solutions, local and global stability analysis of equilibria to examine its epidemiological relevance. To validate the model and estimating the important model parameters and prediction about the disease, we consider the real cases of Italy from $15^{th}$ Feb to $13^{th} $ April 2020. In a series of graphical map, we have presented the comparative study to estimate the current scenarios and to predict the control measurement time boundary of the outbreak.

In this paper we report a new development on the numerical model for aluminum-steam combustion. This model is based on diffusion flame of continuum regime and the thermal equilibrium between the particle and the flow field, which can be used to calculate the aluminum particle combustion model for two phase calculation conditions. The model prediction is in good agreement with the experimental data. A new type of vortex combustor was proposed for the combustion of aluminum and steam, and the mathematical model of the two phase reacting flow with in this combustor was established. The turbulence effects are modeled using the Reynolds Stress Model (RSM) with Linear Pressure-Strain approach, and the Eddy-Dissipation model is used to simulate the gas phase combustion. Aluminum particles are injected into the vortex combustor and form a swirling flow around the chamber and their trajectories are traced using the Discrete Phase Model (DPM). The simulation results show that the vortex combustor can achieve high efficient combustion of aluminum and steam. The influencing factors, such as the eccentric distance of the inlet of aluminum particles, particle size and steam inlet diameter, etc., are studied. The work described in this paper represents an attempt to the design of a vortex combustor in order to increase aluminum combustion efficiency.

Human number understanding is thought to rely on the analogue number system (ANS), working according to Weber’s law. We propose an alternative account, suggesting that symbolic mathematical knowledge is based on a discrete semantic system (DSS), a representation that stores values in a semantic network, similar to the mental lexicon or to a conceptual network. Here, focusing on the phenomena of numerical distance and size effects in comparison tasks, first we discuss how a DSS model could explain these numerical effects. Second, we demonstrate that DSS model can give quantitatively as appropriate a description of the effects as the ANS model. Finally, we show that symbolic numerical size effect is mainly influenced by the frequency of the symbols, and not by the ratios of their values. This last result suggests that numerical distance and size effects cannot be caused by the ANS, while the DSS model might be the alternative approach that can explain the frequency-based size effect.

Local polynomial smoother (LPS) is a weighted local least-squares nonparametric method. It provides a local Taylor series fit of the data at any location and can be directly used in a differential equation to provide a numerical scheme. In this article, we introduce this new nonparametric idea based on local polynomial smoother, for acquiring the numerical solution of the Bagley-Torvik fractional-order differential equations. Furthermore, this paper will present a numerical comparison with some methods, such as legendre operational matrix and pseudo-spectral method. The efficiency and accuracy of the LPS method are demonstrated by two numerical examples.

Recently, there has been a growing interest in the production of electricity from renewable energy sources (RES). The RES investment is characterized by uncertainty, which is long-term, costly, depend on feed-in-tariff and support schemes. In this paper, we address the real option valuation (ROV) of a solar power plant investment. The real option framework is investigated. This framework considers the renewable certificate price, furthermore the cost of delay between establishing and operating the solar power plant. The optimal time of launching the project and assess the value of deferred option are discussed. The new three stage numerical methods are constructed, the Lobatto3C-Milstein (L3CM) methods. The numerical methods are integrated with concept of Black-Scholes option pricing theory, and applied in option valuation for solar energy investment with uncertainty. The numerical results of L3CM, finite difference and Monte Carlo methods are compared to show the efficiency of our methods. Our data set refers to the Arab Republic of Egypt.

The present study aims to develop an efficient numerical method for computing the diffraction and radiation of water waves with horizontal long cylindrical structures, such as floating breakwaters. A higher-order scheme is used to discretize geometry of the structure as well as the relevant physical quantities. As the kernel of this method, Wehausen’s free-surface Green function is calculated by a newly-developed Gauss-Kronrod adaptive quadrature algorithm after elimination of its Cauchy-type singularities. To improve computational efficiency, a Chebyshev approximation approach is applied to a fast calculation of the Green function that needs evaluation thousands of times. In addition, OpenMP parallel technique is used to the formation of influence coefficient matrix, which significantly reduces CPU time. Finally, computations are performed on wave exciting forces and hydrodynamic coefficients for the long cylindrical structures, either floating or submerged. Comparison with other numerical and analytical methods demonstrates good performance of the present method.

The development of computational acoustics allows simulation of sound generation and propagation in complex environment. In particular, meshfree methods are widely used to solve acoustics problems through arbitrarily distributed field points and approximation smoothness flexibility. As a Lagrangian meshfree method, smoothed particle hydrodynamics (SPH) method reduce the difficulty in solving problems with deformable boundaries, complex topologies, or multiphase medium. The traditional SPH method has been applied in acoustic simulation. This study presents the corrective smoothed particle method (CSPM), which is a combination of SPH kernel estimate and Taylor series expansion. The CSPM is introduced as a Lagrangian approach to improve accuracy in solving acoustic wave equations in the time domain. Moreover, a boundary treatment technique based on the hybrid meshfree and finite difference time domain (FDTD) method is proposed to represent different acoustic boundaries with particles. To model sound propagation in pipes with different boundaries, soft, rigid, and absorbing boundary conditions are built with this technique. Numerical results show that the CSPM algorithm is consistent and demonstrates convergence with exact solutions. Main computational parameters are discussed, and different boundary conditions are validated to be effective for benchmark problems in computational acoustics.

In the present paper, laminar forced convection nanofluid flows in a uniformly heated horizontal tube were revisited by direct numerical simulations. Single and two-phase models were employed with constant and temperature-dependent properties. Comparisons with experimental data showed that the mixture model performs better than the single-phase model in the all cases studied. Temperature-dependent fluid properties also resulted in a better prediction of the thermal field. A particular attention was paid to the grid arrangement. The two-phase model was used then confidently to investigate the influence of the nanoparticle size on the heat and fluid flow with a particular emphasis on the sedimentation process. Four nanoparticle diameters were considered: 10, 42, 100 and 200 nm for both copper-water and alumina/water nanofluids. For the largest diameter dnp = 200 nm, the Cu nanoparticles were more sedimented by around 80 %, while the Al₂O₃ nanoparticles sedimented only by 2.5 %. Besides, it was found that increasing the Reynolds number improved the heat transfer rate, while it decreased the friction factor allowing the nanoparticles to stay more dispersed in the base fluid. The effect of nanoparticle type on the heat transfer coefficient was also investigated for six different water-based nanofluids. Results showed that the Cu-water nanofluid achieved the highest heat transfer coefficient, followed by C, Al₂O₃, CuO, TiO₂, and SiO₂, respectively. All results were presented and discussed for four different values of the concentration in nanoparticles, namely φ = 0, 0.6, 1 and 1.6%. Empirical correlations for the friction coefficient and the average Nusselt number were also provided summarizing all the presented results.

From June 18 to 19, 2013, the Ésera river in the Pyrenees, Northern Spain, caused widespread damage due to flooding as a result of torrential rains and sustained snowmelt. We estimate the contribution of snow melt to total discharge applying a snow energy balance to the catchment. Precipitation is derived from sparse local measurements and the WRF-ARW model over three nested domains, down to a grid cell size of 2 km. Temperature profiles, precipitation and precipitation gradient are well simulated, although with a possible displacement regarding the observations. Snowpack melting was correctly reproduced and verified in three instrumented sites, and according to satellite images. We found that the hydrological simulations agree well with measured discharge. Snowmelt represented 33% of total runoff during the main flood event and 23% at peak flow. The snow energy balance model indicates that most of the energy for snow melt during the day of maximum precipitation came from turbulent fluxes. This approach forecast correctly peak flow and discharge during normal conditions at least 24h in advance and could give an early warning of the extreme event 2.5 days before.

The numerical simulation of entropy generation at mixed natural free and forced convection flow in a square lid-driven cavity filled with an eletrically conducting binary fluids saturated porous media in the presence of a magnetic field is performed in this investigation. Both the top and bottom horizontal walls of the cavity are kept at constant and different temperatures while the left and right vertical walls are adiabatic and insulated. The magnetic field is applied in normal direction to the cavity. The Darcy model, including Brinkman term relative to viscous effects and Forchheimer term due to inertial forces, is used for the momentum equations, and the SIMPLER algorithm, based on the control volume finite-element method approach is used to solve the pressure –velocity coupling. The flow pattern and the heat transfer characteristics inside the cavity are presented in the form of isotherms, streamlines and isentropic lines and average Nusselt number trend. The entropy generation rate is determined for various values of Darcy number (10^-3 ≤ Da ≤ 1) and for a range of Hartmann number (0 ≤ Ha ≤ 10²). This study is limited to Pr = 7 which corresponds to water. It is observed that entropy generation is very affected by the variations of dimensionless control parameters considered.

By installing FBG sensors on the geogrids, smart geogrids can both reinforce and monitor the stability for geogrid-reinforced slopes. In this paper, a geogrid-reinforced sand slope model test was conducted in the laboratory and fiber Bragg grating (FBG) sensing technology is used to measure the strain distribution of the geogrid. Based on the model test, the performance of the reinforced soil slope is simulated by finite element software Midas-GTS, and the stability of the reinforced soil slope is analyzed by strength reduction method. The relationship between the geogrid strain and safety factor is set up. The results indicate that the measured strain and calculated results agree very well. The geogrid strain measured by FBG sensor can be applied to evaluate the stability of geogrid-reinforced sand slopes.

Energy dissipation is one of the most important factors in choosing stepped spillways. However, few studies are about energy dissipation with different horizontal face angles. In this paper, the realizable k-ε turbulence model was used to study the flow field and energy dissipation rates of different stepped spillways with five horizontal face angles in the skimming flow regions. It shows that the free water surfaces and energy dissipation rates change with horizontal face angles and the energy dissipation rates initially decrease and then increase as the horizontal face angle increases. And the horizontal face angles cause the great changes of flow field.

A three-dimensional, unstructured grid, hydrodynamic and suspended-sediment transport model (i.e., SELFE-SED) was developed to simulate temporal and spatial variations of suspended sediment and was applied to the subtropical subalpine Tsuei-Feng Lake (TFL) of Taiwan. The model was validated with measured water level and suspended‑sediment concentration in 2009, 2010, and 2011. The overall model simulation results are in quantitative agreement with the observational data. The validated model was then applied to explore the most important parameter that affects the suspended-sediment concentration and to investigate the effect of wind stress on the mean current and suspended‑sediment distribution in this shallow lake. Modeling results of sensitivity analysis reveal that the settling velocity is a crucial parameter and erosion rate is less important in the suspended-sediment transport model. Remarkable lake circulation was found based on the strength of wind speed and wind direction. Strong wind would result in higher mean current in the top layer and suspended-sediment distribution in the top and bottom layers. This study demonstrated that the wind stress played a significant influence on mean circulation and suspended-sediment transport in a shallow lake.

The difficult current environmental situation, caused by construction industry residues containing ceramic materials could be improved by using these materials as recycled aggregates in mortars, with their processing causing a reduction in their use in landfill, contributing to recycling and also minimizing the consumption of virgin materials. Although some research is currently being carried out into recycled mortars, little is known about their stress-strain (σ-ε); therefore this work will provide the experimental results obtained from recycled mortars with recycled ceramic aggregates (with contents of 0, 10, 20, 30, 50 and 100%), such as: the density, the compression strength, as well as the σ-ε curves representative of their behavior. The values obtained from the analysis process of the results are those of: σ (elastic ranges and failure maximum), ε (elastic ranges and failure maximum), and Resilience and Toughness; in order to finally obtain, through numerical analysis, the equations to predict their behavior (related to their recycled content). At the end of the investigation it is established that mortars with recycled ceramic aggregate contents of up to 20% could be assimilated just like mortars with the usual aggregates, and the prediction equations produced could be used in cases of similar applications.

The effects of the surface waves generated by the wind have a significant effect on the currents. A wave current coupled model plays an important role in the design of offshore structures. The interaction between fluids such as incompressible ocean waves and current and offshore structures is significant with many real-time applications in offshore engineering. These coupled models can be applied to Offshore Floating Production Operating and offloading (FPSO), Wind or current turbines and offshore pipelines. The complex issues related to the design are analyzed by using Computational Fluid Dynamics, which requires an investigation of the multiphase flow between wave and current and the structure which is considered restrictive due to the computational cost. If viscous effects are neglected then the single-phase flow models have been recommended, where wave-current interaction have been modelled successfully. Models have been developed where velocities and pressure are computed and the results can be verified with the experimental results available in the literature. In this study the existing numerical methods, mesh types are discussed along with their coupling methods. Here single-phase and multiphase models with small and medium movement are reviewed and their applications are highlighted.

Subsurface temperature is the key parameter in geothermal exploration. An accurate estimation of the reservoir temperature is of high importance and usually done either by interpolation of borehole temperature measurement data or numerical modeling. However, temperature measurements at depths which are of interest for deep geothermal applications (usually deeper than 2 km) are generally sparse. A pure interpolation of such sparse data always involves large uncertainties and usually neglects knowledge of the 3D reservoir geometry or the rock and reservoir properties governing the heat transport. Classical numerical modeling approaches at regional scale usually only include conductive heat transport and do not reflect thermal anomalies along faults created by convective transport. These thermal anomalies however are usually the target of geothermal exploitation. Kriging with trend does allow including secondary data to improve the interpolation of the primary one. Using this approach temperature measurements of depths larger than 1,000 m of the federal state of Hessen/Germany have been interpolated in 3D. A 3D numerical conductive temperature model was used as secondary information. This way the interpolation result reflects thermal anomalies detected by direct temperature measurements as well as the geological structure. This results in a considerable quality increase of the subsurface temperature estimation.

Major marine construction projects, resulting in the release of sediments, are subject to environmental assessment and other regulatory approval processes. An important tool used for this is the development of specialized numerical methods for these marine activities. An integrated set of numerical methods addresses four distinct topics: the near-field release and mixing of suspended sediments into the water column (i.e. the initial dilution zone); the transport of the suspended sediments under the influence of complex ocean currents in the far-field; the settling of the transported suspended sediments onto the seabed; and the potential for resuspension of the deposited sediments due to sporadic occurrences of unusually large near-bottom currents. A review of projects subjected to environmental assessment in the coastal waters of British Columbia from the year 2006 to 2017, is presented to illustrate the numerical models being used and their ongoing development. Improvements include higher resolution model grids to better represent the near-field, the depiction of particle size dependent vertical settling rates and the computation of resuspension of initially deposited sediments, especially in relation to temporary subsea piles of sediments arising from trenching for marine pipelines. The ongoing challenges for this numerical modeling application area are also identified.

Fiber-metal laminates (FMLs) such as Kevlar reinforced aluminum laminate (ARALL), Carbon reinforced aluminum laminate (CARALL), and Glass reinforced aluminum laminate (GLARE) offer great potential for weight reduction applications in automobile and aerospace construction. In order to investigate the feasibility for utilizing such materials in the form of laminates, sheet hydro-bulging tests are studied under the condition of uniform blank holder force for three-layered aluminum and aluminum-composite laminates using orthogonal carbon and Kevlar as well as glass fiber in the middle. The experimental results validate the finite element results and they exhibited that the forming limit of glass fiber in the middle is the highest among the studied materials, while carbon fiber material performs the worst. Furthermore, the crack modes are different for the three kinds of fiber materials investigated in the research. This study provides fundamental guidance for the selection of multi-layer sheet materials in the future manufacturing field.

ABSTRACT. In this paper, we suggest a biomathematical numerical method for analysing mRNA nucleotides sequences based on UA/CG Fibonacci numbers proportions. This method is used to evaluate then compare the spike genes related to the main SARS-CoV2 VARIANTS currently circulating within the world population. The 10 main results proposed to be reproduced by peers are: 1/ SARS-CoV2 genome and spike evolution in one year 2020-2021. 2/ SARS-CoV2 Origins. 3/ Comparing 11 reference variants spikes. 4/ analysing 32 CAL.20C California variant patients spikes. 5/ Toward a meta mRNA Fibonacci gene end message code. 6/ Analysing S501 UK, S484 South Africa and « 2 mutations » INDIA variants. 7/ Suggesting a possible variants spike mRNA palindrome symmetry metastructure improving mRNA stability then infectiousness. 8/ Analysing Fibonacci Metastructures in the mRNA coding for the vaccines PFIZER and MODERNA. 9/ Does the CG-rich modification of the synonymous codons of the spikes of the 2 mRNA vaccines affect the expression and quantity of SARS-CoV2 antibodies? 10/ The exceptional case of the Brazilian variant P.1. Particularly, we suggest the following conjecture at mRNA folding level : CONJECTURE of SARS-CoV2 VARIANTS: The growth of long Fibonacci structures in the shape of "podiums" for almost all of the variants studied (UK, California, South Africa, India, etc.) suggests the probable folding of the Spike mRNA in the form of a "hairpin", which can strengthen the cohesion and the lifespan of this mRNA. Finally, we show that these kinds of Fibonacci matastructures disapear TOTALLY by analysing the published mRNA sequences of PFIZER and MODERNA vaccines. One fact is certain, the two mRNAs of the Moderna and Pfizer vaccines will result in a low functionality of the spike vaccine. This is because their designers by seeking greater stability, have doped to build CG rich sequences which, as soon as they are inserted into the human host, will, paradoxically, seek to mutate, like SARS-CoV2 variants, towards CG ==> UA forms in order to improve their STABILITY and LIFETIME. We conclude using new biomathematics theoretical methods (Master code and numerical standing waves), and comparing the Spikes of the two vaccines Moderna and Pfizer, that there will be very probable differences in stability and shelf life of the two respective mRNAs vaccines. However, “State of the Art” analyzes will disclose that their two protein sequences are strictly identical. By modified their synonymous codons using different strategies, no one can guarantee that the quantity of antibodies generated will be identical in the two cases. We wish to draw attention to the great ADAPTATION power - at the global scale of their genomes - of the most infectious VARIANTS, such as the BRAZIL 20J / 501Y.V3 variant (P.1). This is very worrying for the VACCINES <==> VARIANTS run: We demonstrate how the Brazilian variant P.1 which becomes uncontrollable in Brazil in April 2021 has a level of organization of long metastructures of 17,711 bases covering the genome which is 3.6 more important than that of the 2 reference genomes SARS-CoV2 and worldwide D614G. We suggest that this high level of overall structure of this variant contributes to the stability of this genome and, might explain its greater contagiousness.

Marine cold air outbreaks (MCAOs) are large-scale phenomena in which cold air masses are advected over open ocean. It is well-known that these events are linked to the formation of polar lows and other mesoscale phenomena associated with high wind speeds, and that they therefore in some cases represent a hazard to maritime activities. However, it is still unknown whether MCAOs are generally conducive to higher wind speeds than normal. Here this is investigated by comparing the behaviour of ocean surface wind speeds during MCAOs in three atmospheric reanalysis products with different horizontal grid spacings, along with case studies using a convection-permitting numerical weather prediction model. The study regions are the Labrador Sea and the Greenland–Iceland–Norwegian (GIN) Seas, where MCAOs have been shown to be important for air–sea interaction and deep water formation. The main findings are: 1) Wind speeds during the most extreme MCAO events are stronger than normal and higher than wind speeds during less severe events; 2) The peak times of MCAO usually occur when baroclinic waves pass over the regions; and 3) Reanalyses with grid spacings of more than 50 km appear to underestimate winds driven by the large ocean–atmosphere energy fluxes during MCAOs. It is also shown that while the strong wind episodes during MCAOs generally last for just a few days, MCAOs can persist for up to 50 days. These findings demonstrate that it would be worthwhile to forecast MCAOs, and that it might be possible to do this beyond the standard weather forecasting range of up to 10 days.

Based on the geologic conditions of natural gas hydrate (NGH) accumulations in the Shenhu area, northern slope of the South China Sea, this paper used a method of combining reservoir stimulation technique (RST) with depressurization to investigate the conversion efficiency of marine NGH reservoirs in different intrinsic permeability and initial NGH saturation conditions, and analyze the influence of variably-stimulation effect on marine NGH conversion efficiency in different accumulation conditions, provided a reference scheme for improving the NGH conversion efficiency in the Shenhu area. In this work, we performed calculations for the variations in CH4 production rate and cumulative volume of CH4 in different initial NGH saturation, intrinsic permeability and stimulation effect conditions, variance analysis and range analysis methods were used to analyze the significance of these key factors and their interaction, and investigated the sensitivity of stimulation effect on NGH conversion efficiency, respectively. The simulation results showed that stimulation effect have a significant influence on NGH conversion efficiency, and the influence of interaction between these three factors were not obvious. Possibly most important, we clarified an optimum stimulation effect for higher NGH conversion efficiency under different accumulation conditions, especially in low-permeability and high-saturation, which corresponding stimulation effect were different.

Fractional differential mathematical model unfolding the dynamics of the COVID-19 pandemic in India is presented and explored in this paper. The purpose of this study is to estimate the future outbreak of disease and potential control strategies using mathematical models in India as a whole country as well as in some of the states of the country. This model is calibrated based on reported cases of infections over the month of April 2020 in India. We have used iterative fractional complex transform method to find approximate solutions of the model having modified Riemann Liouville fractional differential operator. We have also carried out a comparative analysis between actual and estimated cumulative cases graphically, moreover, most sensitive parameters for basic reproduction number$(R_0)$ are computed and their effect on transmission dynamics of COVID-19 pandemic is investigated in detail.

In this paper an integrated methodology for the coupling between 1D- and 3D-CFD simulation codes is presented, which has been developed to support the design and calibration of new diesel engines. The aim of the proposed methodology is to couple 1D engine models, which may be available in the early-stage engine development phases, with 3D predictive combustion simulations, in order to obtain reliable estimates of engine performance and emissions for newly designed automotive diesel engines. The coupling procedure features simulations performed in 1D-CFD by means of GT-SUITE and in 3D-CFD by means of Converge, executed within a specifically designed calculation methodology. An assessment of the coupling procedure has been performed by comparing its results with experimental data acquired on an automotive Diesel engine, considering different working points including both part load and full load conditions. Different multiple injection schedules have been evaluated for part-load operation, including pre and post injections. The proposed methodology, featuring detailed 3D chemistry modeling, was proven to be capable to properly assess pollutant formation, specifically to estimate NOx concentrations. Soot formation trend was also well-matched for most of the explored working points. The proposed procedure can therefore be considered as a suitable methodology to support the design and calibration of new Diesel engines, thanks to its ability to provide reliable engine performance and emissions estimations from the early-stage of a new engine development.

Liquid propellants are fast becoming attractive for pulsed plasma thrusters due to their high efficiency and low contamination issues. However, the complete plasma interaction and acceleration processes are still not very clear. Present paper develops a multi-layer numerical model for liquid propellant PPTs. The model proposes a possible acceleration mechanism for liquid fed pulsed plasma thrusters and accurately predicts the propellant utilization capabilities and estimations for the fraction of propellant gas that is completely ionized and accelerated to high exit velocities. Validation of the numerical model and the assumptions on which the model is based on is achieved by comparing the experimental results from two different liquid-fed thrusters developed at the University of Tokyo. Results show 50% of liquid propellant injected is completely ionized and accelerated to high exit velocities (> 50Km/s), whereas, neutral gas contribute to only 7% of the total specific impulse and accelerated to low exit velocity (< 4Km/s). The model shows an accuracy up-to 92%. Optimization methods are briefly discussed to ensure efficient propellant utilization and performance. The model acts as a tool to understand the background physics and to optimize the performance for liquid-fed PPTs.

The United Arab Emirates (UAE) is significantly dependent on desalinated water and groundwater resource, which is expensive and highly energy intensive. Despite the scarce water resource, only 54% of the recycled water was reused in 2015. In this study, an “Oasis” complex comprised of Sustainable Farming Compartments (SFCs) was proposed for reusing treated wastewater to decrease the ambient temperature of the SFC via an evaporative cooling system. A prototype SFC with half the original scale (width = 1.8 m, depth = 1.5 m, front height = 1.2 m back height = 0.9 m) was designed, built, and tested in an environmentally controlled laboratory and field site to evaluate the feasibility and effectiveness of the SFC under the climatic conditions in Abu Dhabi. Based on the experimental results, the temperature drops obtained from the SFC in the laboratory and field site were 5 ̊C at initial relative humidity of 60% and 7- 15 ̊C at initial relative humidity of 50%, respectively. An energy simulation using dynamic numerical simulations was performed in comparison to the results of the experiment. The energy-based dynamic simulation shows good agreement with the experimental results. The total power consumption of the SFC system was approximately three and a half times lower than that of an electrical air conditioner.

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.

Information is the currency of life, but the origin of prebiotic information remains a mystery. We propose transitional pathways from the cosmic building blocks of life to the complex prebiotic organic chemistry that led to the origin of information systems. The prebiotic information system, specifically the genetic code, is segregated, linear, and digital and probably appeared during biogenesis four billion years ago. In the peptide/RNA world, lipid membranes randomly encapsulated amino acids, RNA, and protein molecules, drawn from the prebiotic soup, to initiate a molecular symbiosis inside the protocells. This endosymbiosis led to the hierarchical emergence of several requisite components of the translation machine: tRNAs, aaRS, mRNAs, and ribosomes. When assembled in the right order, the translation machine created biosynthetic polypeptides, a process that transferred information from mRNAs to proteins. This was the beginning of the prebiotic information age. The molecular attraction between tRNA and amino acids led to different stages of the translation machines and the genetic code. tRNA is an ancient molecule that designed and built mRNA for storing the information of its cognate amino acid. Each mRNA strand became the storage device for the genetic information that encoded the amino acid sequences in triplet nucleotides. As information appeared in the digital languages of the codon within mRNA, and the genetic code for protein synthesis evolved, the prebiotic chemistry then became more organized and directional. The origin of the genetic code is enigmatic; herein we propose an evolutionary explanation: the demand for a wide range of specific enzymes in the peptide/RNA world was the main selective pressure for the origin of information-directed protein synthesis. We review three main concepts on the origin and evolution of the genetic code: the stereochemical theory, the coevolution theory, and the adaptive theory. These three theories are compatible with our coevolution model of the translation machines and the genetic code. We suggest biosynthetic pathways as the origin of the specific translation machines which provided the framework for the origin of the genetic code. During translation, the genetic code developed in three stages coincident with the refinement of the translation machines: GNC code developed by the pre-tRNA/pre-aaRS /pre-mRNA machine, SNS code by the tRNA/aaRS/mRNA machine, and finally the universal genetic code by the tRNA/aaRS/mRNA/ribosome machine. Our hypothesis provides the logical and incremental steps for the origin of the programmed protein synthesis. In order to understand the prebiotic information system better, we converted letter codons into numerical codons in the Universal Genetic Code Table. We have developed a software called CATI (Codon-Amino Acid-Translator-Imitator) to translate randomly chosen numerical codons into corresponding amino acids and vice versa. This conversion has granted us insight into how the translation might have worked in the peptide/RNA world. There is great potential in the application of numerical codons to bioinformatics such as barcoding, DNA mining, or DNA fingerprinting. We constructed the likely biochemical pathways for the origin of translation and the genetic code using the Model-View-Controller (MVC) software framework, and the translation machinery step-by-step. Using AnyLogic software we were able to simulate and visualize the entire evolution of the translation machines and the genetic code. The results indicate that the emergence of the information age from the peptide/RNA world was a watershed event in the origin of life about four billion years ago.

Currently, the use of numerical models for reproducing the landscape evolution of a river basin is part of the day-by-day research activities of fluvial engineers and geomorphologists. However, despite landscape modelling is based on a rather long tradition, and scientists and practitioners are trying to schematize the processes involved in the evolution of a landscape since decades, there is still the need for improving both the knowledge of the physical mechanisms and their numerical coding. The present review focuses on the first aspect, discussing the main components of a landscape evolution model and their more common schematizations, presenting possible open questions to be addressed towards an improvement of the reliability of such kind of models in describing the fluvial geomorphology.

Wildland fires are responsible for large socio-economic impacts. Fires affect the environment, damage structures, threaten lives, cause health issues, and involve large suppression costs. These impacts can be mitigated via accurate fire spread forecast to inform the incident management team. We show that a fire forecast system based on a numerical weather prediction (NWP) model coupled with a wildland fire behavior model can provide this forecast. This is illustrated with the Chimney Tops II wildland fire responsible for large socio-economic impacts. The system is run at high horizontal resolution (111 m) over the region affected by the fire to provide a fine representation of the terrain and fuel heterogeneities and explicitly resolve atmospheric turbulence. Our findings suggest that one can use the high spatial resolution winds, fire spread and smoke forecast to minimize the adverse impacts of wildland fires.

The aim of this study is to model hydrodynamic processes of the Istanbul Strait with its stratified flow characteristic and calibrate the most important parameters using local and global search algorithms. For that two open boundary conditions are defined, which are in the North and South part of the Strait. Observed bathymetric, hydrographic, meteorological and water level data are used to set up the Delft3D-FLOW model. First, the sensitivities of model parameters on the numerical model outputs are assessed using PEST toolbox. Then, the model is calibrated based on the objective functions focusing on the flowrates of upper and lower layers. The salinity and temperature profiles of the Strait are only used for model validation. The results show that the calibrated model outputs of Istanbul Strait are reliable and consistent with the in-situ measurements. The sensitivity analysis reveals that the Spatial Low-Pass Filter Coefficient, Horizontal Eddy Viscosity, Prandtl-Schmidt Number, Slope in log-log Spectrum and Manning Roughness Coefficient are most sensitive parameters affecting flowrate performance of the model. The agreement between observed salinity profiles and simulated model outputs is promising whereas the match between observed and simulated temperature profiles is weak showing that the model can be improved particularly for simulating the mixing layer.

Rigid-hulled inflatable boats are extremely practical and popular nowadays. They offer a effective conciliation among usability and costs. Their stable and seaworthy behaviour is guaranteed by performing hydroplaning hulls coupled with unsinkable inflated tubes. At the same time, their design is often based on tradition and preconceptions. Rarely, the design assumptions are validated by the reality or, even, by deeper investigations. In this article, both numerical methods and experimental mechanics techniques are proposed as an essential way for supporting the designers in their decisive tasks. Three different situations are detailed where a numerical or an experimental approach shows its benefit inside the engineering design process: firstly permitting to investigate the behaviour of materials driving the fiberglass selection; then measuring the levels of stress and strain in the hull during sailing; finally, using all available information as a base for developing numerical models of the hull slamming in waves. Even if the discussion is focused on a rigid inflatable boat, large part of its considerations is relevant beyond this special case.

A galvalume corrugated sheet was utilized as formwork for a reinforced concrete beam in flexure. A numerical model was validated to the experimentally obtained data, and further adopted to simulate the behavior of this composite structure under elevated temperatures. The properties and constitutive stress-strain data of the basic materials were obtained from experiments, and superimposed into the finite element model. The study concluded that the load carrying capacity of the member decreased as a direct function on temperature increase, and the cracking moment was very sensitive to the temperature fluctuation. The elevated temperatures also altered the failure mode.

A microstructure-based model is developed to study the effective anisotropic properties (magnetic, dielectric or thermal) of two phase particle-filled composites. The Green’s function technique and the effective field method are used to derive theoretically the homogenized (averaged) properties for a representative volume element containing isolated inclusion and an infinite, chain-structured particles. Those results are compared with the FE approximations conducted for the assumed representative volume element. In addition, as a special case, the Maxwell–Garnett model is retrieved when particle interactions are not taken into consideration. We shall also give some information on the optimal design of the effective anisotropic properties taking into account the shape of magnetic particles.

Cloud-based CNC is an emerging paradigm of Industry 4.0 where computer numerical control (CNC) functionalities are moved to the cloud and provided to manufacturing machines as a service. Among many benefits, C-CNC allows manufacturing machines to leverage advanced control algorithms running on cloud computers to boost their performance at low cost, without need for major hardware upgrades. However, a fundamental challenge of C-CNC is how to guarantee safety and reliability of machine control given variable Internet quality of service, especially on public Internet networks. We propose a three-tier redundant architecture to address this challenge. We then prototype tier one of the architecture on a 3D printer successfully controlled via C-CNC over public Internet connections, and discuss follow-on research opportunities.

Some cognitive functions, shared by humans and certain animals, were acquired early in the course of phylogeny and, in humans, are operational in their primitive form shortly after birth. This is the case for the quantification of discrete objects. The further phylogenetic evolution of the human brain allows such functions to be reconstructed in a much more sophisticated way during child development. Certain functional characteristics of the brain (plasticity, multiple cognitive processes involved in the same response, interactions and substitution relationships between those processes) provide degrees of freedom that open up the possibility of different pathways of reconstruction. The within- and between-individual variability of these developmental pathways offers an original window on the dynamics of development. Here, I will illustrate this theoretical approach to cognitive development—which can be called “reconstructivist” and “pluralistic”—using children's construction of number as an example.

A method for solving the turbulence equations embedded in the sigma ocean general ocean circulation model is proposed. Like the general circulation model, the turbulence equations are solved using the splitting method by physical processes. The turbulence equations are split into two main stages describing transport-diffusion and generation-dissipation processes. Parameterization of turbulence in the framework of equations allows, at the generation-dissipation stage, to use both numerical and analytical solutions and to ensure high efficiency of the algorithm. The results of large-scale ocean dynamics simulation taking into account the parameterization of vertical turbulent exchange are considered. Numerical experiments were carried out using k-omega turbulence model embedded to the Institute of Numerical Mathematics Ocean general circulation Model (INMOM). Both the circulation and turbulence models are solved using the splitting method with respect to physical processes. The coupled model is used to simulate the hydrophysical fields of the North Atlantic and Arctic Oceans for 1948--2009. The model has a horizontal resolution of 0.25 degree and 40 sigma-levels along the vertical. The sensitivity of the solution to the changes in mixing parameterization is studied. Experiments demonstrate that taking into account the climatic annual mean buoyancy frequency improves the reproduction of large-scale ocean characteristics. There is a positive effect of Prandtl number variations for reproducing the upper mixed layer depth. The experiments also demonstrate the computational effectiveness of the proposed approach in solving the turbulence equations.

First of all the paper presents a solid solution of piezoelectric ceramic that was synthesized according to the general formula Pb(1-x)Srx(Ti0.48Zr0.52)(1-y)NbyO3 with x = 0.05 and y = 0.02, using wet ceramic processing technology and using an oxide mix as prime material. The effects of dopants (Sr2+ and Nb5+) on phase constitution, on microstructure and on the dielectric and piezoelectric properties were determinate. The Zr/Ti ratio was chosen near the morphotropic phase boundary of the PZT system in studied composition. The XRD data revealed that the PZT doped composition had tetragonal perovskite structure. Secondly the paper presents the design of the novel piezo-motor based on a surface wave which translates the linear extension of different piezoelectric segments of a piezoelectric cylinder into a rotational bending movement. This rotational bending of the piezoelectric cylinder is then transformed into a continuous rotation of the rotor through a calculated contact. The design of the motor takes advantage of the high piezoelectric constants of the developed material in an optimal way in order to increase the energetic efficiency. A brief mathematical model of electromechanical answer of piezoelectric materials is presented as it was used in the modeling of the material during the numerical simulations.

Wave energy converter (WEC) arrays deployed in coastal regions may create physical disturbances potentially resulting in environmental stresses. Presently, limited information is available on the nature of these physical disturbance or the resultant effects. A quantitative Spatial Environmental Assessment Tool (SEAT) for evaluating potential effects of wave energy converter (WEC) arrays on nearshore hydrodynamics and sediment transport is presented for the central Oregon coast (USA) through coupled numerical model simulations of an array of WECs. Derived climatological wave conditions were used as inputs to the model to allow for the calculation of risk metrics associated with various hydrodynamic and sediment transport variables such as maximum shear stress, bottom velocity, and change in bed elevation. The risk maps provided simple, quantitative, and spatially-resolved means of evaluating physical changes in the vicinity of a hypothetical WEC array in response to varying wave conditions. Near-field risk of sediment mobility was determined to be moderate in the lee of the densely spaced array, where the potential for increased sediment deposition could result in benthic habitat alteration. Modifications to the nearshore sediment deposition and erosion patterns were observed near headlands and topographic features, which could have implications for littoral sediment transport. The results illustrate the benefits of a risk evaluation tool for facilitating coastal resource management at early market marine renewable energy sites.

This paper presents a hydrodynamic numerical model for a pitching wave energy converter (WEC). The model uses potential wave theory and is based on Cummins' equation, with nonlinear hydrostatic restoring stiffness and excitation forces based on instantaneous body position and water surface elevation. The numerical model can include non-linear forces, like quadratic drag, power-take-off and other forces that may account for unknown viscous effects observed in experiments. The paper discusses the applicability and limitations of the code and presents the cases where assumptions and simplifications can be made. The goal is to conclude on the simplest, yet accurate, version of the model by evaluating its accuracy using experimental data. The case study for validation is Floating Power Plant's (FPP's) WEC [1]. In the full-scale commercial project, FPP's device consists of a semisubmersible platform hosting a wind turbine (5-8MW) and 4 WECs, each one connected to the platform by a rotation shaft. Due to the configuration of the platform, strong interactions occur between the WECs and the structure, as they are very closely spaced. In order to validate a numerical model able to simulate these hydrodynamic interactions, wave basin experiments with a similar but simplified setup were performed.

5G heterogeneous network overlaid by millimeter-wave (mmWave) access employs mmWave meshed backhauling as a promising cost-efficient backhaul architecture. Due to the nature of mobile traffic distribution in practice which is both time-variant and spatially non-uniform, dynamic construction of mmWave meshed backhaul is prerequisite to support the varying traffic distribution. Focusing on such scenario of outdoor dynamic crowd (ODC), this paper proposes a novel method to control mmWave meshed backhaul for efficient operation of mmWave overlay 5G HetNet through Software-Defined Network (SDN) technology. Our algorithm is featured by two functionalities, i.e., backhauling route multiplexing for overloaded mmWave small cell base stations (SC-BSs) and mmWave SC-BSs’ ON/OFF status switching for underloaded spot. In this paper, the effectiveness of the proposed meshed network is confirmed by both numerical analyses and experimental results. Simulations are conducted over a practical user distribution modeled from measured data in realistic environments. Numerical results show that the proposed algorithm can cope with the locally intensive traffic and reduce energy consumption. Furthermore, a WiGig (Wireless Gigabit Alliance certified) device based testbed is developed for Proof-of-Concept (PoC) and preliminary measurement results confirm the proposed dynamic formation of the meshed network’s efficiency.

In this paper we introduce the hypo-q-norms on a Cartesian product of algebras of bounded linear operators on Hilbert spaces. A representation of these norms in terms of inner products, the equivalence with the q-norms on a Cartesian product and some reverse inequalities obtained via the scalar reverses of Cauchy-Buniakowski-Schwarz inequality are also given. Several bounds for the norms δ_p, ϑ_p and the real norms η_r,p and θ_r,p are provided as well.

There are many criterion to generalize the concept of numerical radius; one of the most recent interesting generalization is what so called the generalized Euclidean operator radius. Simply, it is the numerical radius of multivariable operators. In this work, several new inequalities, refinements and generalizations are established for this kind of numerical radius.

FRAP technique have been used for decades to measure movements of molecules in 2D. Data obtained by FRAP experiments in cell plasma membranes are assumed to be described through means of two parameters, a diffusion coefficient D (as defined in a pure Brownian model) and a mobile fraction M. Nevertheless, it has also been shown that recoveries can be nicely fit using anomalous sub-diffusion. FRAP at variable radii has been developed using the Brownian diffusion model to access geometrical characteristics of the surrounding landscape of the molecule. Here we performed numerical simulations of continuous time random walk (CTRW) anomalous subdiffusion and interpreted them in the context of variable radii FRAP. These simulations were compared to experimental data obtained at variable radii on living cells using the PH domain of the membrane binding protein EFA6 (exchange factor for ARF6, a small G protein). This protein domain is an excellent candidate to explore the structure of the interface between cytosol and plasma membrane in cells. By direct comparison of our numerical simulations to the experiments, we show that this protein does not exhibit anomalous diffusion in BHK cells. The non Brownian PH-EFA6 dynamics observed here is more related to spatial heterogeneities such as cytoskeleton fences effects.

Three-dimensional (3D) printing technologies are transforming the design and manufacture of components and products across a variety of disciplines, however their application in the construction industry is still limited. Material deposition processes can achieve infinite geometries and have advanced from rapid prototyping and model-scale markets to their application in fabricating functional products, large objects and the construction of full-scale buildings. Many international projects have recently been realized and the construction industry is beginning to utilise these dynamic technologies. The potential advantages for integrating 3D printing into house construction are significant, these include the capacity for mass customization of designs and parameters for functional and aesthetic purposes, reduction in construction waste from highly precise material placement, and the use of recycled waste products in layer deposition materials. With the ultimate goal of improving construction efficiency and decreasing building costs, applying Strand7 Finite Element Analysis software, a numerical model was designed specifically for 3D printing in a cement mix incorporated with recycled waste product High Density Polyethylene (HDPE) and found that construction of an arched truss-like roof was structurally feasible without the need for steel reinforcements. The lab sizes prototypes were manufactured based on the destined numerical model by using a 3D printing technology. Currently available 3D printing technologies can be adopted for building construction and this paper discusses the applications, advantages, limitations and future directions of 3D printing as an innovative and viable solution for affordable house construction.

In this work, some numerical radius inequalities based on the recent Dragomir extension of Furuta's inequality are obtained. Some particular cases are also provided.

This paper presents the results of aerodynamical performance research focused on maintaining the thermal comfort and increasing the energy efficiency of a typical social housing unit located in a high-density urban area. Bioclimatic design strategies are used to develop a sustainable and economic technology in existing housing clusters in Mexico City. A full-scale prototype, built on campus facilities, was used to study the flow conditions around the design. All scaled prototypes implement similar criterion. Furthermore, a scaled prototype is evaluated within a low speed wind tunnel installation. Additionally, numerical simulations were performed at transient state based on previous physical measurements and historical local climatic conditions to determine preferable modifications.

A simple optimization strategy for the computation of 3D finite-differencing kernels on many-cores architectures is proposed. The 3D finite-differencing computation is split direction-by-direction and exploits two level of parallelism: in-core vectorization and multi-threads shared-memory parallelization. The main application of this method is to accelerate the high-order stencil computations in numerical relativity codes.

During the recent decades, global warming by greenhouse gas evolution has attracted worldwide attention and evermore strict regulations thereon have become institutionalized as international policies. In the process, more environment-friendly power generation technologies have been developed utilizing fossil fuels with a view to timely commercialization. As one of such “clean coal” technology, Integrated Gasification Fuel Cell system is a promising power generation means and COS hydrolysis reactor is installed downstream of coal syngas to remove acidic gas constituents such as H2S and COS. The most significant design parameters affecting performance of the COS hydrolysis reactor were selected to be GHSV, (catalytic) reaction temperature, and length ratio and numerical modeling was performed considering heat and fluid flow transfer as well as chemical reaction kinetics. Effect of the selected design parameters on the variation of conversion rate and reactant gas mixture concentration were comprehensively investigated to predict performance of COS hydrolysis reactor. Stochastic modeling of reactor performance was finally performed using Monte Carlo simulation and linear regression fitting.

In the present study, a simple approximation expression is given for the relationship between the period and amplitude of a simple pendulum under magnetic action. The analytical solution presented for the given problem. Two numerical quadrature methods Simpson's and Boole's method were utilized to demonstrate a new approximation of the problem. The results of the numerical quadrature have been compared to the exact solution. Absolute and relative mistakes of the problem have been presented. The Matlab program 2013R has created a numerical method that is used to analyze the outcome, It has been determined that the comparison's outcomes attest to the method's suitability and correctness. Moreover, the results show that numerical solution is suitable for the problem.

Due to the accuracy of numerical calculation of fluid flow inside a hydrocyclone can be obtained using Computational Fluid Dynamics (CFD), highly modified super computers are used to simulate the fluid flow and track particle motion inside a hydrocyclone. This paper deals with the numerical study using three multiphase models viz. Volume of fluid, Mixture and Eulerian model. The dimensions of the hydrocyclone taken into consideration for numerical analysis is same as considered by Rajamani. Validation of axial and tangential velocities at different strategically decided axial stations, RMS axial and tangential velocity profiles of the hydrocyclone is done using Reynolds Stress Model (RSM). The hydrocyclone model has been designed in Creo 3.0 using the same dimensions which later was imported to CFD for meshing. Fine hexagonal mesh numbering up to 5 lacs were constructed to obtain optimum results. Fluid flow was allowed to be developed in ANSYS FLUENT 16.2. Entire simulation took 96 hours to generate results and track particle movements inside the hydrocyclone. The particle tracking has been done using three multiphase model. The first being the volume of fluid was used for validation purposes and the comparison of the Mixture and Eulerian model are the basic focus of this research work. Conclusive results indicate that usage of different multiphase model does not result in variation in particle motion. The slight variation in grade efficiency values are hardly noticeable. The Mixture model and Eulerian model predict lower separation efficiency as compared with Volume of fluid multiphase model.

A cosmological model was developed using the equation of state of photon gas, as well as cosmic time. The primary objective of this model is to see if determining the observed rotation speed of galactic matter is possible, without using dark matter (halo) as a parameter. To do so, a numerical application of the evolution of variables in accordance with cosmic time and a new state equation was developed to determine precise, realistic values for a number of cosmological parameters, such as energy of the universe , cosmological constant Λ, curvature of space k, energy density , age of the universe etc. Some assumptions were put forth in order to solve these equations. The current version of the model partially explains several of the observed phenomena that raise questions. Numerical application of the model has yielded the following results, among others: Initially, during the Planck era, at the very beginning of Planck time, , the universe contained a single photon at Planck temperature Planck energy in the Planck volume. During the photon inflation phase (before characteristic time ~10^-9 ), the number of original photons (alphatons) increased at each unit of Planck time and geometrical progression ~n³, where n is the quotient of cosmic time over Planck time . Then, the primordial number of photons reached a maximum of N~10⁸⁹, where it remained constant. It's primordial photons (alphatons) are still present today and represents the essential of the energy contained in the universe via the cosmological constante . Such geometric growth in the number of photons can bring a solution to the horizon problem through exchange and a photon energy volume that is in phase with that of the the volume energy of the universe. The age of the universe in cosmic time that is in line with positive energy conservation (in terms of conventional thermodynamics) and the creation of proton, neutron, electron, and neutrino masses, is ~76 [Gy] (observed . The predicted total mass (p, n, e and n), based on the Maxwell-Juttner relativistic statistical distribution, is ~7x10⁵⁰ . The predicted cosmic neutrino mass is ≤ 8.69x10^-32 (≤ 48.7 ) if based on observations of SN1987A. The temperature variation of the cosmic microwave background (CMB), as measured by Planck, can be said to be partially due to energy variations in the universe (DU/U) during the primordial baryon synthesis (energy jump from the creation of protons and neutrons). In this model, what is usually referred to as dark energy actually corresponds to the energy of the universe that has not been converted to mass, and which acts on the mass created by the energy-mass equivalence principle and the cosmological gravity field, , associated with the cosmological constant, which is high during primordial formation of the galaxies (<1 [Gy]). A look at the Casimir effect makes it possible to estimate a minimum Casimir pressure and thus determine our possible relative position in the universe at cosmic time 0,1813 (t₀/t_W=13,8 /76,1 ). Therefore, from the observed age of 13,8 [Gy], we can derive a possible cosmic age of 76,1 [Gy]. That energy of the universe, when taken into consideration during the formation of the first galaxies (< 1 [Gy]), provides a relatively adequate explanation of the non-Keplerian rotation of galactic masses. Indeed, such residual, non-baryonic energy, when considered in Newton’s gravity equation, adds the term , which can partially explain, without recourse to dark matter, the rotations of some galaxies, such as M33, UGC12591, UGC2885, NGC3198, NGC253, DDO161, UDG44, the MW and the Coma cluster. Today, in the MW, that cosmological gravity force is in the order of 10²⁶ times smaller than the conventional gravity force. The model predicts an acceleration of the mass in the universe (q~ 0,986); the energy associated with curvature is the driving force behind the expansion of the universe, rather than the energy associated with the cosmological constant . An equation to determine expansion is obtained using the energy form of the Friedmann equation relative to Planck power and cosmic time or Planck force acting at the frontier of the universe moving at c. This constant Planck force, from unknown sources, acts everywhere to the expansion of the universe as a stretching effect on the volume. Finally, the model partly explains the value a₀ of the MOND theory. Indeed, a₀ is not a true constant, but depends on the cosmological constant at the time the great structures were formed (~1 [Gy]), as well as an adjustment of the typical mass and dimension of those great structures, such as galaxies. The constant a₀ is a different expression of the cosmological gravity force as expressed by the cosmological constant, Λ, acting through the energy-mass equivalent during the formation of the structures. It does not put in question the value of G.

The properties of iterative splitting methods with serial versions have been analyzed since recent years, see [1] and [3]. We extend the iterative splitting methods to a class of parallel versions, which allow to reduce the computational time and keep the benet of the higher accuracy with each iterative step. Parallel splitting methods are nowadays important to solve large problems, which can be splitted in subproblems and computed independently with the dierent processors. We present a novel parallel iterative splitting method, which is based on the multi-splitting methods, see [2], [10] and [15]. Such a exibilisation with multisplitting methods allow to decompose large iterative splitting methods and recover the benet of their underlying waveform-relaxation (WR) methods. We discuss the convergence results of the parallel iterative splitting methods, while we could reformulate such an error to a summation of the individual WR methods. We discuss the numerical convergence of the serial and parallel iterative splitting methods and present dierent numerical applications to validate the benet of the parallel versions.

Summary: Studying biological systems generally relies on computational modeling and simulation, e.g., for model-driven discovery and hypothesis testing. Progress in standardization efforts led to the development of interrelated file formats to exchange and reuse models in systems biology, such as SBML, the Simulation Experiment Description Markup Language (SED-ML), or the Open Modeling EXchange format (OMEX). Conducting simulation experiments based on these formats requires efficient and reusable implementations to make them accessible to the broader scientific community and to ensure the reproducibility of the results. The Systems Biology Simulation Core Library (SBSCL) provides interpreters and solvers for these standards as a versatile open-source API in Java™. The library simulates even complex bio-models and supports deterministic Ordinary Differential Equations (ODEs); Stochastic Differential Equations (SDEs); constraint-based analyses; recent SBML and SED-ML versions; exchange of results, and visualization of in silico experiments; open modeling exchange formats (COMBINE archives); hierarchically structured models; and compatibility with standard testing systems, including the Systems Biology Test Suite and published models from the BioModels and BiGG databases. Availability and implementation: SBSCL is freely available at https://draeger-lab.github.io/SBSCL/.

As long as the non-contact digital printing is not a common standard in the corrugated packaging industry, corrugated board crushing is a real issue that affects the load capacity of the boxes. Crushing mainly occurs during the converting of corrugated board (e.g. analog flexographic printing or laminating) and is a process that cannot be avoided. However, as show in this study, it can be controlled. In this work, extended laboratory tests were carried out on the crushing of double-walled corrugated board. The influence of fully controlled crushing (with a precision: ±10 μm) in the range from 10 to 70 % on different laboratory measurements was checked. Most of the typical mechanical tests were performed e.g. edge crush test, four-point bending test, shear stiffness test, torsional stiffness test, etc. on reference and crushed specimens. The residual thickness reduction of the crushed samples was also controlled. All empirical observations and performed measurements were the basis for building an analytical model of crushed corrugated board. The proven and verified model was then used to study the crushing effect of the selected corrugated board on the efficiency of simple packages with various dimensions.

In order to reveal the mechanism and influencing factors of high pressure jet cavitation of nozzle in submerged environment, this study focused on the evolutionary process of cavitation bubbles and combined finite volume method and mixed multi-phase flow model to analyze the cavitation, velocity distribution and experimental cavitation intensity of fishing net cleaning equipment. Results show that the cavitation inception, growth and collapse mainly occurred in the peripheral region of the flow field. Ring-shaped cavitation erosion zone appeared on the test sample target. A lot of small dense erosion pits were densely distributed in the ring-shaped erosion zone, erosion marks were observed in the center. The cavitation erosion intensity was greatly affected by the nozzle structure. As the diameter of nozzle increased from 0.6 mm to 1 mm, the maximum gas volume fraction increased by 8.53%. The nozzle outlet enlargement angle greatly increased the cavitation intensity. The nozzle with an outlet angle of 30° exhibited the optimal cavitation erosion performance . The cavitation volume fraction of the nozzle with short necking structure was slightly larger than that of the nozzle with long necking structure at the same level in the necking length rang of 3mm to 7mm. In terms of the influence of nozzle structure on the cavitation erosion effect, the nozzle diameter D ranked the first, followed by the outlet angle α, and the necking length L was at the last.

The Lower Rio Grande Valley, South Texas is considered one of the more vulnerable coastal areas to flooding related with abrupt climate changes. From 1980-2017, there were 7 flooding events, 57 severe storm events, and 8 tropical cyclone events with losses exceeding $1 billion in the State of Texas, according to NOAA NCEI. Coastal flooding is typically a result of storm surge and heavy rainfall produced by hurricanes and tropical storms. In this study, the two-dimensional hydrodynamic flow circulation model is developed to predict the Lower Rio Grande Valley coastal area inundation due to the hurricane storm surge, especially in the case of Hurricane Beulah, 1967. The tropical cyclone properties and tidal constituents were assigned to the updated watershed geographic information with the bottom bathymetric and roughness data. For model validation, the Hurricane Dolly 2008 storm surge due to Hurricane Beulah at the coast and the storm surge reaches up to approximately 40 kilometers west from the coast through a natural river channel. This model can be used for a reliable engineering tool for the coastal hazard emergency management and disaster mitigation.

Façade elements are a building component that satisfies multiple performance parameters. Among other things, “advanced façades” take advantage of hybrid solutions, such as assembling laminated materials. In addition to the enhanced mechanical properties that are typical of optimally composed hybrid structural components, these systems are energy-efficient, durable, and offer lighting comfort and optimal thermal performance, an example of which is the structural solution developed in collaboration with the University of Zagreb and the University of Ljubljana within the Croatian Science Foundation VETROLIGNUM project. The design concept involves the mechanical interaction of timber and glass load-bearing members without sealing or bonding the glass-to-timber surfaces. Following earlier research efforts devoted to the structural analysis and optimization of thus-assembled hybrid Cross-Laminated Timber (CLT)-glass façade elements, in this paper, special focus is given to a thermal and energy performance investigation under ordinary operational conditions. A simplified numerical model representative of a full-size building is first presented by taking advantage of continuous ambient records from a Live-Lab mock-up facility in Zagreb. Afterwards, a more detailed Finite Element (FE) numerical analysis is carried out at the component level to further explore the potential of CLT–glass façade elements. The collected numerical results show that CLT–glass composite panels can offer stable and promising thermal performance for façades similar to national and European standard requirements.

Recently, Indium Tin Oxide, a highly transparent, well conductive, and CMOS-compatible material, has been paying strong attention to the thermo-optic controlled silicon photonics industry because it allows a miniature of the gap between the core silicon and the heater, thus enabling reducing the electric power consumption and enhancing the switching speed. In this article, we propose an ultralow loss and small-size ITO microheater for the phase shift tuning. The designated microheater is manipulated in realizing a numerical co-design of compact and high bandwidth three-mode converter and three-mode selective router. Simulation results illustrate the 3-dB bandwidth for the three-mode converter and three-mode selective router as much as 100-nm and 40-nm during crosstalk under -25 dB, respectively. Besides, co-designed devices attain relatively large fabrication tolerances corresponding to width and height tolerances of ±50 nm and ±5 nm. In addition, the proposed devices consumed less than 90 mW total power consumption and took a fast switching time below 8 μs. Moreover, both two co-designs can be integrated into an estimated compact footprint of 8 μm2160 μm. Such excellent performances demonstrate the attractive potential of ITO as low-loss thermo-optic phase shifters and open an alternative way for enabling ultrafast and high-speed mode division multiplexing systems and very large-scale photonic integrated circuits.

The dynamic development of electromobility poses challenges to designers regarding not only the efficiency of energy transformation but also the battery life, which is influenced by the stability of its operating temperature. Designing cooling systems is connected not only with the optimization of energy management but also with other environmental parameters, such as noise emission. The paper presents the numerical optimization of an innovative radiator for use in electric buses in terms of energy consumption and noise emission. The results of the numerical studies were verified in laboratory and field conditions, showing a very good convergence of the model with the results of the experiments.

The Gleeble-1500D thermal simulation test machine was used to conduct the isothermal compression test on 21-4N at the strain rate ( ) of 0.01-10s-1, the deformation temperature (T) of 1273-1453K and the maximum deformation is 0.916. The data of the stress-strain ( - )were obtained. Based on the - data, the Johnson-Cook (J-C), modified J-C, Arrhenius and Back-Propagation Artificial Neural Network (BP-ANN) models were established. The accuracy of four models were verified, analyzed and compared. The results show that J-C model has a higher accuracy only under reference deformation conditions. When the deformation condition changes greatly, the accuracy of J-C model is significantly reduced. The coupling effect of T and of modified J-C model is considered, and the prediction accuracy is greatly improved The Arrhenius model introduces Zener-Hollomon (Z) to represent the coupling effect of T and , it has a fairly high prediction accuracy. And it can predict flow stress ( ) accurately at different conditions. The accuracy of BP-ANN model is the highest, but its learning rate is low, the learning and memory are unstable. It has no memory for the weights and thresholds of the completed training. So, there are certain limitations of it in use. Finally, a FEM of the isothermal compression experiment for four models were established, and the distribution of the equivalent stress field, equivalent strain field and temperature field with the deformation degree of 60% were obtained.

Biological diversity (BD) explored by the biological systematics is a complexly organized natural phenomenon and can be partitioned in several aspects defined with references to various causal factors structuring biota. These BD aspects are studied by particular research programs based on specific taxonomic theories (TT). They provide in total a framework for comprehending the structure of the biological systematics and its multi-aspect relations to other fields of biology. General principles of individualizing BD aspects and construing TT as quasi-axiomatics are briefly considered. It stressed that each TT is characterized by a specific combination of interrelated ontological and epistemological premises most adequate to the BD aspect a TT deals with. The following contemporary research programs in systematics are recognized and characterized in brief: phenetic, rational (with several subprograms), numerical, typological (with several subprograms), biosystematic, biomorphic, phylogenetic (with several subprograms), evo-devo. From a scientific pluralism perspective, all these research programs related to particular naturally defined BD aspects are of the same biological and scientific significance and no one of them can pretend to take a privileged position. They elaborate “locally” natural classifications that can be united by a kind of generalized faceted classification.

Tornadoes remain an active subject of observational and numerical research due to the damage and fatalities they cause worldwide as well as poor understanding of their behavior, such as what processes result in their genesis and what determines their longevity and morphology. A numerical model executed on a supercomputer run at high resolution can serve as a powerful tool for exploring the rapidly evolving tornado-scale features within a simulated storm, but saving large amounts data for meaningful analysis can result in unacceptably slow model performance, an unwieldy amount of saved data, and saved data spread across millions of files. In this paper, a system for efficiently saving and managing hundreds of terabytes of compressed model output is described in order to support a supercomputer simulation of a tornadic supercell thunderstorm. The challenges of managing a simulation spanning over a quarter trillion grid volumes across the Blue Waters supercomputer are also described. The simulated supercell produces a long-track EF5 tornado, and the near-tornado environment is described during tornadogenesis, where single upward-growing vortex undergoes several vortex mergers before transitioning into a multiple vortex tornado.

In the present paper four passive micromixers designs (G1, G2, G3 and G4) inspired on distillation columns trays were proposed. The devices performance was assessed by numerical simulations. The mixing performance was investigated for a Reynolds number range from 0.01 to 100 and channel height of 200 µm, 500 µm and 1000 µm for oil/ethanol flow. G1 and G4 designs provided high mixing index (> 0.975). The G1 device achieved superior mixing performance with a moderate pressure drop (about 0.5 MPa) due to the induced flow recirculation pattern for a relative high flow rate of 0.21 L/min, highlighting the potential use of such microdevice for scale-up and numbering-up of microdevices in modular chemical plant processing.

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.

In the steeply inclined and extremely thick coal seams (SIETCS) mined using the horizontal sublevel top-coal caving (HSTCC) method, the uncertainty of gas emission is a safety threat to the mining operations. In order to reduce the occurrence of accidents, the determination of gas emission is crucial. In this paper, we first proposed a prediction model for workers at the floor coalbed to calculate gas emission on site. We then put forward a finite element numerical simulation for researchers to predict gas emission from the floor coalbed. At last, we measured gas emitted from the floor coalbed of SIETCS in Wudong Coal Mine in a specific mining period and used the data to verify the applicability and accuracy of these two gas-emission prediction methods. The results showed that the gas emission from Wudong Coal Mine was 1.08 m³/min calculated based on the prediction model and 1.07 m³/min obtained using the user-defined integration method. Both methods have their own advantages, disadvantages and applicable objects, and are important in predicting gas emission from SIETCS mined using HSTCC method.

For the numerical simulation of components made of short fiber-reinforced composites the correct prediction of the deformation including the elastic and plastic behavior and its spatial distribution is essential. When using purely deterministic modeling approaches the information of the probabilistic microstructure is not included in the simulation process. One possible approach for the integration of stochastic information is the use of random fields. In this study numerical simulations of tensile test specimens are conducted utilizing a finite deformation elastic-ideal plastic material model. A selection of the material parameters covering the elastic and plastic domain are represented by cross-correlated second-order Gaussian random fields to incorporate the probabilistic nature of the material parameters. To validate the modeling approach tensile tests until failure are carried out experimentally, that confirm the assumption of spatially distributed material behavior in both the elastic and plastic domain. Since the correlation lengths of the random fields cannot be determined by pure analytic treatments, additionally numerical simulations are performed for different values of the correlation length. The numerical simulations endorse the influence of the correlation length on the overall behavior. For a correlation length of 5mm a good conformity with the experimental results is obtained. Therefore, it is concluded, that the presented modeling approach is suitable to predict the elastic and plastic deformation of a set of tensile test specimens made of short fiber-reinforced composite sufficiently.

The object of this work is to determine the correlation on the Nusselt number on the individual rows of a four-row tubular finned heat exchanger with continuous fins with a staggered tube arrangement using CFD modelling. Correlations for calculating Darcy-Weisbach friction factors on individual tube rows were also determined. Relationships for the Nusselt number and friction factor derived for the entire exchanger based on CFD modelling were compared with those available in the literature determined using experimental data. The maximum relative differences between the Nusselt number for a four-row exchanger determined experimentally and by CFD modelling are in the range from 22% for a Reynolds number based on a tube's outside diameter of 1,000 to 30% for a Reynolds number of 13,000. The maximum relative differences between the friction factor for a four-row exchanger determined experimentally and by CFD modelling are in the range of 50% for a Reynolds number based on a tube outer diameter of 1,000 to 10% for a Reynolds number of 13,000. The CFD modeling performed shows that in the range of Reynolds numbers based on hydraulic diameters from 150 to 1,400, the Nusselt number for the first row in a four-row finned heat exchanger is about 22% to 15% higher than the average Nusselt number for the entire exchanger. In the range of Reynolds number changes based on hydraulic diameter from 2,800 to 6,000, the Nusselt numbers on the first and second rows of tubes are close to each other. Correlations on Nusselt numbers and friction factors derived for individual tube rows can be used in the design of plate-fin and tube heat exchangers used in equipment such as air-source heat pumps, automotive radiators, air-conditioning systems and in air hot-liquid coolers. In particular, the correlations can be used to select the optimum number of tube rows in the exchanger.

Numerical modelling is becoming a major tool for supporting environmental studies at different scales, thanks to the capability of up-to-date codes in reproducing the natural behaviour in a quite reliable manner. In evaluating the habitat diversity of anthropized rivers, however, many issues are rising because of the intrinsic complexity of the processes involved. Using a reach of the Po River in Italy as a case study, the present works aims to provide an estimate of the changes of the Eco-Environmental Diversity as a response to different constant flow discharges. The goals are achieved by means of two solvers of the iRIC suite, applied in sequence to firstly simulate the fluvial hydrodynamics and subsequently provide an estimate of the habitat conditions. Despite the simplifications intrinsically present in the models and the ones introduced for practical purposes, the results pointed out that the reduction of the flow discharge recently observed can threat the overall biological status of the river. Because of the modelling uncertainties, on the other side, these preliminary outcomes show the need for more research, both in terms of data acquisition and numerical schematization, for adequately evaluate the effects of transient hydrology on the river ecosystems. Moreover, additional field surveys are necessary to calibrate and validate the used model for having sufficiently reliable estimates.

The assimilation of Soil Moisture and Ocean Salinity (SMOS) data into the ECMWF (European Centre for Medium Range Weather Forecasts) H-TESSEL (Hydrology revised - Tiled ECMWF Scheme for Surface Exchanges over Land) model is presented. SMOS soil moisture (SM) estimates have been produced specifically by training a neural network with SMOS brightness temperatures as input and H-TESSEL model SM simulations as reference. This can help the assimilation of SMOS information in several ways: (1) the neural network soil moisture (NNSM) data have a similar climatology to the model, (2) no global bias is present with respect to the model even if regional differences can exist. Experiments performing joint data assimilation (DA) of NNSM, 2 metre air temperature and relative humidity or NNSM-only DA are discussed. The resulting SM was evaluated against a large number of in situ measurements of SM obtaining similar results to those of the model with no assimilation, even if significant differences were found from site to site. In addition, atmospheric forecasts initialized with H-TESSEL runs (without DA) or with the analysed SM were compared to measure of the impact of the satellite information. Although, NNSM DA has an overall neutral impact in the forecast in the Tropics, a significant positive impact was found in other areas and periods, especially in regions with limited in situ information. The joint NNSM, T2m and RH2m DA improves the forecast for all the seasons in the Southern Hemisphere. The impact is mostly due to T2m and RH2m, but SMOS NN DA alone also improves the forecast in July- September. In the Northern Hemisphere, the joint NNSM, T2m and RH2m DA improves the forecast in April-September, while NNSM alone has a significant positive effect in July-September. Furthermore, forecasting skill maps show that SMOS NNSM improves the forecast in North America and in Northern Asia for up to 72 hours lead time.

The present work is focused on the numerical solution of the complete energy equation used in fluid film lubrication. The work was motivated by the fact the complete energy equation has no analytic solution that could be used for validations. Its accuracy and computation time are related to the employed numerical method and to the grid resolution. The natural discretization method (NDM) applied on different grids is systematically compared with the spectral method (the Lobatto Point Colocation Method or LPCM) with different polynomial degrees. A one dimensional inclined slider is used for the numerical tests and the energy equation is artificially decoupled from Reynolds. This approach enables to focus all the attention on the numerical solution of the energy equation. The results show that the LPCM is one or two orders of magnitudes more efficient than the NDM in terms of computation time. The energy equation is then coupled with Reynolds equation in a thermo-hydrodynamic analysis of the same 1D slider; the numerical results confirm again the efficiency of the LPCM. A thermo-hydrodynamic analysis of a two-lobe journal bearing is then presented as a practical application.

A conventional 2D numerical model is improved by incorporating three submodels to consider different effects of secondary flow and a module for cohesive sediment transport. The model is applied to a meandering reach of Yangtze River to investigate secondary flow effects on cohesive sediment deposition, and a preferable submodel is selected based on the flow simulation results. Sediment simulation results indicate that the improved model predictions are in better agreement with the measurements in planar distribution of deposition as the increased sediment deposits caused by secondary current on the convex bank have been well predicted. Secondary flow effects on predicted amount of deposition become more obvious during the period when the sediment load is low and velocity redistribution induced by the bed topography is evident. Such effects vary with the settling velocity and critical shear stress for deposition of cohesive sediment. The bed topography effects can be reflected by the secondary flow submodels.

{\bf Abstract.} \\ {\bf Background:} Bangladesh hosts more than 800,000 Rohingya refugees from Myanmar. The low health immunity, lifestyle, access to good healthcare services, social-security cause this population to be at risk of far more direct effects of COVID-19 than the host population. Therefore, evidence-based forecasting of the COVID-19 burden is vital in this regard. In this study, we aimed to forecast the COVID-19 burden among the Rohingya refugees of Bangladesh to keep up with the disease outbreak's pace, health needs, and disaster preparedness. \noindent{\bf Methodology and Findings:} To estimate the possible consequences of COVID-19 in the Rohingya camps of Bangladesh, we used a modified Susceptible-Exposed-Infectious Recovered (SEIR) transmission model. All the values of different parameters used in this model were from the Bangladesh Government’s database and the relevant emerging literature. We addressed two different scenarios, i.e., the best-fitting model and good fitting model with unique consequences of COVID-19. Our best fitting model suggests that there will be good control over the transmission of the COVID-19 disease. At the end of December 2020, there will be only 169 confirmed COVID-19 cases in the Rohingya refugee camps. The average basic reproduction number ($ \mathcal{R}_{0} $) has been estimated to be 0.7563. \noindent{\bf Conclusion:} Our analysis suggests that, due to the extensive precautions from the Bangladesh government as well as other humanitarian organizations, the coronavirus disease will be under control if the maintenance continues like this. Although detailed and pragmatic preparedness should be adopted for the worst scenario.

Facade elements are known to represent a building component with multiple performance parameters to satisfy. Among others, “advanced facades” take advantage of hybrid solutions, like the assemblage of laminated materials. In addition to enhanced mechanical properties that are typical of optimally composed hybrid structural components, these systems are energy-efficient, durable, and offer lightening comfort and optimal thermal performance. This is the case of the structural solution developed in joint research efforts of the University of Zagreb and the University of Ljubljana, within the Croatian Science Foundation VETROLIGNUM project. The design concept involves the mechanical interaction of timber and glass load-bearing members, without sealing or bonded glass-to-timber surfaces. Laminated glass infilled timber frames are recognized as a new generation of structural members with relevant load-carrying capacity (and especially the enhancement of earthquake resistance of framed systems), but also energy-efficient and cost-effective solutions.

This study identifies new Alkali Activated Composite Mud Mortar to alternate Cement Mortar used for laying masonry works. This study handles both experimental and numerical analysis of strength and durability for mortar, especially bond strength along mortar and brick joints have studied detail and compared with ANSYS. The bond strength has estimated through the triplet brick prism of shear behavior. The strength and deformation characteristics of masonry constituents obtained from these tests depend on the actual composite conduct of masonry. There is always a good reception for environmentally friendly, low energy, new construction materials, and used in preparing mortar/brick unit combinations all around the world. Also, this study fulfills those requirements. Especially, masonry laid with the mud mortar meets all. Fly Ash, Ground Granulated Blast furnace Slag, and Quarry Dust have utilized the new composites' preparation with 8Molority, 10Molority developed excellent strength. The composite specimens are having compression and bond strength, ranging from 4.84-5.14 and 0.12-0.18N/mm2. The mixed ratio of materials has taken (0.5 parts of Fly Ash, 0.5 parts of Steel Slag, 1.75 parts of Soil, and 0.25 parts of Quarry Dust). The samples tested resistance against acidic, sulfate, and alkalinity indicated PH values less than 7 and 9 simultaneously. Compression strength varied from 10-20% and gained more and fit composite used for mortar in masonry works.

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.

We investigated preschool-aged children’s understanding of early fractional tasks and how that performance correlates with fine motor skills and use of gestures while counting. Participants were 33 preschoolers aged 4 to 5 in two Southeastern public elementary schools. Children were tested individually in an interview-like setting. Mathematics tasks were presented in a paper and pencil format and the Grooved Pegboard test assessed fine motor skills. Finally, utilization of gestures was evaluated by taking a behavioral rating of the child’s hand morphology, accuracy of gestures, and synchrony of gestures and spoken word while performing a counting task. Results indicate performance on fractional reasoning tasks significantly predicts both fine motor ability and accuracy of gestures.

Flood routing in flood forecasting issue, calculation the height of flood bands, determining the river boundaries, and estimation of protective facilities for flood –exposed building is applicable. In many cases, due to the lack of measuring stations, the status of the upstream flood generating hydrograph is not known. The purpose of this study is to present an integrated method comprising of an optimization model and a hydrodynamic numerical model for flood modeling to determine the upstream hydrograph using the provided hydrograph at the downstream measuring station of a river. The routing procedure consists of three steps: (1) generating a hypothetical upstream hydrograph using genetic algorithm method; (2) hydrodynamic modeling using a numerical simulation model for flood routing according to the hypothetical hydrograph which is generated in the first step; (3) compare the calculated and observed hydrograph in downstream by using a fitness function. This recommended procedure was named as Reverse Flood Routing Method (RFRM) and was then applied to Karun River, the largest river in Iran. Comparing the generated upstream hydrograph by the RFRM model with the corresponding measured hydrograph at Ahvaz hydrometric station, as an ungauged river location, shows the high accuracy of the recommended model in this study.

Polydimethylsiloxane (PDMS) is one of the most popular elastomers and has been used in different fields, especially in biomechanics research. Among the many interesting features of this material, its hyperelastic behaviour stands out, which allows the use of PDMS in various applications, like the ones that mimic soft tissues. However, the hyperelastic behaviour is not linear and needs detailed analysis, especially the characterization of shear strain. In this work, two approaches, numerical and experimental, were proposed to characterize the effect of shear strain on PDMS. The experimental method was implemented as a simple shear testing associated with 3D digital image correlation and was made using two specimens with two thicknesses of PDMS (2 and 4 mm). A finite element software was used to implement the numerical simulations, in which four different simulations using the Mooney-Rivlin, Yeoh, Gent, and Polynomial hyperelastic constitutive models were performed. These approaches showed that the maximum value of shear strain occurred in the central region of the PDMS, and higher values emerged for the 2 mm PDMS thickness. Qualitatively, in the central area of the specimen, the numerical and experimental results have similar behaviours and the values of shear strain are close. For higher values of displacement and thicknesses, the numerical simulation results move further away from experimental values.

This paper reviews how dynamically adaptive wavelet methods can be designed to simulate atmosphere and ocean dynamics in both flat and spherical geometries. We highlight the special features that these models must have in order to be valid for climate modelling applications. These include exact mass conservation and various mimetic properties that ensure the solutions remain physically realistic, even in the under-resolved conditions typical of climate models. Particular attention is paid to the implementation of complex topography in adaptive models. Using \textsc{wavetrisk} as an example, we explain in detail how to build a semi-realistic global atmosphere or ocean model of interest to the geophysical community. We end with a discussion of the challenges that remain to developing a realistic dynamically adaptive atmosphere or ocean climate models. These include scale-aware subgrid scale parameterizations of physical processes, such as clouds. Although we focus on adaptive wavelet methods, many of the topics we discuss are relevant for adaptive mesh refinement (AMR).

Transient flows result in unbalanced forces and high pressure in pipelines. Under these conditions, the combined effects of flow-induced forces along with sudden pipe displacements can create cracks in the pipeline, especially at the junctions. This situation consequently results in water leakage and reduced operational efficiency of the pipeline. In this study, displacements and stresses in a buried pressurized water transmission pipe installed by pipe jacking method are investigated using numerical modeling and considering interactions between fluid, pipe, and soil. The analyses were performed consecutively under no-flow, steady flow, and transient flow conditions, in order to investigate the effects of flow conditions on displacements and stresses in the system. Analyses of the results show that displacements and stresses in the jointed concrete pipes are significant under transient flow conditions. Moreover, because of pressure transient effects, maximum tensile stresses exceed the tensile strength of concrete at the junctions, leading to cracks and consequent water leakage.

Predictive simulation of concurrent debris flow using only pre-disaster information has proven to be difficult as a result of problems in predicting the location of debris-flow initiation (i.e., slope failure). However, because catchment topography has concave characteristics, with all channels in a catchment joining each other as they flow downstream, it is possible to predict damage to downstream area using relatively inaccurate initiation points. Based on this, this paper presents methodologies employing debris-flow initiation points generated randomly using statistical slope failure prediction. A many-case simulation across numerous initiation points was performed to quantify the effect of slope-failure location in terms of deviations in the predicted water level and terrain deformation. It was found that the relative standard deviation diminished as the points approached the downstream area, indicating a location-based predictability effect.

This article is devoted to Feller's diffusion equation which arises naturally in probabilities and physics (e.g. wave turbulence theory). If discretized naively, this equation may represent serious numerical difficulties since the diffusion coefficient is practically unbounded and most of its solutions are weakly divergent at the origin. In order to overcome these difficulties we reformulate this equation using some ideas from the Lagrangian fluid mechanics. This allows us to obtain a numerical scheme with a rather generous stability condition. Finally, the algorithm admits an elegant implementation and the corresponding Matlab code is provided with this article under an open source license.

The need for quick and easy deflection calculations of various prefabricated slabs causes that simplified procedures and numerical tools are used more and more often. Modelling of full 3D finite element (FE) geometry of such plates is not only uneconomical but often requires the use of complex software and advanced numerical knowledge. Therefore, numerical homogenization is an excellent tool, which can be easily employed to simplify a model, especially when accurate modelling is not necessary. Homogenization allows for simplifying a computational model and replacing a complicated composite structure with a homogeneous plate. Here, a numerical homogenization method based on strain energy equivalence is derived. Using the method proposed, the structure of the prefabricated concrete slabs reinforced with steel spatial trusses is homogenized to a single plate element with an effective stiffness. There is a complete equivalence between the full 3D FE model built with solid elements combined with truss structural elements and the simplified homogenized plate FE model. The method allows for the correct homogenization of any complex composite structures made of both solid and structural elements, without the need to perform advanced numerical analyses. The only requirement is a correctly formulated stiffness matrix of a representative volume element (RVE) and appropriate formulation of the transformation between kinematic constrains on RVE boundary and generalized strains.

Computational Fluid Dynamics (CFD) simulations, based on Reynolds Averaged Navier Stokes (RANS) models, are a useful tool for a wide range of coastal and offshore applications, providing a high fidelity representation of the underlying hydrodynamic processes. Generating input waves in the CFD simulation is performed by a numerical wavemaker (NWM), with a variety of different NWM methods existing for this task. While NWMs, based on impulse source methods, have been widely applied for wave generation in depth averaged, shallow water models, they have not seen the same level of adoption in the more general RANS based CFD simulations, due to difficulties in relating the required impulse source function to the resulting free surface elevation for non-shallow water cases. This paper presents an implementation of an impulse source wavemaker, which is able to self-calibrate the impulse source function to produce a desired wave series in deep or shallow water at a specific point in time and space. Example applications are presented, for a numerical wave tank (NWT), based on the opensource CFD software OpenFOAM, for wave packets in deep and shallow water, highlighting the correct calibration of phase and amplitude. Also, the suitability for cases requiring very low reflection from NWT boundaries is demonstrated. Possible issues in the use of the method are discussed and guidance for good application is given.

This code provides computational facilities to simulate current versus time during the charging of Li-ion cells at desire constant voltage by considering multiscale physical phenomena. This code only considers a powder of active materials (at microscale or nanoscale) and a small part of electrolyte around it as a half cell. Then it is extended to a complete cell by applying correct boundary conditions. This code is very useful by modifying code parameters to understand the effect of the complex shape of active materials powder (surface area and powder size), kind of electrolyte, and the applied voltages on the charging response of Li-ion cell. As a summary, a microscale approach to the design of Li-ion cells has been provided via this code.

Time evolution operators of a strongly ionizing medium are calculated by a time-dependent unitary transformation (TDUT) method. The TDUT method has been employed in quantum mechanical system composed of discrete states. This method is especially helpful for solving molecular rotational dynamics in quasi-adiabatic regimes because the strict unitary nature of the propagation operator allows us to set the temporal step size large; a tight limitation on the temporal step size ($\delta t <<1$) can be circumvented by the strict unitary nature. On the other hand, in a strongly ionizing system where the Hamiltonian is not Hermitian, the same approach cannot be directly applied because it is demanding to define a set of field-dressed eigenstates. In this study, the TDUT method was applied to the ionizing regime using the Kramers-Henneberger frame, in which the strong-field-dressed discrete eigenstates are given by the field-free discrete eigenstates in a moving frame. Although the present work verifies the method for a one-dimensional atom as a prototype, the method can be applied to three-dimensional atoms, and molecules exposed to strong laser fields.

Advancement in physics often results from analyzing numerical data and then creating a theoretical model that can explain and predict those data. In the field of lens design, the reverse is true: longstanding theoretical understanding is being overtaken by more powerful numerical methods.

To improve the effect of multi-point stretch forming of sheet metal, it is proposed in this paper to replace fixed ball head with swinging ball head. According to the multi-point dies with different arrangements, this research establishes the finite element models of the following stretch forming, i.e. fixed ball heads with conventional arrangement, swinging ball heads with conventional arrangement, swinging ball heads with declining staggered arrangement, and swinging ball heads with parallel staggered arrangement and then numerical simulation is performed. The simulation results show that by replacing fixed ball head with swinging ball head, the surface indentation of the formed part was effectively suppressed, the stress and tension strain distribution of the formed part was improved and the forming quality was improved; that the thickness of the elastic pad was reduced, the springback was reduced and the forming accuracy was improved; and that when the ball head was applied to multi-point die with staggered arrangement, better forming result was achieved, where the best forming result was achieved in combining the swinging ball heads with the multi-point die with parallel staggered arrangement. The forming experiments were carried out, and the experimental results were consistent with the trend of numerical simulation results, which verified the correctness of the numerical simulation.

The corrugated board packaging industry is increasingly using advanced numerical tools to design and estimate the load capacity of its products. That is why numerical analyzes are becoming a common standard in this branch of manufacturing. Such trend causes either the use of advanced computational models that take into account the full 3D geometry of the flat and wavy layers of corrugated board, or the use of homogenization techniques to simplify the numerical model. The article presents theoretical considerations that extend the numerical homogenization technique already presented in our previous work. The proposed here homogenization procedure also takes into account the creasing and / or perforation of corrugated board, i.e. processes that undoubtedly weaken the stiffness and strength of the corrugated board locally. However, it is not always easy to estimate how exactly these processes affect the bending or torsional stiffness. What is known for sure is that the degradation of stiffness depends, among other things, on the type of cut, its shape, the depth of creasing, as well as their position or direction in relation to the corrugation direction. The method proposed here can be successfully applied to model smeared degradation in a finite element or to define degraded interface stiffnesses on a crease line or a perforation line.

The previously developed models for fuel droplet heating and evaporation processes, mainly the Discrete Multi Component Model (DMCM), and Multi-Dimensional Quasi-Discrete Model (MDQDM) are investigated for the aerodynamic combustion simulation. The models have been recently improved, and generalised for a broad range of bio-fossil fuel blends so that the application areas are broadened with increased accuracy. The main distinctive features of these models are that they consider the impacts of species thermal conductivities and diffusivities within the droplets to account for the temperature gradient, transient diffusion of species and recirculation. A formulation of fuel surrogates is made, using the recently introduced model, referred to as ‘’Complex Fuel Surrogate Model (CFSM)’’ and analysing their heating, evaporation, and combustion characteristics. The CFSM is aimed to reduce the full composition of fuel to a much smaller number of components based on their mass fractions, and to formulate fuel surrogates. Such approach has provided a proof of concept with the implementation of the developed model into a commercial CFD code ANSYS-Fluent. A case study is made for the CFD modelling of gas-turbine engine using kerosene fuel surrogate. The surrogate is proposed using the CFSM. The model is implemented into ANSYS-Fluent via a user-defined function to provide the first full simulation of the combustion process. Detailed chemical mechanism is also implemented into ANSYS Chemkin for the combustion study.

Several problems of radiative transfer are yet unsolved because of the difficulties of the calculations involved in them, especially if the intervening shapes are geometrically complex. The main goal of our investigation in this domain is to convert the formulas that were previously derived, into a graphical interface based on the projected solid-angle principle. Such procedure is now feasible by virtue of several widely diffused programs for Algorithms Aided Design (AAD). Accuracy and reliability of the process is controlled by means of the analytical software DianaX developed at an earlier stage by the authors. With this new approach the often cumbersome procedure of lighting and thermal exchange calculations can be simplified and made available for the neophyte, with the undeniable advantage of reduced computer time.

Structural vibration induced by dynamic load or natural vibration is a nonnegligible factor in failure analysis. Based on vibrating boundary condition, impact resistance of shape memory alloy reinforced composites is investigated. In this investigation, modified Hashin’s failure criterion, Brinson’s model and visco-hyperelastic model are implemented into the numerical model to charactering the mechanical behavior of glass fiber/epoxy resin laminates, SMAs and interphase, respectively. First, fixed boundary condition is maintained in simulation to verify the accuracy of material parameters and procedures by comparing with experimental data. Then, a series of vibrating boundaries with different frequencies and amplitudes are applied during the simulation process to reveals the effect on impact resistances. The statistics of absorbed energy and contact force indicate that impact resistance of the composite under high frequency and large amplitude is lower than that under low frequency and small amplitude, and summarized by a mathematical expression.

The study of the influence of the nanoparticle volume fraction and aspect ratio of microchannels on the fluid flow and heat transfer characteristics of nanofluids in microchannels is important in the optimal design of heat dissipation systems with high heat flux. In this work, the computational fluid dynamics method was adopted to simulate the flow and heat transfer characteristics of two types of water–Al2O3 nanofluids with two different volume fractions and five types of microchannel heat sinks with different aspect ratios. Results showed that increasing the nanoparticle volume fraction reduced the average temperature of the liquid–solid heat transfer surface and thereby improved the heat transfer capacity of the nanofluids. Meanwhile, the increase of the nanoparticle volume fraction led to a considerable increase in the pumping power of the system. Changing the aspect ratio of the microchannel effectively improved the heat transfer capacity of the heat sink. Moreover, increasing the aspect ratio effectively reduced the average temperature of the heating surface of the heat sink without significantly increasing the flow resistance loss. When the aspect ratio exceeded 30, the heat transfer coefficient did not increase with the increase of the aspect ratio. The results of this work may offer guiding significance for the optimal design of high heat flux microchannel heat sinks.

Air filtration mechanisms in the composite filter media used in practical applications are important and challenging to understand because the component fibers could have various size scales and morphologies. In this work, a three-dimensional digital model of nanofiber-based filter media was reconstructed for the first time based on the X-ray tomography data for the cellulose substrate and the Focused Ion Beam-Scanning Electron Microscope (FIB-SEM) image analysis for the several microns thick (3.82-7.90 μm) electrospun polyvinylidene fluoride (PVDF) nanofiber membrane. Besides the high-resolution model where the details of the fibrous structures were fully resolved, another low-resolution model with approximated unresolved structures was also established. Filtration simulations utilizing these models were conducted considering the drag force, Brownian diffusion and aerodynamic slip. The simulated filtration efficiencies agreed well with the experiments for particles of 70-400 nm, including the most penetrating particle size (MPPS, 100-200 nm). Moreover, the structure-resolved models had higher accuracy but higher computational costs, while the unresolved simulations saved much running time but over-predicted the filtration efficiency, especially for smaller particles (<100 nm). Our study presents a comprehensive strategy for investigating the composite filter media with multiscale complex structures using a combination of advanced characterization technologies and modular simulation models.