Numerical Calculation of Seafloor Synthetic Seismograms Caused by Low Frequency Point Sound Source

Elastic wave on seafloor caused by low frequency noise radiated from ship is called ship seismic wave which can be used to identify ship target. In order to analyze the wave components and the propagating properties of ship seismic wave, the numerical calculation of synthetic seismograms on seafloor aroused by a low frequency point sound source is carried out using a wave number integration technique combined with inverse Fourier transform. According to the numerical example of hard seafloor, the time series of seismic wave on seafloor are mostly composed of interface waves and normal mode waves. Each normal mode wave has a well defined low cut-off frequency, while the interface wave doesn't have. The frequency dispersion of normal mode wave is obvious when frequency is lower than 100Hz, while the interface wave is dispersive only in the infra-sound frequency range. The time series of seismic wave is dominated by the interface wave when the source frequency is less than the minimal cut-off frequency of normal mode wave.

Numerical Calculation for Nonlinear Waves in Water of Arbitrarily Varying Depth with Boussinesq Equations

Based on the high order nonlinear and dispersive wave equation with a dissipative term, a numerical model for nonlinear waves is developed, It is suitable to calculate wave propagation in water areas with an arbitrarily varying bottom slope and a relative depth h/L(0)less than or equal to1. By the application of the completely implicit stagger grid and central difference algorithm, discrete governing equations are obtained. Although the central difference algorithm of second-order accuracy both in time and space domains is used to yield the difference equations, the order of truncation error in the difference equation is the same as that of the third-order derivatives of the Boussinesq equation. In this paper, the correction to the first-order derivative is made, and the accuracy of the difference equation is improved. The verifications of accuracy show that the results of the numerical model are in good agreement with those of analytical Solutions and physical models.

Interior Performances of an Impulse Thruster by Numerical Calculation

Impulse thruster is a kind of actuator used for trajectory correction or attitude control of some in-flight munitions and vehicles.A simple mathematical model was set up to model interior characteristics of an impulse thruster.With this model,effects of some key parameters on interior performances of the impulse thruster were studied.Results show that action time is affected significantly by nozzle throat diameter and volume of combustion chamber,while output impulse is sensitive to charge particle diameter and nozzle throat diameter.Through numerical calculation,the ranges of the optimized values for some key parameters were obtained.

A three-dimensional numerical calculation of the wind-driven thermohaline and tide-induced Lagrangian residual current in the Bohai Sea

On the basis of a three-dimensional weakly nonliear theory of Lagrangian residual current in the Baroclinic shallow seas, a diagnostic numerical calculation of wind-driven, thermohaline and tide-induced Lagrangian residual current in the Bohai Sea is made. The model involves the Richardson number in the eddy viscosity coefficient, wind, thcrmolialine and tidal effects in the focing terms. The runoff of the Huanghe River and a part of the Huanghai Warm Water coming from the Huanghai Sea through the Bohai Sea Strait is also considered. The velocity-splitting method is adopted. The wind-driven circu lation, thermohaline circulation and the tide-induced Lagrangian residual circulation are also obtained individually and analysed. The dynamics of the three main eddies in the Lagrangian mean circulation is discussed. Finally, the numerical result is partly verified with the observed data.

Numerical calculation on temperature field of FGH95 alloy droplet during rapid solidification

The temperature field of FGH95 alloy droplet atomized by plasma rotatingelectrode processing (PREP) during solidification has been calculated through numerical analysisbased on equivalent sensible heat capacity method. And thus the relational curves among temperaturegradient of solid-liquid interface, moving velocity of solid-liquid interface and solid fractionduring solidification have been presented. The results indicate that the relation between averagetemperature gradient of solid-liquid interface and droplet size, and the relation between averagemoving velocity of solid-liquid interface and droplet size can be expressed during solidification.

NUMERICAL CALCULATIONS OF GASEOUS REACTING FLOWS IN A MODEL OF GAS TURBINE COMBUSTORS

This paper describes the numerical calculations of gaseous reaction flows in a model of gas turbine combustors. The profiles of hydrodynamic and thermodynamic patterns in a three-dimensional combustor model are obtained by solving the governing differential transport equations. The well-established numerical prediction algorithm SIMPLE, the modified k-ε turbulence model and k-ε-g turbulent diffusion flame model have been adopted in computations. The β function has been selected as probability density function. The effect of combustion process on flow patterns has been investigated. The calculated results have been verified by experiments. They are in remarkably good agreement.

Aero-thermo-elastic Numerical Calculation and Analysis of Spinning Rocket

The application of computational fluid dynamics/computational solid method(CFD/CSM)on solving the aero-thermo-elastic problem of spinning rocket is introduced.Firstly,the aerodynamic coefficients of a rocket are calculated,and the results are compared with the available experimental data,which verified the accuracy of the CFD output.Then,analysis is carried using ANSYS Workbench multi-physics coupling platform,which includes fluid,thermal,and structural solvers.The results show that spinning causes a significant effect on the deformations and stresses.Furthermore,thermal stresses due to high temperature at the rocket warhead and tail edges have a dominated effect,even more than those produced by aerodynamic forces.Consequently,this important outcome should be taken into consideration during the rocket design stages.

Optimization Method for Indoor Thermal Comfort Based on Interactive Numerical Calculation

In order to implement the optimal design of the indoor thermal comfort based on the numerical modeling method, the numerical calculation platform is combined seamlessly with the data-processing platform, and an interactive numerical calculation platform which includes the functions of numerical simulation and optimization is established. The artificial neural network （ANN） and the greedy strategy are introduced into the hill-climbing pattern heuristic search process, and the optimizing search direction can be predicted by using small samples; when searching along the direction using the greedy strategy, the optimal values can be quickly approached. Therefore, excessive external calling of the numerical modeling process can be avoided, and the optimization time is decreased obviously. The experimental results indicate that the satisfied output parameters of air conditioning can be quickly given out based on the interactive numerical calculation platform and the improved search method, and the optimization for indoor thermal comfort can be completed.

TWO DIMENSIONAL NUMERICAL CALCULATION FOR SALT WATER INTRUSION AT ESTUARIES

Based on 2-dimensional vertically governing equations of salt water intrusion at estuaries, a method is obtained, which can predict quantitatively current velocity and salinity distribution along depth when the intrusion occurs. The volume-controlling method proposed by Patanker and Spalding is used to form the computational pattern and the Power- Law Scheme is used as the diversion pattern of the diffusion term. The comparison between the computational results and the measured ones gives a satisfactory agreement.

Numerical calculation of flow and heat transfer process in the new-type external combustion swirl-flowing hot stove

It is clarified that the important method to improve the blast temperature ofthe small and the middle blast furnaces whose production is about two-thirds of total sum of Chinafrom 1000℃ to 1250-1300℃ is to preheat both their combustion-supporting air and coal gas. The airtemperature of blast furnaces can be reached to 1250-1300℃ by burning single blast furnace coal gasif high speed burner is applied to blast furnaces and new-type external combustion swirl-flowinghot stove is used to preheat their combustion-supporting air. The computational results of the flowand heat transfer processions in the hot stove prove that the surface of the bed of the thermalstorage balls there have not eccentric flow and the flow field and temperature field distribution iseven. The computational results of the blast temperature distribution are similar to thosedetermination experiment data. The numerical results also provide references for developing anddesigning the new-type external combustion swirl-flowing hot stoves.

Numerical calculation of drift current in wind waves

-Drift current induced by wind and waves is investigated with phase-averaged Navier-Stokes equation in which the Reynolds stress is closed by k-ε model. The governing equations are solved by the finite volume method in a system of nonorthogonal coordinates which is fitted to the phase-averaged wave surface. The predicted drift current is fairly reasonable and the drag coefficient of sea-surface predicted with the newly developed interface conditions shows good agreement with previous measurements when breaking waves do not exist.

An improved algorithm for numerical calculation of seismic response spectra

The information of seismic response spectra is key to many problems concerned with aseismic structure and is also helpful for earthquake disaster relief if it is generated in time when earthquake happens. While current numerical calculation methods suffer from poor precision, especially in frequency band near Nyquist frequency, we present a set of improved parameters for precision improvement. It is shown that precision of displacement and velocity response spectra are both further improved compared to current numerical algorithms. A uniform fitting formula is given for computing these parameters for damping ratio range of 0.01-0.9, quite convenient for practical application.

Optimal Design of High-Speed Partial Flow Pumps using Orthogonal Tests and Numerical Simulations

To investigate the influence of structural parameters on the performances and internal flow characteristics of partial flow pumps at a low specific speed of 10000 rpm,special attention was paid to the first and second stage impeller guide vanes.Moreover,the impeller blade outlet width,impeller inlet diameter,blade inclination angle,and number of blades were considered for orthogonal tests.Accordingly,nine groups of design solutions were formed,and then used as a basis for the execution of numerical simulations(CFD)aimed at obtaining the efficiency values and heads for each design solution group.The influence of impeller geometric parameters on the efficiency and head was explored,and the“weight”of each factor was obtained via a range analysis.Optimal structural parameters were finally chosen on the basis of the numerical simulation results,and the performances of the optimized model were verified accordingly(yet by means of CFD).Evidence is provided that the increase in the efficiency and head of the optimized model was 12.11%and 23.5 m,respectively,compared with those of the original model.

A Numerical Study on the Effect of the Backflow Hole Position on the Performances of a Self-Priming Pump

A self-priming pump is a centrifugal pump that has the ability to prime itself. Typically, its performance dependson the configuration of its reflux hole. In this study, the ANSYS FLUENT software is used to investigate the effectsof three different radial positions of the reflux hole on gas-liquid two-phase distribution, pressure pulsation, andimp during self-priming. The research results indicate that: (1) The effective channel size for the reflux liquid toenter the volute varies depending on the location of the reflux hole. The effect of the impeller rotation on thereflux liquid becomes more obvious as the setting distance of the reflux aperture decreases. (2) The position ofthe reflux hole significantly affects the gas phase mass fraction inside the impeller, resulting in a significant reductionin the time it takes for the mass fraction to exceed 80%. (3) The position of the reflux hole significantly affectsthe average pressure on each monitoring surface. (4) Placing the reflux hole at a excessively distant radial distancecan result in an excessive vertical component. (5) The self-priming performance of the pump can be improved tosome extent by placing the return hole at a small radial distance.

Numerical simulation and experimental verification of plasma jet development in gas gap switch

Plasma jet triggered gas gap switch has obvious advantages in fast control switch.The development of the plasma in the ambient medium is the key factor affecting the triggering conduction of the gas switch.However,the plasma jet process and its characteristic parameters are complicated and the existing test methods cannot fully characterize its development laws.In this work,a two-dimensional transient fluid calculation model of the plasma jet process of the gas gap switch is established based on the renormalization-group k-εturbulence equation.The results show that the characteristic parameters and morphological evolution of the plasma jet are basically consistent with the experimental results,which verifies the accuracy of the simulation model calculation.The plasma jet is a long strip with an initial velocity of 1.0 km·s-1and develops in both axial and radial directions.The jet velocity fluctuates significantly with axial height.As the plasma jet enters the main gap,the pressure inside the trigger cavity drops by80%,resulting in a rapid drop in the jet velocity.When the plasma jet head interacts with the atmosphere,the two-phase fluid compresses each other,generating a forward-propelled pressure wave.The plasma jet heads flow at high velocity,a negative pressure zone is formed in the middle part of the jet,and the pressure peak decreases gradually with height.As the value of the inlet pressure increases,the characteristic parameters of the plasma jet increase.The entrainment phenomenon is evident,which leads to an increase in the pressure imbalance of the atmospheric gas medium,leading to a significant Coanda effect.Compared with air,the characteristic parameters of a plasma jet in SF6are lower,and the morphological evolution is significantly suppressed.The results of this study can provide some insight into the mechanism of action of the switch jet plasma development process.

EXPERIMENTAL STUDIES AND NUMERICAL CALCULATION ON TURBULENT BOUNDARY LAYER OF WATER FLOW ALONG CURVED SURFACE

This paper will introduce experimental studies and numerical calculation on turbulent boundary layer of water flow along curved surface in our country in recent years.On the basis of the experimental studies,the effects of curvature and roughness on velocity distribution and pressure distribution and the change of turbulent flow boundary layer on overflow bucket concave surface is discussed.We proposed the empirical formula of velocity, pressure and the change of turbulent flow boundary layer on outlet bucket concave.According to the momentum principle,we have deduced the momentum integral equation full water depth boundary layer and using the element as control unit inside the boundary layer on concave surface of bucket.Combining with continuity equation,we have computed the boundary layer development on the bucket of a spillway.Compared with the field experimental data,the calculation results are available.Under polar coordinates,a mathematical model for simulating time-averaged flow characteristics in concave surface of bucket is established.The turbulent flow field on concave surface of bucket is calculated by SIMPLE method and this mathematical model.The flow velocity field,pressure field,distribution of turbulent kinetic energy, distribution of turbulent energy dissipating rate and distribution of shear stress are available.The calculation value is consistent with measured test data.

NUMERICAL CALCULATION OF DILUTE GAS-PARTICLE FLOWS IN A CYLINDRICAL CONTAINER WITH A ROTATING COVER

A new computing technique is described for the solution'of dilute gas-particle fluid flow problems.The com- mon SIMPLE method is used for the calculation of the gas phase,but for the calculation of the particle phase the MacCormack method(80 style)is used.Using the technique for gas-particle flows in a rotating cylindrical con- tainer,the field distributions of both phases are obtained.An important parameter Sk(Stokes number)is taken re- spectively as 0.01,0.1 and 1.The results show the influence of Sk on the flows.The collisions between particles and the side wall of the container are predicted when Sk=1.

NUMERICAL CALCULATION OF SINGLE BUBBLE MOVEMENT AND ITS CAVITATION BEHAVIOR

This paper presents analyses and numerical calculations of a single bubble movement and its cavitation behav- ior where the influence of air diffusion,pressure fluctuation,relative velocity between the bubble and fluid and the growing bubble burst effect are considered.The calculational resuits provides rather detailed description of a single bubble movement and its cavitation behavior.

Review of Collocation Methods and Applications in Solving Science and Engineering Problems

The collocation method is a widely used numerical method for science and engineering problems governed by partial differential equations.This paper provides a comprehensive review of collocation methods and their applications,focused on elasticity,heat conduction,electromagnetic field analysis,and fluid dynamics.The merits of the collocation method can be attributed to the need for element mesh,simple implementation,high computational efficiency,and ease in handling irregular domain problems since the collocation method is a type of node-based numerical method.Beginning with the fundamental principles of the collocation method,the discretization process in the continuous domain is elucidated,and how the collocation method approximation solutions for solving differential equations are explained.Delving into the historical development of the collocation methods,their earliest applications and key milestones are traced,thereby demonstrating their evolution within the realm of numerical computation.The mathematical foundations of collocation methods,encompassing the selection of interpolation functions,definition of weighting functions,and derivation of integration rules,are examined in detail,emphasizing their significance in comprehending the method’s effectiveness and stability.At last,the practical application of the collocation methods in engineering contexts is emphasized,including heat conduction simulations,electromagnetic coupled field analysis,and fluid dynamics simulations.These specific case studies can underscore collocation method’s broad applicability and effectiveness in addressing complex engineering challenges.In conclusion,this paper puts forward the future development trend of the collocation method through rigorous analysis and discussion,thereby facilitating further advancements in research and practical applications within these fields.

Application of Two-dimensional Viscous CE/SE Method in Calculation of Two-phase Detonation

The method of two-dimensional viscous space-time conservation element and solution element (CE/SE) can be used to calculate the gas-liquid two-phase interior flow field in pulse detonation engine (PDE). In this paper, the evolution of the detonation wave and the distribution of its physical parameters were analyzed. The numerical results show that the change of axial velocity of gas is the same as that of detonation pressure. The larger the liquid droplet radius is, the longer the time to get stable detonation wave is. The calculated results coincide with the experimented results better.