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.

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.

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.

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.

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.

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.

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 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.

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.

This article develops duality principles and numerical results for a large class of non-convex variational models. The main results are based on fundamental tools of convex analysis, duality theory and calculus of variations. More specifically the approach is established for a class of non-convex functionals similar as those found in some models in phase transition. Finally, in the last section we present a concerning numerical example and the respective software.

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.

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.

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.

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.

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 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.

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.

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.

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.

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.

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.

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/.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

In the present trend of constructing taller and longer structures, the application of lightweight aggregate concrete is becoming an increasing important advanced solution in the modern construction industry. In engineering practice, the analysis of lightweight concrete elements is performed using the same algorithms used for normal concrete elements. As an alternative to traditional engineering methods, nonlinear numerical algorithms based on constitutive material models may be used. The paper presents a comparative analysis of curvature calculations for flexural lightweight concrete elements, incorporating analytical code methods EN 1992-1 and ACI 318-14, as well as a numerical analysis using the constitutive model of cracked tensile lightweight concrete recently proposed by the authors. To evaluate the adequacy of the theoretical predictions, experimental data of 51 lightweight concrete beams tested during five different programmes were collected. A comparison of theoretical and experimental results showed that the most accurate predictions are obtained using numerical analysis and the constitutive model proposed by the authors. In the future, the latter algorithm can be used as a reliable tool for improving the design standard methods or numerical modelling of lightweight concrete elements subjected to short-term loading.

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.

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.

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.

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.

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.

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.

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.

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.

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.

The application of instrumented indentation to assess material properties like Young’s modulus and micro-hardness has become a standard method. In recent developments, indentation experiments and simulations have been combined to inverse methods, from which further material parameters as yield strength, work hardening rate, and tensile strength can be determined. In this work, an inverse method is introduced by which material parameters for cyclic plasticity, i.e. kinematic hardening parameters, can be determined. To accomplish this, cyclic Vickers indentation experiments are combined with finite element simulations of the indentation with unknown material properties, which are then determined by inverse analysis. To validate the proposed method, these parameters are subsequently applied to predict the uniaxial stress-strain response of a material with success. The method has been validated successfully for a quenched and tempered martensitic steel and for technically pure copper, where an excellent agreement between measured and predicted cyclic stress-strain-curves has been achieved. Hence, the proposed inverse method based on cyclic nanoindentation, as a quasi-non-destructive method, could complement or even substitute the resource-intensive conventional fatigue testing in the future for some applications.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

{\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.

Non steady applied magnetic field impact on a liquid metals has good prospects for industry. For a better understanding of heat and mass transfer processes under these circumstances, numerical simulations are needed. A combination of finite elements and volumes methods was used to calculate the flow and solidification of liquid metal under electromagnetic influence. Validation of numerical results was carried out by means of measuring with ultrasound Doppler velocimetry technique, as well as with neutron radiography snapshots of the position and shape of the solid/liquid interface. As a result of the first part of the work, a numerical model of electromagnetic stirring and solidification was developed and validated. This model could be an effective tool for analyzing the electromagnetic stirring during the solidification process. In the second part, the dependences of the velocity pulsation amplitude and the melt velocity maximum value on the magnetic field pulsation frequency are obtained. It was found numerically the ability of the pulsating force action to develop higher values of the liquid metal velocity at a frequency close to the MHD resonance. Obtained characteristics give a more detailed description of the electrically conductive liquid behaviour under action of pulsating traveling magnetic field.

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.

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.

One of the most rapidly advancing areas of deep learning research aims at creating models that learn to disentangle the latent factors of variation from a data distribution. However, modeling joint probability mass functions is usually prohibitive, which motivates the use of conditional models assuming that some information is given as input. In the domain of numerical cognition, deep learning architectures have successfully demonstrated that approximate numerosity representations can emerge in multi-layer networks that build latent representations of a set of images with a varying number of items. However, existing models have focused on tasks requiring to conditionally estimate numerosity information from a given image. Here we focus on a set of much more challenging tasks, which require to conditionally generate synthetic images containing a given number of items. We show that attention-based architectures operating at the pixel level can learn to produce well-formed images approximately containing a specific number of items, even when the target numerosity was not present in the training distribution.

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.

Present study offers mathematical calculations of the roll-coating procedure lubricated with upper Convected Maxwell Fluid. An incompressible isothermal viscoelastic fluid is considered with roll and the porous web having uniform velocities. By employing Lubrication Approximation Theory, the desired equations of motion for the fluid concerned over porous web are modelled and analyzed. The suction rate on the web and injection rate at the roll surface are anticipated proportionate. Results for velocity profile and pressure gradient are obtained analytically. Fluid parameters of industrial significance i.e. detachment point, pressure, sheet /roll separating force, power contribution and coating thickness are also calculated numerically. Substantial and monotonic increase is witnessed in these quantities with the increase of flow parameters.

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.

The structure of turbulent flow over non-flat surfaces is a topic of major interest in practical applications in both engineering and geophysical settings. A lot of work has been done in the fully rough regime at high Reynolds numbers where the effect on the outer layer turbulence structure and the resulting friction drag is well documented. It turns out that surface topology plays a significant role on the flow drag especially in the transitional roughness regime and therefore, hard to characterize. Survey of literature shows that roughness function depends on the interaction of roughness height, flow Reynolds number and topology shape. In addition, if the surface topology contains large enough scales then it can impact the outer layer dynamics and in turn modulate the total frictional force. Therefore, it is important to understand the mechanisms underlying drag increase from systematically varied surface undulations in order to better interpret quantifications based on mean statistics such as roughness function. In this study, we explore the mechanisms that modulate the turbulence structure over a two-dimensional (2D) sinusoidal wavy surface with a fixed amplitude, but varying slope. To accomplish this, we model the turbulent flow between two infinitely wide 2D wavy plates at a friction Reynolds number, $Re_{\tau}=180$. We pursue two different but related flavors of analysis. The first one adopts a roughness characterization flavor of such wavy surfaces. The second one focuses on understanding the non-equilibrium near surface turbulence structure and their impact on roughness characterization. Analysis of the different statistical quantifications show strong dependence on wave slope for the roughness function indicating drag increase due to enhanced turbulent stresses resulting from increased production of vertical velocity variance from the surface undulations.

The energy transfer process of laser selective melting is very complex. To study the effect of la-ser selective melting on the microstructure and properties of 18Ni-300 martensitic steel, ABAQUS was used to simulate the temperature of laser cladding 18Ni-300 martensitic steel at different time points and different laser power. The results show that the cross-section shape of the molten pool changes from round to oval With the increase of laser power, the higher the peak value of temperature time curve, the greater the temperature gradient; and the laser clad-ding experiment of 18Ni-300 martensitic steel was carried out, and the microstructure and me-chanical properties of the samples under different laser power were analyzed. The results show that with the increase of laser power, the grain size of the cladding layer becomes smaller and the microstructure becomes more compact; the hardness of the side surface of the sample is higher than that of the upper surface, and the tensile strength and elongation show a trend of first increasing and then decreasing.

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.

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.

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.

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.

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.

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.

The adoption of Pump as Turbines (PATs) both in small scale hydroelectric plants and water supply systems has brought different advantages, the most recognized being cost-effectiveness as compared to other hydroturbines. However, due to their lack of flow control ability, their intolerance to off-design operations constitutes one tough shortfall. Moreover, this technology’s newness leads to PAT flow dynamics still being ununderstood even to date. Therefore, this study intends to numerically investigate the mixed flow PAT’s pump mode flow dynamics for five operating conditions expanding from optimum (1QBEP) to deep part-load (0.41 QBEP) conditions. Moreover the effect of runner blade angle on the same has been investigated where three angles namely -2°, 0°, and 2° have been considered. PAT flow stability was found to deteriorate as the flow decreased, where associated pressure pulsation level worsened at different flow zones. In addition, the blade angle increase led correspondingly increasing flow unsteadiness and pressure pulsation levels, where the pulsation frequencies from rotor-stator interactions were dominant for most of flow zones. This study’s findings are of a crucial importance to both scientific and engineering communities as they contribute to thorough understanding of PAT flow dynamics.

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.

In this work the existence and nonexistence of stationary radial solutions to the elliptic partial differential equation arising in the molecular beam epitaxy are studied. Since we are interested in radial solutions we arrive at the following fourth-order differential equation

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.

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.

The amplified Arctic warming is one of several factors influencing atmospheric dynamics. In this work, we consider a series of numerical experiments to identify the direct role of the Arctic sea ice reduction process in forming climatic trends in the northern hemisphere. Aimed at this, we used two more or less independent mechanisms of ice reduction. The first is traditionally associated with increasing the concentration of carbon dioxide in the atmosphere from the historic level of 360 ppm to 450 ppm and 600 ppm. This growth increases air temperature and decreases the ice volume. The second mechanism is associated with a reduction in the reflectivity of ice and snow. We assume that comparing the results of these two experiments allows us to judge the direct role of ice reduction. The most prominent consequences of ice reduction, as a result, were the weakening of temperature gradient at the tropopause level in mid-latitudes, the slower zonal wind at 50-60∘N, intensification of wave activity in Europe, Western America, and Chukotka, and its weakening in the south of Siberia and Kazakhstan. We also consider how climate change may alter regimes such as blocking and stationary Rossby waves. The study used the INM-CM48 climate system model .

This work is concerned with the existence and uniqueness of boundary value problems defined on semi-infinite intervals. These kinds of problems seldom admit exactly known solutions and, therefore, the theoretical information on their well-posedness is essential before attempting to derive an approximate solution by analytical or numerical means. Our utmost contribution in this context is the definition of a numerical test for investigating the existence and uniqueness of solutions of boundary problems defined on semi-infinite intervals. The main result is given by a theorem relating the existence and uniqueness question to the number of real zeros of a function implicitly defined within the formulation of the iterative transformation method. As a consequence, we can investigate the existence and uniqueness of solutions by studying the behaviour of that function. Within such a context the numerical test is illustrated by two examples where we find meaningful numerical results.

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.

Groundwater models serve as support tools to among others: assess water resources, evaluate management strategies, design remediation systems and optimize monitoring networks. Thus, the assimilation of information from observations into models is crucial to improve forecasts and reduce uncertainty of their results. As more information is collected routinely due to the use of automatic sensors, data loggers and real time transmission systems; groundwater modelers are becoming increasingly aware of the importance of using sophisticated tools to perform model calibration in combination with sensitivity and uncertainty analysis. Despite their usefulness, available approaches to perform this kind of analyses still present some challenges such as non-unique solution for the parameter estimation problem, high computational burden and a need of a deep understanding of the theoretical basis for the correct interpretation and use of their results, in particular the ones related to uncertainty analysis. We present a brief derivation of the main equations that serve as basis for this kind of analysis. We demonstrate how to use them to estimate parameters, assess the sensitivity and quantify the uncertainty of the model results using an example inspired by a real world setting. We analyze some of the main pitfalls that can occur when performing such kind of analyses and comment on practical approaches to overcome them. We also demonstrate that including groundwater flow estimations, although helpful in constraining the solution of the inverse problem as shown previously, may be difficult to apply in practice and, in some cases, may not provide enough information to significantly constrain the set of potential solutions.

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.

Keywords: underground engineering; numerical simulation; excavation length effect, major principal stress; displacement; damage initiation; CPU time

This study develops a predictor-corrector algorithm for the numerical simulation of IVPs involving singular generalized fractional derivatives with Mittag-Leffler kernels. The proposed algorithm converts the considered IVP into a Volterra-type integral equation and then uses Trapezoidal rule to obtain approximate solutions. Numerical approximate solutions of some singular generalized fractional derivative with Mittag-Leffler kernels models have been presented to demonstrate the efficiency and accuracy of the proposed algorithm. The algorithm describes the influence of the fractional derivative parameters on the dynamics of the studied models. The suggested method is expected to be effectively employed in the field of simulating generalized fractional derivative models

Corrugated cardboard is an ecological material, mainly because, in addition to virgin cellulose fibers also the fibers recovered during recycling process are used in its production. However, the use of recycled fibers causes slight deterioration of the mechanical properties of the corrugated board. In addition, converting processes such as printing, die-cutting, lamination, etc. cause micro-damage in the corrugated cardboard layers. In this work, the focus is precisely on the crushing of corrugated cardboard. A series of laboratory experiments were conducted, in which the different types of single-walled corrugated cardboards were pressed in a fully controlled manner to check the impact of the crush on the basic material parameters. The amount of crushing (with a precision of 10 micrometers) was controlled by a precise FEMat device, for crushing the corrugated board in the range from 10 to 70 % of its original thickness. In this study, the influence of crushing on bending, twisting and shear stiffness as well as a residual thickness and edge crush resistance of corrugated board was investigated. Then, a procedure based on a numerical homogenization, taking into account a partial delamination in the corrugated layers to determine the degraded material stiffness was proposed. Finally, using the empirical-numerical method, a simplified calculation model of corrugated cardboard was derived, which satisfactorily reflects the experimental results.

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.

Wave is a common environmental load that often causes serious damages to offshore structures. In addition, the stability for the submarine artificial slope is also affected by the wave loading. Although the landslide of submarine slopes induced by the waves has received wide attention, the research on the influence of solitary-wave is rare. In this study, a 2-D integrated numerical model is developed to investigate the stability of the foundation trench under the solitary wave loading. The Reynolds-Averaged-Stokes (RANS) equations are used to simulate the propagation of a solitary wave, while the current is realized by setting boundary inlet/outlet velocity. The pore pressure induced by the solitary wave is calculated by Darcy’s law and the seabed is characterized by Mohr-Coulomb constitutive model. Firstly, the wave model is validated through the comparison between analytical solution and experimental data. The initial consolidation state of slope under hydrostatic pressure is achieved as the initial state. Then, the factor of stability (FOS) for the slope corresponding to different distance between wave crest and slope top is calculated with the strength reduction method. The minimum of FOS is defined as the stability index for the slope with specific slope ratio during the process of dynamic wave loading. The parametric study is conducted to examine the effects of soil strength parameters, slope ratio and current direction. At last, the influence of upper slope ratio in a two-stage slope is also discussed.

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.

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.

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).

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.

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.

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.

Our aim is to implement and test some less known numerical methods that might be used to solve an advection equation. The description is restricted to information necessary to implement and test the presented methods and it is given only for one-dimensional case.

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.

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.

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.

Studies on the detection of layers with elevated black carbon aerosol (BC) concentrations and the formation conditions of these layers help understand the vertical distribution of BC concentrations, which will provide a basis for the assessment of climate effects and early BC pollution warnings. By using the Weather Research and Forecasting with Chemistry (WRF-Chem) numerical model, we performed a numerical simulation analysis on the authenticity of strong elevated BC concentration layers that were detected by an aircraft in the mixing layer over Harbin, China, which is a high-emission area, on a clear sunny afternoon in the early heating period of 2016. We then discuss possible problems and solutions when non-vertical paths are used to detect the vertical distribution of BC concentrations. Finally, we discuss the favorable conditions for the formation of elevated BC concentration layers by weak vertical flow. The results show that the horizontal variability of BC concentration in the mixing layer in the observation area in Harbin was sufficiently large during the measurement. This produced a false elevated layer, as detected by the aircraft during one round of spiral flight in the mixing layer. The root mean square of the horizontal distribution of BC concentration did not change with height in the mixing layer during the daytime, but it decreased with the thickness of the mixing layer and was higher in the mixing layer than in the free atmosphere. Therefore, the thinner the mixing layer, in which the vertical distribution of the BC concentration is detected in an inclined path, the stronger interference of the horizontal variability on the detected results. When a spiral flight detection path is used, the aircraft should fly at least two rounds in the mixing layer. In the daytime, due to strong turbulence in the mixing layer, weak vertical uplift is not favorable for the occurrence of elevated BC concentration layers in the mixing layer. In the nighttime, if weak vertical uplift is well matched with the BC concentration or its vertical gradient, elevated BC concentration layers can be formed in the atmosphere. Compared with upper layers far from the ground, nighttime elevated layers are easier to form in lower layers near the ground because high BC concentrations or large vertical gradients are more likely to occur in the lower layers. Both cases facilitate the occurrence of large vertical upward transport rates of BC.