We are concerned with the use of some classical spectral collocation methods as well as with the new software system Chebfun in order to compute high order (index) eigenpairs of singular as well as regular Schrodinger eigenproblems. We want to highlight both the qualities as well as the shortcomings of these methods and evaluate them vis-a-vis the usual ones. In order to resolve a boundary singularity we use Chebfun with the simple domain truncation technique. Although this method is equally easy to apply with spectral collocation, things are more nuanced in the case of these methods. A special technique to introduce boundary conditions as well as a coordinate transform which maps an unbounded domain to a nite one are the ingredients. A challenging set of "hard" benchmark problems, for which usual numerical methods (f. d., f. e. m., shooting etc.) fail, are analysed. In order to separate "good"and "bad"eigenvalues we estimate the drift of the set of eigenvalues of interest with respect to the order of approximation and/or scaling of domain parameter. It automatically provides us with a measure of the error within which the eigenvalues are computed and a hint on numerical stability. We pay a particular attention to problems with almost multiple eigenvalues as well as for problems with a mixed (continuous) spectrum. In the latter case we try to numerically highlight its existence. Special attention will be paid to the higher eigenpairs (the pair of eigenvalue and the corresponding eigenfunction approximated by an eigenvector spanning its nodal values).

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

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.

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.

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

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

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

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

The 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 annual energy matrix has been changing in the last years because of the necessity of less dependence on fossil fuels, which are running out on the planet. Therefore, a study was started to simulate the most efficient production of H2 through COMSOL Multiphysics®, as it is a gas that in the future will be essential to supply the world's energy necessities, due to its easy to get and insert in piped gas pipes. From this study, we were able to present an optimization of parameters that presents good indications of how to obtain efficient hydrogen production in an idealized reactor.

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

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.

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.

Transition between free-surface and pressurized flows is an crucial phenomenon in many hydraulic systems, including water distribution systems, urban drainage systems, etc. During the transition, the force exerted on the structures changes drastically, thus it is meaningful to simulate this process. However, severe numerical oscillations are widely observed behind filling-bores, causing unphysical pressure variations and even computation failure. In this paper, some oscillation-suppressing approaches are reviewed and evaluated on a benchmark model. Then a new oscillation-suppressing approach is proposed to admit numerical viscosity when the water surface is at proximity of conduct roof which has first order accuracy. This approach adds numerical viscosity when water surface is at the proximity of conduct roof. It can sufficiently suppress numerical oscillations under an acoustic wave speed of 1000m/s and is simple to apply. In comparison with two experiments, the simulation results of this method show good agreement and little numerical oscillations. The results in this paper can help readers to choose an appropriate oscillation-suppressing method to improve the robustness and accuracy of flow regime transition simulations.

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.

We provide a very general review of the resonant transmission line method for optical fiber problems. The method has been found to work seamlessly for a variety of difficult problems including elliptical and eccentric core fibers as well as “holey”, photonic crystal fibers. A latest version has been shown to offer great versatility with respect to cases of unconventional, inhomogeneous index profiles.

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.

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.

The paper continues the study of efficient algorithms for the computation of zeta functions over the complex plane. We aim to apply the modifications of algorithms to the investigation of underlying fractal structures associated with the Riemann zeta function. We discuss the computational complexity and numerical aspects of the implemented algorithms based on series with binomial-like coefficients.

A multiscale modeling concept for numerical simulation of multiphysics turbulent flow utilizing map-based advection is described. The approach is outlined with emphasis on its theoretical foundations and physical interpretations in order to establish the context for subsequent presentation of the associated numerical algorithms and the results of validation studies. The model formulation is a synthesis of existing methods, modified and extended in order to obtain a qualitatively new capability. The salient feature of the approach is that time advancement of the flow is fully resolved both spatially and temporally, albeit with modeled advancement processes restricted to one spatial dimension. This one-dimensional advancement is the basis of a bottom-up modeling approach in which three-dimensional space is discretized into under-resolved mesh cells, each of which contains an instantiation of the modeled one-dimensional advancement. Filtering is done only to provide inputs to a pressure correction that enforces continuity and to obtain mesh-scale-filtered outputs if desired. The one-dimensional advancement, the pressure correction, and coupling of one-dimensional instantiations using a Lagrangian implementation of mesh-resolved volume fluxes is sufficient to advance the three-dimensional flow without time advancing coarse-grained equations, a feature that motivates the designation of the approach as autonomous microscale evolution (AME). In this sense, the one-dimensional treatment is not a closure because there are no unclosed terms to evaluate. However, the approach is additionally suitable for use as a subgrid-scale closure of existing large-eddy-simulation methods. The potential capabilities and limitations of both of these implementations of the approach are assessed conceptually and with reference to demonstrated capabilities of related methods.

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.

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.

The scope of this work was part of our previous research on Pulsed Power Plasma Stimulation Technique. Electromagnetic fields generated during a pulsed power plasma discharge were simulated using a finite element method in the COMSOL RF module. The field distributions were calculated during and after the pulse. The current paper takes advantage of prior work measuring and modeling the current distribution in the plasma arc. Agreement with laboratory experimental measurements provides support for extension of the model to reservoir scale. The validated model was used to quantify the signal attenuation level in different medium environment.

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

A single fracture is the basic unit of fracture medium, and the roughness of fracture wall surface is an important factor influencing hydraulic characteristics of the flow in bedrock fracture. However, effects of the shape and density of roughness elements (various bulges/pits on rough fracture wall surfaces) on water flow in a single rough fracture have not been thoroughly discovered. Thus the water flow in single fracture with different shapes and densities of roughness elements was simulated by using Fluent software in this study. The results show that in wider fractures the flow rate mainly depends on fracture aperture, while in narrow and close fracture medium the surface roughness of fracture wall is the main factor of head loss of seepage; there is a negative power exponential relation between the hydraulic gradient index and the average fracture aperture, i.e. with the increase of fracture aperture, the relative roughness of fracture and the influence weight of hydraulic gradient both decrease; and in symmetrical-uncoupled fractures there is a super-cubic relation between the discharge per unit width and average aperture. Above results would help to deepen the understanding of rough fracture seepage.

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 this work, some numerical radius inequalities based on the recent Dragomir extension of Furuta's inequality are obtained. Some particular cases are also provided.

In this study, we have considered for numerical solution of a Poisson equation in a domain. we have presented a first and second order finite difference method for solving the system of the intial boundary value. We have presented the numerical solution of the Poisson equations in a two-dimensional finite region. Numerical Solution and exact solution for the presented method has been compared and the figures presented. The present error analysis tables demonstrate the efficiency of the 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.

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

This paper examines the influence of physical parameters on the collapse dynamics of a spherical bubble filled with diatomic gas ($\kappa=7/5$). The present numerical investigation shows that each physical parameter affects the bubble collapse dynamics differently. After comparing the contribution of each physical parameter, it appears that, of all the parameters, the surrounding liquid environment affects the bubble collapse dynamics the most. Meanwhile, surface tension has the weakest influence and can be ignored in the bubble collapse dynamics. However, surface tension must be retained in the initial analysis since this, as well as the pressure difference jointly control initial bubble formation. As an essential part of this study, a general Maple code is provided.

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.

Herein we foremost detailed the numerical modeling of the double absorber layer- methyl ammonium lead iodide– carbon nitride layer solar cell and subsequently provided in-depth insight on the active layer associated with dominant radiative and non-radiative recombination losses limiting the efficiency ( ) of the solar cell. Under recombination kinetics phenomena, we explored the influence of Radiative recombination, Auger recombination, Shockley Read Hall recombination, the energy distribution of defects; Band Tail recombination (Hoping Model), Gaussian distribution, metastable defect states including single donor (0/+), single acceptor (-/0), Double Donor (0/+/2+), double acceptor (2/-/0-), and the interface layer defects on the output characteristics of the solar cell. Setting defect (or trap) density to with uniform energy distribution of defects for all the layers, we achieved the of 24. 16 %. A considerable enhancement in power conversion efficiency was perceived as we reduced the trap density to for the absorber layers. Further, it was observed that for the absorber layer with double donor defect states, the active layer should be carefully synthesized to reduce crystal order defects to keep the total defect density as low as to achieve efficient device characteristics

Dynamic Realistic Failure Process Analysis (DRFPA) is a numerical software based on FEM it is firstly validated by comparing with the experimental test. Such validation is carried out under two consistent concrete specimens under uniaxial compression. After validation, DRFPA is used to study the influence of two parallel pre-existing crack inclination angles on the fracture of concrete under dynamic loads. The evolution of wing cracks and secondary cracks is investigated by observing the initiation, coalescence, and propagation of them. Meanwhile, the relationship between failure mode and inclination angles is discussed via analyzing the acoustic emission events. Finally, it is concluded the relationship between inclination angles of pre-cracks and performance/stability of concrete structure.

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.

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.

This paper revisits a classic problem in physics - Hertz elastic dynamics of two colliding elastic spheres. This study obtains impact period in terms of hypergeometric function and successfully combines Deresiewicz's three segmental solutions into one single solution. Our numerical investigation confirms that Deresiewicz's inversion is a good approximation. As an essential part of this study, a general Maple code is provided.

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

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.

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.

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.

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.

The problems of formation and propagation of surface seismic waves in vertically layered and horizontally inhomogeneous media are studied. A mathematical model has been developed for the propagation of surface seismic waves in an inhomogeneous medium based on the use of the fundamental principles of wave dynamics and Fourier transforms. The existence of various types of groups of these waves is proved, and their main properties are studied. It has been established that the reflection coefficients of different groups of a dispersive wave on models of half-spaces (a layered, inhomogeneous medium) with a sharp velocity differentiation are characterized by the same dependence on the angle of incidence as in the absence of dispersion. In the region of the critical angle, which depends on the period of the incident wave, a local decrease in the intensity of wave reflection is observed, which is associated with the origin of the head dispersive wave. It has been found that the reflection coefficients of Rayleigh waves from low-velocity layers, thin (d / λ2 = 0.205) and thick (d / λ2 = 0.82) are approximately 0.2 ± 0.04 and practically do not change with increasing angle of incidence. On the basis of theoretical studies, it was revealed that reflected surface waves carry information both about the elastic properties of the reflecting medium and about the dimensions of the reflecting horizontal inhomogeneity. An experimental dispersion curve of the group velocity of Rayleigh waves on the Central Asia-Caucasus path has been obtained. The seismic stress state of the "subgrade-base" system was studied. Based on the analysis of the results of numerical experiments, it has been established that in order to obtain a complete picture of the seismically stressed state of the subgrade, it is necessary to take into account the contribution to it of each component of the accelerogram of a real earthquake.

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.

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.

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.

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

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

For embankment slope stability, soil cohesivity is one of the most important shear strength parameters. The effect of cohesiveness on the slope stability of a homogeneous embankment dam under rapid drawdown loading conditions was examined in this study. With the use of numerical modeling in GeoStudio, different situations were explored based on cohesiveness (0 kN/m², 5 kN/m², 10 kN/m^2, and 15 kN/m²) under a 1 m per day drawdown rate. The factor of safety value obtained from the long-term steady-state condition under 15 kN/m² cohesion was equivalent to a 116.8% increase from the one obtained under 0 kN/m² cohesion. The factor of safety values obtained after subjecting the embankment to different soil cohesion levels yielded a p-value of 1.91 x 10^-41, according to the Analysis of Variance. The calculated p-value (alpha value) is less than 0.05, suggesting that the differences between the examined cohesiveness values based on the list of the factor of safety values are statistically significant. The findings derived in this study show that it is significant to capture the effect of material characteristics during the design phase of an embankment dam.

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.

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

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

The conventional explanation of delayed-choice experiments seems to violate our causal intuition. This apparent violation is caused by a misinterpretation of the conventional formulation of quantum mechanics. I reanalyze these experiments using advanced and time-symmetric formulations of quantum mechanics. All three formulations give the same experimental predictions, but the advanced and time-symmetric formulations violate our causal intuition that effects only happen after causes. I explore reasons why our causal intuition may be wrong at the quantum level. I also suggest how conventional causation might be recovered in the classical limit, and speculate on cosmological boundary conditions.

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

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.

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

The steady streaming (SS) phenomenon is gaining increased attention in the microfluidics community, because it can generate net mass flow from the zero-mean vibration. We developed numerical simulation and experimental measurement tools to analyze this vibration induced flow, which has been challenging due to its unsteady nature. Validity of these analysis methods is confirmed by comparing the three-dimensional (3D) flow field induced around a cylindrical micropillar under circular vibration. In the numerical modeling, we directly solved the flow in the Lagrangian frame so that the substrate with a micropillar becomes stationary, and the result was converted to the Eulerian frame to compare them with the experimental results. The present approach enables to avoid the introduction of moving boundary or small perturbation approximation. The flow field obtained by the micro particle image velocimetry (PIV) measurement supported the three-dimensionality observed in the numerical results, which could be important for controlling the mass transport and manipulating particulate objects in the microfluidic systems.

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.

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.

Melt mixing is a convenient method to prepare polymer nanocomposites, but the extent of the dispersion of the solid filler reached is often limited, and may compromise the anticipated performance of these materials during service. Since the efficiency of extensional flows on dispersion is now well recognized, several mixers were designed with the aim of inducing both shear and extensional flow components. This work combines experimental and numerical data to better understand the kinetics of the dispersion of graphite nanoplates in a polypropylene melt, using a mixing device that consists of a series of stacked rings with equal outer diameter and alternating larger and smaller inner diameters, thus creating a series of converging/diverging flows. Numerical simulation of the flow using the opensource OpenFOAM software assuming both inelastic and viscoelastic responses predicted the velocity, streamlines, flow type and shear and normal stress fields for the mixer. Experimental and computed data were combined to determine the trade-off between the local degree of dispersion of the PP/GnP nanocomposite, measured as Area ratio, and the absolute average value of the hydrodynamic stresses multiplied by the local cumulative residence time. From considerations based on a theoretical approach to dispersion, the cohesive strength of the GnP agglomerates studied was estimated to be in the range 3 - 10kPa.

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.

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

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.

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

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

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.

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.

Relativistic atomic and molecular structure calculations start from very different premises to quantum electrodynamics (QED). QED predicts radiative corrections that can nly be approximated crudely in today's atomic structure programs. This paper describes rigorous simple algorithms for including accurate vacuum polarization and self-energy corrections in the GRASP programs.

Motivated by the way animals react in nature, we introduce a fully-connected neural-network model, which incorporates the external pattern presented to the system as a fundamental tool of the recognition process. In this new scenario, in the absence of this external field, memories are not attractors inside basins of attraction, as in usual attractor neural networks, although basins may be created according to the external pattern, thus allowing storing a much larger number of memories. The key point consists in calibrating the influence of the external pattern such as to cancel the noise generated by those memories not correlated with the external pattern.We illustrate how this proposal works by including this new contribution in the standard Hopfield model, showing a significant increase in its recognition capacity (typically by a factor $10^{2}$). As an additional feature of this model, we show its ability to react promptly to changes in the external environment, a crucial attribute of living beings. This procedure can be applied to a wide variety of neural-network models to increase considerably their recognition and reaction capabilities.

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.

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.

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.

This study employs a high-fidelity numerical framework to determine the plastic material flow patterns and temperature distributions that lead to void formation during friction stir welding (FSW), and to relate the void morphologies to the underlying alloy material properties and process conditions. Three aluminum alloys, viz., 6061-T6, 7075-T6, and 5053-H18 were investigated under varying traverse speeds. The choice of aluminum alloys enables investigation of a wide range of thermal and mechanical properties. The numerical simulations were validated using experimental observations of void morphologies in these three alloys. Temperatures, plastic strain rates, and material flow patterns are considered. The key results from this study are: (1) The predicted stir zone and void morphology are in good agreement with the experimental observations, (2) The temperature and plastic strain-rate maps in the steady-state process conditions show a strong dependency on the alloy type and traverse speeds, (3) The material velocity contours provide a good insight into the material flow in the stir zone for the FSW process conditions that result in voids as well as those that do not result in voids. The numerical model and the ensuing parametric studies presented in this work provide a framework for understanding material flow under different process conditions in aluminum alloys, and potentially in other alloys. Furthermore, the utility of the numerical model for making quantitative predictions and investigating different process parameters to reduce void formation is demonstrated.

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.

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 twisted Wang transform distribution family, defined as the composition of parameter shifted inverse CDF function with an original CDF function, is found to be most suitable for matching low shape factor distributions, characterizing hard to fit or to simulate reinsurance portfolio losses for some perils from our previous study. Among them, the best form for matching a hard-to-fit empirical loss distribution for a specific peril, is the Exponential Fractional Extra Power 0 Distribution in (0,1) with CDF:.The simplest yet still a good form of this family is the Transformed Hyperbolic Tangent Distribution with CDF:,which has analytical formulas for the moments. The twisted Wang transform distribution family is compared and confirmed to be superior to all other well-known distribution families through extensive numerical optimization practice, distribution forms guesses, and computer-aided exploration experiments.

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.

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.

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.

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

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

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.

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.

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.

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.

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.

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

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

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

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

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

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.

Recently, there has been increasing attention toward inhaled nanoparticles (NPs) to develop in-halation therapies for diseases associated with the pulmonary system and investigate the toxic ef-fects of hazardous environmental particles on human lung health. Taking advantage of microfluidic technology for cell culture applications, lung-on-a-chip devices with great potential in replicating the lung air-blood barrier (ABB) have opened new research insights in preclinical pathology and therapeutic studies associated with aerosol NPs. Surprisingly, the air interface in such devices, which makes them tremendously unique among the other common liquid-based cellular mi-crosystems and enables exposure of cultivated cells to the drug/toxic aerosols, has been largely disregarded. Accordingly, there exists a significant research gap in the comprehensive under-standing of the NPs’ dynamics in the lung-on-a-chip devices, which is momentous to control the injection of aerosols to the target area and attain a desired physiochemical efficacy. To address this research gap, a numerical parametric study is presented to provide a deep insight into the dynamic behavior of the airborne NPs in a gas-liquid dual-channel lung-on-a-chip device with a porous membrane separating the channels. A finite element multi-physics model is developed to enable a particle tracing investigation simultaneously in both air and medium phases to replicate the in vivo conditions avoiding laborious and costly experimental trials. The results include the impact of fluid flow as well as geometrical properties on the distribution, deposition, and translocation of the NPs with diameters ranging from 10 nm to 900 nm in the lung-on-a-chip. The obtained results strongly suggest the aerosol injection of NPs instead of the aqueous solution for more efficient deposition on the substrate of the air channel and higher translocation to the media channel. The current study proposes to optimize the affecting parameters to control the injection and delivery of aerosol par-ticles into the lung-on-chip device depending on the objectives of biomedical investigations and provides optimized values for some specific cases. Therefore, this study can assist scientists and researchers in complementing their experimental investigation in future preclinical studies on pulmonary pathology associated with inhaled hazardous and toxic environmental particles, as well as therapeutic studies for developing inhalation drug delivery.