Barrier coverage is a fundamental application in wireless sensor networks, which is widely used to protect a region through intruder detection for smart city. In this paper, we study a new branch of barrier coverage, named warning barrier coverage (WBC). Different from the classic barrier coverage, WBC has the inverse protect direction, which protects unexpected visitors by warning them away from the region enclosed by sensor barriers. The WBC holds a promising prospect in many danger keep out applications for smart city. For example, a WBC can enclose the debris area in the sea and alarm any approaching ships in order to avoid their damaging propellers. One special feature of WBC is that the target region is usually dangerous and its boundary is previously unknown. Hence, the scattered mobile nodes need to detect the boundary and form the barrier coverage themselves. It is challenging to form these distributed sensor nodes into a barrier because a node can sense only the local information and there is no global information of the unknown region or other nodes. To this end, we propose a novel solution AutoBar for mobile sensor nodes to automatically form a WBC for smart city. To pursue the high coverage quality, we theoretically derive the optimal distribution pattern of sensor nodes using convex theory. Based on the analysis, we design a fully distributed algorithm that enables nodes to collaboratively move towards the optimal distribution pattern. In addition, AutoBar is able to reorganize the barrier even any node is broken. To validate the feasibility of AutoBar, we develop the prototype of the specialized mobile node, which consists of two kinds of sensors: one for boundary detection and another for visitor detection. Based on the prototype, we conduct extensive real trace driven simulations in various smart city scenarios. Performance results demonstrate that AutoBar outperforms the existing barrier coverage strategies in terms of coverage quality, formation duration, and communication overhead.

Traditional numerical techniques such as sliding mesh, dynamic grid mesh, and others, have many limitations in dealing with flow simulation with large-scale movement of solid boundaries; which is the case for complex-shaped hydraulic turbomachinery such as propellers, pumps and turbines. The immersed boundary (IB) method provides a new approach to solve the above limitations. Therefore, this study proposes an sharp-interface IB method based on the level-set function that is suitable for simulating the flow through turbomachinery with complex geometries. This method is applied to actual three-dimensional numerical simulations of high-Reynolds number propellers using an in-house computational fluid dynamics solver. The results show that the proposed method can provide comparatively accurate predictions of unsteady load coefficients within the propeller flow passage, and capture the correct propeller wake characteristics as well as the interaction between the propeller wake and free surface. This study aims at providing a theoretical basis and engineering reference for the application of the IB method in engineering numerical simulations.

With permafrost degenerating caused by climate change, water accumulation increases in permafrost regions during last decades. Water accumulation will deteriorate existing status of engineering in cold regions. Water accumulation can make a thermal effect on permafrost under construction, even result in failure of the subgrade. Moreover, the thermal effect is related to water temperature. However, temperature variation of water accumulation is complex, and its influence factors include air temperature, environment, scope of water accumulation and so on. In order to make analysis of the damage mechanism of water accumulation on permafrost, it is necessary to explore the internal temperature change of water accumulation. This paper proposes a review about temperature calculation method of water accumulation in cold environment. The thermal calculation method between air and water boundary of water accumulation is summarized. Water temperature change of water accumulation with various type is analyzed. The thermal calculation considering phase transformation in water accumulation is discussed, and heat transfer from the bottom of water accumulation to the underlying soil is further studied. Finally, the key factors which are advantageous to make research about the thermal effect of water accumulation in permafrost is proposed to optimize calculation method.

This study aimed to establish and verify the validity of an acoustic simulation method during sustained phonation of the Japanese vowels /i/ and /u/. The study participants were six healthy adults. First, vocal tract models were constructed based on computed tomography (CT) data, such as the range from the frontal sinus to the glottis, during sustained phonation of /i/ and /u/. Next, cylindrical shapes virtually extended by 12 cm were added to the vocal tract models to imitate the trachea between the tracheal bifurcation and lower part of the glottis. The Kirchhoff–Helmholtz integral equation was formulated as the wave equation for sound propagation, and the boundary element method was used for discretization. As a result, the relative discrimination thresholds of the vowel formant frequencies for /i/ and /u/ against actual voice were 1.1%–10.2% and 0.4%–9.3% for the first formant and 3.9%–7.5% and 5.0%–12.5% for second formant, respectively. In the vocal tract model with nasal coupling, a pole–zero pair was observed at around 500 Hz, and for both /i/ and /u/, a pole–zero pair was observed at around 1000 Hz regardless of the presence or absence of nasal coupling. These findings demonstrated that /i/ and /u/ could be simulated with high validity in a vocal tract model constructed from CT data obtained during sustained phonation using the boundary element method.

The magnetic polarization tensor has a promising capability of determining the geometry and material properties of metallic samples. In this paper, a novel computation method is proposed to estimate the magnetic polarization tensors for the metallic samples using the boundary element method. In this method, the metallic sample is placed in a uniformly distributed magnetic field. Based on assumptions that the excitation frequency and/or the conductivity of the sample is very high, the metallic sample is regarded as a perfect electrical conductor (PEC). Therefore, the scattered field at a certain distance can be simulated. By utilising the boundary element method, the magnetic polarization tensor can be derived from the simulated scattered field. The theoretical calculation is presented and simulations and experiments have been carried out to validate the proposed method. The results from the simulation are matched with the analytical solution for the case of sphere samples. Moreover, there is a good agreement between the simulation results and the experimental results for the copper cylindrical samples.

The famous Toms effect (1948) consists of a substantial increase of the critical Reynolds number when a small amount of soluble polymer is introduced into water. The most noticeable influence of polymer additives is manifested in the boundary layer near solid surfaces. The goal of the present paper is a group analysis of the boundary layer equations in two mathematical models of the flow of aqueous polymer solutions: the second grade fluid (Rivlin and Ericksen, 1955) and the model derived by Pavlovskii (1971). The equations of the unsteady two-dimensional boundary layer in the Pavlovskii and Rivlin-Ericksen models are analyzed for the first time here. These equations have no definite type so that finding their exact solutions is very important in order to understand the mathematical nature of the above mentioned models. The problem of group classification with respect to the arbitrary function of the longitudinal coordinate and time present in the equations, which sets the pressure gradient of the external flow, arises. All functions for which an extension of the admitted Lie group occurs are found. The task includes the ratio of two characteristic length scales. One of them is the Prandtl scale, and another is defined as the square root of the normalized coefficient of relaxation viscosity (Frolovskaya and Pukhnachev, 2018) and does not depend on the characteristics of the motion. The paper contains a number of exact solutions in the Pavlovskii model including a solution describing the flow near a critical point. Among the solutions of the new model of the boundary layer, a special place is taken by the solution of the stationary problem of flow around a rectilinear plate. Within the framework of the Prandtl theory of the boundary layer, such a solution was constructed by Blasius (1908). As is well-known, this solution has a non-removable defect: the transverse velocity near the edge of the plate increases without bound. The introduction of a relaxation term into the model makes it possible to eliminate this singularity.

The three-axes Doppler sodar Latan-3 operated on an oceanographic stationary platform in the coastal zone of the Black Sea in June 2015. The platform is located 450 meters offshore from the southern coast of the Crimea Peninsula in the region of Katsiveli (44.39°N, 33.99°E). The water depth at the site is about 30 meters. The atmospheric boundary layer (ABL) typical for the mediterranean seas was observed when the wind is from the sea. The physical processes typical for the coastal mountain terrain was observed when the wind was from the shore. Complex measurements of the ABL parameters were performed using a sodar. Auxiliary measurements of the ABL parameters were performed using a temperature profiler and an ultrasonic thermometer-anemometer. Observations were made mostly during a fair weather with a pronounced diurnal course of meteorological parameters. Sodar data analysis revealed a strong wave activity in the ABL. Internal gravity waves with amplitudes of up to one hundred meters were regularly observed in a layered turbulence structure under stable conditions. Various forms of Kelvin-Helmholtz billows were observed at the interface between the sea breeze and the return flow aloft and in the low level jets.

Bioconvective flows have attracted attention in recent years due to actual and potential applications. In this paper, we consider a steady and laminar convective MHD flow of a nanofluid with heat, mass and microorganism transfer with a heat source/sink present. In addition, we assume there exists a first order chemical reaction. The governing partial differential equations (PDEs) are reduced to ordinary differential equations (ODEs) using the scaling group transformation and the associated boundary value problem is then solved. The influences of selected governing parameters on the dimensionless velocity, temperature, nanoparticle concentration, density of motile microorganisms, skin friction, heat transfer, mass transfer, and motile microorganism density rates are computed and discussed

We performed for the first time to our knowledge fully ab initio molecular dynamics simulations of additive tribochemistry in boundary lubrication conditions. We consider an organophosphourus additive that has been experimentally shown to reduce friction in steel-on-steel sliding contacts thanks to the tribologically-induced formation of an iron phosphide tribofilm. The simulations allow us to observe in real time the molecular dissociation at the sliding iron interface under pressure and to understand the mechanism of iron phosphide formation. We discuss the role played by the mechanical stress by comparing the activation times for molecular dissociation observed in the tribological simulations at different applied loads with that expected on the basis of the dissociation barrier.

Non-uniform grain boundary sliding can induce strain and rotation incompatibilities at perfectly planar interfaces. Explicit analytic expressions of stress and lattice rotation jumps are thus derived at a planar interface in the general framework of heterogeneous anisotropic thermo-elasticity with plasticity and grain boundary sliding. Both elastic fields are directly dependent on in-plane gradients of grain boundary sliding. It is also shown that grain boundary sliding is a mechanism that may relax incompatibility stresses of elastic, plastic and thermal origin although the latter are not resolved on the grain boundary plane. This relaxation may be a driving force for grain boundary sliding in addition to the traditionally considered local shears on the grain boundary plane. Moreover, the obtained analytic expressions are checked by different kinds of bicrystal shearing finite element simulations allowing grain boundary sliding and where a pinned line in the interface plane aims at representing the effect of a triple junction. A very good agreement is found between the analytic solutions and the finite element results. The performed simulations particularly emphasize the role of grain boundary sliding as a possible strong stress generator around the grain boundary close to the triple line because of the presence of pronounced gradients of sliding.

Large values and gradients of stresses and deformations, triggering concentrations of stresses and deformations, arise in the corner areas of a structure. The action of forced deformations, leading to the finite rupture of the contact between the elements of a structure, also triggers the concentration of stresses, while the rupture reaches an irregular point, a line on the area boundary. The theoretical analysis of the stress-strain state (SSS) of areas with angular cutouts in the boundary under the action of discontinuous forced deformations is reduced to the study of singular solutions to the homogeneous problem of the elasticity theory that has power-related features. The calculation of stress concentration coefficients in the domain of a singular solution to the elastic problem makes no sense. It is experimentally proven that the zone, that is close to the vertex of the angular cutout in the area boundary, has substantial deformations, rotations, and it corresponds to rising values of the first and second derivatives of displacements along the radius in cases of sufficiently small radii in the neighbourhood of the irregular point of the boundary. For such areas, it is necessary to consider the plane problem of the elasticity theory, taking into account the geometric nonlinearity under the action of forced deformations. This will allow analyzing the effect of relations between orders of values of deformations, rotations, and forced deformations on the form of the equation of equilibrium. The purpose of this work is to analyze the effect of relations of deformation orders, rotations, forced deformations on the form of the equilibrium equation in the polar coordinate system for a V-shaped area under the action of forced temperature-induced deformations with regard for the geometrical non-linearity and physical linearity.

The task of weakly supervised point cloud semantic segmentation has received widespread attention and also has been widely used in autonomous driving, robotics, and modern industries. Due to the high dimensionality and large volume of point cloud data, many technical difficulties appeal to weakly supervised semantic processing, especially segmentation. However, weakly-supervised schemes usually provide only partial labelling information of the underlying point cloud data and thus need to effectively extract local features, geometric information, and contextual relationships only using these limited labelled data for supervised learning. To improve weakly-supervised semantic segmentation, we propose a novel segmentation network through the boundary-based feature aggregation based on a K-NN algorithm with down-sampling operation, and also introduce the smoothness loss and Siamese loss for effective segmentation. The experiments on public datasets demonstrate that our presented segmentation method exceeds most of the existing fully supervised and weakly-supervised methods in terms of mIoU. Specifically, our network has high segmentation accuracy on the labels of objects with similar geometrical structures, such as ceiling, wall, floor, chair and table, reaching 91.2%, 98.8%, 83.3%, 75.3% and 80.2%, respectively. Extensive experiments also illustrate its robustness, effectiveness, and generalization of the proposed weakly supervised segmentation network.

Two-dimensional heat diffusion problem with heat source which is quasilinear parabolic problem is examined analytically and numerically. Periodic boundary conditions are used as a boundary conditions. Since the problem is not linear, Picard’s successive approximation theorem is used. Under certain conditions of natural regularity and consistency imposed on the input data, establish the existence, uniqueness and constant dependence of the solution on the data using the generalized Fourier method. As a numerical solution, implicit finite difference scheme is used. The results ob-tained from analytical and the numerical solutions are so close to each other.

Haemodynamic simulations are increasingly used to study vascular diseases like Intracranial Aneurysms (IA) and to further develope treatment options. However, due to limited data, some aspects must rely on heuristics, especially at the simulation’s distal ends. In literature, Murray’s Law is often used to model the outflow split based on vessel cross-section area, but this poses challenges for the communicating arteries in the Circle of Willis (CoW). In this study, we contribute by assessing the impact of Murray’s Law in patient-specific geometries featuring IA at the posterior communication. We simulate different domain extensions, representing common modelling choices. We establish Full CoW simulations as a baseline to evaluate the effect of these modelling assumptions on haemodynamic indicators, focusing on IA growth and rupture-related factors like Wall Shear Stress (WSS) and Oscillatory Shear Index (OSI). Our findings reveal qualitative alterations in haemodynamics when not modeling posterior communication. Comparisons between computing the anterior circulation and computing the whole Circle of Willis reveal quantitative changes in WSS may reach up to 80%, highlighting the significance of modelling choices in assessing IA risks and treatment strategies.

In this paper we overview all analytical engineering solutions for corrosion problems in pressure vessels delivered over the past 60 years. We briefly detail strengths and weaknesses of existing approaches, thus demonstrating the need for novel uniform corrosion analysis methods for both current and new applications. To complement the review, we also present a new analytical model allowing to estimate the lifetime of pressured elastic vessels with mechanically assisted corrosion. Both internal/external pressure and internal/external corrosion are captured and dissolution driven corrosion is also considered. The readily implemented method improves existing analytical approaches and is shown to be effective for thin and thick shells, as well as various loading and corrosion intensity and geometrical conditions.

In the laboratory practice, several standards for testing the bending stiffness of corrugated board are used. There are often cases of tests where the results depend on the way the sample is placed on the supports. The problem arises when the board is five-ply (with two corrugated layers with different corrugation heights) or when the board has asymmetrically selected papers on the flat layers. In this article, we focus our attention on the problem related to boundary conditions, with particular attention to the local effects of the support of the sample. Since the cardboard layers, both flat and corrugated, have a small thickness, a slight deformation of the papers can always be observed at the point of contact between the sample and the support, which affects the readings of the measured stiffness. The paper presents theoretical and numerical analyzes showing how much the method of supporting the sample affects the measured bending stiffness of various samples. Numerical observations were also compared with the results of analyzes presented by other scientists as well as with experimental results.

Mechanochemical technology is developing rapidly, judging by the scientific information in both basic and applied studies. However, many issues and points of view remain to be discussed. This review presents some new key issues for the optimization of mechanochemical processes in theoretical and practical aspects. Emphasis is placed on powder technology aspects, which are not always discussed compared to functional or microscopic viewpoints. The transfer of chemical species across the interparticle interface between organic and inorganic species during the mechanosynthesis of nanocomposite oxides offers many new possibilities, which are discussed here in detail. Since all material transport is preceded by charge transfer, its driving force has been searched for using terminology beyond the well-established electrochemical ones. The role of organic compounds during the whole process is emphasized, regardless of their survival in the final product. The similarity with the pharmaceutical phenomena is pointed out, although its mentality is very different from that of the synthesis of nanocomposites. Not only for the rational amorphization of the active pharmaceutical ingredient, but also for the stabilization of molecular dispersion states with the participation of excipients are discussed. The effects of liquids, either added or formed by mechanochemical auto-liquefaction, are presented with reference to the comparison between wet and dry grinding. Finally, mechanisms of apparent stabilization of mechanically activated states of products are elucidated to investigate the practical applicability of mechanochemically synthesized products

Due to the numerous heat sources in electronic equipment, the periodic nature of stream-wise flow occurs in a cooling channel so often, necessitating the development of an effective technique to improve the heat-cooling convection in such a situation. This paper investigates thermal convection enhancement in a heated-block duct through the periodic boundary conditions using the element-by-element (EBE) treatment in a semi-implicit projection finite element method (FEM) through preconditioned conjugate gradient (PCG) solver. This paper used the time-mean Nusselt number ratio, friction factor, and thermal performance to study how installing rectangular cylinders in the channel works with periodic boundary conditions varying Reynolds numbers (100, 175, and 250). The outcomes show that the rectangular cylinders installed stream-wise above an upstream block promote thermal convection in the heated-block duct. However, increasing the number of rectangular cylinders can raise friction factor. As a result, the case for periodic boundary conditions with a rectangular cylinder above every two blocks has the best thermal performance.

The main goal of this work is to develop a novel boundary element method (BEM) model for analyzing ultrasonic wave propaga-tion in three-temperature anisotropic viscoelastic porous media. Because of strong nonlinearity of ultrasonic wave propagation in three-temperature porous media problems, the analytical or numerical solutions to the problems under consideration are always challenging, necessitating the development of new computational techniques. As a result, we use a new BEM model to solve such problems. A time stepping procedure based on the linear multistep method is obtained after solving the discretized boundary in-tegral equation with the quadrature rule. The calculation of a double integral is required to obtain fundamental solutions, but this increases the total BEM computation time. Our proposed BEM technique is used to solve the current problem and improve the formulation efficiency. The numerical results are graphed to demonstrate the effects of viscosity and anisotropy on the nonlinear ultrasonic stress waves in three-temperature porous media. The validity, accuracy, and efficiency in the proposed methodology were demonstrated by comparing obtained results to the corresponding solution obtained from finite difference method (FDM).

Laminate structures composed of stiff plates and thin soft interlayers are widely used in aerospace, automotive and civil engineering encouraging the development of reliable non-destructive strategies for their condition assessment. In the paper, elastodynamic behaviour of such laminate structures is investigated with emphasis on its application in ultrasonic based NDT and SHM for the identification of interlayer mechanical and interfacial contact properties. A particular attention is given to the practically important frequency range, in which the wavelength considerably exceeds the thickness of the film. Three layer model with spring-type boundary conditions employed for imperfect contact simulation is used for numerical investigation. Novel effective boundary conditions are derived via asymptotic expansion technique and used for analysis of the peculiar properties of elastic guided waves in considered laminates. It is revealed that the thin and soft film influences the behaviour of the laminate mainly via the effective stiffnesses being a combination of the elastic moduli of the film, its thickness and interface stiffnesses. To evaluate each of these parameters separately (or to figure out that the available experimental data are insufficient), a step-wise procedure employing the effective boundary conditions is proposed and tested versus the laser Doppler vibrometry data for Lamb waves in Aluminium/Polymer film/Alumunium structure. The possibility of using film-related thickness resonance frequencies to estimate the film properties and contact quality is also demonstrated. Additionally, the rich family of edge waves is also investigated, and the splitting of fundamental edge waves into pairs is revealed.

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

This work investigates centered polygonal lacunary functions restricted from the unit disk onto symmetry angle space which is defined by the symmetry angles of a given centered polygonal lacunary function. This restriction allows for one to consider only the p-sequences of the centered polygonal lacunary functions which are bounded, but not convergent, at the natural boundary. The periodicity of the $p$-sequences naturally gives rise to a convergent subsequence, which can be used as a grounds for decomposition of the restricted centered polygonal lacunary functions. A mapping of the unit disk to the sphere allows for the study of the line integrals of restricted centered polygonal that includes analytic progress towards closed form representations. Obvious closures of the domain obtained from the spherical map lead to four distinct topological spaces of the "broom topology'' type.

An improved anti-noise morphology vision navigation algorithm is proposed for intelligent tractor tillage in a complex agricultural field environment. At first the two key steps, Guided Filtering and improved anti-noise morphology navigation line extraction were addressed in detail. Then the experiments were carried out in order to verify the effectiveness and advancement of the presented algorithm. Finally, the optimal template and it’s application condition were studied for improving the image processing speed. The comparison experiment results show that the YCbCr color space has minimum time consumption of 0.094 s in comparison with HSV, HIS and 2R-G-B color spaces. The Guided Filtering method can effectively distinguish the boundary between the new and old soil than other competing vanilla methods such as Tarel, Multi-scale Retinex, Wavelet-based Retinex and Homomorphic Filtering inspite of having the fastest processing speed of 0.113 s. The extracted soil boundary line of the improved anti-noise morphology algorithm has best precision and speed compared with other operators such as Sobel, Roberts, Prewitt and Log. After comparing different size of image template, the optimal template with the size of 140×260 pixels can meet high precision vision navigation while the course deviation angle is not more than 7.5°. The maximum tractor speed of the optimal template and global template are 51.41 km/h and 27.47 km/h respectively which can meet real-time vision navigation requirement of the smart tractor tillage operation in the field. The experimental vision navigation results demonstrated the feasibility of the autonomous vision navigation for tractor tillage operation in the field using the new and old soil boundary line extracted by the proposed improved anti-noise morphology algorithm which has broad application prospect.

An improved anti-noise morphology vision navigation algorithm is proposed for intelligent tractor tillage in a complex agricultural field environment. Firstly, the two key steps, Guided Filtering and improved anti-noise morphology navigation line extraction, were addressed in detail. Then the experiments were carried out in order to verify the effectiveness and advancement of the presented algorithm. Finally, the optimal template and its application condition were studied for improving the image processing speed. The comparison experiment results show that the YCbCr color space has minimum time consumption, 0.094 s, compared with HSV, HIS and 2R-G-B color spaces. The Guided Filtering method can enhance the new & old soil boundary effectively than any other methods such as Tarel, Multi-scale Retinex, Wavelet-based Retinex and Homomorphic Filtering, meanwhile, has the fastest processing speed of 0.113 s. The extracted soil boundary line of the improved anti-noise morphology algorithm has best precision and speed compared with other operators such as Sobel, Roberts, Prewitt and Log. After comparing different size of image template, the optimal template with the size of 140×260 pixels can meet high precision vision navigation while the course deviation angle is not more than 7.5°. The maximum tractor speed of the optimal template and global template are 51.41 km/h and 27.47 km/h respectively which can meet real-time vision navigation requirement of the smart tractor tillage operation in the field. The experimental vision navigation results demonstrated the feasibility of the autonomous vision navigation for tractor tillage operation in the field using the new & old soil boundary line extracted by the proposed improved anti-noise morphology algorithm which has broad application prospect.

Abstract: Characterization the 3-D structure of clouds is needed for a more complete understanding of the Earth's radiative and latent heat fluxes. Here we develop and explore a “ray casting” algorithm applied to the Multi-angle Imaging SpectroRadiometer (MISR) on board the Terra satellite, to reconstruct 3-D cloud volumes for observed clouds. The ray casting algorithm is first applied to geometrically simple synthetic clouds to show that, under the assumption of perfect, clear-conservative cloud masks, the reconstruction method yields overestimation whose magnitude depends on the cloud geometry and the resolution of the reconstruction grid relative to the image pixel resolution. The method is then applied to two select MISR scenes, fully accounting for MISR’s viewing geometry for reconstructions over the Earth’s ellipsoidal surface. The MISR Radiometric Camera-by-camera Cloud Masks at 1.1 km resolution and custom cloud masks at 275 m resolution independently derived from MISR RGB channels are used as input cloud masks. A wind correction method, termed “cloud spreading”, is devised and applied to the cloud masks to offset potential cloud movements over short time intervals (around 7 minutes at maximum) between the cameras. The MISR cloud top height product is used as a constraint to reduce the overestimation at the cloud top. The reconstruction results show that their uncertainty is significant when the wind correction is applied, and that they have more refined structures when the input cloud mask has a higher resolution. Recommendations for improving the presented cloud volume reconstructions as well as for future passive remote sensing satellite missions are discussed.

This paper develops feasible control strategies and associated system responses to bring an autonomous front-axle bicycle robot from specified initial conditions to final conditions such that a specific performance index is minimized. To solve the problem, the following approach is used: The feasible controls derived from the normal equations of optimality are substituted into the state and the costate systems and form a combined control-free state-costate system which is vectorized to enable and ease the application of a numerical method. A computer program written in Matlab computer programming language, codes a fourth-order Runge-Kutta numerical method and then solve the combined state-costate system of ordinary differential equations. The obtained results are the feasible bicycle robot trajectory, the feasible state functions, the feasible costate functions and the feasible control functions. Associated Computational Simulations are designed and provided to persuade on the effectiveness and the reliability of the approach.

A family of linear singularly perturbed Cauchy problems is studied. The equations defining the problem combine both partial differential operators together with the action of linear fractional transforms. The exotic geometry of the problem in the Borel plane, involving both sectorial regions and strip-like sets, gives rise to asymptotic results relating the analytic solution and the formal one through Gevrey asymptotic expansions. The main results lean on the appearance of domains in the complex plane which remain intimately related to Lambert W function, which turns out to be crucial in the construction of the analytic solutions. On the way, an accurate description of the deformation of the integration paths defining the analytic solutions and the knowledge of Lambert W function are needed in order to provide the asymptotic behavior of the solution near the origin, regarding the perturbation parameter. Such deformation varies depending on the analytic solution considered, which lies in two families with different geometric features.

The work is devoted to the study of a family of linear initial value problems of partial differential equations in the complex domain, dealing with two complex time variables. The use of a truncated Laplace-like transformation in the construction of the analytic solution allows to overcome a small divisor phenomenon arising from the geometry of the problem and represents an alternative approach to the one proposed in a recent work by the first two authors. The result leans on the application of a fixed point argument and the classical Ramis-Sibuya theorem.

Data obtained by Direct Numerical Simulations (DNS) of pressure-driven turbulent channel flow are studied in the range 180 ≤ Reτ≤ 10,000. Reynolds number effects on the Mean Velocity Profile (MVP) and second order statistics are analyzed with a view of finding logarithmic behavior in the overlap region or even further from the wall, well in the boundary layer’s outer region. Values of Kármán constant for the MVPs and the Townsend-Perry constants for the streamwise and spanwise fluctuation variances are determined for the Reynolds numbers considered. A data-driven model of the MVP, proposed and validated for zero pressure gradient flow over a flat plate, is employed for pressure-driven channel flow by appropriately adjusting Coles’ strength of the wake function parameter, Π. There is excellent agreement between the analytic model predictions of MVP and the DNS-computed MVP as well as of the Reynolds shear stress profile. The skin friction coefficient Cf is calculated analytically. The agreement between the analytical model predictions and the DNS-based computed discrete values of Cf is excellent.

Ship wave is of great interest for wave drag, coastal erosion and ship detection. In this paper, a highly-paralleled numerical scheme is proposed for simulating three-dimensional (3-D) nonlinear Kelvin ship waves effectively. First, a numerical model for nonlinear ship waves is established based on potential flow theory, boundary integral method and Jacobian-free Newton-Krylov (JFNK) method. In order to improve computational efficiency and reduce data storage of JFNK method, a banded preconditioner method is then developed by formulating the optimal bandwidth selection rule. After that, a Graphics Process Unit (GPU) based parallel computing framework is designed, and a GPU solver is developed by using Compute Unified Device Architecture (CUDA) language. Finally, numerical simulations of 3-D nonlinear ship waves under multiple scales are performed by using the GPU and CPU solvers. Simulation results show that the proposed GPU solver is more efficient than the CPU solver with the same accuracy. More than 66% GPU memory can be saved and the computational speed can be accelerated up to 20 times. Hence, the computation time for Kelvin ship waves simulation can be significantly reduced by applying the GPU parallel numerical scheme, which lays a solid foundation for practical ocean engineering.

This study has simulated the typical rainstorm on 20 July 2021 over central east China by using the first-generation Chinese Reanalysis datasets and Global Land Data Assimilation System datasets, and the Noah land surface model coupled with the advanced weather research and forecasting model. Based on this, the gridded planetary boundary layer (PBL) profiles and ensemble states within soil perturbations are collected to investigate the main land-atmosphere coupling characteristics during this modeled rainstorm by using various local coupling metrics and the introduced ensemble statistical metrics. Results have shown that (1) except for the stratospheric thermodynamics and surface thermal over mountain areas, the main characteristics of mid-low layers and surface have been well documented in this modeled rainstorm; (2) the typical coupling intensity is characterized by the dominant morning moistening, early noon weak PBL warming around 2, noontime buoyant mixing temperature deficit around 274 K, daytime PBL and surface latent flux contribution around 100 and 280 W/m2 respectively, and significant afternoon soil-surface latent flux coupling; (3) moist static energy is more significant than PBL height during the relation chains, which is consistent with the significance of surface moistening indicated by local coupling metrics. In general, wet soil contributes greatly to daytime moisture evaporation, which then increases the early noon PBL warming and enhances the noontime buoyant mixing within weak flux contribution. However, this has been suppressed by large-scale forcing such as the upper southwestern inflows of rainstorms, which has further significantly shaped the spatial distribution of statistical metrics in contrast. These quantitatively described local couplings have highlighted both the convection potential diagnoses usage for the local weather application and more applicable coupling threshold diagnoses within the finer spatial investigation.

The wind observations for multiple levels (40–200 m) have been conducted for a long time (2016–2020) on Jeju Island of South Korea. This study aims at understanding the vertical and temporal characteristics of lower atmosphere. Jeju Island is a region located at mid-latitude and is affected by seasonal monsoon wind. The maximum wind speed appears in the lower layer during day time and is delayed in the upper layer during latter time in diurnal cycle. In summer season, the surface layer increases up to 160 m during day time via dominant solar radiation effect, which is higher than those for other seasons. However, the maximum wind speed in winter season appears irregularly among altitudes, and the surface layer is ~100 m, which is lower than that in summer season. It can be attributed to the increase in the mean wind speed in diurnal cycle caused by the strong northwestern wind for winter season. These results imply that the relationship between near-surface and higher altitudes is primarily affected by solar radiation and seasonal monsoon winds. These results are expected to contribute to site selection criteria for wind farms and to the assessment concerning planetary boundary layer modeling.

Magnetorheological fluids (MRFs) are smart materials consisting of micro-scale magnetizable particles suspended in a carrier fluid. The rheological properties of a MRF can be changed from a fluid-state to a solid-state upon the application of an external magnetic field. This study reports the development of a particle-level simulation code for magnetic solid spheres moving through an incompressible Newtonian carrier fluid. The numerical algorithm is implemented within an open-source finite-volume solver coupled with an immersed boundary method (FVM-IBM) to perform fully-resolved simulations. The particulate phase of the MRF is modeled using the discrete element method (DEM). The resultant force acting on the particles due to the external magnetic field (i.e., magnetostatic polarization force) is computed based on the Clausius-Mossotti relationship. The fixed and mutual dipole magnetic models are then used to account for the magnetic (MAG) interactions between particles. Several benchmark flows were simulated using the newly-developed FVM-IBM-DEM-MAG algorithm to assess the accuracy and robustness of the calculations. First, the sedimentation of two spheres in a rectangular duct containing a Newtonian fluid is computed without the presence of an external magnetic field, mimicking the so-called drafting-kissing-tumbling (DKT) phenomenon. The numerical results obtained for the DKT case study are verified against published data from the scientific literature. Second, we activate both the magnetostatic polarization and the dipole-dipole forces and resultant torques between the spheres for the DKT case study. Next, we study the robustness of the FVM-IBM-DEM-MAG solver by computing multi-particle chaining (i.e., particle assembly) in a two-dimensional (2D) domain for area volume fractions of 20% (260 particles) and 30% (390 particles) under vertical and horizontal magnetic fields. Finally, the fourth computational experiment describes the multi-particle chaining in a three-dimensional (3D) domain allowing to study fully-resolved MRF simulations of 580 magnetic particles under vertical and horizontal magnetic fields.

The Oldroyd–B model for a linear viscoelastic fluid is employed to investigate the buoyant flow in a vertical porous layer with permeable boundaries kept at different uniform temperatures. Seepage flow in the viscoelastic fluid saturated porous layer is modelled through an extended version of Darcy’s law taking into account the Oldroyd–B rheology. The basic stationary flow is parallel to the vertical axis and describes a single–cell vertical pattern where the cell has an infinite vertical height. A linear stability analysis of such a basic flow is carried out to determine the onset conditions for a multicellular pattern. The neutral stability curves and the values of the critical Rayleigh number are evaluated numerically for different retardation time and relaxation time characteristics of the fluid.

Haemodynamic indices are widely used in clinical practice for deciding on a particular type of treatment. Low quality of the CT data and tachycardia complicate interpretation of the measured or simulated values. In this work, we present a novel approach for evaluating resistances in terminal coronary arteries. Using 14 measurements from 10 patients, we show that this algorithm retains the accuracy of 1D haemodynamic simulations in less detailed (truncated) geometric models of coronary networks. We also apply the variable systole fraction model to study the effect of elevated heart rate on the values of FFR, CFR and iFR. We conclude that tachycardia may produce both overestimation or underestimation of coronary stenosis significance.

The accuracy of wind field simulation and prediction is one of the most significant parameters in the field of atmospheric science and wind energy. Limited by the observation data, there are few researches on wind energy development. A 3D Doppler wind lidar (DWL) providing the high-vertical-resolution wind data over the urban complex underlying surface in February 2018 was employed to evaluated the accuracy of vertical wind field simulation systematically for the first time. 11 PBL schemes of the Weather Research and Forecasting Model (WRF) were employed in simulation. The model results were evaluated in groups separated by weather (sunny days, haze days and windy days), observation height layers, and various observation wind speeds. The test results presented that the vertical layer altitude of the observation point position was the most important factor. The simulation is fairly well at a height of 1000-2000m, as most of the relative mean bias of wind speed and wind direction are less than 20% and 6% respectively. Below 1000 m, the wind speed and direction biases are about 30%-150% m.s-1 and 6%-30% respectively. Moreover, when the observed wind speed was lower than 5 m.s-1, the bias were usually large, and the wind speed relative mean bias is up to 50-300%. In addition, the accuracy of simulated wind profile is better in 10-15m.s-1 than other speed ranges, and is better up 1000m than below 1000m in the boundary layer. We see that the WRF boundary layer schemes have different applicability to different weather conditions. The WRF boundary layer schemes have significant differences in wind field simulation with larger error under the complex topography. A PBL scheme is not likely to maintain its advantages in the long term under different conditions including altitude and weather conditions.

Starting with unitary quantum dynamics, we investigate how to add quantum measurements. Quantum measurements have four essential components: the furcation, the witness production, an alignment projection, and the actual choice decision. The first two components still lie in the domain of unitary quantum dynamics. The decoherence concept explains the third contribution. It can be based on the requirement that witnesses reaching the end of time on the wave function side and the conjugate one have to be identical. In this way, it also stays within the quantum dynamics domain. The surjection hypothesis explains the actual choice decision. It is based on a two boundary interpretation applied to the complete quantum universe. It offers a simple way to reduce these seemingly random projections to purely deterministic unitary quantum dynamics, eliminating the measurement problem.

Spatiotemporal dynamics of the atmospheric kinetic energy and its components caused by the ordered and turbulent motions of air masses are estimated from minisodar measurements of three velocity vector components and their variances within the lowest 5–200 m layer of the atmosphere, with a particular emphasis on the turbulent kinetic energy. The layered structure of the total atmospheric kinetic energy has been established. From the diurnal hourly dynamics of the altitude profiles of the turbulent kinetic energy (TKE) retrieved from minisodar data, four layers are established by the character of the altitude TKE dependence, namely, the near-ground layer, the surface layer, the layer with a linear TKE increase, and the transitive layer above. In the first layer, the most significant changes of the ТКЕ were observed in the evening hours. In the second layer, no significant changes in the TKE values were observed. A linear increase in the TKE values with altitude was observed in the third layer. In the fourth layer, the TKE slightly increased with altitude and exhibited variations during the entire observation period. The altitudes of the upper boundaries of these layers depended on the time of day. The MKE values were much less than the corresponding TKE values, they did not exceed 50 m²/s². From two to four MKE layers were distinguished based on the character of its altitude dependence. The two-layer structures were observed in the evening and at night (under conditions of the stable atmospheric boundary layer). In the morning and daytime, the four-layer MKE structures with intermediate layers of linear increase and subsequent decrease in the MKE values were observed. Our estimates demonstrated that the ТКЕ contribution to the total atmospheric kinetic energy considerably (by a factor of 2.5–3) exceeded the corresponding МКЕ contribution.

This paper describes the dielectric relaxation spectroscopy (DRS) method for obtaining glass transition temperature and the onset of crystallinity of skim milk powder (SMP). The methodology consists of exposing SMP samples to numerous constant temperatures, T, and water activities, a_w, and measuring the relaxation time. Measurements are based on changes in the loss-tangent height at its peak frequency. Glass transition times are identified at peaks in the loss-tangent height versus time curve and the onset of crystallinity times are identified when the loss-tangent height collapses. The analysis leads to two boundary curves: glass transition T_g, a_w curve and the onset of crystallinity T_c, a_w curve. They separate the SMP stable (long shelf life) and unstable (short shelf life) regions. Two curve fitting algorithms were developed for (a) a_w², t_g curve where t_g is the time to glass transition and a_w², t_c curve where time t_c is the time to onset of crystallinity and (b) T_g, a_w and T_c, a_w boundary curves. Finally, the DRS data are compared to differential scanning calorimetry (DSC) results.

A computational methodology for the hydrodynamic analysis of horizontal axis marine current turbines is presented. The approach is based on a boundary integral equation method for inviscid flows originally developed for marine propellers and adapted here to describe the flow features that characterize hydrokinetic turbines. To this purpose, semi-analytical trailing wake and viscous-flow correction models are introduced. A validation study is performed by comparing hydrodynamic performance predictions with two experimental test cases and with results from other numerical models in the literature. The capability of the proposed methodology to correctly describe turbine thrust and power over a wide range of operating conditions is discussed. Viscosity effects associated to blade flow separation and stall are taken into account and predicted thrust and power are comparable with results of blade element methods that are largely used in the design of marine current turbines. The accuracy of numerical predictions tend to reduce in cases where turbine blades operate in off-design conditions.

This paper aims at computing feasible control strategies and the corresponding feasible state trajectories to drive an autonomous rear-axle bicycle robot from a given initial state to a final state such that the total running cost is minimized. Pontryagin’s Minimum Principle is applied and derives the optimality conditions from which the feasible control functions, expressed as functions of state and costate variables, are substituted into the combined state-costate system to obtain a new free-control state-costate nonlinear system of ordinary differential equations. A computer program was written in Scilab to solve the combined state-costate system and obtain the feasible state functions, the feasible costate functions and the feasible control functions. Associated Computational Simulations were provided to show the effectiveness and the reliability of the approach.

This paper is a slightly modified, abridged version of a previous work ``Parametric Gevrey asymptotics in two complex time variables through truncated Laplace transforms'' motivated by our contribution in the conference ``Formal and Analytic Solutions of Diff. (differential, partial differential, difference, q-difference, q-difference-differential) Equations on the Internet'' (FASnet20). It aims to clarify and give further detail at some crucial points concerning the asymptotic behavior of the solutions of the problems studied in that work.

This paper is a slightly modified, abridged version of a previous work "Parametric Gevrey asymptotics in two complex time variables through truncated Laplace transforms'' motivated by our contribution in the conference "Formal and Analytic Solutions of Diff. (differential, partial differential, difference, q-difference, q-difference-differential) Equations on the Internet'' (FASnet20). It aims to clarify and give further detail at some crucial points concerning the asymptotic behavior of the solutions of the problems studied in that work.

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

Myanmar is located in the Indo-Burma biodiversity hotspot, but lacks understanding of the bio-diversity hotspots and protected areas within its borders. Climate change affects the effective-ness of protected areas and the potential distribution range of species. Based on Myanmar's re-gional characteristics, endangered status, and representativeness, this paper selects species and identifies biodiversity hotspots using the MaxEnt model. We studied the biodiversity within the terrestrial protected area network in Myanmar, as well as the potential for cross-boundary range shifts among protected areas under different climate scenarios. Species populations outside My-anmar's protected area network may not be able to find suitable habitats in the future, further exacerbating climate-related extinction risks. Additionally, some protected areas in Myanmar may not be fully utilizing their biodiversity representation potential, and with scientific spatial planning, Myanmar's protected areas could support a large number of species. In the context of continued warming, identifying new protected areas that encompass the future and current climatic niches of protected areas can increase the capacity for sustained biodiversity conservation. We hope that this study provides insights for increasing the effectiveness of Myanmar's protect-ed area network and emphasizes the urgency of coordinated protection of protected areas and biodiversity hotspots.

This paper introduces an innovative approach to numerically model Structure-Soil-Structure Interaction (SSSI) by integrating the Boundary Element Method (BEM) and the Finite Element Method (FEM) in a coupled manner. To assess the accuracy of the proposed method, a comparative study is undertaken, comparing its outcomes with those generated by the conventional FEM technique. Alongside accuracy, the computational efficiency aspect is crucial for the analysis of large-scale SSSI problems. Hence, the computational performance of the coupled BEM-FEM method undergoes a thorough examination and is compared with that of the standalone FEM method. The results from these comparisons illustrate the superior capabilities of the proposed method in comparison to the FEM method. The novel approach provides more reliable results compared to traditional FEM methods, serving as a valuable tool for engineers and researchers involved in structural analysis and design.

This work investigates a new nanostructured gas-diffusion-layer (nano-GDL) to improve performance of air-cathode Single-Chamber-Microbial-Fuel-Cells (a-SCMFCs). The new nano-GDLs improves the direct oxygen-reduction-reaction by exploiting the best of nanofibers from electrospinning in terms of high surface ratio to volume and high porosity, and laser-based processing to promote adhesion. Nano-GDLs by electrospinning were fabricated directly collecting two nanofibers mats on the same carbon-based electrode, acting as the substrate. Each layer was designed with a specific function: water resistant, oxygen permeable polyvinylidene-difluoride (PVDF) nanofibers served as a barrier to prevent water-based electrolyte leakage, while an inner layer of cellulose nanofibers was added to promote oxygen diffusion towards the catalytic sites. The maximum current density obtained for a-SCMFCs with the new nano-GDLs is (132.2 ± 10.8) mA m^-2, and it doubles the current density obtained with standard PTFE-based GDL (58.5 ± 2.4 mA m^-2), used as reference material. The energy recovery (EF) factor, i.e. the ratio of the power output to the inner volume of the device, was then used to evaluate the overall performance of a-SCMFCs. a-SCMFCs with nano-GDL provided an EF value of 60.83 mJ m^-3: one order of magnitude higher than the value of 3.92 mJ m^-3 obtained with standard GDL

The microstructure of oil pipeline steels is polycrystalline due to the presence of a complex grain boundary (GB) network in them. The performance of them when exposed to stresses in corrosive environments is largely dictated by the interplay between the hydrogen (H) diffusion, the dislocation slip, and their interactions with the GBs. However, the slip-GB interaction in the presence of H and the subsequent cracking along the GBs remains not fully understood up to date. It is challenging to simultaneously resolve the collective dislocation motions away from the GBs together with the H diffusion near the atomically structured GBs using single-scale techniques. To address this challenge, taking bcc iron as a model material, here we present a a concurrent atomistic-continuum (CAC) computational analysis of the interaction between dislocation slip and a H-charged GB. With the large number of dislocations-mediated plastic flow away from the GB together with the atomistic structure evolution at the GB being simultaneously retained, several main findings from our simulations are: (i) the large number of dislocation-mediated pileup-induced internal stress nearby the H-charged GB can be as high as 3GPa; (ii) the more dislocations accumulated nearby the GB, the slower H diffusion ahead of the slip-GB intersection; (iii) H atoms diffuse fast behind the pileup tip, get trapped within the GB, and diffuse slow ahead of the pileup tip; (iv) the local stresses nearby the pileup tip exhibits a strong heterogeneity along the dislocation line direction. This differs from the common wisdom in many existing 2D theories/models. The buildup of such a high local stress heterogeneity leads to inhomogeneous H diffusion within the GB plane. The CAC simulation-predicted local H diffusivity, Dpileup−tip, and local stresses, σ, nearby the GB are then correlated with each other and consolidated into a mechanics model by considering the dislocation pileup as an Eshelby inclusion. These findings will provide researchers with opportunities to: (a) characterize the coupled dynamics between plasticity, H diffusion, and crack initiation underlying the hydrogen-induced cracking (HIC); (b) develop mechanism-based constitutive rules to be used in diffusion-plasticity coupling models for understanding the interplay between mechanical and mass transport in materials at the continuum level; and (c) connect the atomistic deformation physics of polycrystalline materials with their overall performance in aqueous environments, which is currently difficult to achieve in laboratory experiments.

Porous ferroelectric lead zirconate-titanate (PZT) films are a promising material for various electronics applications. This study focuses on understanding how structure-directing agent (polyvinylpyrrolidone) can change structure and electrical properties of porous PZT films prepared by chemical solution deposition. The films with various porosity up to ~ 40 vol.% and pore connectivity changing from 3-0 to 3-3 were prepared and studied by capacitance-voltage, dielectric hysteresis, transient current, photocurrent, and local current techniques. We have found that the linear decrease in the material volume in a porous film is not the only factor governing film properties. The key role in changing electrical properties plays the creation of new internal grain boundaries. This research expands the understanding of physical phenomena in porous ferroelectric films and may enable the development of new materials and devices.

The digital representation of various signals allows, at the subsequent stages of their transmission, to apply correction codes that provide protection against possible errors arising from the action of interference in the communication channel. At the same time, it is important that, with the required correcting ability, these codes have the maximum possible speed. The article presents the results of calculations for linear codes, showing their really achievable limiting capabilities.

The paper describes the Field Boundary Element Method applied to the fracture analysis of a 2D rectangular plate made of Functionally Graded Material to calculate Mode I Stress Intensity Factor. The object of the Field Boundary Element Method is the transversely isotropic plane plate. Its material presents an exponential variation of the elasticity tensor depending on a scalar function of position, i.e., the elastic tensor results from multiplying a scalar function by a constant taken as a reference. Several examples using a parametric representation of the structural response show the suitability of the method that constitutes a sight of Stress Intensity Factor evaluation of Functionally Graded Materials plane plates even in the case of more complex geometries.

The response of the Cold-Electron Bolometers (CEBs), integrated into a 2-D array of dipole antennas, has been measured by irradiation from YBa$_2$Cu$_3$O$_{7-\delta}$ (YBCO) 50 $\mu$m long Josephson junction into the THz region at frequencies from 0.1 to 0.8 THz. The possibility of controlling the amplitude-frequency characteristic is demonstrated by the external magnetic field in the traveling wave regime of a long Josephson junction. YBCO junction has been formed on the bicrystal Zr$_{1-x}$Y$_x$O$_2$ (YSZ) substrate by magnetron sputtering and etching of the film. CEBs have been fabricated using an Al multilayer structure by a self-aligned shadow evaporation technique on Si substrate. Both receiver and oscillator have been located inside the same cryostat at 0.3K and 2.7K plates, respectively.

A forecasting scheme of the thermal updraft velocity based on a theoretical model and data collected from flights records at gliding competitions, is presented. The forecasting scheme was based on the hypothesis that there is linear relationship between the overheat function at ground surface and the temperature difference between soil and air. The proportionality factor of this relationship was determined experimentally using observations recorded during gliding flights. The results showed that based on this simple scheme forecasting thermal convection is possible at any geographical location.

To gain insight into the differences between onshore and offshore atmospheric turbulence, 2 pressure fluctuations were measured for offshore wind under different environmental conditions. 3 A durable piezo-electric sensor was used to sample turbulent pressure data at 50 kHz. Offshore 4 measurements were performed at 100 m height on Germany’s FINO3 offshore platform in the 5 German Bight together with additional meteorological data provided by Deutscher Wetterdienst 6 (DWD). The statistical evaluation revealed that the stability state in the atmospheric boundary has a 7 large impact on turbulent fluctuations. Therefore, we used higher statistical properties (described 8 by so-called shape factors) to the stability state. Data was classified to be either within the unstable, 9 neutral or stable stratification. We found that in case of stable stratification, the shape factor is 10 mostly close to zero, indicating that a thermally stable environment produces closer-to Gaussian 11 distributions. Non-Gaussian distributions were found in unstable and neutral boundary layer states 12 and an occurrence probability was estimated. Possible impact on laminar-turbulent transition on the 13 blade is discussed with application of so-called laminar aerofoils on wind turbine blades. Use of a 14 cut-off frequency to separate load and aerodynamic turbulence is proposed.

Besides the Lagrangian and the Eulerian descriptions, the motion of a body can also be expressed relative to the present configuration of the body, known as the relative motion description. It is interesting to consider such a relative motion description in general to formulate the basic system of field equations for solid bodies. In doing so, when the time increment from the present state is small enough, the nonlinear constitutive equations can be linearized relative to the present state so that the resulting system becomes linear. This will be done for thermoelastic materials with a brief comment on the exploitation of entropy principle in general. Relative Lagrangian formulation is based on the well-known ``small-on-large'' idea, and can be implemented for solving problems with large deformation in successive incremental manner. Some applications of such a formulation in numerical simulations are briefly reviewed.

Keywords: gender equality policy; initial teacher education; teacher educators’ perceptions; gender perspective; university teaching; sustainable development; SDG5; SDG4.

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

Metal-organic frameworks (MOFs) represent the largest class of materials among the crystalline porous materials ever developed and have attracted attention as core materials for separation technology. Their extremely uniform pore aperture and nearly unlimited structural and chemical characteristics have attracted great interest and promise in applying MOFs to adsorptive and membrane–based separations. This paper reviews the recent research and development of MOF membranes for gas separation. Strategies for polycrystalline membranes and mixed matrix membranes are discussed, with a focus on separation systems involving hydrocarbon separation and CO2 capture. Challenges and opportunities for the industrial deployment of MOF membranes are also be discussed, providing guidance for the design and fabrication of future high-performance membranes. The contributions of the underlying mechanism to separation performances, and the adopted strategies and membrane processing technologies for breaking the selectivity/permeability trade–off are discussed.

This study focuses on the formulation and analysis of a numerical method to compute singularly perturbed functional differential equations (SP-FDEs) involving large mixed shifts. SP-FDEs are a class of differential equations with both shifts and small perturbation parameter. Such equations arise in various scientific fields to model phenomena with memory effects and perturbed dynamics. In this particular study, the considered SP-FDEs involve large negative shift (delay) and large positive shift (advance). Dealing with large mixed shifts in SP-FDEs gives rise to challenges in terms of stability analysis, numerical methods, and convergence of a solution. In this paper, the influence of the large mixed shifts is treated by choosing a uniform mesh discretized in such a way that the term with the shifts lie on the nodal points. Then, using the Numerov’s finite difference method, the continuous problem is transformed into a finite difference scheme and the influence of the perturbation parameter is handled by introducing a fitting factor in the difference scheme. Stability estimate and convergence analysis of the developed scheme are investigated and proved. Numerical experiments are carried out to demonstrate the validity and applicability of the developed scheme. It is shown that the numerical scheme obtained in this work is parameter-robust of second order convergence rate.

Transition delaying is of great importance for the drag and heat flux reduction of hypersonic flight vehicles. The first mode within low frequency and the second mode within high frequency exist simultaneously during the transition of hypersonic boundary layer. This paper proposes a novel bi-frequency synthetic jet to suppress low- and high-frequency disturbances at the same time. Orthogonal table and variance analysis are used to compare the control effects of jet with different frequencies, amplitudes and positions. Linear stability analysis results show that, low frequency synthetic jet can suppress the first mode when it is arranged upstream of synchronization point, while the second mode control effect is relatively weak. The higher the high frequency is, the stronger the suppression effect is on the first mode. For the second mode, the suppression effect is only at f2=89.09kHz. The larger the amplitude, the weaker the promoting effect for the first mode and the second mode, and the more obvious the suppressing effect. For the cases with synthetic jet downstream of synchronization point, all levels of the three parameters promote the unstable mode. In terms of the growth rate with the spanwise wave number, the control effect of the same factor and level under different spanwise wave number is different. In order to obtain the optimal control effect on transition, the three factors and the arrangement position of the synthetic jet should be selected as follows: the position is arranged in the upstream, with f1 = 17.82kHz, f2 = 89.9kHz, a =0.007, so that the maximum growth rate of the first mode is reduced by 9.06% and that of the second mode is reduced by 1.28% compared with the uncontrolled state.

In this work, we introduce a scalable and efficient GPU-accelerated methodology for volumetric particle advection and finite-time Lyapunov exponent (FTLE) calculation, focusing on the analysis of Lagrangian Coherent Structures (LCS) in large-scale Direct Numerical Simulation (DNS) datasets across incompressible, supersonic, and hypersonic flow regimes. LCS play a significant role in turbulent boundary layer analysis, and our proposed methodology offers valuable insights into their behavior in various flow conditions. Our novel owning-cell locator method enables efficient, constant-time cell search, and the algorithm draws inspiration from classical search algorithms and modern multi-level approaches in numerical linear algebra. The proposed method is implemented for both multi-core CPUs and Nvidia GPUs, demonstrating strong scaling up to 32,768 CPU cores and up to 62 Nvidia V100 GPUs. By decoupling particle advection from other problems, we achieve modularity and extensibility, resulting in consistent parallel efficiency across different architectures. Our methodology was applied to calculate and visualize the FTLE on four turbulent boundary layers at different Reynolds and Mach numbers, revealing that coherent structures grow more isotropic proportional to the Mach number, and their inclination angle varies along the streamwise direction. We also observed increased anisotropy and FTLE organization at lower Reynolds numbers, with structures retaining coherency along both spanwise and streamwise directions. Additionally, we demonstrated the impact of lower temporal frequency sampling by upscaling with an efficient linear upsampler, preserving general trends with only 10% of the required storage. In summary, we present a particle search scheme for particle advection workloads in the context of visualizing LCS via FTLE that exhibits strong scaling performance and efficiency at scale. Our proposed algorithm is applicable across various domains requiring efficient search algorithms in large structured domains. While this manuscript focuses on the methodology and its application to LCS, an in-depth study of the physics and compressibility effects in LCS candidates will be explored in a future publication.

This study provides, to our knowledge, the first detailed analysis of how surface oil modifies air-sea interaction in a two-way coupled model, i.e., the coupled-ocean-atmosphere-wave-sediment-transport (COAWST) model modified to account for oil-related changes to air-sea fluxes. This study investigates the effects of oil on surface roughness, surface wind, surface and near surface temperature differences, and boundary-layer stability, and how those conditions ultimately affect surface stress. We first conducted twin coupled modeling simulations with and without the influence of oil over the Deepwater Horizon (DWH) oil spill period (April 21 to May 5, 2010) in the Gulf of Mexico. Then we compared the results by using a modularized flux model with parameterizations selected to match those selected in the coupled model adapted to either ignore or account for different atmospheric/oceanic processes in calculating surface stress. When non-oil inputs to the bulk formula treated as unchanged by oil, the surface stress changes are always negative due to oil-related damping of surface roughness alone. However, the oil-related changes to 10-m wind speed and boundary-layer stability are found to play a dominant role in surface stress changes relative to those due to the oil-related surface roughness changes, highlighting that most of the changes in surface stress are due to the oil-related changes in wind speed and boundary-layer stability. Finally, the oil-related changes to surface stress due to the combined oil-related changes in surface roughness, surface wind, and boundary-layer stability are not large enough to have a major impact on the surface current and surface oil transport, indicating that the feedbacks from the surface oil to the surface oil movement itself are minor for forecasting the surface oil transport unless the fractional oil coverage is much larger than found in this study.

new subsonic blower wind tunnel design has been studied both numerically and experimentally, it is also referred to as “blower” wind tunnel. This paper, is initially aimed to address each sequential stage of the wind tunnel design process. Rather than applying the standard method of modelling solely the flow in the test section, a large-scale CFD model of the whole wind tunnel was employed. The loss of every constituent element was calculated and then all the losses are added up to determine the power needed for the wind tunnel operation which is used as “intake fan” boundary conditions in the CFD model. Then, flow uniformity and turbulent intensity measurements in an empty test section using a Pitot-static tube and hot wire anemometer (HWA) were introduced to validate the CFD results. The results showed that flow quality was significantly affected by boundary layer controllers (honeycomb and mesh screens) in the settling chamber and wide-angle diffuser. Investigations were also conducted to evaluate the flow deficit in the wake area behind a convex hump model using HWA. This was additional experimental tests carried out to validate the suitability of the wind tunnel flow aerodynamic research.

This paper considered the issue of agricultural fields boundary recognition in satellite images. A novel algorithm based on the aggregated history of vegetation index data obtained via open satellite data, Sentinel-2, was proposed. The proposed algorithm included several basic steps, namely the detection of parcel regions on aggregated index data; the calculation of aggregated edge maps; the segmentation of parcel regions using the edges obtained; the computation of connected components and their contour extraction. In this paper, we showed that the use of aggregated vegetation index data and boundary maps allow for much more accurate agricultural field segmentation compared to the instant vegetation index approach. The quality of segmentation within regions of Russia and the Ukraine was estimated. The dataset that was used and Python implementation of the proposed algorithm were provided.

We consider a family of nonlinear singularly perturbed PDEs whose coefficients involve a logarithmic dependence in time with confluent Fuchsian singularities that unfold an irregular singularity at the origin and rely on a single perturbation parameter. We exhibit two distinguished finite sets of holomorphic solutions, so-called outer and inner solutions, by means of a Laplace transform with special kernel and Fourier integral. We analyze the asymptotic expansions of these solutions relatively to the perturbation parameter and show that they are (at most) of Gevrey order 1 for the first set of solutions and of some Gevrey order that hinges on the unfolding of the irregular singularity for the second.

The design of thin-walled structures is commonly based on the solutions of linear boundary-valued problems, formulated within well-developed theories for elastic plates and shells. However, in modern appliances, especially in MEMS design, it is necessary to take into account non-linear mechanical effects that become decisive for flexible elements. Among substantial non-linear effects, that significantly change the deformation properties of thin plates, are the effects of residual stresses caused by incompatibility of deformations, which inevitably arise during the manufacture of ultra-thin elements. The development of new methods of mathematical modeling of residual stresses and incompatible finite deformations in plates is the subject of the paper. To this end the local unloading hypothesis is used. It makes possible to define smooth fields of local deformations (inverse implant field) for the mathematical formalization of incompatibility. The main outcomes are field equations, natural boundary conditions and conservation laws, derived from least action principle and variational symmetries with account of implant field. The derivations were carried out in the framework of elasticity theory for simple materials and, over and above, within the Cosserat’s theory of a two-dimensional continuum. As illustrative examples the distributions of incompatible deformations in a circular plate are considered.

In this review we collect some recent achievements in the accurate and efficient solution of the Nonlinear Schrödinger Equation (NLSE), with the preservation of its Hamiltonian structure. This is achieved by using the energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs) after a proper space semi-discretization. The main facts about HBVMs, along with their application for solving the given problem, are here recalled and explained in detail. In particular, their use as spectral methods in time, which allows efficiently solving the problems with spectral space-time accuracy.

The discovery of solitonlike coherent structure (SCS) in boundary layer flows is crucial for understanding turbulence origins and in particular for laminar-turbulence transition. However, the task of finding solutions for solitonlike coherent structure from the Navier-Stokes equations poses a significant challenge. In this paper, based on the author's previous work [B.H. Sun, Exact similarity solutions of unsteady laminar boundary layer flows, Preprints 2023, 2023092117], we are able to study convergent flow boundary layers, whose solution encompasses both shock wave and solitary wave solutions, and their superposition gives rise to solitary-like waves, namely solitonlike coherent structure. It is found that the solitonlike coherent structure can only be obtained by combining the Navier-Stokes equations and mass conservation, since the combined equations will have the third order derivatives respect to coordinate $x$, in other words, without the mass conservation condition, the Navier-Stokes equations does not contain solitonlike coherent structure by itself. Finally proved solitonlike coherent structure is a universal motion form of viscous fluid.

In this paper, existence and stability results of non-local integro-differential equation with Hilfer derivative of order ω∈(1,2) have been investigated. In this derivative is applied in RLC circuit equation of various fractional order. To proposed problem the existence solutions are derived using Schaefer’s fixed point theorem and the Banach contraction principle is using for uniqueness results. To demonstrate the effectiveness and applicability of our theoretical conclusions and two-step Lagrange polynomial interpolation were used to solve four numerical results.

This article addresses the problem of formation control of a quadrotor and one (or more) marine vehicles operating at the surface of the water with the end goal of encircling the boundary of a chemical spill, enabling such vehicles to carry and release chemical dispersants used during ocean cleanup missions to break-up oil molecules. Firstly, the mathematical models of the Medusa class of marine robots and quadrotor aircrafts are introduced, followed by the design of single vehicle motion controllers that allow these vehicles to follow a parameterized path individually using Lyapunov based techniques. At a second stage, a distributed controller using event triggered communications is introduced, enabling the vehicles to perform cooperative path following missions according to a pre-defined geometric formation. In the next step, a real time path planning algorithm is developed that makes use of a camera sensor, installed on-board the quadrotor. This sensor enables the detection in the image of which pixels encode parts of a chemical spill boundary and use them to generate and update in real time a set of smooth B-spline based paths for all the vehicles to follow cooperatively. The performance of the complete system is evaluated by resorting to 3-D simulation software, making it possible to simulate visually a chemical spill. Results from real water trials are also provided for parts of the system, where two Medusa vehicles are required to perform a static lawn-mowing path following mission cooperatively at the surface of the water.

In this work, we studied the Green’s functions of the second order differential operators with involution. Uniform equiconvergence of spectral expansions related to the second-order differential operators with involution is obtained. Basicity of eigenfunctions of the second-order differential operator operator with complex-valued coefficient is established.

The aim of this articles is to investigate the entropy optimization in unsteady MHD flow Darcy-Forchheimer nanofluids towards a stretchable sheet. The surface we tend to think about is porous and stretchy under acceleration. Flow occurs due to the stretching of the surface. Four distinct types of aqueous nanostructures are taken in this examination where copper oxide ( ), copper ( ), titanium dioxide ( ) and aluminum oxide ( ) are the nanoparticles. Irreversibility analysis are discussed through second law of thermodynamics. The expression of energy is mathematically designed and discussed according to heat generation / absorption, dissipation, thermal radiation, and joule heating. The nonlinear PDE (partial differential conditions) is first changed to ODE (normal differential conditions) through appropriate similarity variables. Here we used the numerically embedded solution technique to develop a numerical result for the obtained nonlinear flow expression. Influence of various flow parameter velocity temperature distribution and entropy generation are discussed. Reduction occurs in velocity profile for larger porosity and magnetic parameters. An enhancement in entropy generation and temperature distribution is seen for Brinkman number. An opposite effect is noticed in velocity and temperature through solid volume friction.

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

The primary cilium plays a pivotal role during embryonic development of vertebrates. It acts as a somatic signaling hub for specific pathways, such as sonic hedgehog signaling. In humans, mutations in genes that cause dysregulation of ciliogenesis or ciliary function lead to severe developmental disorders called ciliopathies. Beyond its obvious role in early morphogenesis, growing evidence points towards an essential function of the primary cilium in neural circuit formation in the central nervous system. However, very little is known about a potential role in the formation of the peripheral nervous system. Here, we investigated the presence of the primary cilium in neural crest cells and their derivatives in the trunk of the developing chicken embryo in vivo. We found that neural crest cells, sensory neurons, and boundary cap cells all bear a primary cilium during key stages of early peripheral nervous system formation. Moreover, we described differences in the ciliation of neuronal cultures of different populations from the peripheral and central nervous system. Our results offer a framework for further in vivo and in vitro investigations on specific roles that the primary cilium might play during peripheral nervous system formation.

The boundary layer height (BLH) distinguishes the interface between the lower and free atmosphere, which is a key variable in numerical simulation, aerosol, and environmental pollution studies. Thus, the goal of this work is to propose a novel method conjuncting with numerical regularization to analyze climate characteristics of the marine boundary layer height (MBLH) using 2007-2011 GPS-RO data from the COSMIC mission. While traditionally, the difference method has been used to achieve this aim, herein, we propose an innovative method, in which the bending angle profile gradient was calculated using the numerical regulation method where the regulation parameters are determined by the double-parameter model function method. Then, the MBLH was determined by employing the maximum gradient method to ascertain the height corresponding to the smallest gradient. The results show a correlation between currents and the MBLH—a relationship that has not been previously demonstrated. A low MBLH is associated with seasons and regions where cold ocean currents are prevailing; whereas a high MBLH is observed in the seasons and places where warm currents are prevailing. This correlation was validated by comparing the obtained results with different occultation data—i.e., atmprf and echprf—which also established that atmprf is more sensitive to convective cloud top capture. In seas with active convection, the armprf calculated bending angle is higher than that from echprf. Subsequently, the standard deviation was used to express the MBLH confidence level. The results show that the MBLH standard deviation is highest in low latitudes and lowest in the middle and high latitudes. Furthermore, we analyzed the interannual MBLH variation trend, which displayed a seasonal variation and spatial distribution corresponding with the current and subsolar point. Finally, we conducted a case study in the South China Sea, and identified a distinctive seasonal change and downward trend.

In the framework of 2D circular membrane Micro-Electric-Mechanical-Systems (MEMS), a new non-linear second-order differential model with singularity in the steady-state case is presented in this paper. In particular, starting from the fact that the electric field magnitude is locally proportional to the curvature of the membrane, the problem is formalized in terms of the mean curvature. Then, a result of existence of at least one solution is achieved. Finally, two different approaches prove that the uniqueness of the solutions is not ensured.

The three-dimensional single fin configuration finds application in an intake geometry where the cowl-shock wave interacts with the side-wall boundary-layer. Accurate numerical simulation of such three-dimensional shock/turbulent boundary-layer interaction flows, which are characterized by the appearance of strong crossflow separation, is a challenging task. Reynolds-averaged Navier-Stokes computations using the shock-unsteadiness modified Spalart-Allmaras model is carried out at Mach of 5 at large fin angle of 23^◦. The computed results using the modified model are compared to the standard Spalart-Allmaras model and validated against the experimental data. The focus of work is to implement the modified model and to study the flow physics in detail in the complex region of swept-shock-wave turbulent boundary-layer interaction in terms of the shock structure, expansion fan, shear layer and the surface streamlines. The flow structure is correlated to the wall pressure and skin friction in detail. It is observed that the standard model predicts an initial pressure location downstream of the experiments. The modified model reduces the eddy viscosity at the shock and predicts close to the experiments. Overall, the surface pressure using modified model is predicted accurately at all the locations. The skin friction is under predicted by both the models in the reattachment region and is attributed to the poor performance of turbulence models due to flow laminarization.

Determining precise boundaries and protective zones for historical and cultural objects enables their effective preservation. This is particularly crucial in the face of increasing tourism, urbanization, and natural disasters. Accurately defining boundaries and protective zones can aid in the planning of scientific research, archaeological excavations, and other studies aimed at exploring the history and cultural significance of these objects. The article presents the framework for establishing protective zones around historical and cultural heritage sites existing in Ukraine, using the example of the Citadel defensive complex located in Lviv (Ukraine). It proposes general and detailed conceptual models for the combined application of remote and non-invasive methods for investigating historical and cultural heritage sites. This conceptual article introduces the theory of integrating radar satellite imaging with ground-based georadar imaging. Additionally, it presents a software module that has been developed to analyze collected data on two-dimensional historical and cultural heritage objects, refine their boundaries, and establish protective zones around them.

Fluid mechanics of flow in hydrophobic, rectangular microchannels with finite aspect ratios is of paramount importance. In such microchannels, not only the effect of the side walls should be taken into account, but also the classical assumption of no-slip boundary condition (BC) is no longer valid at the solid-liquid interface. Accordingly, slip flow can occur in microchannels fabricated from surfaces with low wetting conditions, hydrophobic surfaces. Determining the interactions of liquid molecules adjacent to solid surface is still a challenging issue, and it is especially important in small scale domains. Herein, the fluid mechanics of flow through rectangular hydrophobic microchannels has been reconsidered by taking into account the general Navier-slip BCs at the solid-liquid interface. For fully developed incompressible flow in microchannels at low Reynolds number, partial differential equation (PDE) of the momentum equation simplifies to the classical Poisson equation. Accordingly, by analytically solving the Poisson equations with general Navier-slip BCs, the most general forms of velocity distributions, flow rate, friction factor and Poiseuille number have been obtained.

In this study, the multi-relaxation-time lattice Boltzmann method is applied to investigate the oscillatory instability of lid-driven flows in two-dimensional semi-elliptical cavities with different vertical to horizontal aspect ratios K in the range of 1.0--3.0. The program implemented in this study is parallelized using CUDA (compute unified device architecture), a parallel computing platform, and computations are carried out on NVIDIA Tesla K40c GPU. To carry out precise calculations, the CUDA algorithm is extensively investigated, and its parallel efficiency indicates that the maximum speedup is 47.6 times faster. Furthermore, the steady--oscillatory Reynolds numbers are predicted by implementing the CUDA-based programs. The amplitude coefficient is defined to quantify the time-dependent oscillation of the velocity magnitude at the monitoring point. The simulation results indicate that the transition Reynolds numbers correlate negatively with the aspect ratio of the semi-elliptical cavity, and are smaller than those of the rectangular cavity at the same aspect ratio. In addition, the detailed vortex structures of the semi-elliptical cavity within a single period are also investigated when the Reynolds number is larger than the steady--oscillatory value to determine the effects of periodic oscillation of the velocity magnitude.

A number of experimental and numerical studies point out that incorporating a rotating cylinder can superiorly enhance the aerofoil performance, especially for higher velocity ratios. Yet, there have been less or no studies exploring the effects of lower velocity ratio at a higher Reynolds number. In the present study, we investigated the effects of Moving Surface Boundary-layer Control (MSBC) at lower velocity ratios (i.e. cylinder tangential velocity to free stream velocity) and higher Reynolds number on a symmetric aerofoil (e.g. NACA 0021) and an asymmetric aerofoil (e.g. NACA 23018). In particular, the aerodynamic performance with and without rotating cylinder at the leading edge of the NACA 0021 and NACA 23018 aerofoil was studied on the wind tunnel installed at Aerodynamics Laboratory. The aerofoil section was tested in the low subsonic wind tunnel, and the lift coefficient and the drag coefficient were studied for different angles of attack. The experiments were conducted for two Reynolds numbers: 200000 and 250000 corresponding to two free stream velocities: 20 m/s and 25 m/s, respectively, for six different angle of attacks (-5°, 0°, 5°, 10°, 15° and 20°). This study demonstrates that the incorporation of a leading edge rotating cylinder results in an increase of lift coefficient at lower angle of attacks (maximum around 33%) and delay in stall angle (from 10° to 15°) relative to the aerofoil without rotating cylinder.

Global wind observations are fundamental for studying weather and climate dynamics. Most wind measurements come from atmospheric motion vectors (AMVs) by tracking the displacement of cloud or water vapor features. These AMVs generally rely on thermal infrared (IR) techniques for their height assignments, which are subject to large uncertainties in the presence of weak or reversed vertical temperature gradients around the planetary boundary layer (PBL) and with tropopause folding. Stereo imaging can overcome the height assignment problem using geometric parallax for feature height determination. In this study we develop a stereo 3D-Wind algorithm to simultaneously retrieve AMV and height from geostationary (GEO) and low Earth orbit (LEO) satellite imagery and apply it to collocated Geostationary Operational Environmental Satellite (GOES) and Multi-angle Imaging SpectroRadiometer (MISR) imagery. The new algorithm improves AMV and height relative to products from GOES or MISR alone, with an estimated accuracy of <0.5 m/s in AMV and <200 m in height with 2.2 km sampling. The algorithm can be generalized to other LEO-GEO or GEO-GEO combinations for greater spatiotemporal coverage. The technique demonstrated with MISR and GOES has important implications for future high-quality AMV observations, for which a low-cost constellation of CubeSats can play a vital role.

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.

Motivated by diffraction experiments on the $\left(2\sqrt{3}\times \sqrt{3}\right)\mathrm{R}{30}^{\circ }$ reconstructed Si(111) due to deposition of rare earth elements (Dy, Tb) and silicide formation we analyse the splitting and non-splitting of superstructure spots. For this purpose, we model diffraction patterns for one dimensional structures generated by the binary surface technique and use supercell models to keep the analysis simple. Diffraction pattern are calculated in the framework of the kinematical diffraction theory and they are analyzed as a function of the domains and domain boundaries. Basic properties of the diffraction pattern are analyzed for model systems of a two-fold and a three-fold periodicity. The rules derived from these calculations are applied to the "real-world" system of Si(111)-($\left(2\sqrt{3}\times \sqrt{3}\right)\mathrm{R}{30}^{\circ }$)-RESi$2_{}$ (RE = Dy or Tb). Depending on the combination of domains and domain boundaries of different types a plethora of different features are observed in the diffraction patterns. These are analyzed to determine the sizes of both domain boundaries and domains from experimentally observed splitting of specific superstructure spots.

For a linear ordinary differential equation (ODE in short) of the third order, results are presented that supplement the theory of conjugate points and extremal solutions by W. Leighton, Z. Nehari, M. Hanan. It is especially noted the sensitivity of solutions to the initial data, which makes their numerical study difficult. Similar results were obtained for the third-order nonlinear equations of the Emden-Fowler type.

J/ψ, a charmonium bound state made of a charm and an anti-charm quark, has been discovered in the 1970s and confirmed the quark model. Because the mass of charms quarks is significantly above the Quantum Chromodynamics (QCD) scale ΛQCD, charmonia are considered as excellent probes to test perturbative Quantum Chromodynamics (pQCD) calculations. Over the past decades, they have been studied extensively at different high energy colliders. However, their production mechanisms, which involves multiple scales, are still not very well understood. Recently, in high multiplicity p+p collisions at RHIC and at the LHC, a significant enhancement of J/ψ production yield has been observed, which suggests a strong contribution from the Multi-Parton Interaction (MPI). This is different from the traditional pQCD picture where charm quark pairs are produced from a single hard scattering between partons in p+p collisions. In this work, we will report the J/ψ normalized production yield as a function of normalized charged particle multiplicity over a board range of rapidity and event multiplicity in the J/ψ→μ+μ− channel with PHENIX Run 15 p+p data at \sqrt{s}=200 GeV. The results are compared with PYTHIA 8 simulations with the MPI option on and off. Finally, the outlooks of J/ψ in p+Au and Au+p collisions along with Color Glass Condensate (CGC) predictions and of multiplicity dependent ψ(2S)/J/ψ ratio in p+p data will be briefly discussed.

In this paper, a multi-patch weakly singular isogeometric boundary element method (WSIGABEM) for two-dimensional elastostatics is proposed. Since the method based on the weakly singular boundary integral equation, the quadrature techniques, dedicated for the weakly singular and regular integrals, are applied in the method. A new formula for the generation of collocation points is suggested to take the full advantage of the multi-patch technique. The generated collo-cation points are essentially inside the patches without any correction. If the boundary conditions are assumed to be continuous in every patch, no collocation point lies on the discontinuous boundaries, thus simplifying the implementation. The multi-patch WSIGABEM is verified by the simple examples with analytical solutions. The feature of the present multi-patch WSIGABEM is investigated by the comparison with the traditional IGABEM. Furthermore, the combination of the present multi-patch WSIGABEM and the particle swarm optimization algorithm results in a shape optimization method in two-dimensional elastostatics. By changing some specific control points and their weights, the shape optimizations of the fillet corner, the spanner and the arch bridge are verified to be effective.

This study determines the equivalent stress intensity factor (SIF) model that best fits the experimental behavior of low-carbon steel under mixed mode (I and II). The study assessed Tanaka, Richard, and Pook's equivalent SIF models. The theoretical values used for comparison correspond to experimental results in a modified C(T) geometry by machining a hole ahead of the crack tip subjected to fatigue loads with a load ratio of R = 0.1. The comparison involved the SIF for six experimental points and the values computed through the numerical simulation. The Paris, Klesnil, and modified Forman-Newman crack growth models were used with each equivalent SIF to analyze the prediction in the estimated number of cycles. The Klesnil model showed the closest prediction since the error between the calculated and experimentally recorded number of cycles is the lowest. However, the material behavior reflects a reduced crack propagation rate attributed to plasticity in the crack tip. Results suggest that Asaro equivalent SIF conservatively estimates the element lifespan with increasing errors from 2.3% at the start of growth to 27 % at the end of the calculation.

In order to reveal the impact of boundary-spanning activities of cooperative innovation teams on team innovation performance, this paper takes the panel data of 71 cooperative innovation teams from January to November 2022 as the research sample. It introduces intermediary variables (teamwork crafting and individual work crafting) to analyze the impact mechanism of boundary-spanning activities of teams on innovation performance, 71 teams were divided into 41 experimental groups and 30 control groups, and a quasi-natural experiment was conducted on the innovation performance of team boundary-spanning activities using the Double Difference Model (DID).Research has shown that boundary-spanning activities of collaborative innovation teams can promote team innovation performance. Team job crafting has a mediating effect on team innovation performance in boundary-spanning activities of collaborative innovation teams. Team job crafting and individual job crafting mediate between the boundary-spanning activities of collaborative innovation teams and team innovation performance. Further analysis using the double difference model found that compared to teams without boundary-spanning activities, teams with boundary-spanning activities can directly improve team innovation performance. When team reflection is vital, and task interdependence is high, it will promote team innovation performance. This research enriches the research on the effects of boundary-spanning activities of collaborative innovation teams, explores solutions based on quasi-nature, and provides a reference for improving the team innovation performance of collaborative innovation teams.

For the KdV-Burgers equations on cylindrical and spherical waves the development of a regular profile starting from an equilibrium under a periodic perturbation at the boundary is studied. The equations describe a medium which is both dissipative and dispersive. Symmetries, invariant solutions and conservation laws are investigated. For an appropriate combination of dispersion and dissipation the asymptotic profile looks like a periodical chain of shock fronts with a decreasing amplitude (sawtooth waves). The development of such a profile is preceded by a head shock of a constant height and equal velocity which depends on spatial dimension as well as on integral characteristics of boundary condition; an explicit asymptotic for this head shock and a median of the oscillating part is found.

The use of technologies and the Internet poses problems and risks related to digital security. This article presents the results of a study on the evaluation of the digital competence of future teachers in the DigCompEdu European framework. 317 undergraduate students from Spain and Portugal answered a questionnaire with 59 items, validated by experts, in order to assess the level and predominant competence profile in initial training (including knowledge, uses and interactions and attitudinal patterns). The results show that 47% of the participants belong to the profile of teachers at medium digital risk, evidencing habitual practices that involve risks such as sharing information and digital content inappropriately, not using strong passwords, and ignoring concepts such as identity, digital “footprint” and digital reputation. The average valuations of each item in the seven categories show that future teachers have an average competence in the area of digital security. They have good attitudes toward security but less knowledge and fewer skills and practices related to the safe and responsible use of the Internet. Future lines of work are proposed, aimed at responding to the demand for a better prepared and more digitally competent citizenry. The demand for education in security, privacy and digital identity is becoming increasingly important, and these elements form an essential part of initial training.

We show how the Poincar\'e--Riemann--Hilbert boundary-value problem enables us to construct effective estimates of the potential in the Schr\"odinger equation. The apparatus of the three-dimensional inverse problem of quantum scattering theory is developed for this. It is shown that the unitary scattering operator can be studied as a solution of the Poincar\'e--Riemann--Hilbert boundary-value problem. This allows us to go on to study the potential in the Schr\"odinger equation

In this paper we address the approximation of the coupling problem for the wave equation and Maxwell's equations of electromagnetism in the time domain in terms of the electric field, by means of a nodal linear finite-element discretization in space, combined with a classical'explicit finite-difference scheme for the time discretization. Our study applies to the particular case where the dielectric permittivity has a constant value outside a sub-domain, whose closure does not intersect the boundary of the domain where the problem is defined. Inside this sub-domain Maxwell's equations hold. Outside this sub-domain the wave equation holds, which may correspond to the Maxwell's equations with a constant permittivity under certain conditions. We consider as a model the case of first order absorbing boundary conditions. Optimal error estimates that hold in natural norms under reasonable assumptions are given, among which lies a typical CFL condition for hyperbolic equations.

The purpose of the research is to investigate the unique solvability of a nonlocal boundary value problem for a loaded equation of parabolic-hyperbolic type in a special domain. Methods: Using representations of a general regular solution, the existence and uniqueness of the problem posed is proved. Results: A new method for proving the unique solvability of nonlocal boundary value problems for one loaded equation of mixed type in a special domain has been formed and developed. Conclusions: The obtained results make it possible to solve some problems of genetics, immunology, transonic gas dynamics and thermal physics, which are associated with displacement.

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

This paper deals with a non-standard finite difference scheme defined on a quasi-uniform mesh for approximate solutions of the Magneto-Hydro Dynamics (MHD) boundary layer flow of an incompressible fluid past a flat plate for a wide range of the magnetic parameter. We show how to improve the obtained numerical results via a mesh refinement and a Richardson extrapolation. The obtained numerical results are favourably compared with those available in the literature.

Since crystals are made of periodic structures in space, predicting their three period vectors starting from any values based on the inside interactions is a basic theoretical physics problem. For the general situation where crystals are under constant external stress, we derived dynamical equations of the period vectors in the framework of Newtonian dynamics, for pair potentials recently (doi:/10.1139/cjp-2014-0518). The derived dynamical equations show that the period vectors are driven by the imbalance between the internal and external stresses. This presents a physical process where when the external stress changes, the crystal structure changes accordingly, since the original internal stress can not balance the external stress. The internal stress has both a full kinetic energy term and a full interaction term. It is extended to many-body interactions in this paper. As a result, all conclusions in the pair-potential case also apply for many-body potentials.

The nucleation and dynamics of Josephson and Abrikosov vortices determine the critical currents of layered high-T_c superconducting (HTS) thin films, grain boundaries, and coated conductors, so understanding their mechanisms is of crucial importance. Here we treat pair creation of Josephson and Abrikosov vortices in layered superconductors as a secondary Josephson effect, in which each full vortex is viewed as a composite fluid of micro-vortices, such as pancake vortices, which tunnel coherently via a tunneling matrix element. We introduce a two-terminal magnetic (Weber) blockade effect that blocks tunneling below a threshold current, and simulate time-correlated vortex tunneling above threshold. The model shows nearly precise agreement with voltage-current (V-I) characteristics of HTS cuprate grain boundary junctions, which becomes more concave rounded as temperature decreases, and also explains the piecewise linear V-I behavior observed in iron-pnictide bicrystal junctions and other HTS devices. When applied to either Abrikosov or Josephson pair creation, the model explains a plateau seen in plots of critical current vs. thickness of HTS coated conductors. The observed correlation between theory and experiment strongly supports the proposed quantum picture of vortex nucleation and dynamics in layered superconductors.

Different types of resonators are used to create acoustic metameterials and metasurfaces. Recent studies focused on the use of multiple resonators of dipole, quadrupole, octupole and even hexadecapole types. This paper considers the theory of an acoustic metasurface, which is a flat surface with a periodic arrangement of multiple resonators. The sound field reflected by the metasurface is determined. If the distance between the resonators is less than a half the wavelength of the incident plane wave, the far field can be described by a reflection coefficient that depends on the angle of incidence. This allows to characterize the acoustic properties of the metasurface by a homogenized boundary condition, which is a high order tangential impedance boundary condition. The tangential impedance depending on the multiple order of the resonators is introduced. In addition, we analyze sound absorption properties of these metasurfaces, which are a critical factor in determining their performance. The paper presents a theoretical model that accounts for the multiple orders of resonators and their impact on sound absorption. The maximum absorption coefficient for a diffuse sound field, as well as the optimal value for the homogenized impedance, are calculated for arbitrary multipole orders.

Physical processes play important roles in controlling eutrophication and oligotrophication. In stratified lakes, internal waves can cause vertical transport of heat and nutrients without breaking the stratification, through boundary mixing events. Such is the case in tropical Valle de Bravo (VB) lake, where strong diurnal winds drive internal waves, boundary mixing and hypolimnetic warming during stratification periods. We monitored VB during 18 years (2001-2018) when important water-level fluctuations (WLF) occurred, affecting mixing and nutrient flux. Mean hypolimnetic temperature increase (0.06–1.04°C month-1) occurred in all the stratifications monitored. We analyzed temperature distributions and modeled the hypolimnion heat budget to assess vertical mixing between layers (26,618–140,526 m-3h-1), vertical diffusivity coefficient KZ (6.2x10-7–3.3x10-6 m2s-1) and vertical nutrient entrainment to epilimnion on monthly scale. Stability also varied as a function of WLF. Nutrient flux to the epilimnion ranged 0.36–5.99 mg m-2d-1 for soluble reactive phosphorus (SRP) and 5.8–97.1 mg m-2d-1 for dissolved inorganic nitrogen (DIN). During low water-level years, vertical nutrient fluxes increase and can account for up to &gt;40% of the total external nutrients load to the lake. Vertical mixing changes related to WLF affect nutrient recycling, their flux to sediments, ecosystemic metabolic balance and planktonic composition of VB.