THE NAVIER-STOKES EQUATIONS WITH THE KINEMATIC AND VORTICITY BOUNDARY CONDITIONS ON NON-FLAT BOUNDARIES

We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition （without the incompressibility condition）, which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n（n ≥ 3） with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.

How to Trace out Impact-Triggered Effects Globally Scattered around Formation Boundaries: Case Uhry, North Germany (Eocene/Oligocene Boundary)

By focusing on impact-triggered phenomena having occurred synchronously with or shortly prior to formation boundaries, two glass sand pits (Upper Maastrichtian) located near Uhry, North Germany have been studied in regard to the K/T boundary throughout the last 40 years during progressive exploitation of glass sand. However, a clastic sequence of sand, mass flow and pelite deposited in a deep channel of about 10 - 12 m in depth, eroded into the glass sand, surprisingly shows an Upper Eocene/Lower Oligocene age, well defined by a Dinocyst assemblage (Chiripteridium c. galea, Enneado cysta arcuata, Areoligera tauloma = D 12na - D 14na) from a 0.5 meter thick pelite that marks the Rupelian transgression within an estuarian system running northwest/southeastward. The section exposes a high energy mass flow and formerly solid frozen angular glass sand blocks of up to a meter-size embedded in fluvial sand of the channel base. Furthermore, erratic clastics of up to 0.4 meter in diameter appear at the pelite base. The “unusual” Dinocyst assemblage is of autochthonous origin and comprises the fresh water alga Pediastrum Kawraiskyias indicator for cold climate, hitherto only known from Quaternary. Missing pollen indicate a vegetation-less hinterland. Thus, there cannot be any doubt that around the E/O b. at least one “rare event” has happened as verified by short tremendous flooding and significant temperature fall (“cosmic winter”). According to the attitude of the global impact scientific community, these phenomena belong to the spectrum of “indirect effects” of major impacts. Radiometric ages of relevant major impact events underline that both impact craters of Popigai, Russia (100 Kilometer in diameter, 35.7 Ma) and Chesabreake, USA (85 Kilometer in diameter, 35.5 Ma) happened shortly before the E/O b.(33.75 Ma). In addition, a tektite strewn field along the eastern coast of the USA and micro-tektites (Gulf of Mexico, Caribbean Sea, Barbados) yield an age of ~34.4 Ma, close to the E/O b. Consequently, there does exist an extremely high probability that Uhry site hosts impact-triggered products at the E/O b. It should be stressed that the Upper Eocene Epoch comprises an amazingly high number of impact events during the time-span 34.2 - 37.0 Ma.

Theories and Applications of Mixed Quantum-Classical Non-adiabatic Dynamics

Electronically non-adiabatic processes are essential parts of photochemical process, collisions of excited species, electron transfer processes, and quantum information processing. Various non-adiabatic dynamics methods and their numerical implementation have been developed in the last decades. This review summarizes the most significant development of mixed quantum-classical methods and their applications which mainly include the Liouville equa- tion, Ehrenfest mean-field, trajectory surface hopping, and multiple spawning methods. The recently developed quantum trajectory mean-field method that accounts for the decoherence corrections in a parameter-free fashion is discussed in more detail.

The instability of dust acoustic waves in inhomogeneous dusty plasmas with non-adiabatic dust charge fluctuation

This paper investigates the propagation of linear dust acoustic waves in inhomogeneous dusty plasmas due to spatial gradients of dust charge, plasma densities. A linear dispersion relation is obtained with the non-adiabatic dust charge iguctuation and the non-thermally distributed ions. The numerical results show that the inhomogeneity, nonthermal ions and non-adlabatic dust charge iguctuatlon have strong iniguence on the frequency and the damping rate of waves.

New bifurcations of basin boundaries involving Wada and a smooth Wada basin boundary

This paper demonstrates and analyses double heteroclinic tangency in a three-well potential model, which can produce three new types of bifurcations of basin boundaries including from smooth to Wada basin boundaries, from fractal to Wada basin boundaries in which no changes of accessible periodic orbits happen, and from Wada to Wada basin boundaries. In a model of mechanical oscillator, it shows that a Wada basin boundary can be smooth.

Non-adiabatic quantum dynamical studies of Na(3p) + HD(v = 1, j5) = 0)→NaH/NaD + D/H reaction

Non-adiabatic dynamical calculations are carried out for the Na(3 p)+HD(ν = 1, j = 0)→NaH/NaD+D/H reaction on the diabatic potential energy surfaces of Wang et al.(Sci. Rep. 2018, 8, 17960) by using the time-dependent wave packet method. The state-to-state integral cross sections and differential cross sections of two reaction channels(NaH/NaD+D/H)are calculated for collision energy up to 0.4 eV. The cross section branching ratio indicates that the dominant reaction channel changes from NaD+H to NaH+D when the collision energy is larger than 0.227 eV. The products from two reaction channels both prefer to form in vibrationally cold but rotationally hot states, and they both tend to forward scattering.

Elastic one-return boundary element method and hybrid elastic thin-slab propagator for strong-contrast media with sharp boundaries

In this paper, we developed the theory and algorithm of an elastic one-way boundary element method（BEM） and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong contrast media. This approach can takes the advantage of accurate boundary condition of BEM and completely overcomes the weak contrast limitation of the perturbationtheory based one-way operator approach. The one-way BEM is a smooth boundary approximation, which avoids huge matrix operations in exact full BEM. In addition, the one-way BEM can model the primary-only transmitted and reflected waves and therefore is a valuable tool in elastic imaging and inversion. Through numerical tests for some simple models,we proved the validity and efficiency of the proposed method.

Solvent Effects on Excited-State Relaxation Dynamics of Paddle-Wheel BODIPY-Hexaoxatriphenylene Conjugates:Insights from Non-adiabatic Dynamics Simulations

Understanding the excited state dynamics of donor-acceptor(D-A)complexes is of fundamental importance both experimentally and theoretically.Herein,we have first explored the photoinduced dynamics of a recently synthesized paddle-wheel BODIPY-hexaoxatriphenylene(BODIPY is the abbreviation for BF_(2)-chelated dipyrromethenes)conjugates D-A complexes with the combination of both electronic structure calculations and nonadiabatic dynamics simulations.On the basis of computational results,we concluded that the BODIPY-hexaoxatriphenylene(BH)conjugates will be promoted to the local excited(LE)states of the BODIPY fragments upon excitation,which is followed by the ultrafast exciton transfer from LE state to charge transfer(CT).Instead of the photoinduced electron transfer process proposed in previous experimental work,such a exciton transfer process is accompanied with the photoinduced hole transfer from BODIPY to hexaoxatriphenylene.Additionally,solvent effects are found to play an important role in the photoinduced dynamics.Specifically,the hole transfer dynamics is accelerated by the acetonitrile solvent,which can be ascribed to significant influences of the solvents on the charge transfer states,i.e.the energy gaps between LE and CT excitons are reduced greatly and the non-adiabatic couplings are increased in the meantime.Our present work not only provides valuable insights into the underlying photoinduced mechanism of BH,but also can be helpful for the future design of novel donor-acceptor conjugates with better optoelectronic performance.

Non-Adiabatic Molecular Dynamics Simulations of Non-Charge-Transfer and Charge-Transfer Scattering in H^++CO2 at ELab=30 eV

The H^++CO2 reaction at high energies is relevant in atmospheric chemistry,astrophysics,and proton cancer therapy research.Therefore,we present herein a complete investigation of H^++CO2 at ELab=30 eV with the simplest-level electron nuclear dynamics(SLEND)method.SLEND describes nuclei via classical mechanics and electrons with a singledeterminantal Thouless wavefunction.The 3402 SLEND conducted simulations from 42 independent CO2 target orientations provide a full description of all the reactive processes and their mechanisms in this system:non-charge-transfer scattering(NCTS),charge-transfer scattering(CTS),and single C=O bond dissociation;all this valuable information about reactivity is not accessible experimentally.Numerous details of the projectile scattering patterns are provided,including the appearance and coalescence of primary and secondary rainbow angles as a function of the target orientation.SLEND NCTS and CTS differential cross sections(DCSs)are evaluated in conjunction with advanced semi-classical techniques.SLEND NCTS DCS agrees well with its experimental counterpart at all the measured scattering angles,whereas SLEND CTS DCS agrees well at high scattering angles but less satisfactorily at lower ones.Remarkably,both NCTS and CTS SLEND DCSs predict the primary rainbow angle signatures in agreement with the experiment.

A Lattice Boltzmann Boundary Condition at Impermeable Boundaries

A hydrodynamic boundary condition for the lattice Boltzmann model at impermeable boundaries is developed.This boundary condition satisfies both the no-slip condition and the fluid conservation at boundary nodes.Poiseuille flow and Couette flow are calculated with this technique to demonstrate the accuracy of the present boundary condition.

GLOBAL SOLUTIONS TO 1D COMPRESSIBLE NAVIER-STOKES/ALLEN-CAHN SYSTEM WITH DENSITY-DEPENDENT VISCOSITY AND FREE-BOUNDARY

This paper is concerned with the Navier-Stokes/Allen-Cahn system,which is used to model the dynamics of immiscible two-phase flows.We consider a 1D free boundary problem and assume that the viscosity coefficient depends on the density in the form ofη(ρ)=ρ^(α).The existence of unique global H^(2m)-solutions(m∈N)to the free boundary problem is proven for when 0<α<1/4.Furthermore,we obtain the global C^(∞)-solutions if the initial data is smooth.

Delineation of urban growth boundary based on FLUS model under the perspective of land use evaluation in hilly mountainous areas

With rapid economic development,the size of urban land in China is expanding dramatically.The Urban Growth Boundary(UGB)is an expandable spatial boundary for urban construction in a certain period in order to control the urban sprawl.Reasonable delineation of UGB can inhibit the disorderly spread of urban space and guide the normal development of the city.It is of practical significance for the construction of green urban space.The study utilizes GIS technology to establish a land construction suitability evaluation system for Nankang city,which is experiencing rapid urban expansion,and outlines the preliminary UGB under the future land use simulation(FLUS)model.At the same time,considering the coupled coordination of"Production-Living-Ecological Space",and based on the suitability evaluation,we revised the preliminary UGB by combining the advantages of the patch-generating land use simulation(PLUS)model and the convex hull model to delineate the final UGB.The results show that:1)the comprehensive score of the evaluation of the suitability of the construction of land from high to low shows the distribution of the center of the city to the surrounding circle type spread,the center of the city has the highest suitability score.The results of convex hull model show that the urban expansion type of Nankang is epitaxial.In the future,the urban expansion will mainly occur in the northern part of the city.The PLUS model predicts an increase of 3359.97 hm^(2)of construction land in Nankang by 2035,of which 2022.97 hm^(2)is urban construction land.2)The FLUS model has a prediction accuracy of 86.3%and delineates a preliminary UGB area of 9215.07 hm^(2).3)We used the results of the construction suitability evaluation,PLUS model simulation results,and convex hull model predictions to revise the originally delineated UGB.The final delineated UGB area is 8895.67 hm^(2)and it is capable of meeting the future development of the study area.The results of the delineation can promote sustainable urban development,and the delineation methodology can provide a reference basis for the preparation of territorial spatial planning.

Multi-scale physics-informed neural networks for solving high Reynolds number boundary layer flows based on matched asymptotic expansions

Multi-scale system remains a classical scientific problem in fluid dynamics,biology,etc.In the present study,a scheme of multi-scale Physics-informed neural networks is proposed to solve the boundary layer flow at high Reynolds numbers without any data.The flow is divided into several regions with different scales based on Prandtl's boundary theory.Different regions are solved with governing equations in different scales.The method of matched asymptotic expansions is used to make the flow field continuously.A flow on a semi infinite flat plate at a high Reynolds number is considered a multi-scale problem because the boundary layer scale is much smaller than the outer flow scale.The results are compared with the reference numerical solutions,which show that the msPINNs can solve the multi-scale problem of the boundary layer in high Reynolds number flows.This scheme can be developed for more multi-scale problems in the future.

Non-adiabatic State and Heat Energy Release in the Interior of the Earth

The latest results from the earth sciences indicated that the movement of the earth surface substances has a close relationship with the energy exchanges on large scale occurred in the margin between the core and mantle. There is a tendency towards the average temperature evolution of the earth dynamic system from higher to lower on the whole. And more and more evidence showed that the interior of the earth is in a non-adiabatic state, the cooling process occurred in the interior exerted a great influence on the surface environment. The development of the heat energy resources from the interior of the earth would result in the long term accelerating cooling though it may increase the surface temperature temporarily.

Dirac method for nonlinear and non-homogenous boundary value problems of plates

The boundary value problem plays a crucial role in the analytical investigation of continuum dynamics. In this paper, an analytical method based on the Dirac operator to solve the nonlinear and non-homogeneous boundary value problem of rectangular plates is proposed. The key concept behind this method is to transform the nonlinear or non-homogeneous part on the boundary into a lateral force within the governing function by the Dirac operator, which linearizes and homogenizes the original boundary, allowing one to employ the modal superposition method for obtaining solutions to reconstructive governing equations. Once projected into the modal space, the harmonic balance method(HBM) is utilized to solve coupled ordinary differential equations(ODEs)of truncated systems with nonlinearity. To validate the convergence and accuracy of the proposed Dirac method, the results of typical examples, involving nonlinearly restricted boundaries, moment excitation, and displacement excitation, are compared with those of the differential quadrature element method(DQEM). The results demonstrate that when dealing with nonlinear boundaries, the Dirac method exhibits more excellent accuracy and convergence compared with the DQEM. However, when facing displacement excitation, there exist some discrepancies between the proposed approach and simulations;nevertheless, the proposed method still accurately predicts resonant frequencies while being uniquely capable of handling nonuniform displacement excitations. Overall, this methodology offers a convenient way for addressing nonlinear and non-homogenous plate boundaries.

Grain boundary engineering for enhancing intergranular damage resistance of ferritic/martensitic steel P92

Ferritic/martensitic(F/M)steel is widely used as a structural material in thermal and nuclear power plants.However,it is susceptible to intergranular damage,which is a critical issue,under service conditions.In this study,to improve the resistance to intergranular damage of F/M steel,a thermomechanical process(TMP)was employed to achieve a grain boundary engineering(GBE)microstructure in F/M steel P92.The TMP,including cold-rolling thickness reduction of 6%,9%,and 12%,followed by austenitization at 1323 K for 40 min and tempering at 1053 K for 45 min,was applied to the as-received(AR)P92 steel.The prior austenite grain(PAG)size,prior austenite grain boundary character distribution(GBCD),and connectivity of prior austenite grain boundaries(PAGBs)were investigated.Compared to the AR specimen,the PAG size did not change significantly.The fraction of coincident site lattice boundaries(CSLBs,3≤Σ≤29)and Σ3^(n) boundaries along PAGBs decreased with increasing reduction ratio because the recrystallization fraction increased with increasing reduction ratio.The PAGB connectivity of the 6%deformed specimen slightly deteriorated compared with that of the AR specimen.Moreover,potentiodynamic polarization studies revealed that the intergranular damage resistance of the studied steel could be improved by increasing the fraction of CSLBs along the PAGBs,indicating that the TMP,which involves low deformation,could enhance the intergranular damage resistance.

A High-Accuracy Curve Boundary Recognition Method Based on the Lattice Boltzmann Method and Immersed Moving Boundary Method

Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.

Wavelet Multi-Resolution Interpolation Galerkin Method for Linear Singularly Perturbed Boundary Value Problems

In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.

Experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering boundary effect

The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted,where the hydraulic head difference(DH)varies,and the boundary includes non-loss and soil-loss conditions.Soil samples are filled into seven soil storerooms connected in turn.After evaluation,the variation in content of fine sand(ΔR_(f))and the hydraulic conductivity of soils in each storeroom(C_(i))are analyzed.In the non-loss test,the soil sample filling area is divided into runoff,transited,and accumulated areas according to the negative or positive ΔR_(f) values.ΔR_(f) increases from negative to positive along the seepage path,and Ci decreases from runoff area to transited area and then rebounds in accumulated area.In the soil-loss test,all soil sample filling areas belong to the runoff area,where the gentle-loss,strengthened-loss,and alleviated-loss parts are further divided.ΔR_(f) decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part,and C_(i) increases and then decreases along the seepage path.The relationship between ΔR_(f) and Ci is different with the boundary condition.Ci exponentially decreases with ΔR_(f) in the non-loss test and increases with ΔR_(f) generally in the soil-loss test.

Development of D_α band symmetrical visible optical diagnostic for boundary reconstruction on EAST tokamak

To investigate the potential of utilizing visible spectral imaging for controlling the plasma boundary shape during stable operation of plasma in future tokamak, a D_α band symmetric visible light diagnostic system was designed and implemented on the Experimental Advanced Superconducting Tokamak(EAST). This system leverages two symmetric optics for joint plasma imaging. The optical system exhibits a spatial resolution less than 2 mm at the poloidal cross-section, distortion within the field of view below 10%, and relative illumination of 91%.The high-quality images obtained enable clear observation of both the plasma boundary position and the characteristics of components within the vacuum vessel. Following system calibration and coordinate transformation, the image coordinate boundary features are mapped to the tokamak coordinate system. Utilizing this system, the plasma boundary was reconstructed, and the resulting representation showed alignment with the EFIT(Equilibrium Fitting) results. This underscores the system's superior performance in boundary reconstruction applications and provides a diagnostic foundation for boundary shape control based on visible spectral imaging.