Modeling and simulation of solvent behavior and temperature distribution within long stick propellants with large web thickness undergoing drying

Drying is a complicated physical process which involves simultaneous heat and mass transfer in the removal of solvents inside propellants.Inappropriate drying techniques may result in the formation of a hard skin layer near the surface to block the free access of most solvent through for long stick propellants with large web thickness,which lead to lower drying efficiency and worse drying quality.This study aims to gain a comprehensive understanding of drying process and clarify the mechanism of the blocked layer near the propellant surface.A new three-dimensional coupled heat and mass transfer(3D-CHMT)model was successfully developed under transient conditions.The drying experiment results show that the 3DCHMT model could be applied to describe the drying process well since the relative error of the content of solvent between simulation and experiment values is only 5.5%.The solvent behavior simulation demonstrates that the mass transfer process can be divided into super-fast(SF)and subsequent minorfast(MF)stages,and the SF stage is vital to the prevention of the blocked layer against the free access for solvent molecules inside propellant grains.The effective solvent diffusion coefficient(Deff)of the propellant surface initially increases from 3.4×10^(-6)to 5.3×10^(-6)m^(2)/s as the temperature increases,and then decreases to 4.1×10^(-8)m^(2)/s at 60-100 min.The value of Deffof surface between 0-1.4 mm has a unique trend of change compared with other regions,and it is much lower than that of the internal at100 min under simulation conditions.Meanwhile,the temperature of the propellant surface increases rapidly at the SF stage(0-100 min)and then very slowly thereafter.Both the evolution of Deffand temperature distribution demonstrate that the blocked layer near the propellant surface has been formed in the time period of approximately 0-100 min and its thickness is about 1.4 mm.To mitigate the formation of blocked layer and improve its drying quality of finial propellant products effectively,it should be initially dried at lower drying temperature(30-40℃)in 0-100 min and then dried at higher drying temperature(50-60℃)to reduce drying time for later drying process in double base gun propellants.The present results can provide theoretical guidance for drying process and optimization of drying parameters for long stick propellants with large web thickness.

Quasi-3D Numerical Simulation of Tidal Hydrodynamic Field

Based on the 2D horizontal plane numerical model, a quasi-3D numerical model is established for coastal regions of shallow water. The characteristics of this model are that the velocity profiles;can be obtained at the same time when the equations of the value of difference between the horizontal current velocity and its depth-averaged velocity in the vertical direction are solved and the results obtained are consistent with the results of the 2D, model. The circulating flow in the rectangular area induced by wind is simulated and applied to the tidal flow field of the radial sandbanks in the South Yellow Sea. The computational results from this quasi-3D model are in good agreement with analytical results and observed data. The solution of the finite difference equations has been found to be stable, and the model is simple, effective and practical.

A three-dimensional numerical study on the effect of geometric asymmetry on arcjet thruster performance

In an arcjet thruster,the cathode and constrictor degrade with time,and the electrical arc discharge may become unsymmetrical.In this work,a three-dimensional numerical model of a hydrogen plasma arcjet is developed and validated to study the effect of unsymmetrical electric arc discharge on thruster performance.The unsymmetrical arc discharge is realized by introducing a radial shift of the cathode so that the cathode tip offset is 80μm(25%of the constrictor radius).Simulations are conducted for both axially centered cathode(coaxial)and off-centered cathode(non-coaxial)configurations with identical propellant flow rates and input current.Simulations show asymmetrical arc discharge in the non-coaxial cathode configuration,resulting in azimuthally asymmetric Joule heating,species concentrations,and velocity field.This asymmetry continues as the plasma expands in the divergent section of the nozzle.Temperature,species concentrations,and axial velocity exhibit asymmetric radial distribution at the nozzle exit.The computed Joule heating was found to reduce with cathode shift,and consequently,the thrust and specific impulse of the thruster was decreased by about 6.6%.In the case of the non-coaxial cathode,geometric asymmetry also induces a small side thrust.

Meshing effects of the 3-D FEM numerical modeling in seismo-electromagnetics:An application in selectivity of seismic electric signal (SES)

We investigated how density and quality of mesh around interest domain affect electromagnetic （EM） responses of 3D Earth layered media using finite element method （FEM）. Effect of different mesh shapes was also investigated using a method of mixing structured and unstructured mesh. As a case study, we estimated the effects of meshing on selectivity phenomenon of seismic electric signal （SES）. Our results suggest that the relative errors resulting from mesh effects may not be negligible, which may lead to some unconvincing explanation of the SES selectivity based on the numerical modeling results.

Runout prediction of potential landslides based on the multi-source data collaboration analysis on historical cases

Long runout landslides involve a massive amount of energy and can be extremely hazardous owing to their long movement distance,high mobility and strong destructive power.Numerical methods have been widely used to predict the landslide runout but a fundamental problem remained is how to determine the reliable numerical parameters.This study proposes a framework to predict the runout of potential landslides through multi-source data collaboration and numerical analysis of historical landslide events.Specifically,for the historical landslide cases,the landslide-induced seismic signal,geophysical surveys,and possible in-situ drone/phone videos(multi-source data collaboration)can validate the numerical results in terms of landslide dynamics and deposit features and help calibrate the numerical(rheological)parameters.Subsequently,the calibrated numerical parameters can be used to numerically predict the runout of potential landslides in the region with a similar geological setting to the recorded events.Application of the runout prediction approach to the 2020 Jiashanying landslide in Guizhou,China gives reasonable results in comparison to the field observations.The numerical parameters are determined from the multi-source data collaboration analysis of a historical case in the region(2019 Shuicheng landslide).The proposed framework for landslide runout prediction can be of great utility for landslide risk assessment and disaster reduction in mountainous regions worldwide.

Benchmarking of two three-dimensional numerical models in time/space domain to predict railway-induced ground vibrations

In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.

3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing

A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Approach of moving boundary and its application in 3D tidal current simulation of Changjiang (Yangtze) River Estuary based on ECOMSED model

Considering three-dimensional model ECOMSED can not simulate wetting-drying of shoal with its fixed boundary, an approach to represent moving boundary in the model is introduced here. This approach smoothly joints the internal and external mode by making use of wetting and drying technique and is verified by a numerical test which presents a good agreement with the previous test results obtained by other researchers. A three dimensional numerical model is established to simulate the hydrodynamics in spring tide in the Changjiang （Yangtze） River estuary by this modified ECOMSED model which is also validated through the observed field data, the simulation presents a good periodic tidal change. It also successfully simulates the tidal current of computational areas and reproduces the tidal flat intermittent appearance.

Numerical Study on Instantaneous Discharge of Unsorted Particle Cloud in Cross Flow

The mixing characteristics of particles such as dredged sediment of variable size discharged into cross flow are studied by a 3D numerical modal, which is developed to model the particle-fluid two-phase flow. The Eulerian method with the modified k- ε parameterization of turbulence for the fluid phase is used to solve fluid phase, while a Lagrangian method for the sohd phase （particles）, both the processes are coupled through the momentmn sources. In the model the wake turbulence induced by particles has been included as additional source term in the k - ε model; and the variable drift velocities of the particles are treated efficiently by the Lagrangian method in which the particles are tracked explicitly and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumping a cloud of particles is governed by the total buoyancy of the cloud, the drag force on each particle and the velocity of cross-flow. The computed results show a roughly linear relationship between the displacement of the frontal position and the longitudi- nal width of the particle cloud. The particle size in the cloud and the velocity of cross flow dominate the flow behavior. The computed results are compared with the results of laboratory experiments and satisfactory agreement is obtained.

A Higher-Efficient Three-Dimensional Numerical Model for Small Amplitude Free Surface Flows

A higher-efficient three-dimensional non-hydrostatic model is developed to simulate small amplitude free surface flows based on a staggered unstructured grid. In this model, a fractional step algorithm is adopted to solve the Navier-Stokes equations in two major steps. A top-layer pressure method is proposed to minimize the number of vertical layers and subsequently the computational cost. Three classical examples of small amplitude free surface flows are used to demonstrate the capability and efficiency of the model. The satisfactory results demonstrated the capability and efficiency of modelling a range of small amplitude free surface flows with only a small number of vertical layers.

A model study of influence of circulation on the pollutant transport in the Zhujiang River Estuary and adjacent coastal waters

A tracer model with random diffusion coupled to the hydrodynamic model for the Zhujiang River Estuary (Pearl River Estuary, PRE) is to examine the effect of circulations on the transport of completely conservative pollutants. It is focused on answering the following questions: (1) What role does the estuarine plume front in the winter play in affecting the pollutants transport and its distribution in the PRE ? (2) What effect do the coastal currents driven by the monsoon have on the pollutants transport? The tracer experiment results show that: (1) the pollutant transport paths strongly depend on the circulation structures and plume frontal dynamics of the PRE and coastal waters; (2) during the summer when a southwesterly monsoon prevails, the pollutants from the four easterly river inlets and those from the bottom layer of offshore stations will greatly influence the water quality in Hong Kong waters, however, the pollutants released from the four westerly river-inlets will seldom affect the water qual

Contrasted continental rifting via plume-craton interaction： Applications to Central East African Rift

The East African Rift system （EARS） provides a unique system with the juxtaposition of two contrasting yet simultaneously formed rift branches, the eastern, magma-rich, and the western, magma-poor, on either sides of the old thick Tanzanian craton embedded in a younger lithosphere. Data on the pre-rifr, syn-rift and post-rift far-field volcanic and tectonic activity show that the EARS formed in the context of the interaction between a deep mantle plume and a horizontally and vertically heterogeneous lithosphere under far-field tectonic extension. We bring quantitative insights into this evolution by implementing high-resolution 3D thermo-mechanical numerical deformation models of a lithosphere of realistic rheology. The models focus on the central part of the EARS. We explore scenarios of plumelithosphere interaction with plumes of various size and initial position rising beneath a tectonically pre-stretched lithosphere. We test the impact of the inherited rheological discontinuities （suture zones） along the craton borders, of the rheological structure, of lithosphere plate thickness variations, and of physical and mechanical contrasts between the craton and the embedding lithosphere. Our experiments indicate that the ascending plume material is deflected by the cratonic keel and preferentially channeled along one of its sides, leading to the formation of a large rift zone along the eastern side of the craton, with significant magmatic activity and substantial melt amount derived from the mantle plume material. We show that the observed asymmetry of the central EARS, with coeval amagmatic （western） and magmatic （eastern） branches, can be explained by the splitting of warm material rising from a broad plume head whose initial position is slightly shifted to the eastern side of the craton. In that case, neither a mechanical weakness of the contact between the craton and the embedding lithosphere nor the presence of second plume are required to produce simulations that match observations. This result reconciles the passive and active rift models and demonstrates the possibility of development of both magmatic and amagmatic rifts in identical geotectonic environments.

Assessment of an alternative to deep foundations in compressible clays： the structural cell foundation

The new type of deep foundation for buildings on saturated, compressible-low strength clayey soil deposits, branded structural cell essentially consists of a rigid concrete top slab, structurally connected to reinforced concrete peripheral walls （diaphragms） that enclose the natural soil. Accordingly, as the initial volume of the confined soft clays within the lateral stiff diaphragms will remain constant upon loading, the hollowed structural cell will be ＂transformed＂ into a very large cross-section pillar of unit weight slightly higher than that of the natural soft clayey soil. This type of foundation seems to be a highly competitive alternative to the friction pile-box foundations （widely used in Mexico City clays）, due to its economic and environmental advantages. Economies result, for example, from the absence of huge excavations hence sparing the need of earth retaining structures. Further savings result from appreciably smaller concrete volumes required for building the structural cell than the friction pile-box foundation; moreover, the construction time of the former is much shorter than that of the latter. Regarding the impact to the environment, less air contamination follows from the fact that both traffic jams and soil excavation lessen appreciably. Considering these facts and others regarding scheduling, it was decided to replace 48-friction pile-box foundations specified in the master plan project by this new type of foundation. The overall behavior of these cell foundations over a five-year period is fared from close visual observations and their leveling during the first three years after their construction.