Improved Diurnal Cycle of Precipitation on Land in a Global Non-Hydrostatic Model Using a Revised NSAS Deep Convective Scheme

In relatively coarse-resolution atmospheric models,cumulus parameterization helps account for the effect of subgridscale convection,which produces supplemental rainfall to the grid-scale precipitation and impacts the diurnal cycle of precipitation.In this study,the diurnal cycle of precipitation was studied using the new simplified Arakawa-Schubert scheme in a global non-hydrostatic atmospheric model,i.e.,the Yin-Yang-grid Unified Model for the Atmosphere.Two new diagnostic closures and a convective trigger function were suggested to emphasize the job of the cloud work function corresponding to the free tropospheric large-scale forcing.Numerical results of the 0.25-degree model in 3-month batched real-case simulations revealed an improvement in the diurnal precipitation variation by using a revised trigger function with an enhanced dynamical constraint on the convective initiation and a suitable threshold of the trigger.By reducing the occurrence of convection during peak solar radiation hours,the revised scheme was shown to be effective in delaying the appearance of early-afternoon rainfall peaks over most land areas and accentuating the nocturnal peaks that were wrongly concealed by the more substantial afternoon peak.In addition,the revised scheme enhanced the simulation capability of the precipitation probability density function,such as increasing the extremely low-and high-intensity precipitation events and decreasing small and moderate rainfall events,which contributed to the reduction of precipitation bias over mid-latitude and tropical land areas.

Two-Layer Non-Hydrostatic Model for Generation and Propagation of Interfacial Waves

When pycnocline thickness of ocean density is relatively small, density stratification can be well represented as a two-layer system. In this article, a depth integrated model of the two-layer fluid with constant density is considered,and a variant of the edge-based non-hydrostatic numerical scheme is formulated. The resulting scheme is very efficient since it resolves the vertical fluid depth only in two layers. Despite using just two layers, the numerical dispersion is shown to agree with the analytical dispersion curves over a wide range of kd, where k is the wave number and d the water depth. The scheme was tested by simulating an interfacial solitary wave propagating over a flat bottom, as well as over a bottom step. On a laboratory scale, the formation of an interfacial wave is simulated,which also shows the interaction of wave with a triangular bathymetry. Then, a case study using the Lombok Strait topography is discussed, and the results show the development of an interfacial wave due to a strong current passing through a sill.

Non-hydrostatic modelling of regular wave transformation and current circulation in an idealized reef-lagoon-channel system

The wave-induced setup and circulation in a two dimensional horizontal(2DH)reef-lagoon-channel system is investigated by a non-hydrostatic model.The simulated results agree well with observations from the laboratory experiments,revealing that the model is valid in simulating wave transformation and currents over reefs.The effects of incident wave height,period,and reef flat water depth on the mean sea level and wave-driven currents are examined.Results show that the distributions of mean sea level and current velocities on the reef flat adjacent to the channel vary significantly from those in the area close to the side walls.From the wave averaged current field,an obvious alongshore flux flowing from the reef flat to the channel is captured.The flux from the reef flat composes the second source of the offshore rip current,while the first source is from the lagoon.A detailed momentum balance analysis shows that the alongshore current is mainly induced by the pressure gradient between the reef flat and the channel.In the lagoon,the momentum balances are between the pressure and radiation stress gradient,which drives flow towards the channel.Along the channel,the offshore current is mainly driven by the pressure gradient.

Numerical Analysis of the Nonlinear Parameterization of Waves in Currents over a Submerged Sill with a Non-Hydrostatic Model

An investigation of the effects of a uniform current strength direction(following or opposing wave propagation) on the nonlinear transformation of irregular waves over a submerged trapezoidal sill is carried out using SWASH,a non-hydrostatic numerical wave model.The nonlinear parameters(i.e.,asymmetry,skewness,and kurtosis) are calculated,and the empirical formulas for these parameters are presented as a function of the local Ursell number based on the present numerical data measured.In the shoaling area of the submerged sill,the nonlinear characteristics of waves are more obvious when waves propagate in the same direction as the currents than when waves propagate in the opposite direction.Whereas nonlinear parameters grow with the strengthening of the following currents over the crest,they tend to decrease as the adverse current velocity increases over the crest area of the submerged sill.

Simulation of wave scattering over a floating platform in the ocean with a coupled CFD-IBM model

A numerical study of linear wave scattering over a floating platform has been simulated by an efficient numericalmodel in this letter.The non-hydrostatic model is used to simulate the free surface and the uneven bottom.For thesolid body modelling,the immersed boundary method(IBM)is implemented by introducing a virtual boundaryforce into the momentum equations to emulate the boundary conditions.This implementation enhances theability of the model to simulate interactions between waves and floating structures.A numerical case involvingwave interactions with a floating platform is studied to validate the numerical model.By simulating the wavepropagation,the numerical model captures the variation of the wave scattering very well,which verifies theperformance of the numerical model and the robust strategy of the IBM.

A Hybrid Finite-Volume and Finite Difference Scheme for Depth- Integrated Non-Hydrostatic Model

A depth-integrated, non-hydrostatic model with hybrid finite difference and finite volume numerical algorithm is proposed in this paper. By utilizing a fraction step method, the governing equations are decomposed into hydrostatic and non-hydrostatic parts. The first part is solved by using the finite volume conservative discretization method, whilst the latter is considered by solving discretized Poisson-type equations with the finite difference method. The second-order accuracy, both in time and space, of the finite volume scheme is achieved by using an explicit predictor-correction step and linear construction of variable state in cells. The fluxes across the cell faces are computed in a Godunov-based manner by using MUSTA scheme. Slope and flux limiting technique is used to equip the algorithm with total variation dimensioning property for shock capturing purpose. Wave breaking is treated as a shock by switching off the non-hydrostatic pressure in the steep wave front locally. The model deals with moving wet/dry front in a simple way. Numerical experiments are conducted to verify the proposed model.

Maximum initial primary wave model for low-Froude-number reservoir landslides based on wave theory

The impulse waves induced by large-reservoir landslides can be characterized by a low Froude number.However,systematic research on predictive models specifically targeting the initial primary wave is lacking.Taking the Shuipingzi 1#landslide that occurred in the Baihetan Reservoir area of the Jinsha River in China as an engineering example,this study established a large-scale physical model(with dimensions of 30 m×29 m×3.5 m at a scale of 1:150)and conducted scaled experiments on 3D landslide-induced impulse waves.During the process in which a sliding mass displaced and compressed a body of water to generate waves,the maximum initial wave amplitude was found to be positively correlated with the sliding velocity and the volume of the landslide.With the increase in the water depth,the wave amplitude initially increased and then decreased.The duration of pressure exertion by the sliding mass at its maximum velocity directly correlated with an elevated wave amplitude.Based on the theories of low-amplitude waves and energy conservation,while considering the energy conversion efficiency,a predictive model for the initial wave amplitude was derived.This model could fit and validate the functions of wavelength and wave velocity.The accuracy of the initial wave amplitude was verified using physical experiment data,with a prediction accuracy for the maximum initial wave amplitude reaching 90%.The conversion efficiency(η)directly determined the accuracy of the estimation formula.Under clear conditions for landslide-induced impulse wave generation,estimating the value ofηthrough analogy cases was feasible.This study has derived the landslide-induced impulse waves amplitude prediction formula from the standpoints of wave theory and energy conservation,with greater consideration given to the intrinsic characteristics in the formation process of landslide-induced impulse waves,thereby enhancing the applicability and extensibility of the formula.This can facilitate the development of empirical estimation methods for landslide-induced impulse waves toward universality.

Contrasts of bimodal tropical instability waves(TIWs)-induced wind stress perturbations in the Pacific Ocean among observations,ocean models,and coupled climate models

The coupling between wind stress perturbations and sea surface temperature(SST)perturbations induced by tropical instability waves(TIWs)in the Pacific Ocean has been revealed previously and proven crucial to both the atmosphere and ocean.However,an overlooked fact by previous studies is that the loosely defined“TIWs”actually consist of two modes,including the Yanai wave-based TIW on the equator(hereafter eTIW)and the Rossby wave-based TIW off the equator(hereafter vTIW).Hence,the individual feedbacks of the wind stress to the bimodal TIWs remain unexplored.In this study,individual coupling relationships are established for both eTIW and v TIW,including the relationship between the TIW-induced SST perturbations and two components of wind stress perturbations,and the relationship between the TIW-induced wind stress perturbation divergence(curl)and the downwind(crosswind)TIW-induced SST gradients.Results show that,due to different distributions of eTIW and vTIW,the coupling strength induced by the eTIW is stronger on the equator,and that by the vTIW is stronger off the equator.The results of any of eTIW and vTIW are higher than those of the loosely defined TIWs.We further investigated how well the coupling relationships remained in several widely recognized oceanic general circulation models and fully coupled climate models.However,the coupling relationships cannot be well represented in most numerical models.Finally,we confirmed that higher resolution usually corresponds to more accurate simulation.Therefore,the coupling models established in this study are complementary to previous research and can be used to refine the oceanic and coupled climate models.

Analysis and Evaluation of the Global Aerosol Optical Properties Simulated by an Online Aerosol-coupled Non-hydrostatic Icosahedral Atmospheric Model

Aerosol optical properties are simulated using the Spectral Radiation Transport Model I~）r Aerosol Species （SPRINTARS） coupled with the Non-hydrostatic ICosahedral Atmospheric Model （NICAM）. The 3-year global mean all-sky aerosol optical thickness （AOT） at 550 nm, theAngstr/Sm Exponent （AE） based on AOTs at 440 and 870 nm, and the single scattering albedo （SSA） at 550 nm are estimated at 0.123, 0.657 and 0.944, respectively. For each aerosol species, the mean AOT is within the range of the AeroCom models. Both the modeled all-sky and clear-sky results are compared with observations from the Moderate Resolution Imaging Spectroradiometer （MODIS） and the Aerosol Robotic Network （AERONET）. The simulated spatiotemporal distributions of all-sky AOTs can generally reproduce the MODIS retrievals, and the correlation and model skill can be slightly improved using the clear-sky results over most land regions. The differences between clear-sky and all-sky AOTs are larger over polluted regions. Compared with observations from AERONET, the modeled and observed all-sky AOTs and AEs are generally in reasonable agreement, whereas the SSA variation is not well captured. Although the spatiotemporal distributions of all-sky and clear-sky results are similar, the clear-sky results are generally better correlated with the observations. The clear-sky AOT and SSA are generally lower than the all-sky results, especially in those regions where the aerosol chemical composition is contributed to mostly by sulfate aerosol. The modeled clear-sky AE is larger than the all-sky AE over those regions dominated by hydrophilic aerosol, while the＇opposite is found over regions dominated by hydrophobic aerosol.

Comparisons of Wave Force Model Effects on the Structural Responses and Fatigue Loads of a Semi-Submersible Floating Wind Turbine

The selection of wave force models will significantly impact the structural responses of floating wind turbines.In this study,comparisons of wave force model effects on the structural responses and fatigue loads of a semi-submersible floating wind turbine(SFWT)were conducted.Simulations were performed by employing the Morison equation(ME)with linear or second-order wave kinematics and potential flow theory(PFT)with first-or second-order wave forces.A comparison of regular waves,irregular waves,and coupled wind/waves analyses with the experimental data showed that many of the simulation results and experimental data are relatively consistent.However,notable discrepancies are found in the response amplitude operators for platform heave,tower base bending moment,and tension in mooring lines.PFT models give more satisfactory results of heave but more significant discrepan-cies in tower base bending moment than the ME models.In irregular wave analyses,low-frequency resonances were captured by PFT models with second-order difference-frequency terms,and high-frequency resonances were captured by the ME models or PFT models with second-order sum-frequency terms.These force models capture the response frequencies but do not reasonably predict the response amplitudes.The coupled wind/waves analyses showed more satisfactory results than the wave-only analyses.However,an important detail to note is that this satisfactory result is based on the overprediction of wind-induced responses.

Simulating the evolution of focused waves by a two-layer Boussinesq-type model

Accurate simulation of the evolution of freak waves by the wave phase focusing method requires accurate linear and nonlinear properties,especially in deep-water conditions.In this paper,we analyze the ability to simulate deep-water focused waves of a two-layer Boussinesq-type(BT)model,which has been shown to have excellent linear and nonlinear performance.To further improve the numerical accuracy and stability,the internal wavegenerated method is introduced into the two-layer Boussinesq-type model.Firstly,the sensitivity of the numerical results to the grid resolution is analyzed to verify the convergence of the model;secondly,the focused wave propagating in two opposite directions is simulated to prove the symmetry of the numerical results and the feasibility of the internal wave-generated method;thirdly,the limiting focused wave condition is simulated to compare and analyze the wave surface and the horizontal velocity of the profile at the focusing position,which is in good agreement with the measured values.Meanwhile the simulation of focused waves in very deep waters agrees well with the measured values,which further demonstrates the capability of the two-layer BT model in simulating focused waves in deep waters.

A Hybrid Dung Beetle Optimization Algorithm with Simulated Annealing for the Numerical Modeling of Asymmetric Wave Equations

In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.

Predicting hepatocellular carcinoma: A new non-invasive model based on shear wave elastography

BACKGROUND Integrating conventional ultrasound features with 2D shear wave elastography(2D-SWE)can potentially enhance preoperative hepatocellular carcinoma(HCC)predictions.AIM To develop a 2D-SWE-based predictive model for preoperative identification of HCC.METHODS A retrospective analysis of 884 patients who underwent liver resection and pathology evaluation from February 2021 to August 2023 was conducted at the Oriental Hepatobiliary Surgery Hospital.The patients were divided into the modeling group(n=720)and the control group(n=164).The study included conventional ultrasound,2D-SWE,and preoperative laboratory tests.Multiple logistic regression was used to identify independent predictive factors for RESULTS In the modeling group analysis,maximal elasticity(Emax)of tumors and their peripheries,platelet count,cirrhosis,and blood flow were independent risk indicators for malignancies.These factors yielded an area under the curve of 0.77(95%confidence interval:0.73-0.81)with 84%sensitivity and 61%specificity.The model demonstrated good calibration in both the construction and validation cohorts,as shown by the calibration graph and Hosmer-Lemeshow test(P=0.683 and P=0.658,respectively).Additionally,the mean elasticity(Emean)of the tumor periphery was identified as a risk factor for microvascular invasion(MVI)in malignant liver tumors(P=0.003).Patients receiving antiviral treatment differed significantly in platelet count(P=0.002),Emax of tumors(P=0.033),Emean of tumors(P=0.042),Emax at tumor periphery(P<0.001),and Emean at tumor periphery(P=0.003).CONCLUSION 2D-SWE’s hardness value serves as a valuable marker for enhancing the preoperative diagnosis of malignant liver lesions,correlating significantly with MVI and antiviral treatment efficacy.

Evaluation of Nonbreaking Wave-Induced Mixing Parameterization Schemes Based on a One-Dimensional Ocean Model

Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.

A Higher-Efficient Non-Hydrostatic Finite Volume Model for Strong Three-Dimensional Free Surface Flows and Sediment Transport

In order to accurately simulate strong three-dimensional （3-D） free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional （2-D） grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.

Modeling and Implementation of a Data Security and Protection Medium Using the Generated Key Based on Electromagnetic Wave Propagation Theories

Today, the advent of quantum computers and algorithms is calling into question the semantic security of symmetrical and asymmetrical cryptosystems. The security of objects connected to the network, which must provide a security service and protect the privacy of users by providing protection against attacks such as identity theft, denial of service, eavesdropping and unauthorised access to personal and sensitive data. It is therefore necessary to find a robust method of using the key that is effective in protecting and preventing data tampering. In this paper, we design and implement a security and data protection method using a key generated on the basis of electromagnetic wave propagation theories. Modelling and implementation of a data security and protection method using a key generated on the basis of electromagnetic wave propagation theories.

Electrostatic Attraction and Repulsion Explained and Modelled Mathematically Using Classical Physics—A Detailed Mechanism Based on Particle Wave Functions

The phenomenon of electrical attraction and repulsion between charged particles is well known, and described mathematically by Coulomb’s Law, yet until now there has been no explanation for why this occurs. There has been no mechanistic explanation that reveals what causes the charged particles to accelerate, either towards or away from each other. This paper gives a detailed explanation of the phenomena of electrical attraction and repulsion based on my previous work that determined the exact wave-function solutions for both the Electron and the Positron. It is revealed that the effects are caused by wave interactions between the wave functions that result in Electromagnetic reflections of parts of the particle’s wave functions, causing a change in their momenta.

Traveling Wave Solutions of a SIR Epidemic Model with Spatio-Temporal Delay

In this paper, we studied the traveling wave solutions of a SIR epidemic model with spatial-temporal delay. We proved that this result is determined by the basic reproduction number R0and the minimum wave speed c*of the corresponding ordinary differential equations. The methods used in this paper are primarily the Schauder fixed point theorem and comparison principle. We have proved that when R0>1and c>c*, the model has a non-negative and non-trivial traveling wave solution. However, for R01and c≥0or R0>1and 0cc*, the model does not have a traveling wave solution.

Atlantic blocking events in a simplified nonlinear baroclinic model for local finite-amplitude wave activity

为研究北大西洋阻高的形成机制,本文在局部有限振幅波活动(LWA)框架下进行了一系列数值实验.采用的数值模型能显式地描绘出两种重要的大气内部动力过程,即非线性纬向位涡通量和Rossby波包传播.模拟结果显示,这两种动力学过程均是形成大西洋阻高的重要机理.具体地,非线性纬向位涡通量和Rossby波包传播,分别是大西洋阻高南部和北部LWA形成的主导因子.因此,本研究综合了前人关于大西洋阻高的研究成果,为其形成机理提供了新的认识.

Vertical two-dimensional non-hydrostatic pressure model with single layer

The vertical two-dimensional non-hydrostatic pressure models with multiple layers can make prediction more accurate than those obtained by the hydrostatic pres- sure assumption. However, they are time-consuming and unstable, which makes them unsuitable for wider application. In this study, an efficient model with a single layer is developed. Decomposing the pressure into the hydrostatic and dynamic components and integrating the x-momentum equation from the bottom to the free surface can yield a horizontal momentum equation, in which the terms relevant to the dynamic pressure are discretized semi-implicitly. The convective terms in the vertical momentum equation are ignored, and the rest of the equation is approximated with the Keller-box scheme. The velocities expressed as the unknown dynamic pressure are substituted into the continuity equation, resulting in a tri-diagonal linear system solved by the Thomas algorithm. The validation of solitary and sinusoidal waves indicates that the present model can provide comparable results to the models with multiple layers but at much lower computation cost.