The exact similarity solutions (also called as special exact solutions) of two dimensional laminar boundary layer were obtained by Blasius in 1908, however, no similarity solutions for two dimensional turbulent boundary layers have been reported in the literature. With the help of dimensional analysis and invariance principle, Prandtl mixing length $\ell=\kappa y$ for one dimensional turbulent boundary layer is extended to $\ell(x,y)=\kappa y (1-\frac{y}{\delta})\sqrt{\frac{\nu}{U\delta}} $ for the two dimensional turbulent boundary layers, furthermore, with a similarity transformation, we successfully transform the two dimensional turbulent boundary layers partial differential equations into a single ordinary differential equation $f'''+ ff''+\beta(1-f'^2)+\kappa^2[\eta^2(1-\eta)^2f''|f''|]'=0$. As an application, similarity solutions of the two dimensional turbulent boundary layer on a flat plate at zero incidence have been studied in detail. To solve the ordinary differential equation numerically, a complete Maple code is provided.

The exact similarity solutions of two dimensional laminar boundary layer were obtained by Blasius in 1908, however, for two dimensional turbulent boundary layers, no similarity solutions (special exact solutions) have ever been found. In the light of Blasius' pioneer works, we extend Blasius similarity transformation to the two dimensional turbulent boundary layers, and successfully transform the two dimensional turbulent boundary layers partial differential equations into a single ordinary differential equation. By author's Maple code, we numerically solve the ordinary differential equation and produce some useful quantities.

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.

The identification 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 [Sun (2023) in Ref.31], 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 do exit in all kind of flows.

Quinacridone (QA) and 2,9-dimethylquinacridone (DQA) are synthetic substances suitable as a hard, abrasion-resistant, self-organizing gliding layer on ice, snow, and water. For sustainable use, a large number of parameters must be considered to demonstrate that these non-biogenic substances and their by-products and degradation products are harmless to humans and the environment in the quantities released. For this task, available experimental data are used and supplemented for all tautomers by numerous relevant physical, chemical, toxicological and ecotoxicological estimated values based on various Quantitative Structure Activity Relationship (QSAR) methods. On the one hand, the low solubility of QA and DQA leads to stable gliding layers and thus low abrasion and uptake by plants, animals, and humans. On the other hand, the four hydrogen bond forming functional groups per molecule allow nanoparticle decomposition and enzymatic degradation in natural environments. All available data justify a sustainable use of QA and DQA as a gliding layer. The assessment of the toxicological properties is complemented with an investigation of the size and morphology of DQA particles as well as field tests indicating an excellent performance as a gliding layer on snow.

Panel data of our interest consist of a moderate number of panels, while the panels contain a small number of observations. An estimator of common breaks in panel means without a boundary issue for this kind of scenario is proposed. In particular, the novel estimator is able to detect a common break point even when the change happens immediately after the first time point or just before the last observation period. Another advantage of the elaborated change point estimator is that it results in the last observation in situations with no structural breaks. The consistency of the change point estimator in panel data is established. The results are illustrated through a simulation study. As a by-product of the developed estimation technique, a theoretical utilization for correlation structure estimation, hypothesis testing, and bootstrapping in panel data is demonstrated. A practical application to non-life insurance is presented as well.

The shear stress generated by wall turbulence is the main cause of wall friction resistance in turbulent flow through pipes. This paper investigates the impact of inserting rest periods (regions of constant Reynolds number) within the pulsating operating cycle of velocity on the fully turbulent flow at large time-averaged Reynolds numbers, using the Large Eddy Simulation (LES) method. The study aims to explore the effect of increasing rest periods within pulsations on drag resistance. The dimensionless shear stress and drag reduction rates during different time periods of rest were analyzed. Numerical simulation results indicate that the pulsating velocity operation mode does not necessarily lead to drag reduction; it may even result in increased resistance. Inserting rest periods within the pulsation cycle can achieve drag reduction effects, with the maximum drag reduction rate reaching 21.8\%. Comparisons with experimental data and Direct Numerical Simulation (DNS) from the literature validate the feasibility of using the LES method for pipe pulsating operation modes.

NongXiangXing Baijiu (NBJ) is produced using a solid-state fermentation system, but the impact of stratified fermented grains on fermentation remains unclear. Therefore, in a typical distillery plant, we evaluated the physical and chemical composition, volatile compound profile, and microbial community of fermented grains in NBJ both above (FG-A) and below (FG-B) the Huangshui line. Significant differences in fermentation parameters between FG-A and FG-B were observed after 30 days of fermentation (P < 0.05). Additionally, the partial least squares projection to latent structure discriminant analysis (PLS-DA) revealed distinct differences in volatile compounds between FG-B and FG-A, identifying 26 discriminant markers. The diversity of short-chain fatty acids (SCFAs) and their esters was higher in FG-B compared to FG-A. Furthermore, microbial diversity and abundance differed significantly between the two layers of fermented grains (P < 0.05). Correlation and pathway enrichment analyses indicated that the higher SCFA content in FG-B could be attributed to the greater abundance of acid-producing microorganisms compared to FG-A. This study highlights the differences between the two layers of fermented grains in NBJ fermentation, offering new insights into solid-state fermentation and expanding the current understanding of the traditional NBJ fermentation process.

Nowadays, rotor-stator mixers are widely used in the chemical process industry to carry out the dispersion and de-agglomeration processes for manufacturing the homogenous electrode slurries. Development of the process simulation models for the turbulent mixing of these slurries helps to optimize the mixing parameters such as rotational speed and dispersion time. In this work, computational fluid dynamics (CFD) approach through ANSYS Fluent are explored to determine and analyze the flow patterns, velocity fields and viscosity distribution during the dispersion of anode slurries. It utilizes a non-homogeneous Eulerian multiphase approach and the typical k-ε model to create turbulence fields and compute the mixing rate. This work significantly improves our understanding of the relationship between rotor speed, dispersion time, viscosity distribution, and mixer geometry in rotor-stator mixers.

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.

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.

As waterflooding multi-layer reservoirs get into the high-water-cut stage, inter-layer conflicts become increasingly serious, leading to a worsening development effect over time. Production layers regroup is an effective approach for resolving inter-layer conflicts and improving waterflooding efficiency. At the current stage, there are limitations for most of the methods of production layers regroup. This article proposes a smart method for optimizing the layer regroup plan based on genetic algorithm. Comprehensively considering various factors that affect the regroup of layers, based on the combination principle of "smaller intra-group variance and larger inter-group variance of each influencing factor are expected ", genetic algorithm is used to calculate the fitness value of the initial combination schemes, the advantageous schemes with higher fitness values are selected as the basis of the next generation. Then crossover and mutation operations are performed to those advantageous schemes to generate new schemes. Through continuous selection and evolution, until the global optimal solution with the highest fitness value is found, the optimal combination scheme is determined. Comparative analysis with numerical simulation results demonstrates the reliability of this intelligent method, with an increased oil recovery of 4.34% for the sample reservoir. Unlike selecting a preferable plan from a limited number of predefined combination schemes, this method is an automatic optimization to solve the optimal solution of the problem. It improves both efficiency and accuracy as compared to conventional reservoir engineering methods, numerical simulation methods and most of the mathematical methods, thus could provide effective guidance for EOR strategies of waterflooding reservoirs in high-water-cut stage.

Biomaterials are the basis for the development of medicine because they allow safe contact with a living organism. The aim of this work was to produce innovative oxide layers with a porous microstructure on the surface of commercially pure titanium Grade 4 (CpTi G4), and to characterize their properties as drug carriers. Anodization of the CpTi G4 subjected to mechanical grinding and electrochemical polishing was carried out in a solution of 1M ethylene glycol with the addition of 40 g of ammonium fluoride at a voltage of 20 V for 2, 18, 24 and 48 hours at room temperature. It was found that the longer anodization time, the greater the number of pores formed on the CpTi G4 surface as revealed by FE-SEM method, and the greater the surface roughness determined in profilometric tests. As the anodizing time increases, the amount of the drug in the form of gentamicin sulfate incorporated into the resulting pores decreases. The most favorable drug release kinetics profile determined by UV-VIS absorption spectroscopy was found for the CpTi G4 anodized for 2 h making them the best choice when it comes to the drug carrier allowing full use of the therapeutic dose of the drug.

High mountain rainforests are vital in the global energy and carbon cycle. Understanding the exchange of energy and carbon plays an important role to reflect responses to climate change. In this study, an eddy-covariance (EC) measurement system installed in the high Andean Mountains of South Ecuador was used. As EC measurements are affected by heterogeneous topography and vegetation height, the main objective was to estimate the effect of the sloped terrain and the forest on the turbulent energy and carbon fluxes considering the energy balance closure (EBC) and the heat storage. The results showed that the performance of the EBC was generally good and estimated to be 79.5%. This could be improved when the heat storage effect was considered. Based on the variability of the residuals in the diel modifications in the imbalances were highlighted. Particularly, during daytime residuals were largest (56.9 W/m² on average) with a clear overestimation. At nighttime mean imbalances were rather weak (6.5 W/m²) and mostly positive, while strongest underestimations developed in the transition period to morning hours (down to -100 W/m²). With respect to the Monin-Obukhov stability parameter ((z-d)/L) and the friction velocity (u*) it was revealed that largest overestimations evolved in weak unstable and very stable conditions associated with large u* values. In contrast, underestimation was related to very unstable conditions. The estimated carbon fluxes were independently modelled with a non-linear regression using a light-response relationship and reached a good performance (R²=0.52). All fluxes were additionally examined in the annual course to estimate whether both, energy and carbon fluxes, resembled the microclimatological conditions of the study site. This unique study demonstrated that EC measurements provide valuable insights into land surface -atmosphere interactions and contribute to our understanding of energy and carbon exchanges. Moreover, the flux data provide an important basis to validate coupled atmosphere ecosystem models.

Delaying the onset of laminar-turbulent transition is an attractive method in reducing skin friction drag, especially on streamlined bodies where Tollmien-Schlichting instabilities are the dominating mechanism for transition. Miniature vortex generators (MVGs) offer an effective approach to attenuate these instabilities by generating counter-rotating vortex pairs. They are placed in pairs within an array and resemble small winglet-type elements. The conventional methodology involves adjusting MVG parameters and conducting computationally expensive DNS and/or downstream stability analyses to assess their effectiveness. However, analyzing the vortex parameters of MVG-generated vortices can potentially guide a more targeted approach in modifying MVG parameters and identifying critical factors for transition delay. Therefore, this study investigates changes in three primary MVG parameters, namely inner distance, periodicity and height, and utilizes computational fluid dynamics (CFD) analysis to create a dataset that examines the characteristics of the generated counter-rotating vortex pairs and their potential in drag reduction. The objective is to establish correlations among these parameters and their influence on delaying transition.

We primarily investigate the existence of solutions for fractional neutral integro-differential equations subjected to Neumann-type boundary conditions, which is crucial for understanding natural phenomena. Taking into account factors such as neutral type, fractional-order integrals, and fractional-order derivatives, we employ probability density functions, Laplace transforms, and resolvent operators to formulate a well-defined concept of a mild solution for the specified equation. Following this, by integrating fixed point theorems, we establish the existence of mild solutions under more relaxed conditions.

The dynamics of the topological charge of a vortex optical beam propagating in a turbulent air with accounting for the cubic nonlinearity is theoretically considered. On a number of examples, we show that the optical beam self-focusing manifests itself ambiguously depending on optical wave power. At near-critical values of beam power, self-focusing leads to enhanced spatial localization of optical vortices and substantial suppression of vortex walk-off relative to the beam axis caused by air turbulence. However, with increasing of optical intensity the modulation instability imposed by cubic nonlinearity becomes significant and contributes jointly with medium turbulence to faster divergence of vortex beams.

A single mode waveguide grating coupler based on multiple Si-SiO₂ pairs onto Si substrate has been designed. Numerical analysis has been carried out to calculate optimum thickness of the layers of Si-SiO₂ that ensures the constructive interference between reflected waves and actual guided wave for high coupling efficiency. Based on the results, an optimal design is developed and modeled by using a 2-D finite difference time domain (2-D FDTD) simulator that dictates a coupling efficiency of as much as 78% (-1.07 dB) at the wavelength of 1550 nm, and a 1-dB bandwidth of 77 nm. The numerical method will be useful to calculate the optimum thicknesses of the layers for any reflector based grating coupler of different materials.

This review is a first attempt at bringing together various concepts from research on wall- and magnetically-bounded turbulent flows. Brief reviews of both fields are provided: The main similarities identified are coherent (turbulent) structures, flow generation and transport barriers. Examples are provided and discussed.

Electrochromic measurements have been performed throughout the full composition range of reactive magnetron-sputtered mixed Titanium oxide and Molybdenum oxide (TiO2-MoO3) to search for enhanced Coloration Efficiency (CE). Optical parameters and composition have been determined and mapped by using Spectroscopic Ellipsometry (SE). Scanning Electron Microscopy (SEM) with Energy-Dispersive X-ray Spectroscopy (EDS). EDS has been used to check the SE results as it was used in our earlier paper [N. T. Ismaeel et al, Materials 2023, 16(12), 4204] for (TiO2-SnO2). Ti and Mo targets were put separately from each other, and the Indium-Tin-Oxide (ITO) covered glass and Si-probes on a glass substrate (30 cm × 30 cm) were moved under the two separated targets (Ti and Mo) in a reactive Argon-Oxygen (Ar-O2) gas mixture. By using this combinatorial process, all the compositions (from 0 to 100%) were achieved in the same sputtering chamber after one sputtering process. The CE (the change of light transmission for the unit electric charge) of mixed metal oxides (TiO2-MoO3) that are deposited by reactive sputtering has been studied in an electrochemical cell. This paper aims to assess the results of investigations of such materials showing the enhanced electrochromic behavior compared to the pure materials and to enhance the Coloration Efficiency.

A Cu₂ZnSnS₄ (CZTS) film with a thickness of approximately 1.5 μm was fabricated on a Mo-coated glass substrate by annealing a CZTS precursor fabricated from nanoparticle ink. The chemical states of the elements in the CZTS thin film in the depth direction were studied to identify the presence of secondary phases, which are detrimental to the performance of solar cells containing CZTS. X-ray diffraction was unable to detect any secondary phases in CZTS because of their small relative amount. Instead, X-ray photoelectron spectroscopy (XPS), which is highly sensitive to chemical state, was conducted at different depths in the CZTS film to further check the presence of secondary phases. XPS analysis revealed peaks shift consistent with the presence of secondary phases. For the CZTS film annealed in a S atmosphere at 575 °C for 3 h, the film surface consisted of a secondary-phase layer composed of CuS, ZnS, and SnSx (x=1 or 2) originating from the decomposition of CZTS. At depths below 80 nm, the film was pure CZTS. Formation of MoS2 at the CZTS–Mo interface was confirmed by XPS analysis of Mo and S.

Quantum turbulence deals with the phenomenon of turbulence in quantum fluids, such as superfluid helium and trapped Bose-Einstein condensates (BECs). Although much progress has been made in understanding quantum turbulence, several fundamental questions remain to be answered. In this work, we investigated the entropy of a trapped BEC in several regimes, including equilibrium, small excitations, the onset of turbulence, and a turbulent state. We considered the time evolution when the system is perturbed and let to evolve after the external excitation is turned off. We derived an expression for the entropy consistent with the accessible experimental data, that is, using the assumption that the momentum distribution is well-known. We related the excitation amplitude to different stages of the perturbed system, and we found distinct features of the entropy in each of them. In particular, we observed a sudden increase in the entropy following the establishment of a particle cascade. We argue that entropy and related quantities can be used to investigate and characterize quantum turbulence.

This paper attempts to clarify an issue regarding the lasting unsolved problem of turbulence, namely the closure problem. This study shows that all perspectives regarding the numbers of un- known quantities in the Reynolds turbulence equations stem from the misunderstandings of physics of the Reynolds stress tensor. The current literatures have a consensus that the Reynolds stress tensor has six unknowns; however, this study shows that the Reynolds stress tensor actually has only three ones, namely the three components of fluctuation velocity. With this new understanding, the closed turbulence equations for incompressible flows are proposed.

Turbulence is a fluid dynamic problem refractory to mathematical treatment. Examining a theoretical model of liquid water flowing in a cylinder at different Raleigh numbers, we propose a novel approach to elucidate the first stages of turbulent flows. The weakly bonded molecular assemblies of liquid water form a fluctuating branched polymer in which every micro-cluster displays different density. Against the common view of liquid water as an incompressible and continuous fluid, we suggest that the occurrence of transient local aggregates could be able to generate the vortices and eddies that are the hallmarks of turbulence. We quantify the local changes in velocity, diameter and density required to engender “obstacles” to the average flow. Then, we show how these microstructures, equipped with different Raleigh numbers and characterized by high percolation index, could generate boundary layers that contribute to micro-vortices production. We conclude that the genesis of turbulence cannot be assessed in terms of collective phenomena, rather is sustained, among many other factors, by the underrated microscopic inhomogeneities of fluids like liquid water.

Under suitable continuity and uniqueness conditions, solutions of an α order Caputo fractional boundary value problem are differentiated with respect to boundary values and boundary points. This extends well-known results for nth order boundary value problems. The approach used is a standard technique and makes heavy use of recent results for differentiation of solutions of Caputo fractional intial value problems with respect to initial conditions and continuous dependence for Caputo fractional boundary value problems.

Localized turbulence is common in geophysical flows, where the roles of rotation and stratification are paramount. In this study, we investigate the evolution of a stratified turbulent cloud under rotation. Recognizing that a turbulent cloud is composed of vortices of varying scales and shapes, we start our investigation with a single eddy. Compared to an eddy under pure rotation, the stratified eddy shows the physical manifestation of a known potential vorticity mode, appearing as a static stable vortex. In addition, the expected shift from inertial waves to inertial-gravity waves is observed. In our numerical simulations of the turbulent cloud, carried out at a constant Rossby number over a range of Froude numbers, stratification causes columnar structures to deviate from vertical alignment. This deviation increases with increasing stratification, slowing the expansion rate of the cloud. The observed characteristics of these columnar structures are consistent with the predictions of linear theory, particularly in their tilt angles and vertical growth rates, suggesting a significant influence of inertial-gravity waves. Using Lagrangian particle tracking, we have identified regions where wave activity dominates over turbulence. In scenarios of milder stratification, these inertial-gravity waves are responsible for a significant energy transfer away from the turbulent cloud, a phenomenon that attenuates with increasing stratification.

Today, Colebrook’s equation is mostly accepted as an informal standard for modeling of turbulent flow in hydraulically smooth and rough pipes including transient zone in between. The empirical Colebrook’s equation relates the unknown flow friction factor (λ) with the known Reynolds number (R) and the known relative roughness of inner pipe surface (ε/D). It is implicit in unknown friction factor (λ). Implicit Colebrook’s equation cannot be rearranged to derive friction factor (λ) directly and therefore it can be solved only iteratively [λ=f(λ, R, ε/D)] or using its explicit approximations [λ≈f(R, ε/D)]. Of course, approximations carry in certain error compared with the iterative solution where the highest level of accuracy can be reached after enough number of iterations. The explicit approximations give a relatively good prediction of the friction factor (λ) and can reproduce accurately Colebrook’s equation and its Moody’s plot. Usually, more complex models of approximations are more accurate and vice versa. In this paper, numerical values of parameters in various existing approximations are changed (optimized) using genetic algorithms to reduce maximal relative error. After this improvement computational burden stays unchanged while accuracy of approximations increases in some of the cases very significantly.

This study examines the differential impact of proactive and reactive job crafting behaviors on job well-being, drawing from the perspective of boundary theory. The increasing popularity of work connectivity behavior after-hours (WCBA) after work has attracted widespread attention from scholars on the relationship between WCBA and employee well-being. One view suggests that the impact of WCBA is negative, while another view suggests it is positive. Obviously, the impact of (WCBA) on well-being at work is still contradictory. To clarify the complexity of the impact of WCBA on well-being at work, based on boundary theory, we divided WCBA into proactive WCBA (PC) and reactive WCBA (RC), and examined the double-edged sword effect of WCBA on well-being at work, as well as the mediating mechanisms of job control and work-to-home conflict, and the moderated effects of boundary segmentation preferences. This study uses empirical sampling method to collect data from 125 employees for a period of five days for quantitative research. The results show that, first, PC has an inverted U-shaped effect on job control, and job control is negatively related to well-being at work. Thus the mediating effect of job control is significant. Second, RC has a negative impact on job control, and there is also a negative relationship between job control and well-being at work. Therefore, the mediating effect of job control is significant. Third, PC and RC are positively correlated with work-to-home conflict, and work-to-home conflict has a significant positive impact on well-being at work. Therefore, the mediating effect of work-family conflict is significant. Fourth, the study also found that integration preference moderates the relationship between work-to-home conflict and PC on well-being at work, that is, the mediating effect of work-to-home conflict is stronger for employees with low integration preference. This study enriches our understanding of WCBA, PC and RC from the perspective of proactive and passive employee behaviors. The study also provides a new interpretation of the impact of WCBA on well-being at work and offers valuable insights for the sustainable development in digital social transformation and the application of boundary theory and the theory of empowerment-subjugation in achieving the United Nations Sustainable Development Goals. In addition, the study deepens our understanding of the heterogeneous regulatory role played by work-family integration preferences in influencing well-being at work under different types of WCBA.

In the present paper we investigate the conservative conditions in Quadratic Gravity. It is shown explicitly that the Bianchi identities lead to the conservative condition of the left-hand-side of the (gravitational) field equation. Therefore, the total energy-momentum tensor is conservative in the bulk (like in General Relativity). However, in Quadratic Gravity it is possible to have singular hupersurfaces separating the bulk regions with different behavior of the matter energy-momentum tensor or different vacua. They require special consideration. We derived the conservative conditions on such singular hypersurfaces and demonstrated the very possibility of the matter creation. In the remaining part of the paper we considered some applications illustrating the obtained results.

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.

Atmospheric turbulence, which produces chaotic motions in the planetary boundary layer, can inhibit mixing between fast-reacting species produced or released at different locations. This segregation process modifies the effective rate at which reactions occur between these species, and are not appropriately accounted for in coarse resolution models since these models assume complete mixing of tracers within each gridbox. Here we present a few examples of LES-based simulations applied to chemically reactive species in a forested area with high emissions of biogenic hydrocarbons, an urban area rich in anthropogenic emissions; and a maritime area with high emissions of reduced sulfur species.

This paper proposed an explicit and simple representation of velocity fluctuation and the Reynolds stress tensor in terms of the mean velocity field. The proposed turbulence equations are closed. The proposed formulations reveal that the mean vorticity is the key source of producing turbulence. It is found that there are no velocity fluctuation and turbulence if there were no vorticity. As a natural consequence, the laminar- turbulence transition condition was obtained in a rational way.

The synchronization of complex networks is believed to be closely related to network topology. How to study synchronization of real network without clear topology and parameters? In this Letter, a synchronized stability boundary equation is derived that is system- and disturbance-independent and applicable to arbitrarily coupled grids. The results imply that spontaneous synchronization in a network may be network independent. These conclusions provide new research paths for network synchronization and analyze the synchronization stability of grids in a unified way.

The diffraction-like process displayed by a spatially localized matter wave is here analyzed in a case where the free evolution is frustrated by the presence of hard-wall-type boundaries (beyond the initial localization region). The phenomenon is investigated in the context of a nonrelativistic, spinless particle with mass m confined in a one-dimensional box, combining the spectral decomposition of the initially localized wave function (treated as a coherent superposition of energy eigenfunctions) with a dynamical analysis based on the hydrodynamic or Bohmian formulation of quantum mechanics. Actually, such a decomposition has been used to devise a simple and efficient analytical algorithm that simplifies the computation of velocity fields (flows) and trajectories. As it is shown, the development of space-time patters inside the cavity depends on three key elements: the shape of the initial wave function, the mass of the particle considered, and the relative extension of the initial state with respect to the total length spanned by the cavity. From the spectral decomposition it is possible to identify how each one of these elements contribute to the localized matter wave and its evolution; the Bohmian analysis, on the other hand, reveals aspects connected to the diffraction dynamics and the subsequent appearance of interference traits, particularly recurrences and full revivals of the initial state, which constitute the source of the characteristic symmetries displayed by these patterns. It is also found that, because of the presence of confining boundaries, even in cases of increasingly large box lengths, no Fraunhofer-like diffraction features can be observed, as happens when the same wave evolves in free space. Although the analysis here is applied to matter waves, its methodology and conclusions are also applicable to confined modes of electromagnetic radiation (e.g., light propagating through optical fibers).

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

One route to reducing CO2 emissions is to improve the energy efficiency of machines. For example, conventional combustion engines are being downsized (and down-speeded), are now running on lower viscosity engine lubricants (such as SAE 0W-20 or lower viscosity grades) and often also have stop-start systems fitted (to prevent engine idling when the vehicle is stopped). Some of these changes result in higher levels of mixed and boundary friction, and so accurate estimation of mixed/boundary friction losses is becoming of increased importance, both for estimating friction losses, and also for estimating wear volumes. Traditional approaches to estimating mixed/boundary friction, which employ real area of contact modelling such as the Greenwood-Tripp model [1] are widely used, but recent experimental data suggests that these approaches underestimate mixed/boundary friction losses [2-5]. A new model has been developed, based on experimental data, which estimates the proportion of mixed/boundary lubrication, X, as a function of the value (where is the ratio of the oil film thickness separating the surfaces to the combined root mean square surface roughness of the interface surfaces). The precise equation that describes the way in which X varies with takes the form of a “reverse S-curve”, which makes physical sense since S-curves arise naturally in growth processes, and in mixed/boundary lubrication, the real area of contact of rough surfaces grows as the load (and 1/) increases. Numerical estimates of mixed/boundary lubrication losses in internal combustion engines are presented here and compared with published data [2,6]. In addition, the improved method described here is used to estimate both the financial cost of mixed/boundary lubrication for today’s passenger car fleet, and the CO2 emissions associated with these friction losses.

In this study, the Bohler K 190 steel was used. The steel was manufactured by the powder metallurgy (PM) process. The boronizing process was carried out in the range of 1173 to 1323 K, for 1-10 h. The samples were boronized in solid medium, called the Durborid powder mixture. For the microstructural observations, the scanning electron microscopy was utilized for determining the morphology of interfaces and measuring the layers’ thicknesses. The phase composition of boride layers was also determined with X-ray diffraction analysis. To investigate the redistribution of chemical elements redistribution during the boronizing process, the EDS mapping and EDS point analysis were used. The boride layers were constituted by FeB and Fe2B phases except for 1173 K for 1 h. The values of Vickers microhardness of Fe2B, FeB and transition zone were estimated. Finally, to assess the boron activation energies in FeB and Fe2B, the so-called integral method was applied and the results in terms of activation energies were confronted with the literature data.

The involvement of macular function in its preganglionic elements, during the neurodegenerative process of multiple sclerosis (MS), is controversial. In this case-control observational and retrospective study, we assessed multifocal electroretinogram (mfERG) responses from 41 healthy Controls, 41 relapsing-remitting MS patients without optic neuritis (ON) (MS-noON Group), 47 MS patients with ON: 27 with full recovery of high-contrast best corrected visual acuity (BCVA) (MS-ON-G Group) and 20 with poor recovery of BCVA (MS-ON-P Group). MfERG N1 and P1 implicit times (ITs), and N1-P1 response amplitude densities (RADs) were measured from concentric rings (R) with increasing foveal eccentricity: 0-5° (R1), 5-10° (R2), 10-15° (R3), 15-20° (R4), 20- 25° (R5), and from retinal sectors [superior, nasal, inferior and temporal] between 0-15° and 0- 25°. In MS-ON-P Group, mean mfERG RADs detected from R1 (0-5°) and from the central nasal quadrant (0-15°) were significantly reduced (p<0.01) with respect to those of Control, MS-noON and MS-ON-G Groups. No other significant differences between Groups for any mfERG parameters were found. Our results suggest that in MS, exclusively after ON with poor recovery of BCVA, the neurodegenerative process can induce dysfunctional mechanisms involving photoreceptors and bipolar cells of the fovea and of the more central nasal macular area.

In recent years, the field of deep learning achieved considerable success in pattern recognition, image segmentation and may other classification fields. There are a lot of studies and practical applications of deep learning on images, video or text classification. In this study, we suggest a method for changing the architecture of the most performing CNN models with the aim of designing new models to be used as stand-alone networks or as a component of an ensemble. We propose to replace each activation layer of a CNN (usually a ReLu layer) by a different activation function stochastically drawn from a set of activation functions: in this way the resulting CNN has a different set of activation function layers.

The Cobb angle is critical in assessing adolescent idiopathic scoliosis (AIS) patients. This study aimed to evaluate the error in selecting the cranial and caudal limit vertebrae on AIS digital X-rays of experienced and novice observers and its correlation with the error in measuring the Cobb angle and determining the length of the scoliotic curves. Using the TraumaMeter v.873 software, eight raters (four experienced and four novices) independently evaluated 68 scoliotic curves on 42 AIS X-rays across three sessions one month apart according to Hopkins’ criteria of absolute reliability. Two spine specialist experts using the Cobb method established the gold standard measurements. The percentage of error in the selection of the cranial vertebrae (the most-tilted vertebra above the apex) is higher than in the caudal vertebrae (the most-tilted vertebrae below the apex) (44.7%, CI95% 41.05-48.3 compared to 35%, CI95% 29.7-40.4). The error in the choice of the most tilted vertebrae, in number of vertebrae per scoliotic curve, was MBE (mean bias error) =0.45, CI95% 0.38-0.52 in the cranial and MBE=0.35, CI% 0.69-0.91 in the caudal. The percentage of errors in the choice of boundary vertebrae is lower for the experienced than for the novices, both in cranial and caudal. The experienced chose at least one boundary vertebrae incorrectly in 52% of the meas-urements and the novices in 65.1%. The error in the Cobb angle did not depend on the number of boundary vertebrae chosen incorrectly (r=0.198, p=0.111). Of the 50% of Cobb angle measurements closest to the actual value, 55.4% (454) were obtained by incorrectly selecting at least one boundary vertebra. There was a positive correlation (r=0.673, p=0.000) between the error in selecting the boundary vertebrae and the error in determining the length of the scoliotic curves. The error in quantifying the length of scoliotic curves was MBE=0.43, CI95% 0.35-0.51 vertebrae per curve. We can conclude that errors in selecting boundary vertebrae are common among experienced and novice observers, with a greater frequency of error in the cranial vertebrae. There is no correlation between the error in selecting the boundary vertebrae and the error in measuring the Cobb angle. Contrary to the consensus, the accuracy of determining the length of the scoliotic curve is limited by the Cobb method’s reliance on the correct selection of boundary vertebrae. Inappropriate choice of boundary vertebrae leads to an estimated error in determining the length of the scoliotic curves of 0.805 vertebrae for each boundary vertebra chosen incorrectly.

In this study, the generality and prediction accuracy of a generalised series model for large eddy simulation of premixed and non-premixed turbulent combustion is explored. The model is based on Taylor series expansion of the chemical source term in scalar space and implemented into OpenFOAM. The mathematical model does not depend on the combustion regimes, and has the correct limiting behaviour. The numerical error sources are outlined and analysed. The model is first applied to a piloted methane/air non-premixed jet flame (Sandia Flame D). The statistical (time-averaged and RMS) results agree well with experimental measurements, particularly with regard to the mixture fraction, velocity, temperature, and concentrations of major species CH4, CO2, H2O, and O2. The concentrations of intermediates CO and H2 are over-predicted, due to the limitations of the reduced reaction mechanism employed. Then, a Bunsen piloted flame is simulated. Most of the statistical properties of both reactive species and progress variable are well reproduced. The only major discrepancy evident is in the temperature, which is attributed to the experimental uncertainties of temperature fields in the pilot stream. These findings have demonstrated the model`s generality for both premixed and non-premixed combustion simulation, as well as the accuracy of prediction of reactive species distribution.

RANS simulations have been broadly used to investigate turbulence in the oceans and atmosphere. Within these environments there are a multitude of tracers undergoing reactions (eg. phytoplankton growth, chemical reactions). The distribution of these reactive tracers is strongly influenced by turbulent mixing. With a 50 member ensemble of two dimensional Rayleigh-Taylor induced turbulent mixing, we show that the dynamics of a reactive tracer growing according to Fisher’s equation are poorly captured by the mean. A fluctuation-dependent sink introduced by Reynolds averaging Fisher’s equation transfers concentration from the mean to the fluctuations. We compare the dynamics of the reactive tracer with those of a passive tracer. The reaction causes the concentration of reactive tracer to increase thereby increasing Fickian diffusion and allowing the reactive tracer to diffuse into turbulent structures that the passive tracer cannot reach. A positive feedback between turbulent mixing and fluctuation growth is identified. We show that eddy viscosity and diffusivity parameterizations fail to capture the bulk trends of the system and identify a need for negative eddy diffusivities. One must therefore be cautious when interpreting RANS results for reactive tracers.

Wind loads have become one of the key influence factors for the running safety of vehicles and comfort of passengers. Investigation on the wind speed spectrum characteristics of a moving vehicle under turbulent crosswinds is of great influence. Expressions of the wind speed spectrum of a moving vehicle was obtained from the von Kármán spectrum based on Taylor’s frozen flow hypothesis. The influence factors, including the ratio of the vehicle speed to the wind speed and the wind yaw angle from 15° to 175°, were analyzed. The maximum value of the wind speed spectrum and the corresponding frequency were studied as well. The results show that the maximum values of the wind speed spectrum of the moving vehicle were larger than those of the static vehicle. The maximum value of the wind speed spectrum corresponding to the moving vehicle first increased and then decreased as the wind yaw angle increased. Some of the frequencies corresponding to the longitudinal wind speed spectrum values of the moving vehicles were smaller than those of the static vehicle. For moving vehicles, the frequency values corresponding to the maximum values of the longitudinal wind speed spectrum first increased and then decreased as the ratio of the vehicle speed to the wind speed and the wind yaw angle increased.

This paper presents a novel theoretical and computational methodology for the generation of an onset turbulent field with prescribed properties in the numerical simulation of an arbitrary viscous flow. The methodology is based on the definition of a suitable distribution of volume force terms in the right–hand side of the Navier–Stokes equations. The distribution is represented by harmonic functions that are randomly variable in time and space. The intensity of the distribution is controlled by a simple PID strategy in order to obtain that the generated turbulent flow matches a prescribed turbulence intensity. A further condition is that a homogeneous isotropic flow is es- tablished downstream of the region where volume force terms are imposed. Although it is general, the proposed methodology is primarily intended for the computational modelling of hydrokinetic turbines in turbulent flows representative of tidal or riverine installations. A first numerical applica- tion is presented by considering the injection of homogeneous and isotropic turbulence with 16% intensity into a uniform unbounded flow. The analysis of statistical properties as auto-correlation, power spectral density, probability density functions, demonstrates that the generated flow tends to achieve satisfactory levels of stationarity and isotropy, whereas the simple control strategy used determines overestimated turbulent intensity levels.

This paper gives a new unified formula for the Newtonian fluids valid for all pipe flow regimes from laminar to the fully rough turbulent. It includes laminar, unstable sharp jump from laminar to turbulent, and all types of the turbulent regimes: smooth turbulent regime, partial non-fully developed turbulent and fully developed rough turbulent regime. The formula follows the inflectional form of curves as suggested in Nikuradse’s experiment rather than monotonic shape proposed by Colebrook and White. The composition of the proposed unified formula consists of switching functions and of the interchangeable formulas for laminar, smooth turbulent and fully rough turbulent flow. The proposed switching functions provide a smooth and a computationally cheap transition among hydraulic regimes. Thus, the here presented formulation represents a coherent hydraulic model suitable for engineering use. The model is compared to existing literature models, and shows smooth and computationally cheap transitions among hydraulic regimes.

In this paper, we propose an efficient numerical method to solve the diffusive logistic model with a free boundary which is often used to simulate the spreading of the new or invasive species. The boundary movement is tracked by the level-set method, where the Hamilton-Jacobi weighted essentially nonoscillatory (HJ-WENO) scheme is utilized to capture the boundary curve embedded by the Cartesian grids via the embedded boundary method. Then the radial basis function-finite difference (RBF-FD) method is adopted to the spatial discretization and the IMEX scheme is considered for time integration. A variety of numerical examples are carried out to demonstrate the evolution of diffusive logistic model in different initial boundaries.

This article is the fourth part of a scientific project under the general title “Geometrized vacuum physics based on the Algebra of Signatures". In the first three articles [1,2,3], the foundations of the Algebra of Stignatures were laid and the main aspects of the kinematics of vacuum layers were considered. This article continues the development of the mathematical apparatus of the proposed project, in particular, the dynamics of vacuum layers is developed based on the Algebra of Signatures. The development of this direction of research (with simplifications related to Riemann's differential geometry) led to the possibility of a geometrized representation of the electric field strength and magnetic field induction. This geometrized mathematical apparatus allows one to interpret the electromagnetic field as an interweaving of accelerated and rotational flows of the adjacent layers of vacuum. The proposed dynamic models of accelerated movements and rotations of vacuum layers can provide a theoretical basis for the development of “zero” (i.e. vacuum) technologies.

We fabricate large-area atomically thin MoS₂ layers through the direct transformation of crystalline molybdenum trioxide (MoO₃) by sulfurization at relatively low temperatures. The obtained MoS₂ sheets are polycrystalline (~10-20 nm single-crystal domain size) with areas of up to 300×300 µm² with 2-4 layers in thickness and show a marked p-type behaviour. The synthesized films are characterized by a combination of complementary techniques: Raman spectroscopy, X-ray diffraction, transmission electron microscopy and electronic transport measurements.

The turbulent Poiseuille flow between two parallel plates is one of the simplest possible physical situations, and it has been studied intensively. In this paper, we propose a modified Prandtl-van Driest mixing length that satisfies both boundary conditions and wall damping effects. With our new formulations, we numerically solve the problem and, moreover, propose an approximate analytical solution of mean velocity. As applications of our solution, an approximate analytical friction coefficient of turbulent Poiseuille flow is proposed.

The mean flow and turbulence statistics of the flow through a simplified urban environment, which is an active research area in order to improve the knowledge of turbulent flow in cities, is investigated. This is useful for civil engineering, pedestrian comfort and for health concerns caused by pollutant spreading. In this work, we provide analysis of the turbulence statistics obtained from well-resolved large-eddy simulations (LES). A detailed analysis of this database reveals the impact of the geometry of the urban array on the flow characteristics and provides for a good description of the turbulent features of the flow within a simplified urban environment. The most prominent features of this complex flow include coherent vortical structures such as the so-called arch vortex, the horseshoe vortex and the roof vortex. These structures of the flow have been identified by an analysis of the turbulence statistics. The influence of the geometry of the urban environment (and particularly the street width and the building height) on the overall flow behavior have also been studied. Finally, the well-resolved LES results were compared with the experimental database from Monnier et al. to discuss differences and similarities between the respective urban configurations.

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

Distribution systems can form islands when faults occur. Each island represents a subsection with variable boundaries subject to the location of a fault(s) in the system. A subsection with variable boundaries is referred to as island in this paper. For operation in autonomous mode, it is imperative to detect the island configurations and stabilize these subsections. This paper presents a novel scheme for the detection of islanding boundaries and stabilizing the system during autonomous operation. In the first stage, a boundary detection method is proposed to detect the configuration of the island. In the second stage a dynamic voltage sensitivity factor (DVSF) is proposed to assess the dynamic performance of the system. In the third stage, a wide area load shedding program is adopted based on DVSF to shed the load in weak busbars and stabilize the system. The proposed scheme is validated and tested on a generic 18-bus system using a combination of EMTDC/PSCAD and MATLAB softwares.

MPC (Magneto-Photonic Crystal) Optimisation is a feature-rich Windows software application designed to enable researchers to analyze the optical and magneto-optical spectral properties of multilayers containing gyrotropic constituents. A set of computational algorithms aimed at enabling the design optimization and optical or magneto-optical (MO) spectral analysis of 1D magnetic photonic crystals (MPC) is reported, together with its Windows software implementation. Relevant material property datasets (e.g., the optical spectra of refractive index, absorption, and gyration) of several important optical and MO materials are included, enabling easy reproduction of the previously published results from the field of MPC-based Faraday rotator development, and an effective demonstration-quality introduction of future users to the multiple features of this package. We also report on the methods and algorithms used to obtain the absorption coefficient spectral dispersion datasets for new materials, for which the film thickness, transmission spectrum, and refractive index dispersion function are known.

Quasi-interpolation is a powerful tool for approximating functions using radial basis functions (RBFs) such as Gaussian kernels. This avoids solving large systems of equations as in RBF interpolation. However, quasi-interpolation with Gaussian kernels on compact intervals can have significant errors near the boundaries. This paper proposes a quasi-interpolation method with Gaussian kernels using Chebyshev points and boundary corrections to improve the approximation near the boundaries. The boundary corrections use a linear approximation of the function beyond the interval to estimate the truncation error and add correction terms. Numerical studies on test functions show that the proposed method reduces errors significantly near boundaries compared to quasi-interpolation without corrections, for both equally spaced and Chebyshev points. The convergence and accuracy with the boundary corrections are generally better with Chebyshev points compared to equally spaced points. The proposed method provides an efficient way to perform quasi-interpolation on compact intervals while controlling the boundary errors.

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.

This paper proposes a data-driven approach to the calculation of ship overtaking boundary conditions in narrow channels based on ship overtaking characteristics by the Automatic Identification System (AIS) data. First, we collect the ship tracking process data in a narrow channel by processing the original AIS data, followed by a normality test and correlation analysis, and finally, the boundary conditions are obtained using a weighted K-means clustering algorithm. to identify the relevant data. We applied the proposed approach to estimate the overtaking boundary conditions of ships in the straight channel section of Shanghai South Passage. The method yields boundary conditions that significantly reduce collision risk when a ship engages in pursuit behavior in narrow waterways, providing a scientific basis for the overtaking of Marine Autonomous Surface Ships (MASS) in autonomous navigation within such environments, further promoting the automation level of USVs and autonomous vessels.

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.

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.

The efficiency of organic solar cells (OSCs) is influenced by various factors, among which environmental temperature plays a significant role. Previous studies have shown that the thermal stability of these cells can be enhanced by incorporating a third component into their structure. Ternary organic solar cells, in particular, have shown promising results in improving thermal stability. A well-designed electron transport layer (ETL) can significantly bolster thermal stability by facilitating efficient charge transport and reducing charge recombination. In this study, we investigated the effect of temperature on the efficiency of inverted ternary structures, ranging from 300 K to 400 K. The structures examined include FTO/SnO2/PM6:D18:L8-BO/PEDOT: PSS/Ag, FTO/Spiro-OMeTAD/PM6:D18:L8-BO/PEDOT: PSS/Ag, and FTO/PC60BM/PM6:D18:L8-BO/PEDOT: PSS/Ag. Through simulations employing three different electron transport layers—SnO2, Spiro-OMeTAD, and PC60BM—we observed that at 340 K (65°C) these structures maintained 92% of their original high efficiency of ~20% at 300 K, which represents a high level of thermal stability under the high PCE condition.

The study of coherent vortices remains an urgent field of singular optics of vortex beams. In this paper, coherent properties of partially coherent vortex Bessel—Gaussian optical beams propagating through a uniform medium (free space) or a turbulent atmosphere are examined theoretically. The consideration is based on analytical solution of the equation for the transverse second-order mutual coherence function of the field of a partially coherent optical radiation in a turbulent atmosphere. For the partially coherent Bessel—Gaussian beam, the second-order mutual coherence function of the source field is taken as a Gaussian Shell-model. In this approximation, we analyze the behavior of the coherence degree and the integral coherence scale of these beams as a function of the propagation pathlength, propagation conditions, and beam parameters, such as the radius of the Gauss factor of the beam, parameter of the Bessel factor of the beam, topological charge, and correlation width of the source field of partially coherent radiation. In addition, the integral coherence scale of the partially coherent vortex Bessel—Gaussian beam is compared with that of the partially coherent Gaussian beam. It is found that as a partially coherent vortex Bessel—Gaussian beam propagates through a turbulent atmosphere, there appear not two (as might be expected: one due to atmospheric turbulence and another due to partial coherence of the source field), but only one ring dislocation of the coherence degree (due to the simultaneous effect of both these factors on the optical radiation). In addition, it is shown that the dislocation of the coherence degree that affects significantly the beam coherence level is formed only for beams, for which the coherence width of the source field is larger than the diameter of the first Fresnel zone.

This research provides a computationally efficient numerical tool for predicting 3-D turbulent flows over hydrofoils. It addresses a gap in existing numerical methods by extending the capabilities of the laminar VIScous Vorticity Equation (VISVE) solver to handle 3-D turbulent flow scenarios. The implementation integrates the k − ω SST model into the 3-D VISVE solver using the Finite Volume Method (FVM), extending its capability to handle turbulent flows. The turbulence model is also parallelized to improve computational efficiency. The application of the newly developed model to 3-D hydrofoil cases expands its scope and serves as preparatory work for future rotational propeller applications. Rigorous testing includes a convergence analysis concerning grid numbers and time step sizes, and validation against a Reynolds-Averaged Navier-Stokes (RANS) solver ensures the reliability and accuracy of the enhanced VISVE solver for turbulent flow simulations. The turbulent VISVE solver offers advantages such as reduced computational domain and costs through a vorticity-based approach. Turbulence concentration within boundary layers and free shear flows does not compromise the advantages of the small computational domain. The simplified meshing process, automatically generated by denoting the number of panels on the hydrofoil, makes the solver accessible for researchers and engineers.

Swirling coal flames are encountered in swirl burners of boiler furnaces and cyclone coal combustors, and were studied using large-eddy simulation (LES) in reported references. However, the structures of swirling coal flames remain to be further studied. In the present paper LES of swirling coal flames is conducted for two combustors using two gas combustion models. The statistical simulation results are accessed by experiments and the instantaneous simulation results are used to analyze the specific features of flame structures of four cases. It is seen that different inlet flow conditions and different secondary-air flow rate lead to different flame structures

A formula for the turbulent Prandtl number is derived in terms of local variables available from two-equations turbulence models. The derivation is a direct consequence of the expected square additivity of the molecular and turbulent parameters defining the effective viscosity and the effective conductivity. The formula does not degenerate and leads to a Kays like formulation if approximated. Its utility and its limitation are discussed.

Based on the insight gained from single jet analysis performed earlier [1], CFD analysis on turbulent jets impinging on one another at an angle was performed. Multiple impingement angles were considered for this study to gain better understanding of the parameters affecting resultant jet growth and velocity distribution. From the study of single jet, it was concluded that the SST k-ω turbulence model was the ideal turbulence model capable of accurately predicting the flow physics of the jets exiting a fully developed pipe at low Reynolds number. Hence, for the study of impinging jets, SST k-ω turbulence model was used to study the velocity and jet growth characteristics. It became evident that the mesh alignment with the velocity vector at exit of the pipe domain plays a crucial role in the accuracy of the results. The parameter used to evaluate this condition was identified as False Diffusion and was observed to affect the TKE parameter significantly. Methods to reduce false diffusion are also discussed in this article. Based on the mesh obtained from the grid sensitivity study, jets impinging at 30, 45 and 60 degrees at Reynolds number of 7500 were numerically analyzed. It was observed that the profile of the resultant jet closely matched with the prediction of elliptical profile predicted by past researchers, [2] and [3]. Also, it was seen that higher jet growth was predicted in case of jets impinging at a higher impingement angle.

We present information theoretic search strategies for single and multi-robot teams to localize the source of biochemical contaminants in turbulent flows. The robots synthesize the information provided by sporadic and intermittent sensor readings to optimize their exploration strategy. By leveraging the spatio-temporal sensing capabilities of a mobile sensing network, our strategies result in control actions that maximize the information gained by the team while optimizing the time spent localizing the position of the biochemical source. By leveraging the team's ability to obtain simultaneous measurements at different locations, we show how a multi-robot team is able to speed up the search process resulting in a collaborative information theoretic search strategy. We validate our proposed strategies in both simulations and experiments.

Research subjects of this study are four representative locations in the industrial complex, in the city of Banja Luka, Republic of Srpska, Bosnia and Herzegovina. 16 polycyclic aromatic hydrocarbons (PAHs) (∑16PAHs), humus and pH were determined. The main objective of the paper is to determine the concentration levels, to assess the probable sources of PAHs contamination in soil and groundwater and to determine the ecological risk. The ∑16PAHs in soil (at depths of 30 cm, 100 cm, 200 cm, 300 cm and 400 cm) ranged from 0.99 to 2.24 mg/kg, from 0.34 to 0.46, from 0.24 to 0.32, from 0.13 to 0.27 and from 0.13 to 0.47, with mean values of 1.70 mg/kg, 0.40 mg/kg, 0.28 mg/kg, 0.20 mg/kg and 0.26 mg/kg, respectively. The ∑16PAHs in groundwater ranged from 0.23 to 4.50 mg/m³, with a mean value of 1.42 mg/m³. Surface soil and groundwater are heavily contaminated. All values of ∑PAHs in soil layers were lower in the depths of the soil. Factor analysis indicates three sources of contamination, i.e. principal component (PC) PC1 (pyrogenic), PC2 (petrogenic) and PC3 (biomass), with 52.39%, 26.14% and 8.46% of the total variance, respectively. ∑PAH and PAHs indicate high ecological risk for most PAHs, which decreases with soil depth.

Starting from the analysis of the lack of positivity of the Cattaneo heat equation, the article addresses the thermodynamic relevance of the positivity constraint in irreversible thermodynamics, than is at least as important as the entropic constraints. The fulfillment of this condition in hyperbolic models leads to the parametrization of the concentration fields with respect to internal variables associated with the microscopic dynamics. Using Brownian motion theory as a landmarking example for deriving macroscopic transport equations from the equations of motion at the particle/molecular level, we discuss two typical problems involving hydrodynamic interactions at microscale: surface chemical reactions at a solid interface of a diffusing reactant, and mass-balance equations in a complex viscoelastic fluid, in which the physics of the interaction leads either to the overcome the parabolic diffusion model or to consider the parametrization of the concentration with respect to the degrees of freedom associated with the relaxation dynamics of the solvent fluid.

Introduction: The concept of hypersexuality belongs to modern parlance, according to a predominantly clinical meaning, understood as a psychological and behavioural alteration as a result of which sexually motivated stimuli are sought in inappropriate ways, often experienced in a way that is not completely satisfactory. Methods: Literature up to February 2023 was reviewed. Results: Forty-two articles were included in the review. Conclusion: Hypersexuality is a potentially clinically relevant condition, consisting of one or more dysfunctional and pathological behaviours of one's sexual sphere and graded according to the severity of impairment of subjective acting out; for this reason, the Perrotta Hypersexuality Global Scale (PH-GS) is suggested, which distinguishes high-functioning forms (pro-active and dynamic hypersexuality) from those of attenuated and corrupted functioning (dysfunctional and pathological hypersexuality of grade I and II). Future research is hoped to address the practical needs of this condition, such as the exact etiopathology, the role of oxytocin in dopaminergic hypotheses (and its ability to attenuate the symptomatology suffered by the patient in terms of manic drive), the best structural and functional personality framing of the subject, and the appropriate therapy to pursue.

Fibreglass composite plates are used for applications ranging from aerospace and automotive industries to construction industry. Vibroacoustic properties of the fibreglass plate has been investigated for different boundary conditions. The theory of the plate vibration and sound radiation has been used to determine the deflection of the composite plate at different locations on its surface. 3-Dimentional (3D) deformations and their corresponding contour velocities of plate have been computed for a range of frequencies. It has been shown that clamping the plate at four edges delayed first resonance of the plate by 15 Hz and second resonance of the plate by 30 Hz. Vibroacoustic indicators and the radiation impedance matrix of the plate have also been calculated. The results of this study show that fibreglass plate could be used in building industry to dampen the vibration of machine and structural vibration of buildings, and to sandwich building structures to form panels in order to attenuate airborne sound and to lower noise transmission of structural borne sound

In this work, meshless methods are applied for the solution of two-dimensional steady-state heat conduction problems with nonlocal multi-point boundary conditions (NMBC). These meshless procedures are based on multiquadric radial basis function (MQ RBF) and its modified version (i.e. integrated MQ RBF). The proposed meshless methods which were recently published in \cite{Reutskiy2016} is compared with standard collocation method (i.e. Kansa's method). Three different sorts of solution domain are considered in which Dirichlet boundary condition is specified on some part of the boundary and is related to the unknown function values at a discrete set of interior points. The influence of NMBC on the accuracy and condition number of the system matrix associated to the proposed methods is investigated. The sensitivity of the shape parameter is also analyzed in the proposed methods. Performance of the proposed approaches in terms of accuracy and efficiency is confirmed on the benchmark problems.

Recently, the conformable derivative and its properties have been introduced. In this paper, we propose and prove some new results on conformable Boundary Value Problems. First, we introduce a conformable version of classical Sturm´s separation, and comparison theorems. For a conformable Sturm-Liouville problem, Green's function is constructed, and its properties are also studied. In addition, we propose the applicability of the Green´s Function in solving conformable inhomogeneous linear differential equations with homogeneous boundary conditions, whose associated homogeneous boundary value problem has only trivial solution. Finally, we prove the generalized Hyers-Ulam stability of the conformable inhomogeneous boundary value problem.

Crop type and field boundary mapping enable cost-efficient crop management on the field scale and serve as the basis for yield forecasts. Our study uses a data set with crop types and corresponding field borders from the federal state of Bavaria, Germany, as documented by farmers from 2016 to 2018. The study classified corn, winter wheat, barley, sugar beet, potato, and rapeseed as the main crops grown in Upper Bavaria. Corresponding Sentinel-2 data sets include the normalised difference vegetation index (NDVI) and raw band data from 2016 to 2018 for each selected field. The influences of clouds, raw bands, and NDVI on crop type classification are analysed, and the classification algorithms, i.e., support vector machine (SVM) and random forest (RF), are compared. Field boundary detection and extraction are based on non-iterative clustering and a newly developed procedure based on Canny edge detection. The results emphasise the application of Sentinel’s raw bands (B1–B12) and RF, which outperforms SVM with an accuracy of up to 94%. Furthermore, we forecast data for an unknown year, which slightly reduces the classification accuracy. The results demonstrate the usefulness of the proof-of-concept and its readiness for use in real applications.

Carbon fibers (CFs) cannot be directly used for the preparation of CF paper because of their chemically inert nature. Herein, the surface of CF was modified via the spontaneous oxidative self-polymerization of dopamine. By taking full advantage of the spontaneous oxidation and self-polymerization properties of PD to maintain the maximum strength of CFs, a polydopamine-modified CF paper (PDA-CFP) with excellent performance was prepared using PD-modified CFs (PDA-CFs), which increased the proportion of hydrophilic functional groups on the surface of carbon fibers, increased the O/C ratio on the CF surface by 6 times, and improved the bond strength between the modified CF and the adhesive by making full use of the interaction force between polydopamine and PVA fibers. In this way, the primary properties of the CF paper was improved. Overall, the results showed that the dispersion of CF was considerably improved by the dopamine modification. In addition, the primary physical properties of PDA-CFP were better than those of virgin CF paper (CFP-0). PDA-CFP exhibited a maximum tensile strength of 2.04 kN/m2, a minimum resistivity of 0.06055 Ω·cm and a minimum porosity of 72.4%. The tightness could be increased by up to 12.1%.

The Sun is a huge gaseous body. However, we cannot observe events in the inner Sun due to the convection zone opacity according to previous models. Therefore, the flares originate from the front surface of the Sun. But the current study relied on the distance distribution of X-Ray solar flares, which concluded that the inner layers have much lower opacity than expected. It is even less than what was expected by the latest models based on helioseismology. This means that the flares may originate from the solar interior or solar core, and perhaps from the backside surface, and even appear to us from the frontside surface. Which the re-estimate and correct the currently listed solar flare’s location is needed. Additionally, the flare’s distance illustrations the solar interior layers and appears their boundaries from the core to the photosphere. This method allows us to monitor the variation of the core’s radius with time. The model of the flare’s distance has been developed in current study. But this needs to redevelopment after re-estimating the solar flares locations.

Mixed cadmium tellurides – cadmium sulfides thin layers were formed on the polyamide PA 6. Monotelluropentathionic acid (H2TeS4O6) was used as precursor of tellurium and sulfur. We applied a low-temperature, nontoxic and cost-effective SILAR method. Cadmium telluride (CdTe) and sulfide (CdS) layers were formed through consecutive reactions of sorbed/diffused chalcogens species from telluropentathionate anion (TeS4O62–) with functional groups of polyamide and alkaline cadmium sulfate. The pseudo-second-order rate and Elovich kinetic models best fitted to quantify an uptake of chalcogens and cadmium on PA6. The effects of chalcogens and Cd on the structure and optical properties of PA6 were characterized using UV/Vis and IR spectra. The clear changes of these properties found in dependence on the concentration and exposure time in precursor's solutions. Fourier transform infrared spectroscopy and an ultraviolet-visible spectroscopy were applied in order to evaluate the effect of the chalcogens having particles on the changes in structure of polyamide 6 films depending on the exposure time in the solution of the chalcogens precursor and its concentration. The optical band gap energy of formed layers was found to be in the order of 1.52–2.36 eV. Studies of scanning electron microscopy and atomic force microscopy reveal that the diameter of the average grain is about 30 nm. The grains are conical in shape and unevenly distributed all over the surface of the substrate.

The studies of laminar unsteady boundary layer flows is crucial for understanding laminar-turbulence transition and origins of turbulence. However, the task of finding its solutions poses a significant challenge. In this paper, we propose a novel approach by introducing a similar transformation to convert the 2D unsteady laminar boundary layer equations into a single partial differential equation. By applying this transformation, we are able to obtain the "exact" similarity solutions for the 2D unsteady laminar boundary layer equations, specifically for the case of flat plate boundary flow. Notably, this is the first time that such an exact solution has been obtained. Application of this exact solution is being used to solve Stokes' first problem in 2D boundary layer. Perspectives on the transition fron Rayleigh solution to Blasius solution is provided.

The objective to fast-track the mapping and registration of large numbers of unrecorded land rights globally, leads to the experimental application of Artificial Intelligence (AI) in the domain of land administration, and specifically the application of automated visual cognition techniques for cadastral mapping tasks. In this research, we applied and compared the ability of rule-based systems within Object Based Image Analysis (OBIA), as opposed to human analysis, to extract visible cadastral boundaries from Very high resolution (VHR) World View-2 image, in both rural and urban settings. From our experiments, machine-based techniques were able to automatically delineate a good proportion of rural parcels with explicit polygons where the correctness of the automatically extracted boundaries was 47.4% against 74.24% for humans and the completeness of 45% for machine, as against 70.4% for humans. On the contrary, in the urban area, automatic results were counterintuitive: even though urban plots and buildings are clearly marked with visible features such as fences, roads and tacitly perceptible to eyes, automation resulted in geometrically and topologically poorly structured data, that could neither be geometrically compared with human digitised, nor actual cadastral data from the field. These results provide an updated snapshot with regards to the performance of contemporary machine-drive feature extraction techniques compared to conventional manual digitising.

Laser heat treatment is applied after coating. Evaluation of the results was performed by studying the microstructures by metallographic, SEM/EDX microscopy and the mechanical properties were obtained by microscopic hardness and abrasion resistance. The objective of this study was to investigate the effect of laser heat treatment on the wear resistance of metal coatings. The results show the influence of the microstructure and chemical composition of the electrodes used on the micro-hardness and wear resistance of the metal coatings.

The color of a thin copper oxide layer formed on a copper plate was transformed from reddish-brown into blue-purple by irradiation with 5 keV Ar+ ions to a fluence as low as 1 1015 Ar+ cm–2. In the unirradiated copper oxide layer, the copper valence state of Cu2+ as well as Cu+ and/or Cu0 was included as indicated by the presence of a shake-up satellite line in a photoemission spectrum. While for the irradiated one, the satellite line decreased in intensity, indicating that irradiation resulted in the reduction from Cu2+ to Cu+ and/or Cu0. Furthermore, nuclear reaction analysis using a 16O(d, p)17O reaction with 0.85 MeV deuterons revealed a significant loss of oxygen (51015 O atoms cm–2) in the irradiated layer. Thus, the chromatic change observed in the present work originated in the irradiation-induced reduction of a copper oxide.

This paper presents the results of the investigation of composite sinters W-TiB2 which were used as an electrode in the process of electro-spark deposition (ESD) and the examination of the deposited layers. The scope of the study includes detailed characteristics of powder mixtures, composite sinters made using the spark plasma sintering method (SPS) and layers deposited in the electro-spark process. The ESD process, using the W+30 vol.% TiB2 electrode, was carried out using an automated device. The substrates were made of copper and aluminium. The topography analysis of the surfaces of the composite layer and the evaluation of their wear resistance properties are also presented. The analysis of the results of research showed the possibility to obtain, using the SPS method, composite materials of good quality, which can be used as electrodes in the ESD process. The obtained layers had increased wear resistance in relation to the substrate material.

The effect of growth temperature and precursor flows on the doping level and surface morphology of Ge-doped GaN layers was researched. The results show that germanium is more readily incorporated at low temperature, high growth rate and high V/III ratio, thus revealing a similar behavior to what was previously observed for indium. V-pit formation can be blocked at high temperature but also at low V/III ratio, the latter of which however causing step bunching.

A century-old scientific conundrum is solved in this paper. The Prandtl mixing length modelled plane boundary turbulent flow is described by: $\frac{du^+}{d{y^+} }+\kappa^2{y^+} ^2\left( \frac{du^+}{d{y^+} }\right)^2=1$, together with boundary condition $ {y^+} =0:\, u^+=0$. Only an approximate solution to this nonlinear ordinary differential equation (ODE) has been sought so far, however, the exact solution to this ODE has not been obtained. By introducing a transformation, $2\kappa y^+=\sinh \xi$, I successfully find the exact solution of the ODE as follows: $u^+=\frac{1}{\kappa}\ln(2\kappa {y^+} +\sqrt{1+4\kappa^2{y^+} ^2})-\frac{2{y^+} }{1+\sqrt{1+4\kappa^2{y^+} ^2}}$.

AAs the global energy market undergoes a wholesale transformation accelerated by the need to decarbonise, a rapid transition to renewable energy and the mass deployment of distributed energy resources, autonomous energy networks or microgrids are emerging as an attractive mechanism for the delivery of electricity to end users. Yet in Australia, at least, relatively little is known about key aspects of microgrids that are fundamental to their successful deployment, not least the more commercial and economic elements rather than the purely technical. Drawing on the extant global literature on microgrids, this paper explores the most important of these aspects including business models, ownership and investment. Identifying the ambiguity, inconsistency and uncertainty evident in many of the feasibility studies currently in train across Australia, this paper highlights specific areas for future research that need to be addressed if the full potential of microgrids is to be realised in the context of a global energy transition both domestically and internationally.

Quantum turbulence is characterized by many degrees of freedom interacting non-linearly to produce disordered states, both in space and time. The advances in trapping, cooling, and tuning the interparticle interactions in atomic Bose-Einstein condensates (BECs) make them excellent candidates for studying quantum turbulence. In this work, we investigate the decaying regime of quantum turbulence in a trapped BEC. Although much progress has been made in understanding quantum turbulence, other strategies are needed to overcome some intrinsic difficulties. We present an alternative way of investigating this phenomenon by defining and computing a characteristic length scale, which possesses relevant characteristics to study the establishment of the quantum turbulent regime. One intrinsic difficulty related to these systems is that absorption images of BECs are projected to a plane, thus eliminating some of the information present in the original momentum distribution. We overcome this difficulty by exploring the symmetry of the cloud, which allows us to reconstruct the three-dimensional momentum distributions with the inverse Abel transform. We present our analysis with both the two- and three-dimensional momentum distributions, discussing their similarities and differences. We argue that the characteristic length allows us to visualize the time evolution of the turbulent state intuitively.

In this paper, we introduce a comprehensive solution aimed at enhancing the understanding of three-dimensional atmospheric pollutant dispersion. This innovative solution develops a generalized model that extends previous research and is applicable to all parameterization schemes of the equation, including wind speed profiles and turbulent diffusion coefficients, while incorporating the dry deposition criterion. Our methodology involves subdividing the atmospheric boundary layer into distinct sublayers, which facilitates a detailed examination of pollutant dispersion dynamics. Extensive validation with data from the Hanford experiment has demonstrated the accuracy of this solution in simulating pollutant concentrations. The results exhibit a strong correlation between projected and observed concentrations, underscoring the statistical reliability of our approach. This validation situates the statistical indices of our solution within an acceptable range, confirming its accuracy in predicting atmospheric pollutant dispersion. These findings thus establish our solution as a valid and effective method for studying complex environmental phenomena.

The Lorentz boost and rapidity for a radially falling frame are derived for the Schwarzschild metric in Kruskal-Szekeres coordinates. It is found through this analysis that the spacetime is radially length contracted toward the metric source in the falling frame as a result of the Lorentz boosting relative to the rest frame (which is shown to be consistent with the rest frame). When the falling frame reaches the event horizon, the Lorentz boost goes to infinity, causing the falling frame to become light-like, but trapped at the horizon. It is shown that in the falling frame at the horizon, the density of the source goes to infinity as a result of the length contraction.

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.

The Chebyshev collocation method (ChC) implemented as Chebfun is used in order to solve a class of second order one-dimensional singular and genuinely nonlinear boundary value problems. Efforts to solve these problems with conventional ChC have generally failed, and the outcomes obtained by finite differences or finite elements are seldom satisfactory. We try to fix this situation using the new Chebfun programming environment. However, for toughest problems we have to loosen the default Chebfun tolerance in Newton's solver as the ChC runs into trouble with ill-conditioning of the spectral differentiation matrices. Although in such cases the convergence is not quadratic the Newton updates decrease monotonically. This fact, along with the decreasing behaviour of Chebyshev coefficients of solutions, suggests that the outcomes are trustworthy, i.e., the collocation method has exponential (geometric) rate of convergence or at least an algebraic rate. We consider first a set of problems that have exact solutions or prime integrals and then another set of benchmark problems that do not possess these properties. Actually, for each test problem carried out we have determined how Chebfun solution converges, its length, the accuracy of the Newton method and especially how well the numerical results overlap with the analytical ones (existence and uniqueness).

Active flow control of canonical laminar separation bubbles by steady and harmonic vortex generator jets (VGJs) was investigated using direct numerical simulations. Both control strategies were found to be effective in controlling the laminar boundary-layer separation. However, the present results indicate that using the same blowing amplitude, harmonic VGJs were more effective and efficient in reducing the separated region than the steady VGJs considering the fact that the harmonic VGJs use less momentum than the steady case. For steady VGJs, longitudinal structures formed immediately downstream of injection location led to formation of hairpin-type vortices causing an earlier transition to turbulence. Symmetric hairpin vortices were shown to develop downstream of the forcing location for the harmonic VGJs as well. However, the increased control effectiveness for harmonic VGJs flow control strategy is attributed to the fact that shear-layer instability mechanism was exploited. As a result, disturbances introduced by VGJs were strongly amplified leading to development of large-scale coherent structures, which are very effective in increasing the momentum exchange, thus, limiting the separated region.

Amorphous ferrite-type rare-earth (RE) substituted garnets and garnet-oxide nanocomposite layers are prepared on clear glass substrates by using RF magnetron sputter-deposition process. By using a combination approach employing custom-built spectrum-fitting software in conjunction with Swanepoel’s envelope method, the spectral dispersion function of optical constants and the layer thicknesses are derived accurately from the transmission spectra of the as-deposited samples. The effects of excess metal-oxides added to the base material systems during the co-deposition process are found to affect the refractive index and the optical absorption coefficients of garnet-oxide composites. A number of optical constant datasets are presented, enabling the experimentalists to design nanophotonic or integrated-optics devices employing these functional materials.

In this paper, we present a study of various ZnO/SiO2-stacked thin film structures for flexible micro-energy harvesting devices. Two groups of micro-energy harvesting devices, SiO2/ZnO/SiO2 micro-energy generators (SZS-MGs) and ZnO/SiO2/ZnO micro-energy generators (ZSZ-MGs), were fabricated by stacking both SiO2 and ZnO thin films, and the resulting devices were characterized. With a particular interest in the fabrication of flexible devices, all the ZnO and SiO2 thin films were deposited on ITO-coated PEN substrates using an RF magnetron sputtering technique. The effects of the thickness and/or position of the SiO2 films on the device performance were investigated by observing the variations of output voltage in comparison with that of a control sample. As a result, compared to the ZnO single-layer device, all the ZSZ-MGs showed much better output voltages, while all the SZS-MG showed only slightly better output voltages. Among the ZSZ-MGs, the highest output voltages were obtained from the ZSZ-MGs where the SiO2 thin films were deposited using a deposition power of 150 W. Overall, the device performance seems to depend significantly on the position as well as the thickness of the SiO2 thin films in the ZnO/ SiO2-stacked multilayer structures, in addition to the processing conditions.

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

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.

It is shown that, for a de Sitter Universe, the Hartle-Hawking (HH) wave function can be obtained in a simple way starting from the Friedmann- Lemaitre-Robertson-Walker (FLRW) line element of cosmological equations. An oscillator having imaginary time is indeed derived starting from the Hamiltonian obtaining the HH condition. This proposes again some crucial matter on the meaning of complex time in cosmology. In order to overcome such difficulties, we propose an interpretation of the HH framework based on de Sitter Projective Holography.

In this article, we establish the existence of and then compare smallest eigenvalues for nabla fractional boundary value problems involving a fractional difference boundary condition, using the theory of u0-positive operators with respect to a cone in a Banach space.

The effects of the solvent used for the active layer materials of an OLED based on TADF emitters plays a fundamental role in solution-deposited devices. This work focuses the effects on the performance of different solvents employed to fabricate a very simple two organic layer OLED based on a green TADF emitter, under the concept of host:guest matrix. From the different results of the main figures of merit, it was possible to conclude that the OLED that used toluene as solvent for the active layer reached a maximum EQE of 14%, almost the maximum already obtained for this emitter in more complex device structures. With the analysis of the charge transport processes, it was possible to establish an explanatory model for the obtained results. Through impedance spectroscopy, additional characterization about the nature of charge transport processes was carried out. With these results, it was possible to corelate the relaxation times, with the electrical properties of the active layer, and infer about the interaction between the electrical charges and the defect levels opening new possibilities to further development in the printed OLEDs.

Quasistatic nanoindentation is a proven tool that provides information on the micromechanical behavior of wood cell walls. However, quasistatic tests cannot probe the time-dependent mechanical behavior shown by wood. Nanoindentation dynamic mechanical analysis (nanoDMA) can mesaure the viscoelastic properties of wood cell walls. This research aimed to study the quasistatic and viscoelastic properties of individual radiata pine wood (Pinus radiata D. Don) cell wall layers. To minimize variability and retrieve both properties at the same locations, a load function composed of a multiload-quasistatic function followed by dynamic reference frequency segments was developed. Nanoindentations were then performed on the S2 layer and compound corner middle lamella (CCML) of unembedded latewood cells. Because the S2 layer is anisotropic, both transverse and longitudinal-tangential wood planes were studied. In the transverse plane, the average results of the quasistatic elastic moduli (Es) for the S2 layer and CCML were 15.7 GPa and 4.6 GPa, respectively. In the longitudinal-tangential plane, the E_s was 3.9 GPa. In the transverse section, the hardness (H) of the S2 layer and CCML were 331 MPa and 277 MPa, respectively, and in the longitudinal-tangential section H was 244 MPa. To acquire the viscoelastic properties, measurements were made over more than three decades of frequency. An increase of the storage modulus (E'), and a reduction of the loss modulus (E'') and loss factor (tan⁡ δ ) as frequency increased were found in both wood orientations. The quasi-static and dynamic indentations equivalent at 0.1 Hz showed similar values for Es and E’. This study contributes to our knowledge of wood cell wall micromechanical properties.

A case study is presented which demonstrates the value and validity of a novel approach to the use of consolidated amateur (‘ham’) radio reception reports as indicators of the presence of intense ionospheric sporadic E (Es). It is shown that the use of amateur data can provide an important supplement to other techniques, allowing the detection and tracking of Es where no suitable ionosonde or other measurements are available. The effectiveness of the approach is demonstrated by reference to ionosonde data, and the advantages and limitations of the technique are discussed.

Archie model is the basis of calculating oil saturation, but there are some limitations when using this model to calculate oil saturation in water-flooded layer. The main reason is that the main parameters, such as rock resistivity and formation water resistivity, are constantly changing dynamic parameters in the underground with the influence of injected water. Considering that water-flooded layers changes with injection waters, influence factors of rock resistivity and formation water resistivity of primary parameters are analyzed. Considering the dynamic data of water cut is the most reactive underground fluid characteristics of real information, combined with dynamic and static, and the “double ratio” model of later development by the same sedimentary layer is established, which realizes the inversion of rock resistivity and formation water resistivity, then modified Archie model relating to water-flooded layers of the same sedimentary explaining formation. The explanation of actual data indicates that the “double ratio” model well considers the dynamic variation of production data, which makes the inversion of the flooded rock resistivity relatively accurate, besides, the modified Archie model can accurately calculate the oil saturation of water-flooded layers with a reasonable result, which offers scientific basis for the predicting of remaining oil distribution rules.

In this work, we study the dynamical behaviors of the electromagnetic fields and material responses in the hyperbolic metamaterial consisting of periodically arranged metallic and dielectric layers. The thickness of each unit cell is assumed to be much smaller than the wavelength of the electromagnetic waves, so the effective medium concept can be applied. When electromagnetic (EM) fields are present, the responses of the medium in the directions parallel to and perpendicular to the layers are like that of Drude and Lorentz media, respectively. We derive the energy density of the EM fields and the power loss in the effective medium based on Poynting theorem and the dynamical equations of the polarization field. We also show that the Lagrangian density of the system can be constructed. The Euler-Lagrangian equations yield the correct dynamical equations of the electromagnetic fields and the polarization field in the medium. The canonical momentum conjugates to every dynamical field can be derived from the Lagrangian density via differentiation or variation with respect to that field. We apply Legendre transformation to this system, and find that the resultant Hamitonian density is identical to the energy density, up to an irrelevant divergence term.