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

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

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

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

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

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

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

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

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

This paper deals with the study of the existence and non existence of solutions of a three parameter's family of nonlinear fractional differential equation with mixed-integral boundary value conditions. We consider the $\alpha$-Riemann-Liouville fractional derivative, with $\alpha \in (1,2]$. In order to deduce the existence and non existence results, we first study the linear equation, by deducing the main properties of the related Green's functions. We obtain the optimal set of parameters where the Green's function has constant sign. After that, by means of the index theory, the nonlinear boundary value problem is studied. Some examples, at the end of the paper, are showed to illustrate the applicability of the obtained results.

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

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

We introduce a class of logarithmic wave equation. We study the global existence of week solution for this class of equation. We deal with the initial boundary value problem of this class. Using the Galerkin method and the Gross logarithmic Sobolev inequality we establish the main theorem of existence of week solution for this class of equation arising from Q-Ball Dynamic in particular.

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

This study investigates the solution of a boundary value problem using two numerical methods: the Rayleigh-Ritz method and the finite element method (FEM). The aim is to compare the performance of these methods and assess the reliability of FEM as a generalization of the Rayleigh-Ritz approach for more complex problems. The Rayleigh-Ritz method and the linear element formulation of the finite element method were employed to solve the boundary value problem. A detailed comparison of the results obtained from both methods was performed. Graphical illustrations were used to present the solutions, and potential sources of error were analyzed, including element and domain approximation errors, round-off errors, and the impact of using linear rather than quadratic elements in FEM. The solutions generated by both methods were found to be in close agreement, demonstrating that FEM is a viable alternative to the Rayleigh-Ritz method for solving boundary value problems. The minor discrepancies observed can be attributed to approximation errors and the choice of linear elements in the finite element analysis. This work highlights the applicability and effectiveness of both the Rayleigh-Ritz method and FEM in solving boundary value problems. It underscores the finite element method’s flexibility, especially in handling more complex boundary conditions and geometries, and contributes to the understanding of the factors influencing the accuracy of numerical methods in structural analysis.

Event detection provides a crucial foundation for intelligent text analysis and supports numerous downstream applications. However, existing detection models often fail to meet the high accuracy requirements and typically require a human–machine collaborative detection mode involving “machine detection followed by manual audit and correction.” Analysis reveals that the recall rate of machine detection is a major factor limiting the efficiency of human–machine collaboration. To address this at its root and ensure high recall, an event trigger recommendation task is introduced, and a span-based recommendation model for Chinese, SpETR (Span-based Event Trigger Recommender), is developed. This model integrates confidence information from both trigger identification and event classification stages to recommend candidate triggers and their event types for each event. To enhance recall, a semantic boundary smoothing method is proposed, which assigns soft labels to spans around gold triggers based on semantic completeness and part-of-speech, combined with a modified focal loss function to make the model predictions smoother. Experiments on DuEE and a self-built dataset show that SpETR is most effective when the number of recommendations is set to 3, achieving recall rates of 97.76% and 98.70%, respectively, comprising improvements of 4.30% and 5.22% over the optimal baseline models. In terms of human–machine collaborative detection, SpETR saved 30.6% and 25.8% of time compared to traditional detection models on two datasets, and the average F1 score of the annotation results was also improved by 5.35% and 3.48%, respectively.

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

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

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

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

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

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

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

In this paper, the existence of two nontrivial stationary solutions for the nonlinear Gilpin Ayala two species competition model is given by using the mountain pass lemma, and the local stability criterion of the trivial solution is given by using Lyapunov function method. Based on the local stability criterion, we give some suggestions on how to avoid the population extinction. This is, when the population is on the verge of extinction, we should try our best to avoid the diffusion behavior and reduce the diffusion coefficient, otherwise the species are easy to go extinct. Numerical example shows the effectiveness of the proposed method.

In this work, we study a boundary control problem for the evolutionary Stokes equations, under mixed boundary conditions, in two and three dimensions. The Stokes equation is a valuable tool in cardiovascular biomechanics, especially for modeling blood flow in low-velocity and high-viscosity conditions. Its application in the microvasculature, the design of medical devices, and the study of cardiovascular diseases helps improve the understanding and treatment of various pathological conditions. We provide a comprehensive theoretical framework to address the analysis and the derivation of a system of first-order optimality conditions that characterizes the solution of the control problem. Finally, solution-finding algorithm is proposed and illustrated with some simulations.

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

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

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

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

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

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

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

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

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

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

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

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

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

We consider boundary value problems for a nonlinear mass transfer model, which generalizes the classical Boussinesq approximation, under inhomogeneous Dirichlet boundary conditions for the velocity and the substance’s concentration. It is assumed that the leading coefficients of viscosity and diffusion and the buoyancy force in the model equations depend on concentration. We develop mathematical apparatus of studying inhomogeneous boundary value problems under consideration. It is based on using of a weak solution of the boundary value problem and construction of liftings of the inhomogeneous boundary data. They remove the inhomogeneity of the data and reduce original problems to equivalent homogeneous boundary value problems. Based on this apparatus we will prove the theorem of the global existence of a weak solution to boundary value problem under study and establish important properties of the solution. In particular, we will prove the validity of the maximum principle for the substance’s concentration. We will also establish sufficient conditions for the problem data, ensuring the local uniqueness of weak solutions.

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

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

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

This paper studies an effective finite difference scheme for solving two-dimensional Heston stochastic volatility option pricing model problems. A dynamically balanced up-downwind strategy for approximating the cross-derivative is implemented and analyzed. Semi-discretized and spatially nonuniform platforms are utilized. The numerical method comprised is simple, straightforward with reliable first order overall approximations. The spectral norm is used throughout the investigation and numerical stability is proven. Simulation experiments are given to illustrate our results.

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

Meshfree particle method, which is always regarded as a pure Lagrangian approach, is easily represented complicated domain topologies, moving boundaries, and multiphase media. Solving acoustic problems with the mesfree particle method forms a branch of the acoustic wave modeling field, namely, particle-based computational acoustics (PCA). The aim of this paper is to improve the accuracy of using the PCA method to solve two-dimensional acoustic problems, and realize the particle representation with a hybrid meshfree and finite-difference time-domain (FDTD) method for acoustic boundary conditions at both the plane and curved surface. As a widely used Lagrangian meshfree method, the smoothed particle hydrodynamics (SPH) based on the support domain and the kernel function has developed rapidly in recent years. The traditional SPH method is easily implements parallel processing and has been applied in sound wave simulation. As a corrective method with higher accuracy than SPH, the acoustic propagation and scattering in the time domain is simulated with the corrective smoothed particle method (CSPM). Moreover, a hybrid meshfree-FDTD boundary treatment technique is utilized to represent different acoustic boundaries in the Lagrangian approach. In this boundary treatment technique, the parameter value of virtual particles is obtained with the FDTD method, which concerns truncation errors based on the Tayler series expansion. Soft, rigid, and Mur’s absorbing boundary conditions are developed to simulate sound waves in finite and infinite domain. Results of modeling acoustic propagation and scattering show that CSPM is accurate and convergence with exact solutions, and different acoustic boundaries are validated to be effective in the computation.

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

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

Transient conjugate heat transfer measurements under varying temperature and velocity inlet boundary conditions at incompressible flow conditions were performed for flat plate and ribbed channel geometries. Therefrom local adiabatic wall temperatures and heat transfer coefficients are determined. Those data are analyzed using typical heat transfer correlations, e.g. Nu=CRe^mPr^n determining the local distribution of C and m. It is shown that they are closely linked. A relationship lnC=A−mB is observed with A and B as modeling parameters. They could be related to parameters in log law or power law representations for turbulent boundary layer flows. The parameter m is shown to have a close link to local pressure gradients and therewith near wall streamlines as well as friction factor distributions. A normalization of the C parameter allows to derive a Reynolds analogy factor and therefrom local wall shear stresses.

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.

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

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

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

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

Modeling fluid flow in three-dimensional (3D) Discrete Fracture Networks (DFNs) is of relevance in many engineering applications, such as hydraulic fracturing, oil/gas production, geothermal energy extraction, nuclear waste disposal and CO₂ sequestration. A new Boundary Element Method (BEM) technique with discontinuous quadratic elements and a parallel Domain Decomposition Method (DDM) is presented herein for the simulation of the steady-state fluid flow in 3D DFN systems with wellbores, consisting of planar fractures having arbitrary properties and wellbore trajectories. Numerical examples characterized by DFNs of increasing complexity are investigated to evaluate the accuracy and efficiency of the presented technique. The results show that accurate solutions can be obtained with less nodes than with mesh-based methods (e.g. Finite Element Method). In addition, the DDM algorithm used provides a quite fast convergence. The simulation results of the fluid flow around intersections among traces (linear intersections between fractures), intersections between traces and a fracture boundaries, and wellbore intersections is accurate. Source code is available at : https://github.com/BinWang0213/PyDFN3D.

Soil-structure interaction (SSI) leads to a modification of dynamic properties of a structure, but due to the complexity on analysis, it is traditionally assumed in seismic designs that the structure is fixed-supported at the ground, which unavoidably gives potential dangers to seismic performances of the structure. This study focuses on the possibility of the detrimental effects of SSI by comparing the dynamic responses of SSI system to fixed-base structure. The SSI model consists of an oscillator, a three-dimensional foundation and a multilayered half-space. In indirect boundary element method is employed, combined with non-singular Green’s functions of distributed loads to calculate the foundation impedances for substructure method. Although the SSI mechanism dampens peak amplitude of structure response in frequency domain, case studies on four buildings to 42 earthquakes in time history show a possibility of 15%-20% that SSI amplifies the dynamic responses, depending on the frequency component of incident waves. Also, a neural network algorithm is proposed to provide a new idea of involving SSI analysis into the traditional seismic designs, to eliminate potential hazards of SSI.

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

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

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

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

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

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

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

For studying of the finite-size behavior of the Ising model under different boundary conditions we propose an alternative to the standard transfer matrix technique approach based on Abel\`{e}s theorem and Chebyshev polynomials. Using it we easily reproduce the known for periodic boundary conditions results for Lee-Yang zeros, the exact position space renormalization group transformation, etc., and extend them deriving new results for antiperiodic boundary conditions. Note that in the latter case one has a nontrivial order-parameter profile, which we also calculate, where the average value of a given spin depends on the distance from the seam with opposite bond in the system. It is interesting to stress that under both boundary conditions, the one-dimensional case exhibits Schottky anomaly.

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

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

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 aim of this paper is the analysis of arbitrarily loaded isotropic rectangular Kirchhoff plates supported at all corner points, whereby the edges are arbitrarily supported (simply supported, clamped, or free). Analytical methods commonly used include the single trigonometric series of Lévy, the double trigonometric series of Navier, the power series expansions, etc. In this paper the deflection surface was approximated with the sum of a particular solution to the governing differential equation (GDE) and two single series. The terms of the single series were the product of an unknown function of an independent variable and a trigonometric function of the other independent variable, whereby the trigonometric functions were chosen to vanish at the corner points. On the one hand the terms of the series were required to satisfy the homogeneous GDE, leading to two uncoupled differential equations, one for each unknown function, and so the approximate solution satisfied exactly the GDE. On the other hand the boundary conditions were satisfied only at selected collocation points along the boundary, the number of collocation points in each direction corresponding to the number of terms of the associated series. Results for several plates supported at all corner points with arbitrary edge support conditions were presented and showed good agreement with the exact results, the accuracy increasing the more terms of series were considered. Cantilevered plates and plates having two adjacent edges free will be analyzed using this approach in future research.

A diffusion welding of coarse-grained and ultrafine-grained (UFG) specimens of titanium near-α alloy Ti-5Al-2V used in nuclear power engineering was made by Spark Plasma Sintering. The failure of the welded specimens in the conditions of hot salt corrosion and of electrochemical corrosion were shown to have preferentially intercrystalline character. In the case of presence of macrodefects, crevice corrosion of the welded joints was observed. The resistance of the alloys against the intercrystalline corrosion was found to be determined by the concentration of vanadium at the titanium grain boundaries, by the size and volume fraction of the β-phase particles and by the presence of micro- and macropores in the welded joints. The specimens of the welded joints of the UFG alloy have higher hardness, hot salt corrosion resistance and the electrochemical corrosion.

The planetary boundary layer height (PBLH) is one of the key factors in influencing the dispersion of the air pollutants in the troposphere and hence the air pollutant concentration on ground level. For this reason, accurate air pollutant concentration depends on the performance of PBLH prediction. Recently, ceilometer, a lidar instrument to measure cloud base height, has been used by atmospheric scientists and air pollution control authorities to determine the mixing level height (MLH) in improving forecasting and understanding the evolution of aerosol layers above ground at a site. In this study, ceilometer data at an urban (Lidcombe) and a rural (Merriwa) location in the New South Wales, Australia was used to validate the PBLH prediction from two air quality models (CCAM-CTM and WRF-CMAQ) as well as to understand the aerosol transport from sources to receptor point at Merriwa for the three case studies where high PM10 concentration was detected in each of the three days. The results show that the PBLH prediction by the two air quality models corresponds reasonably well with observed ceilometer data and the cause and source of high PM10 concentration at Merriwa can be found by using ceilometer MLH data to corroborate with back trajectory analysis of transport of aerosols to the receptor point at Merriwa. Of the three case studies, one had aerosol source from north and north west of Merriwa in remote NSW where windblown dust is the main source, and the other two had sources from south and south east of Merriwa where anthropogenic sources dominate,

This work implements a stable algorithm in Field-Programmable Gate Arrays (FPGAs) for which two architectures (unrolling the loop) were developed and analyzed. The algorithm recovers sources located at the boundary separating two homogeneous media that make up a two-dimensional non-homogeneous region from measurements on the boundary of such region. The problem of recovering these sources is an ill-posed inverse problem as small errors in the measurement can produce important changes in the source location. Inverse source problems have many applications in different areas, such as engineering and medicine, making the proposed implementation important. The first architecture (mode one) allows considering different operating speeds, which is an advantage depending on whether we work with fast or slow signals. The second one (mode two) reduces resource consumption by exploiting the characteristics of the source identification algorithm, which is an advantage for multichannel problems such as inverse electrocardiography or electroencephalography. The architectures were tested on four FPGAs of the 7 Series of Xilinx: Spartan-7 xc7s100fgga484, Virtex-7 xc7v585tffg1157, Kintex-7 xc7k70tfbg484, and Artix-7 xc7a35tcpg236. The two FPGA implementations of the algorithm were validated using synthetic examples implemented in MATLAB. The results presented here can be extended to concentric spheres and complex geometries.

In this article, the multitouch option, also called the n  touch option (or the “baseball” option when n  3 ) is analyzed and valued in closed form. This is a kind of barrier option that sets a gradual knockout / knock-in mechanism based on the number of times the underlying asset has crossed a predefined barrier in various time intervals before expiry. The higher the number of predefined time intervals during which the barrier has been touched, the lower the value of a knock-out contract at expiry, and conversely for a knock-in one. Multitouch options can be viewed as an extension of step barrier options, preserving the ability of the latter to adjust the exposure to risk over time, while eliminating the notorious danger of “sudden death” that holders of step barrier options are faced with. Unlike occupation time derivatives, the payoff at expiry does not depend on the amount of time spent outside the authorized range, but on the number of passages beyond the authorized range.

The study of long-term average weather patterns is known as climatology. It is a distinct field of study from meteorology and can be broken down into several subdivision. In order to predict the future, the knowledge of climatology is essential. In other words, with the help of climatology, we can figure out how likely it is that snow and hail will fall to the ground, and how much solar thermal radiation can reach a certain location etc. Climatology often focuses on how the climate has changed over time and how it has affected people and events. Both meteorology and climatology fall under the general term "meteorology", in particular, they are subdivision of research in the same field. In case of predicting the weather, meteorologists use variables such as humidity, air pressure, and temperature. This article's primary objective is to familiarize engineers with the fundamentals of climate and its processes so that they can effectively apply this knowledge to comprehend the climatic impact on water resources systems.

A numerical study for the effect of crest width, breaking parameter and trunk permeability on hydrodynamics and flow behavior in the vicinity of rubble-mound, permeable, zero-freeboard breakwaters (ZFBs) is presented. The modified two dimensional Navier-Stokes equations for two-phase flows in porous media with a Smagorinsky model for the subgrid scale stresses were solved numerically. An immersed-boundary/level-set method was used. The numerical model was validated for the cases of wave propagation over a submerged impermeable trapezoidal bar and over a low-crested permeable breakwater. Five cases of breakwaters were examined, and the main results are: (a) The size of the crest width, B, does not notably affect the wave reflection, vorticity and currents in the seaward region of ZFBs, while wave transmission, currents in the leeward side, and mean overtopping discharge, all decrease with increasing B. A non-monotonic behavior of the wave setup is also observed. (b) As the breaking parameter decreases, wave reflection, transmission, currents, mean overtopping discharge, and wave setup decrease. This observation is also verified by relevant empirical formulas. (c) As the ZFB trunk permeability decreases, an increase of the wave reflection, currents, wave setup, and a decrease of wave transmission and mean overtopping discharge is observed.

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

The article proposes various new approximate analytical solutions of the boundary value prob-lem for the non-stationary system of Nernst-Planck-Poisson (NPP) equations in the diffusion layer of an ideally selective ion-exchange membrane at overlimiting current densities. As is known, the diffusion layer in the general case consists of a space charge region and a region of local electroneutrality. The proposed analytical solutions of the boundary value problems for the non-stationary system of Nernst-Planck-Poisson equations are based on the derivation of a new singularly perturbed nonlinear partial differential equation for the potential in the space charge region (SCR). This equation can be reduced to a singularly perturbed inhomogeneous Burgers equation, which, by the Hopf-Cole transformation, is reduced to an inhomogeneous singularly perturbed linear equation of parabolic type. Inside the extended SCR, there is a suffi-ciently accurate analytical approximation to the solution of the original boundary value prob-lem. The electroneutrality region has a curvilinear boundary with the SCR, and with an un-known boundary condition on it. The article proposes a solution to this problem. The new ana-lytical solution methods developed in the article can be used to study non-stationary boundary value problems of salt ion transfer in membrane systems.

Ever since Varela and Maturana proposed the concept of autopoiesis as the minimal requirement for life, there has been a focus on cellular systems that erect topological boundaries to separate themselves from their surrounding environment. Here we reconsider whether the existence of such a spatial boundary is strictly necessary for self-producing entities. This work presents a novel computational model of a minimal autopoietic system inspired by dendrites and molecular dynamics simulations in three-dimensional space. A series of simulation experiments where the metabolic pathways of a particular autocatalytic set are successively inhibited until autocatalytic entities that could be considered autopoietic are produced. These entities maintain their distinctness in an environment containing multiple identical instances of the entities without the existence of a topological boundary. This gives rise to the concept of a metabolic boundary which manifests as emergent self-selection criteria for the processes of self-production without any need for unique identifiers. However, the adoption of such a boundary comes at a cost, as these autopoietic entities are less suited to their simulated environment than their autocatalytic counterparts. Finally, this work showcases a generalized metabolism-centered approach to the study of autopoiesis that can be applied to both physical and abstract systems alike.

This research program investigated the effects of layer thickness (50 and 100 microns) on the microstructure and mechanical properties of electron beam powder bed fusion (EBPBF) additive manufacturing of Inconel 625 alloy. The as-built 50 and 100 micron layer thickness components were also heat treated at temperatures above 1100 oC, which produced a recrystallized grain structure containing annealing twins in the 50 micron layer thickness components, and a duplex grain structure consisting of islands of very small equiaxed grains dispersed in a recrystallized, large-grain structure containing annealing twins. The heat treated component microstructures and mechanical properties were compared with the as-built components in both the build direction (vertical) and perpendicular (horizontal) to the build direction. Vickers microindentation hardness (HV) values for the vertical and horizontal geometries averaged 227 and 220 for the as-built 50 and 100 micron layer components, and 185 and 282 for the corresponding heat treated components. The yield stress values were 387 MPa and 365 MPa for the as-built layer horizontal and vertical 50 micron layer geometries, and 330 MPa and 340 MPa for the as-built 100 micron layer components. For the heat treated 50 micron components, the yield stress values were 340 and 321 MPa for the horizontal and vertical geometries, and 581 and 489 MPa for the 100 micron layer components, respectively. The elongation for the 100 micron layer as-built horizontal components was 28% in contrast to 65% for the corresponding 100 micron heat treated layer components, an increase of 132% for the duplex grain structure. However, the coarse grains containing annealing twins and the equiaxed small grain islands in the duplex structure for the heat treated components contained continuous carbides in the grain boundaries, and this may indicate sensitization and a reduction in corrosion resistance. These findings point to the potential mechanical property advantage for heat treatment of Inconel 625 alloy 100 micron layer thickness components fabricated by EBPBF.

Height assignment is an important problem for satellite measurements of Atmospheric Motion Vectors (AMVs) that are interpreted as winds by forecast and assimilation systems. Stereo methods assign heights to AMVs from the parallax observed between observations from different vantage points in orbit while tracking cloud or moisture features. In this paper, we fully develop the stereo method to jointly retrieve wind vectors with their geometric heights from geostationary satellite pairs. Synchronization of observations between observing systems is not required. NASA and NOAA stereo-winds codes have implemented this method and we have processed large datasets from GOES-16, -17, and Himawari-8. Our retrievals are validated against rawinsonde observations and demonstrate the potential to improve forecast skill. Stereo winds also offer an important mitigation for the loop heat pipe anomaly on GOES-17 during times when warm focal plane temperatures cause infra-red channels that are needed for operational height assignments to fail. We also examine several application areas, including deep convection in tropical cyclones, planetary boundary layer dynamics, and fire smoke plumes, where stereo methods provide insights into atmospheric processes. The stereo method is broadly applicable across the geostationary ring where systems offering similar Image Navigation and Registration (INR) performance as GOES-R are deployed.

We extensively studied serigraphic screen-printed commercial silver flake inks loaded with silicon inclusions in order to achieve pinning at the grain boundaries. Based on grain size measurements using electron backscattered diffraction (EBSD), commercial silver ink with silicon microparticle content of 5 wt.% shows a significant grain growth retardation compared to the pristine silver ink, which stabilizes electrical conductivity up to 700 ◦C via Zener pinning mechanism. The modified silicon-loaded silver ink experiences a two times increase in grain size when heated up to 700 ◦C, compared to a seven times increase for the pristine silver ink. In turn, this enables operation temperatures significantly higher than the conventional operational window of microparticle-based silver inks, which are usually limited to 400 ◦C. Using isothermal exposures of 10 minutes up to 4 hours this phenomenon is observed at temperatures ranging from 250 ◦C to 900 ◦C. The electrical conductivity stability, grain size evolution and oxide contents were studied up to 4 hours. The activation energy of silver ink with silicon inclusions is 54% lower than for the pristine silver ink due to pining effect, which retards grain growth via the Zener mechanism. Most importantly, the electrical resistivity becomes stable up to 700 ◦C, which is more than twice the operation limit for off-the-shelf screen-printable silver flake inks. Hence, we demonstrate that adding controlled amounts of silicon particles to silver inks towards grain structure engineering can open new vistas of possibilities for screen-printed metallic inks.

In this paper, we study nonlocal problems for a fractional diffusion equation and a degenerate hyperbolic equation with singular coefficients in the lower terms. The uniqueness of the solution to the problem is proved by the method of energy integrals. The existence of a solution is equivalently reduced to the question of the solvability of Volterra integral equations of the second kind and a fractional differential equation. For a particular solution of the proposed problem, its visualization is carried out for various values of the order of the fractional derivative. It is shown that the order of the derivative affects the intensity of the diffusion process (subdiffusion), as well as the shape of the wave front.

This article provides exact analytical formulae for various kinds of rainbow step barrier options. These are highly flexible and sophisticated multi-asset barrier options based on the following principle : the option life is divided into several time intervals on which different barriers are monitored w.r.t. different underlying assets. From a mathematical point of view, new formulae are provided for the first passage time of a multidimensional geometric Brownian motion to a boundary defined as a step function. The article shows how to implement the obtained valuation formulae in a simple and very efficient manner. Numerical results highlight a strong sensitivity of rainbow step barrier options to the correlations between the underlying assets.

The reason for excessive multidecade deflections is still unclear after 70 years of the cantilever method that has been applied in building large-span segmental concrete bridges. The failure of the Koror-Babeldaob bridge gives the author a clue that the reason might be the boundary condition. The asymmetric properties of the bearing and prestress tendon draw the author’s attention. To understand the process of stress redistribution, the virtual stiffness method which is derived from the principle of minimum potential energy was created, and the phenomenon the author suggested named Material Inelasticity and Prestress Loss Coupling Effect is discovered. The unexplainable excessive deflection and prestress loss are found in the new load model. The calculation of excessive deflection is discussed. The new model fits the experience well. The optimization design methods and construction key points are proposed.

It is known that the elastic characteristics of polyamides change with increasing temperature, in particular, the Young's modulus decreases significantly. This fact is practically not taken into account in design calculations of metal-polymer plain (MP) bearings, operating under conditions of the boundary and dry friction. The purpose of the study is the analysis of the effect of temperature on the change of the Young's modulus and, accordingly, the contact strength and triboresource according to the developed method of calculating MP bearings. MP bearings with a bushing made of polyamide PA6 reinforced with glass or carbon dispersed fibers were investi-gated. Quantitative and qualitative regularities of change of the maximum contact pressures and resource of the bearings at temperature increase under conditions of boundary and dry friction are established. The pressures in the bearing bushing made of PA6+30GF will be lower than for the bushing made of PA6+30CF. The resource of the bushing made of PA6+30CF will be significantly greater than for PA6+30GF. For thermoplastic polymers, the increase in temperature will have a useful practical effect due to the decrease in the rigidity of the polymer composites of the bearing bushing.

The stele concept is one of the oldest enduring concepts in plant biology. This paper reviews the concept and its foundations, and builds an argument for an updated view of steles and their evolution. The history of studies of stelar organization has generated a widely ranging array of definitions of the stele that determine the way we classify steles and construct scenarios about the evolution of stelar architecture. Because at the level of the organism biological evolution proceeds by, and is reflected in, changes in development, concepts of structure need to be grounded in development in order to be relevant in an evolutionary perspective. For the stele, most of the traditional definitions that incorporate development have viewed it as the totality of tissues that either originate from procambium – currently the prevailing view – or are bordered by a boundary layer (e.g., endodermis). A definition of the stele that would bring consensus between these perspectives recasts the stele as a structural entity of dual nature. Here, I review briefly the history of the stele concept, basic terminology related to stelar organization, and traditional classifications of the steles. I then revisit boundary layers from the perspective of histogenesis as a dynamic mosaic of developmental domains. I use classic and recent anatomical and molecular data to reaffirm and explore the importance of boundary layers for stelar organization. Drawing on data from comparative anatomy, developmental regulation, and the fossil record, I offer a model for a stele concept that integrates both the boundary layer and the procambial perspective, consistent with a dual nature of the stele. The dual stele model posits that stelar architecture is determined in the apical meristem by two major cell fate specification events: a first one that specifies a provascular domain and its boundaries, and a second event that specifies a procambial domain (which will mature into conducting tissues) from cell subpopulations of the provascular domain. If the position and extent of the developmental domains defined by the two events are determined by different concentrations of the same morphogen (most likely auxin), then the distribution of this organizer factor in the shoot apical meristem, as modulated by changes in axis size and the effect of lateral organs, can explain the different stelar configurations documented among tracheophytes. This model provides a set of working hypotheses that incorporate assumptions and generate implications that can be tested empirically. The model also offers criteria for an updated classification of steles that is in line with current understanding of plant development. In this classification, steles fall into two major categories determined by the configuration of boundary layers – boundary protosteles and boundary siphonosteles, each with subtypes defined by the architecture of the vascular tissues. Validation the dual stele model and, more generally, in-depth understanding of the regulation of stelar architecture, will necessitate targeted efforts in two areas: (i) the regulation of procambium, vascular tissue, and boundary layer specification in all extant vascular plants, considering that most of the diversity in stelar architecture is hosted by seed-free plants, which are the least explored in terms of developmental regulation; (ii) the configuration of vascular tissues and, especially, boundary layers, in as many extinct species and lineages as possible.

With the growth of air traffic demand in busy airspace, there is an urgent need for airspace sectorization to increase air traffic throughput and ease the pressure on controllers. The purpose of this paper is to develop a method framework that can perform airspace sectorization automatically, reasonably, which can be used as an advisory tool for controllers as an automatic system, especially for eliminating irregular sector shapes generated by simulated annealing algorithm (SAA) based on region growth method. Two graph cutting method, dynamic Monte Carlo method by changing location of flexible vertices (MC-CLFV) and Monte Carlo method by radius changing (MC-RC) were developed to eliminating irregular sector shapes generated by SAA in post-processing. The experimental results show that the proposed method framework of AS can automatically and reasonably generate sector design schemes that meet the design criteria. Our methodology framework and software can provide assistant design and analysis tools for airspace planners to design airspace, improve the reliability and efficiency of airspace design, and reduce the burden of airspace planners. In addition, this lays the foundation for reconstructing airspace with more intelligent method.

Goals: To use visual logistics for interpreting COVID-19 molecular and rapid antigen test (RAgT) performance, determine prevalence boundaries where risk exceeds expectations, and evaluate benefits of recursive testing along home, community, and emergency spatial care paths. Methods: Mathematica/open access software helped graph relationships, compare performance patterns, and perform recursive computations. Results: Tiered sensitivity/specificity comprise: T1) 90%/95%; T2) 95%/97.5%; and T3) 100%/≥99%, respectively. In emergency medicine, median RAgT performance peaks at 13.2% prevalence, then falls below T1, generating risky prevalence boundaries. RAgTs in pediatric ERs/EDs parallel this pattern with asymptomatic worse than symptomatic performance. In communities, RAgTs display large uncertainty with median prevalence boundary of 14.8% for 1/20 missed diagnoses, and at prevalence >33.3-36.9% risk 10% false omissions for symptomatic subjects. Recursive testing improves home RAgT performance. Home molecular tests elevate performance above T1, but lack adequate validation. Conclusions: Widespread RAgT availability encourages self-testing. Asymptomatic RAgT and PCR-based saliva testing present the highest chance of missed diagnoses. Home testing twice, once just before mingling, and molecular-based self-testing help avoid false omissions. Community and ER/ED RAgTs can identify contagiousness in low prevalence (<22%). Real-world trials of performance, cost-effectiveness, and public health impact could identify home molecular diagnostics as the optimal diagnostic portal.

The Review examines studies aimed at the organization of energy (non-mechanical) control of high-speed flow/flight using spatially multi-component plasma structures and combined energy deposition. The Review covers selected works on experimental obtaining and numerical modeling multi-component plasma structures and use of sets of actuators based on the plasma of such a spatially type for the purposes of control of shock wave/bow shock wave – energy source interaction, as well as control of shock wave – boundary layer interaction. A series of works on repetitive multiple laser pulse plasma structures is also analyzed from the point of view of examining shock wave/bow shock wave boundary layer interaction. Self-sustained theoretical models for laser dual-pulse, multi-mode laser pulses and self-sustained glow discharge are also considered. Separate sections are devoted to high-speed flow control using combined physical phenomena and numerical prediction of flow control possibilities using thermal longitudinally layered plasma structures. The wide possibilities for organization and applying spatially multi-component structured plasma for the purposes of high-speed flow control are demonstrated.

Abstract: A prevalence boundary (PB) marks the point in prevalence where the false omission rate, RFO=FN/(TN+FN), exceeds the tolerance limit for missed diagnoses. Objectives were to mathematically analyze rapid antigen test (RAgT) performance, determine why PBs are breeched, and evaluate the merits of testing three times over five days, now required by the US Food and Drug Administration for asymptomatic persons. Equations were derived to compare test performance patterns, calculate PBs, and perform recursive computations. An independent July 2023 FDA-university-commercial evaluation of RAgTs provided performance data used in theoretical calculations. Tiered sensitivity/specificity comprise: Tier-1) 90%, 95%; Tier-2) 95%, 97.5%; and Tier-3) 100%, ≥99%, respectively. Repeating a T2 test improves the PB from 44.6% to 95.2% (RFO 5%). In the FDA-university-commercial evaluation, RAgTs generated sensitivity of 34.4%, which improved to 55.3% when repeated, then 68.5% with the third test. With RFO=5%, PBs were 7.37/10.46/14.22%, respectively. PB analysis suggests RAgTs should achieve clinically proven sensitivity of 91.0-91.4%. When prevalence exceeds PBs, missed diagnoses can perpetuate virus transmission. Repeating low-sensitivity RAgTs delays diagnosis. In homes, high-risk settings, and hotspots, PB breaches may prolong contagion, defeat mitigation, facilitate new variants, and transform outbreaks into endemic disease. Molecular diagnostics can help avoid these potential vicious cycles