Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose...Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.展开更多
The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interf...The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the platecovered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presented to investigate the effects of different physical quantities, such as the thickness of the plate, Young’s modulus, the ratios of the densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid. Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.展开更多
A weak nonlinear model of a two-layer barotropic ocean with Rayleigh dissipation is built.The analytic asymptotic solution is derived in the mid-latitude stationary wind field,and the physical meaning of the correspon...A weak nonlinear model of a two-layer barotropic ocean with Rayleigh dissipation is built.The analytic asymptotic solution is derived in the mid-latitude stationary wind field,and the physical meaning of the corresponding problem is discussed.展开更多
In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of th...In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness dl, and lower layer thick-ness d2, instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehaute＇s plot for free surface waves if water depth ratio r= d1/d2 approaches to infinity and the upper layer water density p1 to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of σ=（P2 - Pl）/P2 → 1.0 and r 〉 1.0. In the end, several figures of the validity ranges for various interfacial wavetheories in the two-layer fluid are given and compared with the results for surface waves.展开更多
In this paper, the diffraction of water waves by a vertically floating cylinder in a two-layer fluid of a finite depth is studied. Analytical expressions for the hydrodynamic loads on the vertically floating cylinder ...In this paper, the diffraction of water waves by a vertically floating cylinder in a two-layer fluid of a finite depth is studied. Analytical expressions for the hydrodynamic loads on the vertically floating cylinder are obtained by use of the method of eigenfunction expansions. The hydrodynamic loads on the vertically floating cylinder in a two-layer fluid inelude not only the surge, heave and pitch exciting forces due to the incident wave of the surface-wave mode, but also those due to the incident wave of the internal-wave mode. This is different from the case of a homogenous fluid. Some given examples show that, for a two-layer fluid system with a small density difference, the hydrodynamic loads for the surface-wave mode do not differ significantly from those due to surface waves in a single-layer fluid, but the hydrodynamic loads for the internal-wave mode are important over a wide range of frequencies. Moreover, also considered are the free surface and interface elevations generated by the diffraction wave due to the incident wave of the surface-wave and interhal-wave modes, and transfer of energy between modes.展开更多
Many new forms of Boussinesq-type equations have been developed to extend the range of applicability of the classical Boussinesq equations to deeper water in the Study of the surface waves. One approach was used by Nw...Many new forms of Boussinesq-type equations have been developed to extend the range of applicability of the classical Boussinesq equations to deeper water in the Study of the surface waves. One approach was used by Nwogu （1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） to improve the linear dispersion characteristics of the classical Boussinesq equations by using the velocity at an arbitrary level as the velocity variable in derived equations and obtain a new form of Boussinesq-type equations, in which the dispersion property can be optimized by choosing the velocity variable at an adequate level. In this paper, a set of Boussinesq-type equations describing the motions of the interracial waves propagating alone the interface between two homogeneous incompressible and inviscid fluids of different densities with a free surface and a variable water depth were derived using a method similar to that used by Nwogu （1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） for surface waves. The equations were expressed in terms of the displacements of free surface and density-interface, and the velocity vectors at arbitrary vertical locations in the upper layer and the lower layer （or depth-averaged velocity vector across each layer） of a two-layer fluid. As expected, the equations derived in the present work include as special cases those obtained by Nwogu （1993, J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） and Peregrine （1967, J. Fluid Mech. 27, 815-827） for surface waves when the density of the upper fluid is taken as zero.展开更多
A quite general coupled variable coefficient modified KdV (VCmKdV) equation in a two-layer fluid systemis derived by means of the reductive perturbation method.Making use of the CK's direct method,some similarityr...A quite general coupled variable coefficient modified KdV (VCmKdV) equation in a two-layer fluid systemis derived by means of the reductive perturbation method.Making use of the CK's direct method,some similarityreductions of the coupled VCmKdV equation are obtained and their corresponding group explanations are discussed.Some exact solutions of the coupled equations are also presented.展开更多
The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an...The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an approximation for the free surface, and the lower one was bounded below by an impermeable bottom surface having a small deformation; the channel was unbounded in the horizontal directions. Assuming irrotational motion, the perturbation technique was employed to calculate the first-order corrections of the velocity potential in the two fluids by using Green＇s integral theorem suitably with the introduction of appropriate Green＇s functions. Those functions help in calculating the reflection and transmission coefficients in terms of integrals involving the shape ftmction c（x） representing the bottom deformation. Three-dimensional linear water wave theory was utilized for formulating the relevant boundary value problem. Two special examples of bottom deformation were considered to validate the results. Consideration of a patch of sinusoidal ripples （having the same wave number） shows that the reflection coefficient is an oscillatory function of the ratio of twice the x-component of the wave number to the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large. Similar results were observed for a patch of sinusoidal ripples having different wave numbers. It was also observed that for small angles of incidence, the reflected energy is greater compared to other angles of incidence up to π/ 4. These theoretical observations are supported by graphical results.展开更多
A previous study （Song. 2004. Geophys Res Lett, 31 （15）：L15302） of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into...A previous study （Song. 2004. Geophys Res Lett, 31 （15）：L15302） of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into a more general case of two-layer fluid with a top free surface. The rigid boundary condition on the upper surface is replaced by the kinematical and dynamical boundary conditions of a free surface, and the equations describing the random displacements of free surface, density-interface and the associated velocity potentials in the two-layer fluid are solved to the second order using the same expansion technology as that of Song （2004. Geophys Res Lett, 31 （15）：L15302）. The results show that the interface and the surface will oscillate synchronously, and the wave fields to the first-order both at the free surface and at the density-interface are made up of a linear superposition of many waves with different amplitudes, wave numbers and frequencies. The second-order solutions describe the second-order wave-wave interactions of the surface wave components, the interface wave components and among the surface and the interface wave components. The extended solutions also include special cases obtained by Thorpe for progressive interracial waves （Thorpe. 1968a.Trans R Soc London, 263A：563～614） and standing interracial waves （Thorpe. 1968b. J Fluid Mech, 32：489-528） for the two-layer fluid with a top free surface. Moreover, the solutions reduce to those derived for random surface waves by Sharma and Dean （1979.Ocean Engineering Rep 20） if the density of the upper layer is much smaller than that of the lower layer.展开更多
The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the fr...The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the free surface and the interface. This paper is concerned with the derivation of Green functions in the three dimensional case of a stationary source oscillating. The source point is located either in the upper or lower part of a two-layer fluid of finite depth. The derivation is carried out by the method of singularities. This method has an advantage in that it involves representing the potential as a sum of singularities or multipoles placed within any structures being present. Furthermore, experience shows that the systems of equations resulted from using a singularity method possess excellent convergence characteristics and only a few equations are needed to obtain accurate numerical results. Validation is done by showing that the derived two-layer Green function can be reduced to that of a single layer of finite depth or that the upper Green function coincides with that of the lower, for each case. The effect of the density on the internal waves is demonstrated. Also, it is shown how the surface and internal wave amplitudes are compared for both the wave modes. The fluid in this case is considered to be inviscid and incompressible and the flow is irrotational.展开更多
The coefficients embodied in a Boussinesq-type model are very important since they are determined to optimize the linear and nonlinear properties.In most conventional Boussinesq-type models,these coefficients are assi...The coefficients embodied in a Boussinesq-type model are very important since they are determined to optimize the linear and nonlinear properties.In most conventional Boussinesq-type models,these coefficients are assigned the specific values.As for the multi-layer Boussinesq-type models with the inclusion of the vertical velocity,however,the effect of the different values of these coefficients on linear and nonlinear performances has never been investigated yet.The present study focuses on a two-layer Boussinesq-type model with the highest spatial derivatives being 2 and theoretically and numerically examines the effect of the coefficient on model performance.Theoretical analysis show that different values for(0.13≤α≤0.25)do not have great effects on the high accuracy of the linear shoaling,linear phase celerity and even third-order nonlinearity for water depth range of 0<kh≤10(k is wave number and h is water depth).The corresponding errors using different values are restricted within 0.1%,0.1%and 1%for the linear shoaling amplitude,dispersion and nonlinear harmonics,respectively.Numerical tests including regular wave shoaling over mildly varying slope from deep to shallow water,regular wave propagation over submerged bar,bichromatic wave group and focusing wave propagation over deep water are conducted.The comparison between numerical results using different values of,experimental data and analytical solutions confirm the theoretical analysis.The flexibility and consistency of the two-layer Boussinesq-type model is therefore demonstrated theoretically and numerically.展开更多
Internal waves are one of the various phenomena that occur at sea,and they affect acoustic equipment and sea density measurement equipment.In this study,internal waves are simulated using computational fluid dynamics ...Internal waves are one of the various phenomena that occur at sea,and they affect acoustic equipment and sea density measurement equipment.In this study,internal waves are simulated using computational fluid dynamics method in the presence of a submarine in a pre-stratified fluid.Several scenarios were implemented by Froude number changes and submersible velocity by using the Navier-Stokes k-εturbulence model.Results indicate that the realizable k-εturbulence model gives better results than the RNG k-εmodel and the internal waves flow in this model are well represented,which increases the wavelength of the internal waves by increasing the Froude number and floating velocity,while the internal angle of the Kelvin waves is decreased.We also observe that increasing the floating velocity causes the turbulent velocity contours to increase due to the drag coefficient and its relationship with the Reynolds number.The Reynolds number increases with the increasing velocity of the float motion.The results indicate the efficiency of this method in the discovery of subsurface objects.展开更多
Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer fluid with a top free surface and a fiat bottom. The solutions were deduced from the general form of ...Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer fluid with a top free surface and a fiat bottom. The solutions were deduced from the general form of linear fluid dynamic equations of two-layer fluid under the f-plane approximation, and wave-induced tangential stress were estimated based on the solutions obtained. As expected, the solutions derived from the present work include as special cases those obtained by Sun et al. （2004. Science in China, Ser. D, 47（12）： 1147-1154） for geostrophic small amplitude surface wave solutions and wave-induced tangential stress if the density of the upper layer is much smaller than that of the lower layer. The results show that the interface and the surface will oscillate synchronously, and the influence of the earth＇s rotation both on the surface wave solutions and the interfacial wave solutions should be considered.展开更多
In this paper,the existence and propagation characteristics of Rossby waves in a two-layer cylindrical fluid are studied.Firstly,based on the dimensionless baroclinic quasi-geostrophic vortex equations including exoge...In this paper,the existence and propagation characteristics of Rossby waves in a two-layer cylindrical fluid are studied.Firstly,based on the dimensionless baroclinic quasi-geostrophic vortex equations including exogenous and dissipative,we derive new(2+1)-dimensional coupled Boussinesq equations describing wave propagation in polar coordinates by employing a multiscale analysis and perturbation method.Then,the Lie symmetries and conservation laws of the coupled Boussinesq equations are analyzed.Subsequently,by using the(G’/G)-expansion method,the exact solutions of the(2+1)-dimensional coupled Boussinesq equations are obtained.Finally,the effects of coupling term coefficients on the propagation characteristics of Rossby waves are analyzed.展开更多
There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle th...There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle these kind of mixed boundary-value problems associated with the Laplace’s equation （or Helmholtz equation） arising in the study of waves propagating through solids or fluids. One of the widely used methods in wave structure interaction is the multipole expansion method. This expansion involves a general combination of a regular wave, a wave source, a wave dipole and a regular wave-free part. The wave-free part can be further expanded in terms of wave-free multipoles which are termed as wave-free potentials. These are singular solutions of Laplace’s equation or two-dimensional Helmholz equation. Construction of these wave-free potentials and multipoles are presented here in a systematic manner for a number of situations such as two-dimensional non-oblique and oblique waves, three dimensional waves in two-layer fluid with free surface condition with higher order partial derivative are considered. In particular, these are obtained taking into account of the effect of the presence of surface tension at the free surface and also in the presence of an ice-cover modelled as a thin elastic plate. Also for limiting case, it can be shown that the multipoles and wave-free potential functions go over to the single layer multipoles and wave-free potential.展开更多
The kinetic characteristics of plasma-assisted oxidative pyrolysis of ammonia are studied by using the global/fluid models hybrid solution method.Firstly,the stable products of plasma-assisted oxidative pyrolysis of a...The kinetic characteristics of plasma-assisted oxidative pyrolysis of ammonia are studied by using the global/fluid models hybrid solution method.Firstly,the stable products of plasma-assisted oxidative pyrolysis of ammonia are measured.The results show that the consumption of NH_(3)/O_(2)and the production of N_(2)/H_(2)change linearly with the increase of voltage,which indicates the decoupling of nonequilibrium molecular excitation and oxidative pyrolysis of ammonia at low temperatures.Secondly,the detailed reaction kinetics mechanism of ammonia oxidative pyrolysis stimulated by a nanosecond pulse voltage at low pressure and room temperature is established.Based on the reaction path analysis,the simplified mechanism is obtained.The detailed and simplified mechanism simulation results are compared with experimental data to verify the accuracy of the simplified mechanism.Finally,based on the simplified mechanism,the fluid model of ammonia oxidative pyrolysis stimulated by the nanosecond pulse plasma is established to study the pre-sheath/sheath behavior and the resultant consumption and formation of key species.The results show that the generation,development,and propagation of the pre-sheath have a great influence on the formation and consumption of species.The consumption of NH_(3)by the cathode pre-sheath is greater than that by the anode pre-sheath,but the opposite is true for OH and O(1S).However,within the sheath,almost all reactions do not occur.Further,by changing the parameters of nanosecond pulse power supply voltage,it is found that the electron number density,electron current density,and applied peak voltages are not the direct reasons for the structural changes of the sheath and pre-sheath.Furthermore,the discharge interval has little effect on the sheath structure and gas mixture breakdown.The research results of this paper not only help to understand the kinetic promotion of non-equilibrium excitation in the process of oxidative pyrolysis but also help to explore the influence of transport and chemical reaction kinetics on the oxidative pyrolysis of ammonia.展开更多
Configuring computational fluid dynamics(CFD)simulations typically demands extensive domain expertise,limiting broader access.Although large language models(LLMs)have advanced scientific computing,their use in automat...Configuring computational fluid dynamics(CFD)simulations typically demands extensive domain expertise,limiting broader access.Although large language models(LLMs)have advanced scientific computing,their use in automating CFD workflows is underdeveloped.We introduce a novel approach centered on domain-specific LLM adaptation.By fine-tuning Qwen2.5-7B-Instruct on NL2FOAM,our custom dataset of 28,716 natural language-to-OpenFOAM configuration pairs with chain-of-thought(CoT)annotations enables direct translation from natural language descriptions to executable CFD setups.A multi-agent system orchestrates the process,autonomously verifying inputs,generating configurations,running simulations,and correcting errors.Evaluation on a benchmark of 21 diverse flow cases demonstrates state-of-the-art performance,achieving 88.7%solution accuracy and 82.6%first-attempt success rate.This significantly outperforms larger general-purpose models such as Qwen2.5-72B-Instruct,DeepSeek-R1,and Llama3.3-70B-Instruct,while also requiring fewer correction iterations and maintaining high computational efficiency.The results highlight the critical role of domain-specific adaptation in deploying LLM assistants for complex engineering workflows.Our code and fine-tuned model have been deposited at https://github.com/YYgroup/AutoCFD.展开更多
This work investigates water-based micropolar hybrid nanofluid(MHNF) flow on an elongating variable porous sheet.Nanoparticles of diamond and copper have been used in the water to boost its thermal conductivity. The m...This work investigates water-based micropolar hybrid nanofluid(MHNF) flow on an elongating variable porous sheet.Nanoparticles of diamond and copper have been used in the water to boost its thermal conductivity. The motion of the fluid is taken as two-dimensional with the impact of a magnetic field in the normal direction. The variable, permeable, and stretching nature of sheet's surface sets the fluid into motion. Thermal and mass diffusions are controlled through the use of the Cattaneo–Christov flux model. A dataset is generated using MATLAB bvp4c package solver and employed to train an artificial neural network(ANN) based on the Levenberg–Marquardt back-propagation(LMBP) algorithm. It has been observed as an outcome of this study that the modeled problem achieves peak performance at epochs 637, 112, 4848, and 344 using ANN-LMBP. The linear velocity of the fluid weakens with progression in variable porous and magnetic factors.With an augmentation in magnetic factor, the micro-rotational velocity profiles are augmented on the domain 0 ≤ η < 1.5 due to the support of micro-rotations by Lorentz forces close to the sheet's surface, while they are suppressed on the domain 1.5 ≤ η < 6.0 due to opposing micro-rotations away from the sheet's surface. Thermal distributions are augmented with an upsurge in thermophoresis, Brownian motion, magnetic, and radiation factors, while they are suppressed with an upsurge in thermal relaxation parameter. Concentration profiles increase with an expansion in thermophoresis factor and are suppressed with an intensification of Brownian motion factor and solute relaxation factor. The absolute errors(AEs) are evaluated for all the four scenarios that fall within the range 10^(-3)–10^(-8) and are associated with the corresponding ANN configuration that demonstrates a fine degree of accuracy.展开更多
This paper investigates the capabilities of large language models(LLMs)to leverage,learn and create knowledge in solving computational fluid dynamics(CFD)problems through three categories of baseline problems.These ca...This paper investigates the capabilities of large language models(LLMs)to leverage,learn and create knowledge in solving computational fluid dynamics(CFD)problems through three categories of baseline problems.These categories include(1)conventional CFD problems that can be solved using existing numerical methods in LLMs,such as lid-driven cavity flow and the Sod shock tube problem;(2)problems that require new numerical methods beyond those available in LLMs,such as the recently developed Chien-physics-informed neural networks for singularly perturbed convection-diffusion equations;and(3)problems that cannot be solved using existing numerical methods in LLMs,such as the ill-conditioned Hilbert linear algebraic systems.The evaluations indicate that reasoning LLMs overall outperform non-reasoning models in four test cases.Reasoning LLMs show excellent performance for CFD problems according to the tailored prompts,but their current capability in autonomous knowledge exploration and creation needs to be enhanced.展开更多
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 dens...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.展开更多
基金Supported by the National Natural Science Foundation of China(Grant No.40674063)National Hi-tech Research and Development Program of China(863Program)(Grant No.2006AA09Z311)
文摘Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.
基金sponsored by the National Basic Research Program of China(973 Program,Grant No.2014CB046203)the National Natural Science Foundation of China(Grant No.11072140)
文摘The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the platecovered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presented to investigate the effects of different physical quantities, such as the thickness of the plate, Young’s modulus, the ratios of the densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid. Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.
基金Project supported by the National Basic Research Program of China (Grant No. 2011CB403501)the National Natural Science Foundation of China (GrantNos. 41175058,41275062,and 11202106)
文摘A weak nonlinear model of a two-layer barotropic ocean with Rayleigh dissipation is built.The analytic asymptotic solution is derived in the mid-latitude stationary wind field,and the physical meaning of the corresponding problem is discussed.
基金the Knowledge Innovation Project of CAS(KJCX-YW-L02)the National 863 Project of China(2006AAO9A103-4)+1 种基金China National Oil Corporation in Beijing(CNOOC)the National Natural Science Foundation of China(10672056).
文摘In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness dl, and lower layer thick-ness d2, instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehaute＇s plot for free surface waves if water depth ratio r= d1/d2 approaches to infinity and the upper layer water density p1 to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of σ=（P2 - Pl）/P2 → 1.0 and r 〉 1.0. In the end, several figures of the validity ranges for various interfacial wavetheories in the two-layer fluid are given and compared with the results for surface waves.
基金This project was financially supported by the National Natural Science Foundation of China(Grant No.10572092)the High Technology Research and Development Programof China(863Program,Grant Nos.2006AA09Z352 and 2006AA09196-6)
文摘In this paper, the diffraction of water waves by a vertically floating cylinder in a two-layer fluid of a finite depth is studied. Analytical expressions for the hydrodynamic loads on the vertically floating cylinder are obtained by use of the method of eigenfunction expansions. The hydrodynamic loads on the vertically floating cylinder in a two-layer fluid inelude not only the surge, heave and pitch exciting forces due to the incident wave of the surface-wave mode, but also those due to the incident wave of the internal-wave mode. This is different from the case of a homogenous fluid. Some given examples show that, for a two-layer fluid system with a small density difference, the hydrodynamic loads for the surface-wave mode do not differ significantly from those due to surface waves in a single-layer fluid, but the hydrodynamic loads for the internal-wave mode are important over a wide range of frequencies. Moreover, also considered are the free surface and interface elevations generated by the diffraction wave due to the incident wave of the surface-wave and interhal-wave modes, and transfer of energy between modes.
文摘Many new forms of Boussinesq-type equations have been developed to extend the range of applicability of the classical Boussinesq equations to deeper water in the Study of the surface waves. One approach was used by Nwogu （1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） to improve the linear dispersion characteristics of the classical Boussinesq equations by using the velocity at an arbitrary level as the velocity variable in derived equations and obtain a new form of Boussinesq-type equations, in which the dispersion property can be optimized by choosing the velocity variable at an adequate level. In this paper, a set of Boussinesq-type equations describing the motions of the interracial waves propagating alone the interface between two homogeneous incompressible and inviscid fluids of different densities with a free surface and a variable water depth were derived using a method similar to that used by Nwogu （1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） for surface waves. The equations were expressed in terms of the displacements of free surface and density-interface, and the velocity vectors at arbitrary vertical locations in the upper layer and the lower layer （or depth-averaged velocity vector across each layer） of a two-layer fluid. As expected, the equations derived in the present work include as special cases those obtained by Nwogu （1993, J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） and Peregrine （1967, J. Fluid Mech. 27, 815-827） for surface waves when the density of the upper fluid is taken as zero.
基金The project supported by National Natural Science Foundation of China under Grant Nos.10547124,10475055,and 90503006the Youth Foundation of Shanghai Jiao Tong University
文摘A quite general coupled variable coefficient modified KdV (VCmKdV) equation in a two-layer fluid systemis derived by means of the reductive perturbation method.Making use of the CK's direct method,some similarityreductions of the coupled VCmKdV equation are obtained and their corresponding group explanations are discussed.Some exact solutions of the coupled equations are also presented.
文摘The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an approximation for the free surface, and the lower one was bounded below by an impermeable bottom surface having a small deformation; the channel was unbounded in the horizontal directions. Assuming irrotational motion, the perturbation technique was employed to calculate the first-order corrections of the velocity potential in the two fluids by using Green＇s integral theorem suitably with the introduction of appropriate Green＇s functions. Those functions help in calculating the reflection and transmission coefficients in terms of integrals involving the shape ftmction c（x） representing the bottom deformation. Three-dimensional linear water wave theory was utilized for formulating the relevant boundary value problem. Two special examples of bottom deformation were considered to validate the results. Consideration of a patch of sinusoidal ripples （having the same wave number） shows that the reflection coefficient is an oscillatory function of the ratio of twice the x-component of the wave number to the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large. Similar results were observed for a patch of sinusoidal ripples having different wave numbers. It was also observed that for small angles of incidence, the reflected energy is greater compared to other angles of incidence up to π/ 4. These theoretical observations are supported by graphical results.
基金supported by the National Science Foundation for Distinguished Young Scholars of China under contract No.40425015the Cooperative Project of Chinese Academy Sciences and the China National 0ffshore 0il Corporation("Behaviours of internal waves and their roles on the marine stuctures").
文摘A previous study （Song. 2004. Geophys Res Lett, 31 （15）：L15302） of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into a more general case of two-layer fluid with a top free surface. The rigid boundary condition on the upper surface is replaced by the kinematical and dynamical boundary conditions of a free surface, and the equations describing the random displacements of free surface, density-interface and the associated velocity potentials in the two-layer fluid are solved to the second order using the same expansion technology as that of Song （2004. Geophys Res Lett, 31 （15）：L15302）. The results show that the interface and the surface will oscillate synchronously, and the wave fields to the first-order both at the free surface and at the density-interface are made up of a linear superposition of many waves with different amplitudes, wave numbers and frequencies. The second-order solutions describe the second-order wave-wave interactions of the surface wave components, the interface wave components and among the surface and the interface wave components. The extended solutions also include special cases obtained by Thorpe for progressive interracial waves （Thorpe. 1968a.Trans R Soc London, 263A：563～614） and standing interracial waves （Thorpe. 1968b. J Fluid Mech, 32：489-528） for the two-layer fluid with a top free surface. Moreover, the solutions reduce to those derived for random surface waves by Sharma and Dean （1979.Ocean Engineering Rep 20） if the density of the upper layer is much smaller than that of the lower layer.
基金supported by the National Natural Science Foundation of China (Grant No. 50779008)
文摘The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the free surface and the interface. This paper is concerned with the derivation of Green functions in the three dimensional case of a stationary source oscillating. The source point is located either in the upper or lower part of a two-layer fluid of finite depth. The derivation is carried out by the method of singularities. This method has an advantage in that it involves representing the potential as a sum of singularities or multipoles placed within any structures being present. Furthermore, experience shows that the systems of equations resulted from using a singularity method possess excellent convergence characteristics and only a few equations are needed to obtain accurate numerical results. Validation is done by showing that the derived two-layer Green function can be reduced to that of a single layer of finite depth or that the upper Green function coincides with that of the lower, for each case. The effect of the density on the internal waves is demonstrated. Also, it is shown how the surface and internal wave amplitudes are compared for both the wave modes. The fluid in this case is considered to be inviscid and incompressible and the flow is irrotational.
基金supported by the National Natural Science Foundation of China(Grant Nos.51779022,51809053,and 51579034)the Innovation Team Project of Estuary and Coast Protection and Management(Grant No.Y220013)the Open Project Fund of State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology(Grant No.LP19015).
文摘The coefficients embodied in a Boussinesq-type model are very important since they are determined to optimize the linear and nonlinear properties.In most conventional Boussinesq-type models,these coefficients are assigned the specific values.As for the multi-layer Boussinesq-type models with the inclusion of the vertical velocity,however,the effect of the different values of these coefficients on linear and nonlinear performances has never been investigated yet.The present study focuses on a two-layer Boussinesq-type model with the highest spatial derivatives being 2 and theoretically and numerically examines the effect of the coefficient on model performance.Theoretical analysis show that different values for(0.13≤α≤0.25)do not have great effects on the high accuracy of the linear shoaling,linear phase celerity and even third-order nonlinearity for water depth range of 0<kh≤10(k is wave number and h is water depth).The corresponding errors using different values are restricted within 0.1%,0.1%and 1%for the linear shoaling amplitude,dispersion and nonlinear harmonics,respectively.Numerical tests including regular wave shoaling over mildly varying slope from deep to shallow water,regular wave propagation over submerged bar,bichromatic wave group and focusing wave propagation over deep water are conducted.The comparison between numerical results using different values of,experimental data and analytical solutions confirm the theoretical analysis.The flexibility and consistency of the two-layer Boussinesq-type model is therefore demonstrated theoretically and numerically.
文摘Internal waves are one of the various phenomena that occur at sea,and they affect acoustic equipment and sea density measurement equipment.In this study,internal waves are simulated using computational fluid dynamics method in the presence of a submarine in a pre-stratified fluid.Several scenarios were implemented by Froude number changes and submersible velocity by using the Navier-Stokes k-εturbulence model.Results indicate that the realizable k-εturbulence model gives better results than the RNG k-εmodel and the internal waves flow in this model are well represented,which increases the wavelength of the internal waves by increasing the Froude number and floating velocity,while the internal angle of the Kelvin waves is decreased.We also observe that increasing the floating velocity causes the turbulent velocity contours to increase due to the drag coefficient and its relationship with the Reynolds number.The Reynolds number increases with the increasing velocity of the float motion.The results indicate the efficiency of this method in the discovery of subsurface objects.
基金supported by the National Science Fund for Distinguished Young Scholars of China undercontract No 40425015 the Knowledge Innovation Programs of the Chinese Academy of Sciences under contract Nos KZCX1-YW-12and KZCX2-YW-201
文摘Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer fluid with a top free surface and a fiat bottom. The solutions were deduced from the general form of linear fluid dynamic equations of two-layer fluid under the f-plane approximation, and wave-induced tangential stress were estimated based on the solutions obtained. As expected, the solutions derived from the present work include as special cases those obtained by Sun et al. （2004. Science in China, Ser. D, 47（12）： 1147-1154） for geostrophic small amplitude surface wave solutions and wave-induced tangential stress if the density of the upper layer is much smaller than that of the lower layer. The results show that the interface and the surface will oscillate synchronously, and the influence of the earth＇s rotation both on the surface wave solutions and the interfacial wave solutions should be considered.
基金Supported by the Nature Science Foundation of Shandong Province of China(ZR2018MA017)。
文摘In this paper,the existence and propagation characteristics of Rossby waves in a two-layer cylindrical fluid are studied.Firstly,based on the dimensionless baroclinic quasi-geostrophic vortex equations including exogenous and dissipative,we derive new(2+1)-dimensional coupled Boussinesq equations describing wave propagation in polar coordinates by employing a multiscale analysis and perturbation method.Then,the Lie symmetries and conservation laws of the coupled Boussinesq equations are analyzed.Subsequently,by using the(G’/G)-expansion method,the exact solutions of the(2+1)-dimensional coupled Boussinesq equations are obtained.Finally,the effects of coupling term coefficients on the propagation characteristics of Rossby waves are analyzed.
文摘There is a large class of problems in the field of fluid structure interaction where higher-order boundary conditions arise for a second-order partial differential equation. Various methods are being used to tackle these kind of mixed boundary-value problems associated with the Laplace’s equation （or Helmholtz equation） arising in the study of waves propagating through solids or fluids. One of the widely used methods in wave structure interaction is the multipole expansion method. This expansion involves a general combination of a regular wave, a wave source, a wave dipole and a regular wave-free part. The wave-free part can be further expanded in terms of wave-free multipoles which are termed as wave-free potentials. These are singular solutions of Laplace’s equation or two-dimensional Helmholz equation. Construction of these wave-free potentials and multipoles are presented here in a systematic manner for a number of situations such as two-dimensional non-oblique and oblique waves, three dimensional waves in two-layer fluid with free surface condition with higher order partial derivative are considered. In particular, these are obtained taking into account of the effect of the presence of surface tension at the free surface and also in the presence of an ice-cover modelled as a thin elastic plate. Also for limiting case, it can be shown that the multipoles and wave-free potential functions go over to the single layer multipoles and wave-free potential.
基金Fundamental Research Funds for the Central Universities(M23JBZY00050)National Natural Science Foundation of China(22278032)。
文摘The kinetic characteristics of plasma-assisted oxidative pyrolysis of ammonia are studied by using the global/fluid models hybrid solution method.Firstly,the stable products of plasma-assisted oxidative pyrolysis of ammonia are measured.The results show that the consumption of NH_(3)/O_(2)and the production of N_(2)/H_(2)change linearly with the increase of voltage,which indicates the decoupling of nonequilibrium molecular excitation and oxidative pyrolysis of ammonia at low temperatures.Secondly,the detailed reaction kinetics mechanism of ammonia oxidative pyrolysis stimulated by a nanosecond pulse voltage at low pressure and room temperature is established.Based on the reaction path analysis,the simplified mechanism is obtained.The detailed and simplified mechanism simulation results are compared with experimental data to verify the accuracy of the simplified mechanism.Finally,based on the simplified mechanism,the fluid model of ammonia oxidative pyrolysis stimulated by the nanosecond pulse plasma is established to study the pre-sheath/sheath behavior and the resultant consumption and formation of key species.The results show that the generation,development,and propagation of the pre-sheath have a great influence on the formation and consumption of species.The consumption of NH_(3)by the cathode pre-sheath is greater than that by the anode pre-sheath,but the opposite is true for OH and O(1S).However,within the sheath,almost all reactions do not occur.Further,by changing the parameters of nanosecond pulse power supply voltage,it is found that the electron number density,electron current density,and applied peak voltages are not the direct reasons for the structural changes of the sheath and pre-sheath.Furthermore,the discharge interval has little effect on the sheath structure and gas mixture breakdown.The research results of this paper not only help to understand the kinetic promotion of non-equilibrium excitation in the process of oxidative pyrolysis but also help to explore the influence of transport and chemical reaction kinetics on the oxidative pyrolysis of ammonia.
基金supported by the National Natural Science Foundation of China(Grant Nos.52306126,22350710788,12432010,11988102,92270203)the Xplore Prize.
文摘Configuring computational fluid dynamics(CFD)simulations typically demands extensive domain expertise,limiting broader access.Although large language models(LLMs)have advanced scientific computing,their use in automating CFD workflows is underdeveloped.We introduce a novel approach centered on domain-specific LLM adaptation.By fine-tuning Qwen2.5-7B-Instruct on NL2FOAM,our custom dataset of 28,716 natural language-to-OpenFOAM configuration pairs with chain-of-thought(CoT)annotations enables direct translation from natural language descriptions to executable CFD setups.A multi-agent system orchestrates the process,autonomously verifying inputs,generating configurations,running simulations,and correcting errors.Evaluation on a benchmark of 21 diverse flow cases demonstrates state-of-the-art performance,achieving 88.7%solution accuracy and 82.6%first-attempt success rate.This significantly outperforms larger general-purpose models such as Qwen2.5-72B-Instruct,DeepSeek-R1,and Llama3.3-70B-Instruct,while also requiring fewer correction iterations and maintaining high computational efficiency.The results highlight the critical role of domain-specific adaptation in deploying LLM assistants for complex engineering workflows.Our code and fine-tuned model have been deposited at https://github.com/YYgroup/AutoCFD.
基金the Deanship of Research and Graduate Studies at King Khalid University for funding this work through large Research Group Project (Grant No. RGP2/198/45)Project supported by Prince Sattam bin Abdulaziz University (Grant No. PSAU/2025/R/1446)。
文摘This work investigates water-based micropolar hybrid nanofluid(MHNF) flow on an elongating variable porous sheet.Nanoparticles of diamond and copper have been used in the water to boost its thermal conductivity. The motion of the fluid is taken as two-dimensional with the impact of a magnetic field in the normal direction. The variable, permeable, and stretching nature of sheet's surface sets the fluid into motion. Thermal and mass diffusions are controlled through the use of the Cattaneo–Christov flux model. A dataset is generated using MATLAB bvp4c package solver and employed to train an artificial neural network(ANN) based on the Levenberg–Marquardt back-propagation(LMBP) algorithm. It has been observed as an outcome of this study that the modeled problem achieves peak performance at epochs 637, 112, 4848, and 344 using ANN-LMBP. The linear velocity of the fluid weakens with progression in variable porous and magnetic factors.With an augmentation in magnetic factor, the micro-rotational velocity profiles are augmented on the domain 0 ≤ η < 1.5 due to the support of micro-rotations by Lorentz forces close to the sheet's surface, while they are suppressed on the domain 1.5 ≤ η < 6.0 due to opposing micro-rotations away from the sheet's surface. Thermal distributions are augmented with an upsurge in thermophoresis, Brownian motion, magnetic, and radiation factors, while they are suppressed with an upsurge in thermal relaxation parameter. Concentration profiles increase with an expansion in thermophoresis factor and are suppressed with an intensification of Brownian motion factor and solute relaxation factor. The absolute errors(AEs) are evaluated for all the four scenarios that fall within the range 10^(-3)–10^(-8) and are associated with the corresponding ANN configuration that demonstrates a fine degree of accuracy.
基金supported by the National Natural Science Foundation of China Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”(Grant No.11988102)the National Natural Science Foundation of China(Grant No.12202451).
文摘This paper investigates the capabilities of large language models(LLMs)to leverage,learn and create knowledge in solving computational fluid dynamics(CFD)problems through three categories of baseline problems.These categories include(1)conventional CFD problems that can be solved using existing numerical methods in LLMs,such as lid-driven cavity flow and the Sod shock tube problem;(2)problems that require new numerical methods beyond those available in LLMs,such as the recently developed Chien-physics-informed neural networks for singularly perturbed convection-diffusion equations;and(3)problems that cannot be solved using existing numerical methods in LLMs,such as the ill-conditioned Hilbert linear algebraic systems.The evaluations indicate that reasoning LLMs overall outperform non-reasoning models in four test cases.Reasoning LLMs show excellent performance for CFD problems according to the tailored prompts,but their current capability in autonomous knowledge exploration and creation needs to be enhanced.
基金financially supported by the Institut Teknologi Bandung Research(Grant No.107a/I1.C01/PL/2017)
文摘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.
