Surface Boundary Research Articles

Nanometer-Resolved Mapping of Organic Cation Migration Behavior in Methylammonium Lead Halide Perovskites.

The performance and stability of organic metal halide perovskite (OMHP) optoelectronic devices have been associated with ion migration. Understanding of nanoscale resolved organic cation migration mechanism would facilitate structure engineering and commercialization of OMHP. Here, we report a three-dimensional approach for in situ nanoscale infrared imaging of organic ion migration behavior in OMHPs, enabling to distinguish migrations along grain boundary and in crystal lattice. Our results reveal that organic cation migration along OMHP film surface and grain boundaries (GBs) occurs at lower biases than in crystal lattice. We visualize the transition of organic cation migration channels from GBs to volume upon increasing electric field. The temporal resolved results demonstrate the fast MA+ migration kinetics at GBs, which is comparable to diffusivity of halide ions. Our findings help understand the role of organic cations in the performance of OMHP devices, and the proposed approach holds broad applicability for revealing migration mechanisms of organic ions in OMHPs based optoelectronic devices.

Improving Equatorial Upper Ocean Vertical Mixing in the NOAA/GFDL OM4 Model

AbstractDeficiencies in upper ocean vertical mixing parameterizations contribute to tropical upper ocean biases in global coupled general circulation models, affecting their simulated ocean heat uptake and ENSO variability. To better understand these deficiencies, we develop a suite of ocean model experiments including both idealized single column models and realistic global simulations. The vertical mixing parameterizations are first evaluated using large eddy simulations as a baseline to assess uncertainties and evaluate their implied turbulent mixing. Global models are then developed following NOAA/GFDL's 0.25° nominal horizontal grid spacing OM4 (uncoupled) configuration of the MOM6 ocean model, with various modifications that target biases in the original model. We test several enhancements to the existing mixing schemes and evaluate them against observational constraints from Tropical Atmosphere Ocean moorings and Argo floats. In particular, we find that we can improve the diurnal variability of mixing in OM4 via modifications to its surface boundary layer mixing scheme, and can improve the net mixing in the upper thermocline by reducing the background vertical viscosity, allowing for more realistic, less diffuse currents. The improved OM4 model better represents the mixing, leading to improved diurnal deep‐cycle variability, a more realistic time‐mean tropical thermocline structure, and a better Pacific Equatorial Undercurrent.

Impacts of building modifications on the turbulent flow and heat transfer in urban surface boundary layers

This study examines turbulent airflow and upward heat transport in real urban environments using a building-resolving large-eddy simulation model to understand the characteristics of turbulent airflow and upward heat transport when geometrical distributions of buildings are modified. The target areas were two real urban districts within Osaka City, Japan, having different morphological features. In the numerical experiments, the initial condition was set to a neutral condition in which temperature is uniformly distributed vertically, and buildings emitted heat at a constant rate. The results in the two districts indicated that the features of turbulence and heat transport distinctly differed with different building arrangement. Specifically, taller buildings significantly decelerated airflows and induced warming behind buildings. More high-rise buildings (which resulted in a larger building variability) in a district with a larger building density caused a large heat flux and warming at higher levels. The sensitivity experiments in which a density and height variability of buildings were modified showed that a building density at higher levels and a building height variability significantly influenced warming at upper levels. An increased building height variability weakened wind speed and disturbed horizontal heat advection, whereas a large building density caused numerous heat sources.

Skin-friction from temperature and velocity data around a wall-mounted cube

This paper reports an algorithm for measuring the time-averaged skin friction vector field τ¯(X)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\overline{\\pmb {\ au }}(\\pmb {X})$$\\end{document} starting from time-resolved temperature maps, acquired by a functional coating of temperature-sensitive paint. The algorithm is applied to a large area around a wall-mounted cube, immersed in the turbulent boundary layer over a flat plate. The method adopts a relaxed version of the Taylor Hypothesis operating on time-resolved maps of temperature fluctuations T′\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T'$$\\end{document} measured on the slightly warmer bounding surface. The procedure extracts U¯T(X)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\overline{U}}_T(\\pmb {X})$$\\end{document}, the celerity of displacement of T′\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T'$$\\end{document}, as the best approximation of the forecasting provided by the frozen turbulence assumption near the wall, where its rigorous application is inappropriate. The τ¯(X)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\overline{\\pmb {\ au }}(\\pmb {X})$$\\end{document} estimation is based on the hypothesis of a linear relationship between U¯T(X)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\overline{U}}_T(\\pmb {X})$$\\end{document} and U¯U(X)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\overline{U}}_U(\\pmb {X})$$\\end{document}, chained to the one between U¯U(X)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\overline{U}}_U(\\pmb {X})$$\\end{document} and U¯τ(X)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\overline{U}}_\ au (\\pmb {X})$$\\end{document}. We assess the outcomes of the proposed algorithm against those derived by the 2D and 3D Lagrangian particle tracking (LPT) methodology ’Shake-The-Box’, whose advent has made available high-quality near-wall flow field information. Furthermore, data from high-density 2D time-resolved LPT allows exploring the suitability of the linear relationships chain between U¯T(X)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\overline{U}}_T(\\pmb {X})$$\\end{document} and U¯τ(X)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\overline{U}}_\ au (\\pmb {X})$$\\end{document} in the proposed context.

Study of different theories of thermoelasticity under the propagation of Rayleigh waves in thermoelastic medium

The present research is on the propagation of Rayleigh waves in a homogenous thermoelastic solid half-space by considering the compact form of six different theories of thermoelasticity. The medium is subjected to an insulated boundary surface that is free from normal stress, tangential stress, and a temperature gradient normal to the surface. After developing a mathematical model, a dispersion equation is obtained with irrational terms. To apply the algebraic method, this equation must be converted into a rational polynomial equation. From this, only those roots are filtered out, which has satisfied both of the above equations for the propagation of waves decaying with depth. With the help of these roots, different characteristics are computed numerically, like phase velocity, attenuation coefficient, and path of particles. Various particular cases are compared graphically by using phase velocity and attenuation coefficient. The elliptic path of surface particles in Rayleigh wave propagation is also presented for the different theories using physical constants of copper material for different depths and thermal conductivity.

A Flexible Hierarchical Framework for Implicit 3D Characterization of Bionic Devices.

In practical applications, integrating three-dimensional models of bionic devices with simulation systems can predict their behavior and performance under various operating conditions, providing a basis for subsequent engineering optimization and improvements. This study proposes a framework for characterizing three-dimensional models of objects, focusing on extracting 3D structures and generating high-quality 3D models. The core concept involves obtaining the density output of the model from multiple images to enable adaptive boundary surface detection. The framework employs a hierarchical octree structure to partition the 3D space based on surface and geometric complexity. This approach includes recursive encoding and decoding of the octree structure and surface geometry, ultimately leading to the reconstruction of the 3D model. The framework has been validated through a series of experiments, yielding positive results.

New charged anisotropic solution in f(Q)-gravity and effect of non-metricity and electric charge parameters on constraining maximum mass of self-gravitating objects

In the present article, A new class of singularity-free charged anisotropic stars is derived in f(Q)-gravity regime. To solve the field equations, we assume a particular form of anisotropy along with an electric field and obtain a new exact solution in f(Q)-gravity. The explicit mathematical expression for the model parameters is derived by the smooth joining of the obtained solutions with the exterior Reissner–Nordstrom de-Sitter solution across the bounding surface of a compact star along with the requirement that the radial pressure vanishes at the boundary. We have modeled four self-gravitating pulsar objects such as LMC X-4, PSR J1903+327, PSR J1614-2230, and GW190814 in our current study and predict the radii of these objects that fall between 8 and 10 km. Furthermore, the physical validity of the solution is performed for self-gravitating object PSR J1614-2230 with mass 1.97±0.04M⊙\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1.97\\pm 0.04~M_{\\odot }$$\\end{document} with radius 10 km. The solution successfully fulfills all the physical requirements along with the stability and hydrostatic equilibrium conditions for a well-behaved model. The non-metricity f(Q)-parameter χ1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _{_1}$$\\end{document} and electric charge parameter η\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\eta $$\\end{document} play an important role in the maximum mass of the objects. The maximum mass increases when χ1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _{_1}$$\\end{document} and η\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\eta $$\\end{document} increase but a non-collapsing stable object can be obtained when χ1≤0.0205\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _{_1}\\le 0.0205$$\\end{document} and η≤0.0006\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\eta \\le 0.0006$$\\end{document}.

An efficient pseudoelastic pure P‐mode wave equation and the implementation of the free surface boundary condition

AbstractBased on the elastic wave equation, a pseudoelastic pure P‐mode wave equation has been recently derived by projecting the wavefield along the wavefront normal direction. This pseudoelastic pure P‐mode wave equation offers an accurate simulation of P‐wave fields with accurate elastic phase and amplitude characteristics. Moreover, considering no S‐waves are involved, it is computationally more efficient than the elastic wave equation, making it an excellent choice as a forward simulation engine for P‐wave exploration. Here, we propose a new pseudoelastic pure P‐mode wave equation and apply the stress image method to it to implement the free surface boundary condition. The new pseudoelastic wave equation offers significantly improved computational efficiency compared to the previous pseudoelastic wave equation. Additionally, the wavefields simulated by this new pseudoelastic wave equation exhibit clear physical interpretations. We evaluate the accuracy of the new wave equation in simulating elastic P‐waves by employing a model with high‐velocity contrasts. We find that this new equation, which purely admits P‐waves, though having exact amplitude and phase behaviour as the elastic waves for transmission components, the amplitudes slightly suffer in the scattering scenario. The difference in amplitude between the elastic and our pseudoelastic increases as the contrast in velocity at the interface (interlayer velocity ratio) increases, especially the S‐wave velocities. This has negative implications on scattering from the free surface boundary condition or the sea bottom interface, especially if the shear wave velocity below the surface or the sea bottom is high. However, in cases where, like for land data in the Middle East, the transition to a free surface is smoother, the accuracy of the pseudoelastic equation is high. In all cases, regardless of the interlayer velocity ratio, the accuracy of the pseudoelastic wave equation in simulating the elastic case, for scattered waves, exceeds that of the acoustic wave equation in phase and amplitude.

Investigating surface and internal waves and viscosity effects in two-layer fluids with a streamlined moving body using fully nonlinear simulation and experiment

Waves generated by a moving submerged moving body were analyzed in a two-layer density fluid. The experiments were performed in a two-dimensional towing tank and compared to the numerical results from a two-dimensional, fully nonlinear numerical towing tank (NTT) based on potential flow. The experimental conditions were determined by the relationship between the depth-dependent Froude number and the critical Froude number in baroclinic mode. The study was performed with various body velocities, changing the upper and lower fluid depth while maintaining the total water depth. An artificial damping coefficient was applied to the entire free surface and interface boundary condition to simulate the viscosity effect of the fluid interface. When the baroclinic mode was dominant, internal precursor soliton waves and depression waves were similar to the experimental results with strong nonlinearity. In contrast, the trailing waves were significantly different from the experimental results without applying an artificial damping scheme. A comparison with the experimental data showed that the effects of the fluid viscosity can be simulated with an appropriate damping coefficient. As the body velocity increases, the depressed or rising waves associated with the baroclinic mode are mixed with the trailing waves of the barotropic mode. As the depth of the upper fluid decreases, the characteristics of the baroclinic mode decrease, and the barotropic mode becomes dominant under the same velocity conditions.

The influence of the bounding surface on the structural ordering of short chains of oligoetherimides.

In this study, we have conducted a comparative analysis of the structural ordering of short oligoetherimide chains (dimers) near the bounding surface, depending on the structure of that surface. In order to clarify the possibility of oligoetherimide ordering along the symmetry axes of graphene, two types of bounding surfaces were considered: graphene, with a regular discrete position of interaction centers (carbon atoms), and a smooth, structureless impermeable wall. The chemical structures of the considered dimers consist of two repeating units of BPDA-P3, ODPA-P3, or aBPDA-P3 thermoplastic polyetherimides. Using all-atom molecular dynamics simulations, the process of structural ordering of the dimers near the surface of the graphene or wall was established. The ODPA-P3 and BPDA-P3 dimers form an ordered state near the graphene surface, while the aBPDA-P3 dimers do not demonstrate structural ordering. The simulation results confirmed that the ordering direction of the BPDA-P3 and ODPA-P3 dimers near the graphene surface is chosen randomly. Comparison of the oligoetherimide structure formed near the attracting wall without a symmetrical location of the interaction centers shows the similarity of the ordering of dimers near the graphene surface and the wall. As in the case of the graphene surface, the ordering of oligoetherimide molecules near the structureless wall demonstrates one direction of ordering. Therefore, we confirmed that the key factor for the onset of ordering is the presence of a confining surface, rather than the symmetrical arrangement of interaction centers in the substrate structure.

A Surface Determination Technique for Dimensional and Geometrical Analysis in Industrial X-ray Computed Tomography

Industrial X-ray computed tomography (XCT) is a nondestructive technique that can measure workpieces with non-accessible internal features or multimaterial components and assess the dimensional properties of assemblies in assembled states. Surface determination is one of its most crucial steps, which consists of determining boundary surfaces between a solid material and the surrounding air or between different solid materials. It allows for extracting surface points and assessing different features of the object from the data acquired through XCT scans. This task is particularly complex because of challenges associated with material properties, artefacts and noise. The main objective of this work is to assess not just the dimensional but also the geometric characteristics of industrial parts, which requires a more accurate definition of surface points. Therefore, we propose a new surface determination technique (SDT) in XCT to achieve subvoxel accuracy in determining surface points. We demonstrated the effectiveness and stability of our method by comparing it with other SDTs documented in the literature or with results from commercial software. The validation involved measuring various attributes, such as diameter, cylindricity and flatness, of a multi-stepped aluminium part calibrated by a coordinate measuring machine.

Anisotropic Durgapal-Fuloria compact stars in f(R) gravity

This study presented a new exact solution for anisotropic compact stellar objects within the framework of f(R)=R+αR2 gravity. In this context, the Durgapal-Fuloria metric potential has been employed to solve the field equation derived for f(R) theory. Furthermore, we have derived the generalized Darmois-Israel junction condition necessary for seamlessly connecting the interior region to the Schwarzschild exterior metric across the boundary hypersurface of the star in the context of f(R) gravity, and the interior solution is matched with the Schwarzschild exterior metric over the bounding surface of a compact star. These junction conditions stipulate that the pressure must not be zero at the boundary and should be proportional to the non-linear terms of f(R) gravity, a crucial aspect often overlooked by many researchers when investigating compact stellar models. Additionally, we derived the values of these parameters by using observational data of various compact stars (CSs), namely Her X-1, SAX J1808.4-3658, SMC X-1, LMC X-4, Cen X-3, 4U 1820-30, PSR J1903+327, 4U 1608-52, Vela X-1, and PSR J1416-2230. This approach enables us to investigate the comprehensive analysis of solutions numerically and graphically. We conducted various physical tests, including gradient of energy density and pressures, anisotropy, stability, equilibrium conditions, energy-density constraints, mass function, compactness, redshift, and adiabatic index, to assess the feasibility of our models. Our findings demonstrate the consistent behavior of our models provides a satisfactory physical situation as far as the observational results are confirmed.

Drift Control of High-Dimensional Reflected Brownian Motion: A Computational Method Based on Neural Networks

Motivated by applications in queueing theory, we consider a stochastic control problem whose state space is the d-dimensional positive orthant. The controlled process Z evolves as a reflected Brownian motion whose covariance matrix is exogenously specified, as are its directions of reflection from the orthant’s boundary surfaces. A system manager chooses a drift vector [Formula: see text] at each time t based on the history of Z, and the cost rate at time t depends on both [Formula: see text] and [Formula: see text]. In our initial problem formulation, the objective is to minimize expected discounted cost over an infinite planning horizon, after which we treat the corresponding ergodic control problem. Extending the earlier work by Han et al. [Han J, Jentzen A, Weinan E (2018) Solving high-dimensional partial differential equations using deep learning. Proc. Natl. Acad. Sci. USA 115(34):8505–8510], we develop and illustrate a simulation-based computational method that relies heavily on deep neural network technology. For the test problems studied thus far, our method is accurate to within a fraction of 1% and is computationally feasible in dimensions up to at least [Formula: see text].

A modeling method for thermal steady-state simulation of the four-layer printed circuit board.

Thermal simulation of a Printed Circuit Board (PCB) can help identify potential overheating risks in the circuit. The proposed modeling method combines analytical temperature solutions and numerical approximations. Only Fourier-series analytical solutions related to the prepreg-layer surfaces need to be calculated, rather than the entire structure. Heat transfer through the lateral sides of a PCB is approximately considered as part of the compensated heat flux of the insulating-layer surface boundaries. Heat diffusion within or between metal layers is numerically approximated using the finite volume method. The core layer is treated as "thermally-thick". Temperature-dependent boundary conditions are considered through iterations. A test solver was developed based on the method. The modeling accuracy was validated by comparison with COMSOL Multiphysics for a four-layer structure with a moderate degree of discretization. Additionally, a PCB for generating DC 3.3V was designed, tested, and modeled, with the modeling results confirmed by the thermal images. The electro-thermal analysis of the distribution of electric potential and Joule heating in traces and vias was integrated into the PCB model. The layout maps of the PCB were further adjusted to reduce Joule heating in the output circuit, and the improvement on reducing the IR drop and hotspot temperature was examined.

Transformation of surface [formula omitted]-forces in high grade elasticity

The expression of the elastic energy, postulated for a given material, specifies through the principle of virtual work which external loadings turn out to be admissible. For the elastic bodies, the energy of which depends on the gradients of the placement higher than (H+1) (or, equivalently, on the derivatives of the deformation gradient higher than H), external force densities can be prescribed over the boundary surface, which do work on the H-th derivative of the virtual displacements along the normal direction. Such generalized forces with H=1,2 have been dealt with in the literature, referred to as double and triple forces: when H is increased, however, the analytical evaluation of these work terms may become prohibitive. Each Eulerian H-force density of this kind, acting on the current boundary surface, turns out to be equivalent to a set of diverse actions defined in the Lagrangian configuration. These actions include: (i) surface terms doing work on the virtual displacement vector; (ii) surface terms doing work on the higher order derivatives of the virtual displacement along the normal direction; (iii) edge terms, doing work on the virtual displacement or on its higher order gradients, which may be further reduced. In this study, a recursive algorithm is utilized to generate all these Lagrangian contributions for any H. Such an approach rests on the binomial expansion of the tensor products of complementary surface projectors, symmetrized in the work duality, and on the integration by parts combined with the generalized divergence theorem.

Intelligent computing applications to study the tri-hybrid nanofluid past over the stretched surface

In this study, a recurrent neural network with the Levenberg-Marquardt backpropagation algorithm is used to examine the effects of various parameters on cooling mechanisms involving tri-hybrid nanofluid consisting of Al2O3, Cu, TiO2 and magnetohydrodynamic stagnation point flow. This research focus on nanofluid in the existence of an asymmetrically stretched disc in water, which plays an imprortant role for understanding heat transfer in addition to the dynamics of momentum and thermal boundary surface. Transform the complex partial differential system through suitable similarity transformations into the ordinary differential equations for numerical computing. These equations have been solved numerically by using the Adams method in Mathematica to generate datasets for various parameters.The current study demonstrates that the velocity profile rises with higher magnetic, velocity slip, and mass suction parameters. The fluid thermal level tends to decrease with an increase in mass suction parameters, whereas an increase in the magnetic parameter and Eckert number raises the fluid’s thermal level. Furthermore, it is noticed that the flow resistance and internal friction of the fluid are proportional to the viscosity parameter. The proposed scheme has shown its effectiveness in a wide range of scenarios and cases. This was achieved by comparing the mean square error metric for the squeezing flow tri-hybrid nanofluid problem with the obtained results. In addition, various statistical measures like state transitions, fitness assessments, error histograms, and regression analyses, provide further evidence of the solver’s robustness and reliability. The excellent agreement between the reference datasets and refined numerical solutions through the scheme showcases the method’s effectiveness, achieving remarkable accuracy levels ranging from 10-6 to 10-12.

Determination of five-parameter grain boundary characteristics in nanocrystalline Ni-W by scanning precession electron diffraction tomography

Determining the full five-parameter grain boundary characteristics from experiments is essential for understanding grain boundaries impact on material properties, improving related models, and designing advanced alloys. However, achieving this is generally challenging, in particular at nanoscale, due to their 3D nature. In our study, we successfully determined the grain boundary characteristics of an annealed nickel-tungsten alloy (NiW) nanocrystalline needle-shaped specimen (tip) containing twins using Scanning Precession Electron Diffraction (SPED) Tomography. The presence of annealing twins in this face-centered cubic (fcc) material gives rise to common reflections in the SPED diffraction patterns, which challenges the reconstruction of orientation-specific virtual dark field (VDF) images required for tomographic reconstruction of the 3D grain shapes. To address this, an automated post-processing step identifies and deselects these shared reflections prior to the reconstruction of the VDF images. Combined with appropriate intensity normalization and projection alignment procedures, this approach enables high-fidelity 3D reconstruction of the individual grains contained in the needle-shaped sample volume. To probe the accuracy of the resulting boundary characteristics, the twin boundary surface normal directions were extracted from the 3D voxelated grain boundary map using a 3D Hough transform. For the sub-set of coherent Σ3 boundaries, the expected {111} grain boundary plane normals were obtained with an angular error of <3° for boundary sizes down to 400 nm². This work advances our ability to precisely characterize and understand the complex grain boundaries that govern material properties.

Exploring the effects of the alkaline earth metal cations on the electronic structure, photostability and radiation hardness of lead halide perovskites

The growing interest to the application of perovskite solar cells (PSC) in aerospace requires great attention to the design of new absorber materials with enhanced photostability and radiation hardness. We addressed this challenge through the B-site engineering by partial replacing of Pb2+ in the complex halides with Ca2+, Sr2+, and Ba2+. Change in the material electronic properties suggests the incorporation of the substituent cations in the perovskite lattice. However, when the loading of alkaline earth halides is high, they tend to segregate to the film surface and grain boundaries. The modification of MAPbI3 and Cs0.12FA0.88PbI3 by replacing a 1% of Pb2+ with alkaline earth metal cations resulted in improvements in the photostability and radiation hardness under exposure to gamma rays and electrons beam. Depending on the material formulation, it was possible to mitigate the photochemical and radiation-induced Pb0 formation and the univalent cation phase segregation. The best of the designed perovskite absorbers could successfully tolerate high doses of gamma rays (5.5 MGy) and electron fluences (3 × 1016 e/cm2), thus outperforming significantly the non-modified material. These findings feature the incorporation of alkaline earth metal cations as a promising approach for the development of PSC with enhanced radiation hardness for aerospace applications.

Effect of Rarefactions and Convective Heat Change on Free Convective Unsteady MHD Flow in a Slip-Flow Regime Past a Vertical Wall with Convective Surface Boundary Condition

The study investigates the unsteady free convective two-dimensional MHD flow past a vertical porous plate with convective surface boundary condition in porous medium in slip flow regime under the action of variable suction velocity. Analytical solutions are obtained for the system using perturbation technique that converts non-linear coupled governing partial differential equations into non-dimensional form of ordinary differential equations. Effects of variable suction velocity, rarefactions parameter and heat change parameter are analysed and discussed graphically for various values of effective physical parameter such as Grasshof number, Magnetic field parameter, Prandtl number, Permeability parameter on fluid velocity and temperature, skin friction, and heat transfer.

SnSe2 Quantum Dots and Chlorhexidine Acetate Suppress Synergistically Non-radiative Recombination Loss for High Efficiency and Stability Perovskite Solar Cells.

Non-radiative recombination losses limit the property of perovskite solar cells (PSCs). Here, a synergistic strategy of SnSe2QDs doping into SnO2 and chlorhexidine acetate (CA) coating on the surface of perovskite is proposed. The introduction of 2D SnSe2QDs reduces the oxygen vacancy defects and increases the carrier mobility of SnO2. The optimized SnO2 as a buried interface obviously improves the crystallization quality of perovskite. The CA containing abundant active sites of ─NH2/─NH─, ─C═N, CO, ─Cl groups passivate the defects on the surface and grain boundary of perovskite. The alkyl chain of CA also improves the hydrophobicity of perovskite. Moreover, the synergism of SnSe2QDs and CA releases the residual stress and regulates the energy level arrangement at the top and bottom interface of perovskite. Benefiting from these advantages, the bulk and interface non-radiative recombination loss is greatly suppressed and thereby increases the carrier transport and extraction in devices. As a result, the best power conversion efficiency (PCE) of 23.41% for rigid PSCs and the best PCE of 21.84% for flexible PSCs are reached. The rigid PSC maintains 89% of initial efficiency after storing nitrogen for 3100 h. The flexible PSCs retain 87% of the initial PCE after 5000 bending cycles at a bending radius of 5 mm.