We firstly present two popular space discretization models of the Nonlinear Schrödinger Equation in dark solitons motion : the Direct-Discrete model and the Ablowitz-Ladik model. Applying the midpoint scheme to the space discretization models, we get two time-space discretization models: the Crank-Nicolson method and the New-Difference method. Secondly, we demonstrate that the solutions of the two space discretization models converge to the solution of the Nonlinear Schrödinger Equation. Also, we prove that the convergence order of the two time-space discretization models are O(h2+τ2) in discrete L2-norm error estimates. Finally, the numerical experiments agree well with the proven theoretical results.

Globally, saline lakes occupying 23% by area 44% by volume among all the lakes might desiccate by 2025 due to agricultural diversion, illegal encroachment, pollution, and invasive species. India’s largest saline lake, Sambhar is currently shrinking at the rate of 4.23% due to illegal saltpan en-croachment. This research article aims to identify the trend of migratory birds and monthly wetland status. Birds survey was conducted for 2019, 2020 and 2021 and combined with literature data of 1994, 2003, and 2013 for visiting trend, feeding habit, migratory and resident ratio, and ecological diversity index analysis. Normalized Difference Water Index was scripted in Google Earth Engine. Results state that it has been suitable for 97 species. Highest NDWI values for the was whole study period was 0.71 in 2021 and lowest 0.008 in 2019 which is highly fluctuating. The decreasing trend of migratory birds coupled with decreasing water level indicates the dubious status for its existence. If the causal factors are not checked, it might completely desiccate by 2059 as per its future prediction. Certain steps are suggested that might help conservation. Least, the cost of restoration might exceed the revenue generation.

We study the center of energy (CE) before and after the separation of superposed wave from a moving medium (MM). It is assumed that two out-of-phase mechanical transverse waves propagating from the opposite directions on a medium moving at non-relativistic speeds are superposed and the superposed portion (SP) is separated from the MM at that moment. We consider the CE of the SP before and after the separation from the MM. The location of CE (LCE) of the SP seems to be at the center of it at the moment of superposition. The SP rotates due to the separation from the MM since the velocity of each portion symmetric with respect to the center of the SP is equal in magnitude and opposite in direction. The magnitudes of velocities of the symmetric portions become different as soon as the SP begins to rotate with the separation from the MM. Then the energies of their symmetric portions are not the same, so the LCE of the SP is not at the center of it. As a result, the LCE of it looks different before and after the separation from the MM. We must find a solution to keep the LCE of the SP constant. We propose that two out-of-phase mechanical transverse waves (MTWs) propagating from the opposite directions on a MM originally have hidden difference in relativistic energy and it suddenly appears within the range observable in Newtonian mechanics when the SP starts rotating. This means that the point of view of hidden difference in relativistic energy is necessary to keep the LCE constant.

In this paper, initially a mathematical model is formulated for transient frequency of power system considering time delays which occur while transmitting the control signals in open communication infrastructure. Time delay negligence in a power system leads to improper measurement of frequency variation in power system. The study of impact of time delays on the stability of power system is performed by estimating the decay rate of frequency wave form using Kalman Filter (KF). In power system, there is a possibility of multiple time delays. This paper also focusses on developing Interacting Multiple Model (IMM) Algorithm with multiple model space using Kalman Filter (KF) as state estimator tool. The multiple time delays in power system is considered as multiple model space. The result shows that KF provides better estimate of correct model for a particular input-set. The qualitative properties of Riccati difference equation (RDE) in terms of state error covariance of IMMKF are also analyzed and presented.

The surface enhanced fluorescence（SEF）detection bases by plasmonic nanopillars array with nanoparticles has opened up a new gate in the application of biological imaging and sensing. The fluorescence enhancement of the probe molecule depends on its position in equilibrium, which is close to the hot spot leading to the electromagnetic field enhancement, but not too close to the metal surface resulting in quenching. Here, a large scale SiO₂-Ag-cicada wing SEF substrate was fabricated by magnetron sputtering with correction enhancement factor of 797.6. Thereinto the cicada wing provides the skeleton of the nanopillars array structure, the deposited Ag constructs two kinds of hot spots, and SiO₂ forms a separation layer to prevent quenching. Moreover, the substrate exhibited good reproducibility, high sensitivity with low limits of detection (LOD) and high stability for oxidation resistance. We propose that SEF substrate with modification of SiO₂ can not only improve the enhancement performance, but also expanding its application in the biological investigations.

Giant plasmonic surface enhancement has been observed in gold coated micron sized inverse pyramids entrapping a gold nanoparticle. The amplification of both surface enhanced Raman and photoluminescence signals was found to be dependent on the diameter of trapped gold nanoparticle and around 50-fold enhancement was detected for 250nm diameter sample relatively to the 50nm one. Finite differential time domain simulations, performed to determine the near-field distribution in the structure, showed that when the nanoparticle protrudes into the hotspot zone of the void, coupling of electromagnetic field occurs and the plasmon-related near-field enhancement is concentrated into the close vicinity of the nanoparticle, mainly into the close gaps around the tangential points of the curved sphere and the flat pyramid surface. This results in a more than 15 times increase of the near-field intensity, compared to the empty void.

The fractional telegraph partial differential equation with fractional Atangana-Baleanu-Caputo (ABC) derivative is studied. Laplace method is used to find the exact solution of this equation. Stability inequalities are proved for the exact solution. Difference schemes for the implicit finite method are constructed. The implicit finite method is used to deal with modelling the fractional telegraph differential equation defined by Caputo fractional of Atangana-Baleanu (AB) derivative for different interval. Stability of difference schemes for this problem is proved by the matrix method. Numerical results with respect to the exact solution confirm the accuracy and effectiveness of the proposed method.

Seaweeds are regarded as one of the valuable coastal resources because of their usage in human food, cosmetics, and other industrial items. They also play a significant role in providing nourishment, shelter, and breeding grounds for fish and many other sea species. This study introduces a newly developed seaweed enhancing index (SEI) using spectral bands of near-infrared (NIR) and shortwave infrared (SWIR) of Landsat 8 satellite data. The seaweed patches in the coastal waters of Karachi, Pakistan were mapped using SEI, and its performance was compared with other commonly used indices - Normalized Difference Vegetation Index (NDVI) and Floating Algae Index (FAI). The accuracy of the mapping results obtained from SEI, NDVI, and FAI was checked with field verified seaweed locations. The purpose of the field surveys was to validate the results of this study and to evaluate the performance of SEI with NDVI and FAI. The performance of SEI was found better than NDVI and FAI in enhancing submerged patches of the seaweed pixels what other indices failed to do.

Earthquake-triggered landslides represent a significant seismic-related disaster, posing threats to both the lives and properties of individuals in affected areas. Furthermore, they can result in road and river blockages, as well as other secondary disasters, significantly impacting post-earthquake rescue efforts. Efficient and accurate rapid assessment of high-risk landslide zones carries important implications for decision-making in disaster response and for mitigating potential secondary disasters. This study uses the high-intensity Zone VII to IX of the 5 September 2022 Luding Ms6.8 earthquake as a case study. The "difference method" and the "cumulative displacement method," based on the simple Newmark model, are employed to assess earthquake-triggered landslides. The evaluation results from both methods demonstrate that the areas posing an ex-tremely high risk of earthquake-triggered landslides are predominantly situated on the western side of the Xianshuihe Fault. Actual landslide data verification shows that both models have high predictive accuracy, with the difference method model slightly outper-forming the cumulative displacement method model. Moreover, the study recommends determining threshold values for each landslide hazard interval with significant physical meaning using past strong earthquake data as soon as possible when utilizing the difference method to assess the risk of earthquake-triggered landslides.

Under the action of cement hydration heat, construction environment, thermal insulation measures, and pipe cooling system, the mass concrete pile cap produces a complex temperature field inside the structure, which makes it difficult to control the internal surface temperature difference, the maximum adiabatic temperature rise and the surface temperature of the structure. In this paper, the mass concrete pile cap of a super large bridge (the length, width, and height are 26.40m，20.90m, and 5.00m, respectively, and the central pier pile cap is constructed with concrete with a strength of C40) is taken as the research support project. The controlling factors affecting the temperature field of the mass concrete pile cap are found by comparing the temperature monitoring of the pile cap with the calculated value of finite element software simulation analysis. Through the finite element software Midas Civil and Midas FEA, the control factors affecting the temperature field of the main pier cap (concrete mold temperature, concrete surface convection coefficient, ambient temperature, and pipe cooling system parameters, etc.) are numerically analyzed, and the influence law and influence degree of the control factors are found out.

Spatial and temporal variation of the free electron concentration in the ionosphere affects SAR interferograms, in particular at low radar frequencies. In this work the identification, estimation and compensation of ionospheric path delay phases in PALSAR-3 and NISAR-L interferograms is discussed. Both these L-band sensors acquire simultaneously SAR data in a main spectral band and in an additional, spectrally separated, narrower second band to support the mitigation of ionospheric path delays. The methods presented permit separating the dispersive and the non-dispersive phase terms based on the double-difference interferogram between the two available spectral bands and the differential interferogram of the main band. The applicability of the proposed methods is demonstrated using PALSAR-3-like data that was simulated based on PALSAR-2 SM1 mode data.

The Black population, often treated as ethnically homogenous, faces a constant challenge in accessing and utilizing healthcare services. This study examines the intragroup differences in health-seeking behavior among the multiple ethnic subgroups that comprise the Black community. A cross-sectional study was conducted among 239 adults ≥18 years who self-identified as Black in the United States and Canada. A multiple logistic regression model was fitted to evaluate the relationship between health-seeking behaviors and ethnic origin controlling for selected social and health related factors. The mean age of participants was 38.6 years, 31% were male, and 20% were unemployed. Sixty-one percent reported very good or excellent health status, and 59.7% reported not being treated for chronic disease at the time of the survey. Aging (OR=1.05, CI:1.01–1.09) and being female (OR= 0.30, CI: 0.14–0.65) were associated with favorable health-seeking behaviors. Compared with participants who had earned high school diplomas only, those who reported having earned graduate degrees were 75% less likely to have engaged in favorable health-seeking behavior (OR=0.25, 95% CI: 0.07–0.86). In the Black community of our sample, regardless of ethnic origins or country of birth, the factors associated with health-seeking behavior were age, gender, and educational attainment. Keywords: African American; African; Caribbean; Healthcare Utilization; Ethnic Disparities; Health-seeking behavior; intragroup difference, black populations, ethnic subgroups

Bhutan located in the Hindu Kush Himalayan (HKH) region consists of several glaciers and glacial lakes at higher elevations.  With the rapid change in global temperature, glaciers are found to melt at an accelerated rate. This rapid melt of glaciers gets accumulated in weak moraine walls forming a glacial lake, posing a major threat to the downstream communities. As per the Bhutan Glacial Lake Inventory 2021, Bhutan has 567 glacial lakes. Furthermore, the Phochhu basin has the maximum glacial lakes (0.05%) and of which 9 are PDGL. Hence a need for a time monitoring system is imminent. With the availability of free High-resolution Satellite Imageries and Advanced Remote Sensing tools, it has been a sine qua no for monitoring glacial lakes in high areas. Therefore, using the Google Earth Engine and Qgis Platform a semi-automated technique was used to generate glacial lake inventories for Phochhu Sub-basin for the year 2021.  We found out that there are 166 glacial lakes covering an area of 24.051 km2.

Considering the significance of thermodynamic cycles in the global energy system, it is necessary to develop new general classes of thermodynamic cycles to relieve current energy and environmental problems. Inspired by the relationship between power cycles and refrigeration cycles, we realize that general classes of thermodynamic cycles should occur in pairs with opposite functions. Here we reverse class 1 heating cycles to obtain another new general class of thermodynamic cycles named class 2 heating cycles (HC-2s). HC-2s have two basic forms, and each contains six thermodynamic processes. HC-2s present the simplest and most general approach to utilizing the temperature difference between a medium-temperature heat source and a low-temperature heat sink to achieve efficient high-temperature heating. HC-2s fill the gaps that have existed since the origin of thermal science, and they will play significant roles in the global sustainable energy system.

A transform approach based on a variable initial time (VIT) formulation is developed for discrete-time signals and linear time-varying discrete-time systems or digital filters. The VIT transform is a formal power series in z^(-1) which converts functions given by linear time-varying difference equations into left polynomial fractions with variable coefficients, and with initial conditions incorporated into the framework. It is shown that the transform satisfies a number of properties that are analogous to those of the ordinary z-transform, and that it is possible to do scaling of z^(- i) by time functions which results in left-fraction forms for the transform of a large class of functions including sinusoids with general time-varying amplitudes and frequencies. Using the extended right Euclidean algorithm in a skew polynomial ring with time-varying coefficients, it is shown that a sum of left polynomial fractions can be written as a single fraction, which results in linear time-varying recursions for the inverse transform of the combined fraction. The extraction of a first-order term from a given polynomial fraction is carried out in terms of the evaluation of z^(i) at time functions. In the application to linear time-varying systems, it is proved that the VIT transform of the system output is equal to the product of the VIT transform of the input and the VIT transform of the unit-pulse response function. For systems given by a time-varying moving average or an autoregressive model, the transform framework is used to determine the steady-state output response resulting from various signal inputs such as the step and cosine functions.

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.

In this paper we discuss a problem of generalization of binomial distributed triangle, that is sequence A287326 in OEIS. The main property of A287326 that it returns a perfect cube n as sum of n-th row terms over k; 0 ≤ k ≤ n−1 or 1 ≤ k ≤ n, by means of its symmetry. In this paper we have derived a similar triangles in order to receive powers m = 5; 7 as row items sum and generalized obtained results in order to receive every odd-powered monomial n^2m+1; m ≥ 0 as sum of row terms of corresponding triangle. In other words, in this manuscript are found and discussed the polynomials D_m(n,k) and U_m(n,k), such that, when being summed up over k in some range with respect to m and n returns the monomial n^2m+1.

Although Eulerian approaches are standard in computational acoustics, they are less effective for certain classes of problems like bubble acoustics and combustion noise. A different approach for solving acoustic problems is to compute with individual particles following particle motion. In this paper, a Lagrangian approach to model sound propagation in moving fluid is presented and implemented numerically, using three meshfree methods to solve the Lagrangian acoustic perturbation equations (LAPE) in the time domain. The LAPE split the fluid dynamic equations into a set of hydrodynamic equations for the motion of fluid particles and perturbation equations for the acoustic quantities corresponding to each fluid particle. Then, three meshfree methods, the smoothed particle hydrodynamics (SPH) method, the corrective smoothed particle (CSP) method, and the generalized finite difference (GFD) method, are introduced to solve the LAPE and the linearized LAPE (LLAPE). The SPH and CSP methods are widely used meshfree methods, while the GFD method based on the Taylor series expansion can be easily extended to higher orders. Applications to modeling sound propagation in steady or unsteady fluids in motion are outlined, treating a number of different cases in one and two space dimensions. A comparison of the LAPE and the LLAPE using the three meshfree methods is also presented. The Lagrangian approach shows good agreement with exact solutions. The comparison indicates that the CSP and GFD method exhibit convergence in cases with different background flow. The GFD method is more accurate, while the CSP method can handle higher Courant numbers.

We performed the accuracy assessment of three different Normalized Difference Water Indices (NDWIs) in water bodies during April 2019, a period in which floods occurred in a large proportion of the Southwest of the Buenos Aires Province (Argentina). The accuracy of the estimations using spaceborne medium-resolution multi-spectral imaging, and the reliability of three NDWIs to highlight shallow water features in satellite images, was evaluated using a high resolution airbone imagery as ground-truth. It is shown that these indices computed using Landsat 8 and Sentinel-2 imagery are only loosely correlated to the actual flooded area in shallow waters. Indeed, NDWI values vary significantly depending on the satellite mission used and the type of index computed.

This paper applied the space-time generalized finite difference scheme for solving the nonlinear dispersive shallow water waves as the modified Camassa−Holm equation, modified Degasperis-Procesi equation, Fornberg-Whitham equation, and its modified form. The proposed meshless numerical scheme was composed of the space-time generalized finite difference method, the two-step Newton-Raphson method, and the time-marching method. The space-time approach can treat the temporal derivative as one of the spatial derivatives. This numerical technique enables all the partial derivatives in the governing equation can be discretized by a spatial discretization method, and the mixed derivative can efficiently deal with using the proposed meshless numerical scheme. The space-time generalized finite difference method is advanced from the Taylor series expansion and the moving-least square method. The numerical discretization process is only related to functional data and weighting coefficients on the central node and its nearby nodes. Thus, the matrix system composed of nonlinear algebraic equations will be a sparse matrix and can be efficiently solved by the two-step Newton-Raphson method. Furthermore, the time-marching method was utilized to proceed with the space-time domain along the time axis. In this paper, several numerical examples were tested to verify the capability of the proposed space-time generalized finite difference scheme.

This paper presents an approach to the vibration analysis of axially functionally graded non-prismatic Timoshenko beams (AFGNPTB) using the finite difference method (FDM). The characteristics (cross-sectional area, moment of inertia, elastic moduli, shear moduli, and mass density) of axially functionally graded beams vary along the longitudinal axis. The Timoshenko beam theory covers cases associated with small deflections based on shear deformation and rotary inertia considerations. The FDM is an approximate method for solving problems described with differential equations. It does not involve solving differential equations; equations are formulated with values at selected points of the structure. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. In this paper, differential equations were formulated with finite differences, and additional points were introduced at the beam’s ends and at positions of discontinuity (supports, hinges, springs, concentrated mass, spring-mass system, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends and to satisfy the boundary and continuity conditions. Moreover, grid points with variable spacing were also considered, the grid being uniform within beam segments. Vibration analysis of AFGNPTB was conducted with this model, and natural frequencies were determined. Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of AFGNPTB, considering the damping. The results obtained in this study showed good agreement with those of other studies, and the accuracy was always increased through a grid refinement.

Green innovation has become the mainstream of the era, and New Urbanization is an inevitable choice in the process of urbanization in China. Focusing on the topics of green innovation and new urbanization, much work has been done to find their factors respectively while the relationship between the two remains to be explored. Hence, in this article, representative indicators of new urbanization and green innovation are selected to study the Yangtze River Delta City Group from the perspective of both the entire urban agglomeration and a single city, in terms of time and space, using the entropy method and the coupling model. The results show that (1). Green innovation promotes the new urbanization development and there is a synergistic effect between the two systems. (2). The level of economic development is an important factor that affects the degree of coupling degree and coordination degree between the two interactions, and its influence is better than the spatial effect. (3). Green innovation and new urbanization have positive spatial autocorrelation and regional agglomeration (there are High-High, Low-Low, and High-Low collections).

Solar energy is a renewable resource of energy that is broadly utilized and has the least emissions among the renewable energies. In this study, machine learning methods of artificial neural networks (ANNs), least squares support vector machines (LSSVM), and neuro-fuzzy are used for advancing prediction models for thermal performance of a photovoltaic-thermal solar collector (PV/T). In the proposed models, the inlet temperature, flow rate, heat, solar radiation, and the sun heat have been considered as the inputs variables. Data set has been extracted through experimental measurements from a novel solar collector system. Different analyses are performed to examine the credibility of the introduced approaches and evaluate their performance. The proposed LSSVM model outperformed ANFIS and ANNs models. LSSVM model is reported suitable when the laboratory measurements are costly and time-consuming, or achieving such values requires sophisticated interpretations.

This paper presents an approach to the vibration analysis of axially functionally graded (AFG) non-prismatic Euler-Bernoulli beams using the finite difference method (FDM). The characteristics (cross-sectional area, moment of inertia, elastic moduli, and mass density) of AFG beams vary along the longitudinal axis. The FDM is an approximate method for solving problems described with differential equations. It does not involve solving differential equations; equations are formulated with values at selected points of the structure. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. In this paper, differential equations were formulated with finite differences, and additional points were introduced at the beam’s ends and at positions of discontinuity (supports, hinges, springs, concentrated mass, spring-mass system, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends and to satisfy the boundary and continuity conditions. Moreover, grid points with variable spacing were also considered, the grid being uniform within beam segments. Vibration analysis of AFG non-prismatic Euler-Bernoulli beams was conducted with this model, and natural frequencies were determined. Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of AFG non-prismatic Euler-Bernoulli beams, considering the damping. The results obtained in this paper showed good agreement with those of other studies, and the accuracy was always increased through a grid refinement.

This article is a big picture of how electrical noise or conducted Electromagnetic Interference (EMI) is generated and mitigated in power converters. It gives an overview of what EMI in power converters is – from generation through to conduction and mitigation. It is meant to cover the complete subject as a summary so that the reader will have an outline of how to control conducted EMI by design (where possible) and how to mitigate by filtering. A clear distinction is made between Differential Mode (DM) and Common Mode (CM) EMI generation and mitigation. By using a boost converter as an example the trade-offs for DM noise control are discussed. It is shown how CM EMI is generated in a boost converter using the concept of the “Imbalance Difference Model” (IDM). Practical measurements for an in-line power filter is given showing the effect of the filter on the total EMI of a boost converter. Measurements for the CM current produced due to the imbalance difference for different values of the boost conductor are also shown. The main contribution of this study is linking CM noise generation to DM EMI. It is shown that CM noise is a direct consequence of DM noise (although circuit imbalance and coupling to a common ground also play a role). This paper will be useful to designers seeking the “bigger picture” of how EMI is generated in power converters and what can be done to mitigate the noise.

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.

How perceptual limits can be overcome has long been examined by psychologists. This study investigated whether visual cues, blindfolding, visual-auditory synesthetic experience and music training could facilitate a smaller frequency difference limen (FDL) in a gliding frequency discrimination test. It was hoped that the auditory limits could be overcome through visual facilitation, visual deprivation, involuntary cross-modal sensory experience or music practice. Ninety university students, with no visual or auditory impairment, were recruited for this one-between (blindfold/visual cue) and one-within (control/experimental session) designed study. A MATLAB program was prepared to test their FDL by an alternative forced-choice task (gliding upwards/gliding downwards/no change) and two questionnaires (Vividness of Mental Imagery Questionnaire & Projector-Associator Test) were used to assess their tendency to synesthesia. Participants with music training showed a significantly smaller FDL; on the other hand, being blindfolded, being provided with visual cues or having synesthetic experience before could not significantly reduce the FDL. However, the result showed a trend of reduced FDLs through blindfolding. This indicated that visual deprivation might slightly expand the limits in auditory perception. Overall, current study suggests that the inter-sensory perception can be enhanced through training but not though reallocating cognitive resources to certain modalities. Future studies are recommended to verify the effects of music practice on other perceptual limits.

Pavement friction plays a crucial role in ensuring the safety of road networks. Accurately assessing friction levels is vital for effective pavement maintenance and management strategies employed by state highway agencies. Traditionally, friction evaluations have been conducted on a case-by-case basis, focusing on specific road sections. However, this approach fails to provide a comprehensive assessment of friction conditions across the entire road network. This paper introduces a hybrid clustering algorithm, namely the combination of density-based spatial clustering of applications with noise (DBSCAN) and Gaussian mixture model (GMM), to perform pavement friction performance rating across a statewide road network. A large, safety-oriented dataset is first generated by integrating network friction and vehicle crash data based on the attributes contributing possibly to friction related crashes. One-, two-, and multi-dimensional clustering analyses, respectively, are then performed to rate pavement friction. The Chi-square test is further employed to validate and identify the practical ratings. It is shown that by effectively capturing the hidden, intricate patterns within the integrated, complex dataset and prioritizing friction-related safety attributes, the hybrid clustering algorithm can produce pavement friction ratings that align effectively with the current practices of the Indiana Department of Transportation (INDOT) in friction management.

This study examines the optimization of ground source heat pump (GSHP) systems and energy piles for sustainable urban development, with a focus on Southeast Asia. GSHPs, which utilize geothermal energy for indoor HVAC needs, offer a sustainable alternative to traditional systems by utilizing consistent subsurface temperatures for heating and cooling. The study highlights the importance of understanding thermal movement within the soil, especially in soft marine clays prevalent in Southeast Asia, to improve GSHP system efficiency. Using a one-dimensional finite-difference model, the study examines the effects of soil thermal conductivity and density on system performance. The results show that GSHP systems, especially when integrated with energy piles, significantly reduce electricity consumption and greenhouse gas emissions, underscoring their potential to mitigate the urban heat island effect in densely populated areas. Despite challenges posed by the region's hot and humid climate, which could affect long-term effectiveness, the study highlights the need for further study, including field experiments and advanced modeling techniques, to optimize GSHP configurations and fully exploit geothermal energy in urban environments. The study's insights into soil thermal dynamics and system design optimization contribute to the advancement of sustainable urban infrastructure development.

A recent study demonstrated that noninvasive measurements of cortical hemodynamics and metabolism in the resting human prefrontal cortex can facilitate quantitative metrics of unilateral mitochondrial-hemodynamic coupling and bilateral connectivity in infraslow oscillation (ISO) frequencies in young adults (YA). The ISO includes three distinct vasomotions with endogenic (E), neurogenic (N), and myogenic (M) frequency bands. The goal of this study was to prove the hypothesis that there are significant differences between young and older adults (OA) in the unilateral coupling (uCOP) and bilateral connectivity (bCON) in the prefrontal cortex. Accordingly, we performed measurements from 24 OA (67.2 ± 5.9 years of age) using the same 2-channel broadband near-infrared spectroscopy (bbNIRS) setup and resting-state experimental protocol as those in the recent study (Shahdadian et al, Cerebral Cortex Communications, 2022). After quantification of uCOP and bCON in three E/N/M frequencies and statistical analysis, we demonstrated that OA had significantly weaker bilateral hemodynamic connectivity but significantly stronger bilateral metabolic connectivity than YA in the M band. Furthermore, OA exhibited significantly stronger unilateral coupling on both prefrontal sides in all E/N/M bands, particularly with a very large effect size in the M band (> 1.9). These age-related results clearly support our hypothesis and were well interpreted following neurophysiological principles. The key finding of this paper is that the neurophysiological metrics of uCOP and bCON are highly associated with age and may have the potential to become meaningful features for human brain health and be translatable for future clinical applications, such as early detection of Alzheimer’s disease.

Taste signaling is a complex process that is linked to obesity and its associated metabolic syndromes. The sweet taste is mediated through a heterodimeric G protein coupled receptor (GPRC) in a species-specific manner and at multi-tissue specific levels. The sweet receptor recognizes a large number of ligands with structural and functional diversities to modulate different amplitudes of downstream signaling pathway(s). The human sweet-taste receptor has been extremely difficult to study by biophysical methods due to inadequate methods for producing large homogeneous quantities of the taste-receptor protein and a lack of reliable in vitro assays to precisely measure productive ligand binding modes leading to activity upon their interactions with the receptor protein. We report a multimodal high throughput assays to monitor ligand binding, receptor stability and conformational changes to model the molecular interactions between ligand-receptor. We applied saturation transfer difference nuclear magnetic resonance spectroscopy (STD-NMR) complemented by differential scanning calorimetry (DSC), circular dichroism (CD) spectroscopy, and intrinsic fluorescence spectroscopy (IF) to characterize binding interactions. Our method using complementary NMR and biophysical analysis is advantageous to study the mechanism of ligand binding and signaling processes in other GPCRs.

Carbon fiber reinforced polymer (CFRP) bending structure is widely used in the aviation field. The occurrence and propagation of delamination damage will reduce the safety of in-service bending structures. As a detection tool with high damage sensitivity, Lamb wave can effectively identify delamination damage. In this study, the signal difference coefficient (SDC) was introduced to quantify delamination damage and evaluate the sensitivity of A0 and S0 mode lamb waves to delamination damage. The simulation results show that compared with the S0 mode Lamb wave, the A0 mode Lamb wave exhibits higher delamination damage sensitivity. The delamination damage can be quantified based on the strong correlation between the SDC and the delamination damage size. The control effect of the linear array PZT phase time-delay method on Lamb wave mode was investigated by simulation. The phase time-delay method realizes the generation of a single-mode Lamb wave, which can separately excite the A0 mode and S0 mode Lamb wave to identify delamination damage of different sizes. The A0 mode Lamb wave was excited by the developed one-dimensional linear comb transducer (LCT), which was used to conduct the detection experiment on the CFRP bending plate with delamination damage sizes of Φ6.0 mm, Φ10.0 mm, and Φ15.0 mm. The experimental results verify the correctness of the simulation. According to the Hermite interpolation results of the finite-element simulation data, the relationship between the delamination damage size and the SDC was fitted by the Gaussian function, which can accurately quantify the delamination damage. The absolute error of the delamination damage quantification with fitting expression is not larger than 0.8 mm, and the percentage error is not more than 8%.

This paper presents an approach to the elastic analysis of beams subjected to Saint-Venant torsion using Green’s theorem and the finite difference method (FDM). The Saint-Venant torsion of beams, also called free torsion or unrestrained torsion, is characterized by the absence of axial stresses due to torsion; only shear stresses are developed. A solution to this torsion problem consists of finding a stress function that satisfies the governing equation and the boundary conditions. The FDM is an approximate method for solving problems described with differential equations; it does not involve solving differential equations, equations are formulated with values at selected nodes of the structure. In this paper, the beam’s cross-section was discretized using a two-dimensional grid and additional nodes were introduced on the boundaries. The introduction of additional nodes allowed us to apply the governing equations at boundary nodes and satisfy the boundary conditions. Beams with solid sections as well as multiply connected cross-sections were analyzed using this model; shear stresses and localized stresses at reentrant corners, torsion constant, and warping displacements were determined. Furthermore, beams with thin-walled closed sections, single-cell or multiple-cell, were analyzed using the Prandtl stress function whereby a linear distribution of the shear stresses over the thickness was considered; closed-form solutions for shear stresses and torsion constant were derived. The results obtained in this study showed good agreement with the exact results for rectangular cross-sections, and the accuracy was increased through a grid refinement. For thin-walled closed sections, the shear stresses obtained at the centerline using the closed-form solutions were in agreement with the values using Bredt’s analysis but the maximal values in the cross-section, which did not necessarily occur at the position with the smallest thickness, were higher; in addition, the results using the closed-form solutions were in good agreement with those using FDM.

This paper presents an approach to the Timoshenko beam theory (TBT) using the finite difference method (FDM). The Timoshenko beam theory covers cases associated with small deflections based on shear deformation and rotary inertia considerations. The FDM is an approximate method for solving problems described with differential equations. It does not involve solving differential equations; equations are formulated with values at selected points of the structure. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. The model developed in this paper consisted of formulating differential equations with finite differences and introducing additional points at the beam’s ends and at positions of discontinuity (concentrated loads or moments, supports, hinges, springs, brutal change of stiffness, spring-mass system, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends. Moreover, grid points with variable spacing were considered, the grid being uniform within beam segments. First-order, second-order, and vibration analyses of structures were conducted with this model. Furthermore, tapered beams were analyzed (element stiffness matrix, second-order analysis, vibration analysis). Finally, a direct time integration method (DTIM) was presented; the FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. The results obtained in this paper showed good agreement with those of other studies, and the accuracy was increased through a grid refinement. Especially in the first-order analysis of uniform beams, the results were exact for uniformly distributed and concentrated loads regardless of the grid.

This research delves into the intricate challenges confronting the agricultural sector, with a specialized focus on mitigating infections in tomato crops, particularly powdery mildew induced by the Leveillula Taurica pathogen. Tomatoes, renowned for their nutritional richness, are vital to global food security. However, conventional methodologies for disease detection exhibit both laborious processes and limited accuracy. In response to these challenges, this study advocates for an innovative fusion of hyperspectral imaging and deep learning methodologies to detect crop disease and its severity. The systematic workflow commenced with the curation of a dataset, involving the acquisition of live images through OpenCV, followed by conversion to RGB format and subsequent feature extraction utilizing a pre-trained visual geometry group (VGG-16) model for enhanced analysis. Sequentially, RGB images were transformed into simulated hyperspectral images (SHSI) leveraging a Neural Network generator model, offering a distinctive viewpoint on spectral information. This novel approach transcends conventional constraints by delivering a three-dimensional perspective, seamlessly integrating spatial and spectral dimensions for holistic data acquisition. The SHSI is further transmuted into a 3D visualization cube comprehensive grasp of spatial and spectral aspects encompassing spectral, spatial, and Haralick features. The research concludes with severity detection, categorized as low, moderate, or high, employing a Gaussian Mixture Model (GMM) and K-means for visualization.

Material identification based on R_value (transparency natural logarithm ratio of low-energy to high-energy) of line-scan dual-energy X-ray transmission (DE-XRT) has a good prospect for industrial application. Unfortunately, the DE-XRT signals before attenuation cannot be directly measured, whereas their precision is very important to R_value. Therefore, a context-based signal synthesis method which takes the filtered signals that remove high-frequency noises and retain low-frequency actual fluctuations as the reference value, and takes into account the vertical (forward/column) direction and horizontal (scanning/row) direction anisotropy of line-scan images was proposed. The vertical is a time series with continuity of signal trend; the horizontal is a spatial characteristic with the fluctuation synchronization within the same row signals. The special synthesis evaluations of curve synthesis difference and surface synthesis difference were also proposed. Experimental results show that the tow evaluations are both only about 0.0007, and it only takes 35 ms to complete the surface synthesis of 119×119 pixels on the CPU with 3.4 GHz main frequency. The presented method can achieve high signal synthesis precision and calculation real-time, so as to facilitate DE-XRT material identification. It can be extended to improve the precision of numerical synthesis in other line-scan uncomplete signals.

In the paper, by virtue of convolution theorem for the Laplace transforms, Bernstein's theorem for completely monotonic functions, and other techniques, the author finds necessary and sufficient conditions for a difference constituted by four derivatives of a function involving trigamma function to be completely monotonic.

Computation of Common Middle Point seismic sections and their subsequent time migration and diffraction imaging provides very important knowledge about the internal structure of 3D heterogeneous geological media and are key elements for successive geological interpretation. Full-scale numerical simulation, that computes all single shot seismograms, provides a full understanding of how the features of the image reflect the properties of the subsurface prototype. Unfortunately, this kind of simulations of 3D seismic surveys for realistic geological media needs huge computer resources, especially for simulation of seismic waves’ propagation through multiscale media like cavernous fractured reservoirs. In order to significantly reduce the query of computer resources we propose to model these 3D seismic cubes directly rather than the shot-by-shot simulation with subsequent CMP stacking. To do that we modify the well known "exploding reflectors principle" for 3D heterogeneous multiscale media by use of the finite-difference technique on the base of grids locally refined in time and space. To be able to simulate realistic models and acquisition we develop scalable parallel software, which needs reasonable computational costs. Numerical results for simulation of Common Middle Points sections and their time migration are presented and discussed.

This paper presents an approach to the Euler-Bernoulli beam theory (EBBT) using the finite difference method (FDM). The EBBT covers the case of small deflections, and shear deformations are not considered. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected points of the structure. Generally, the finite difference approximations are derived based on fourth-order polynomial hypothesis (FOPH) and second-order polynomial hypothesis (SOPH) for the deflection curve; the FOPH is made for the fourth and third derivative of the deflection curve while the SOPH is made for its second and first derivative. In addition, the boundary conditions and not the governing equations are applied at the beam’s ends. In this paper, the FOPH was made for all of the derivatives of the deflection curve, and additional points were introduced at the beam’s ends and positions of discontinuity (concentrated loads or moments, supports, hinges, springs, etc.). The introduction of additional points allowed us to apply the governing equations at the beam’s ends and to satisfy the boundary and continuity conditions. Moreover, grid points with variable spacing were also considered, the grid being uniform within beam segments. First-order analysis, second-order analysis, and vibration analysis of structures were conducted with this model. Furthermore, tapered beams were analyzed (element stiffness matrix, second-order analysis). Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. The results obtained in this paper showed good agreement with those of other studies, and the accuracy was increased through a grid refinement. Especially in the first-order analysis of uniform beams, the results were exact for uniformly distributed and concentrated loads regardless of the grid. Further research will be needed to investigate polynomial refinements (higher-order polynomials such as fifth-order, sixth-order…) of the deflection curve; the polynomial refinements aimed to increase the accuracy, whereby non-centered finite difference approximations at beam’s ends and positions of discontinuity would be used.

Throughout the world, salinity is a major environmental issue that limits agricultural productivity, particularly in arid and semi-arid regions. In addition, climate change is the most important reason for the salinization of agricultural soils in the world, so it is now essential to find solutions to increase salinity tolerance in plants. This study investigated the potential of arbuscular mycorrhizal fungi (AMF) inoculation to enhance the growth and yield performances of flax under different salinity levels by conducting a pot experiment. The experiment was laid out in a two-factor completely randomized design including AMF inoculation (AMF+: with and AMF-: without inoculation) and irrigation water salinity (0, 50, 100, and 150 mM NaCl). According to the results, it is evident that salt stress caused negative physiological effects, including limited growth, reduced photosynthesis, and decreased nitrogen (N) and phosphorus (P) content in shoots and roots of flax plants. Moreover, mycorrhizal association improved the salt tolerance of the plants by increasing chlorophyll content, enhancing N and P shoot and root contents, and consequently yield parameters, such as seed, and stem fiber yield, particularly at moderate salt concentrations (50 and 100 mM NaCl). As a result of using AMF, flax plants grown under salt stress exhibited tolerance, suggesting that AMF could be applied in saline environments to maintain ecological stability.

Alcohol use disorder (AUD) is a chronically relapsing disease characterized by loss of control in seeking and consuming alcohol (ethanol) driven by recruitment of brain stress systems. However, AUD differs among the sexes: men are more likely to develop AUD, but women progress from casual to binge drinking and heavy alcohol use more quickly. The central amygdala (CeA) is a hub of stress and anxiety, with corticotropin releasing factor (CRF)-CRF1 receptor and GABAergic signaling dysregulation occurring in alcohol dependent male rodents. However, we recently showed that GABAergic synapses in female rats are less sensitive to the acute effects of ethanol. Here, we used patch clamp electrophysiology to examine the effects of alcohol dependence on the CRF-modulation of rat CeA GABAergic transmission of both sexes. We found that GABAergic synapses of naïve female rats were unresponsive to CRF application compared males, although alcohol dependence induced a similar CRF responsivity in both sexes. In situ hybridization revealed that females had less CeA neurons containing mRNA for the CRF1 receptor (Crhr1) than males, but in dependence, the percentage of Crhr1-expressing neurons in females increased, unlike males. Overall, our data provide evidence for sexually dimorphic CeA CRF system effects on GABAergic synapses in dependence.

This paper presents an approach to the two-dimensional analysis of elastic isotropic deep beams using the finite difference method (FDM). Deep beams are subjected to in-plane loading and present a shear span to height ratio of less than 2.50; consequently, Euler-Bernoulli beam theory and Timoshenko beam theory do not apply. Deep beams analysis is generally conducted using numerical methods such as the finite element method and to a lesser extent the FDM; the strut-and-tie model and the stress field method are also widely utilized. Analytical approaches usually make use of the Airy stress function, where stresses are formulated in terms of the stress function; however, the exact solution of this function satisfying all of the boundary conditions can hardly be found, even for simple cases. In this paper, deep beams were analyzed using the FDM. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected nodes of the structure. Therefore, the deep beam was discretized with a two-dimensional grid, and additional nodes were introduced at the boundaries and at positions of discontinuity (openings, brutal change of material properties, non-uniform grid spacing), the number of additional nodes corresponding to the number of boundary conditions at the node of interest. The introduction of additional nodes allowed us to apply the governing equations at boundary nodes and satisfy the boundary and continuity conditions. An Airy stress function approach and a displacement potential function approach were considered in this study whereby strong formulations of equations (equilibrium, kinematic, and constitutive) were set. Stress and stability analyses were carried out with this model; furthermore, deep beams of varying stiffness, layered beams, and beams having openings were analyzed. For slender beams, the results obtained with the Airy stress function approach showed good agreement with those of the Euler-Bernoulli beam theory, and for deep beams, the shapes of stress distributions were in good agreement with a proper understanding of the behavior of structures. On the other hand, the displacement potential function approach delivered unsatisfactory results, probably due to the use of an inefficient equation solver; a more powerful tool will be needed in future research for this purpose.

The combination of interstitial hyperthermia treatment (IHT) with high dose rate brachytherapy (HDR-BT) for prostate cancer treatment and has the potential to improve clinical outcome, since it highly enhances the efficiency of cell kill, especially when applied simultaneously. Therefore, we have developed the ThermoBrachy applicators. To effectively apply optimal targeted IHT, treatment planning is considered essential. However, treatment planning in IHT is rarely applied since it is regarded difficult to accurately calculate the deposited energy in the tissue in a short enough time for clinical practice. In this study, we investigated various time-efficient methods for fast computation of the electromagnetic (EM) energy deposition resulting from the ThermoBrachy applicators. Initially, we investigated the use of an electro-quasistatic solver. Next, we extended our investigation to the application of geometric simplifications. Furthermore, we investigated the validity of the superpositioning principle, which can enable adaptive treatment plan optimization without the need for continuous recomputation of the EM field. Finally, we evaluated the accuracy of the methods by comparing them to the golden standard Finite-Difference Time-Domain calculation method using gamma-index analysis. The simplifications considerably reduced the computation time needed, improving from >12 h to a few seconds. All investigated methods showed excellent agreement with the golden standard by showing a >99% passing rate with 1%/0.5 mm Dose Difference and Distance-to-Agreement criteria. These results allow the proposed electromagnetic simulation method to be used for fast and accurate adaptive treatment planning.

This paper presents an approach to the Kirchhoff-Love plate theory (KLPT) using the finite difference method (FDM). The KLPT covers the case of small deflections, and shear deformations are not considered. The FDM is an approximate method for solving problems described with differential equations. The FDM does not involve solving differential equations; equations are formulated with values at selected points of the structure. Generally in the case of KLPT, the finite difference approximations are derived based on the fourth-order polynomial hypothesis (FOPH) and second-order polynomial hypothesis (SOPH) for the deflection surface. The FOPH is made for the fourth and third derivative of the deflection surface while the SOPH is made for its second and first derivative; this leads to a 13-point stencil for the governing equation. In addition, the boundary conditions and not the governing equations are applied at the plate edges. In this paper, the FOPH was made for all of the derivatives of the deflection surface; this led to a 25-point stencil for the governing equation. Furthermore, additional nodes were introduced at plate edges and at positions of discontinuity (continuous supports/hinges, incorporated beams, stiffeners, brutal change of stiffness, etc.), the number of additional nodes corresponding to the number of boundary conditions at the node of interest. The introduction of additional nodes allowed us to apply the governing equations at the plate edges and to satisfy the boundary and continuity conditions. First-order analysis, second-order analysis, buckling analysis, and vibration analysis of plates were conducted with this model. Moreover, plates of varying thickness and plates with stiffeners were analyzed. Finally, a direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of structures, with damping taken into account. In first-order, second-order, buckling, and vibration analyses of rectangular plates, the results obtained in this paper were in good agreement with those of well-established methods, and the accuracy was increased through a grid refinement.

Consciousness is a special type of interaction between subjects exchanging by lingua quanta (phonemes). A set of lingua quanta composes a thesaurus placed on the edge of chaos. Its library is a memory, modification of which is due to tuning of memristive neural elements scattered in the brain volume. The memristive neural model takes into consideration two types of neurons - excitatory and inhibitory and the current leakage, all at the body temperature, T = 309 K. At such temperatures only heavy ions like hydrogen ions - protons can pass robustly through the water medium of the brain. The robust ion transport is proton water wires supported by the Grotthuss mechanism. Final aims of the ions are the gap junctions (electric synapses) linking nearest neurons. Following these observations, a model of excitable nervous tissue is constructed. Namely, a difference mapping on the basis of sigmoid curves can reproduce chaotic modes of neural activity proved by positive values of the Lyapunov exponent. The edge of chaos is placed near the bifurcation boundary dividing the chaos and the periodic convulsive activity typical at the epileptic discharges. In this region the self-sustained spiral waves occur. An intermittent activity of the competing excitatory and inhibitory neurons is observed on the edge of chaos. The intermittent electrical activity of neural tissues is demonstrated by records both from different literature issues and records made by the author together with DR. A. Dudkin on slices of CA1 field of the hippocampus.

During middle age, women are less susceptible to nonalcoholic steatohepatitis (NASH) than men. Thus, we investigated the underlying molecular mechanisms behind these sexual differences using an established rat model of NASH. Mature female and male stroke-prone spontaneously hypertensive 5/Dmcr rats were fed control or high-fat-cholesterol (HFC) diets for 2, 8, and 14 weeks. Although HFC-induced hepatic fibrosis was markedly less severe in females than in males, only minor gender differences were observed in expression levels of cytochrome P450 enzymes (CYP)7A1, CYP8B1 CYP27A1, and CYP7B1, and multidrug resistance-associated protein 3, and bile salt export pump, which are involved in fibrosis-related bile acid (BA) kinetics. However, the BA detoxification-related enzymes UDP-glucuronosyltransferase (UGT) and sulfotransferase (SULT) 2A1, and the nuclear receptors constitutive androstene receptor (CAR) and pregnane X receptor (PXR), were strongly suppressed in HFC fed males, and were only slightly changed in HFC-diet fed females. Expression levels of the farnesoid X receptor and its small heterodimer partner were similarly regulated in a gender-dependent fashion following HFC feeding. Hence, the pronounced female resistance to HFC-induced liver damage likely reflects sustained expression of the nuclear receptors CAR and PXR and the BA detoxification enzymes UGT and SULT.

This paper first explains the theory of the origin of force from a new perspective, explaining that the origin of all forces is caused by the "missing of energy in the time dimension" (non-conservation in a certain period of time), and constructs a unified theory of force on this basis. Then, through the unified theory of force, the following phenomena are reinterpreted: (1) The common origin of the inverse square law of gravitation and electromagnetic force; (2) The origin of the repulsive force of electrons and protons at small distances; (3) The mystery of the rotation curve of galaxies and the new explanation of "dark matter"; (4) explain the origin of the magnetic field; (5) Predict that moving magnets have the same radiation effect as charged particles (electron magnetic moments also have radiation effects); (6) Reinterpret Maxwell's displacement current; (7) Reinterpret synchrotron radiation and bremsstrahlung; (8) Explain that the strong and weak forces are short-range forces; (9) Predict the true structure of black holes and the origin of gravitational waves;(10) Possible explanation for the origin of the speed of light.

We point out an inconsistency in Newton's equations of celestial mechanics. A set of differential equations implied by Newton's equations are shown to be free of this inconsistency. We then investigate the integrals of motion associated with this relative difference system.

With reference to electromagnetic interaction and Abdus Salam’s strong (nuclear) gravity, 1) Square root of ‘reciprocal’ of the strong coupling constant can be considered as the strength of nuclear elementary charge. 2) ‘Reciprocal’ of the strong coupling constant can be considered as the maximum strength of nuclear binding energy. 3) In deuteron, strength of nuclear binding energy is around unity and there exists no strong interaction in between neutron and proton. $G_s\cong 3.32688\times {10}^{28}{\mathrm{m}}^3{\mathrm{kg}}^{-1}{\sec }^{-2}$ being the nuclear gravitational constant, nuclear charge radius can be shown to be, $R_0\cong \frac{2G_s{m}_p}{c^2}\cong 1.24\mathrm{fm}.$ $e_s\cong \left(\frac{G_s{m}_p^2}{\hslash c}\right)e\cong 4.716785\times {10}^{-19}\mathrm{C}$ being the nuclear elementary charge, proton magnetic moment can be shown to be, ${\mu }_p\cong \frac{e_s\hslash }{2m_p}\cong \frac{e{G}_s{m}_p}{2c}\cong 1.48694\times {10}^{-26}\mathrm{J}{.\mathrm{T}}^{-1}.$ ${\alpha }_s\cong {\left(\frac{\hslash c}{G_s{m}_p^2}\right)}^2\cong 0.1153795$ being the strong coupling constant, strong interaction range can be shown to be proportional to $\exp \left(\frac{1}{{\alpha }_s^2}\right).$ Interesting points to be noted are: An increase in the value of ${\alpha }_s$ helps in decreasing the interaction range indicating a more strongly bound nuclear system. A decrease in the value of ${\alpha }_s$ helps in increasing the interaction range indicating a more weakly bound nuclear system. From $Z\cong 30$ onwards, close to stable mass numbers, nuclear binding energy can be addressed with, ${\left(B\right)}_{A_s}\cong Z\times \left\{\left(\frac{1}{{\alpha }_s}+1\right)+\sqrt{\sqrt{30\times 31}}\right\}\left(m_n-m_p\right)c^2\approx Z\times 19.66\mathrm{MeV}.$ With further study, magnitude of the Newtonian gravitational constant can be estimated with nuclear elementary physical constants. One sample relation is, $\left(\frac{G_N}{G_s}\right)\cong \frac{1}{2}{\left(\frac{m_e}{m_p}\right)}^{10}\left[\sqrt{\frac{G_F}{\hslash c}}/\left(\frac{\hslash }{m_{e}c}\right)\right]$ where $G_N$ represents the Newtonian gravitational constant and $G_F$ represents the Fermi’s weak coupling constant. Two interesting coincidences are, ${\left(m_p/m_e\right)}^{10}\cong \exp \left(1/{\alpha }_s^2\right)$ and $2G_s{m}_e/c^2\cong \sqrt{G_F/\hslash c}.$