Acoustic viscoelastic modeling by frequency-domain boundary element method

Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green＇s function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.

Wind Structure in an Intermediate Boundary Layer Model Based on Ekman Momentum Approximation

A quasi three–dimensional, intermediate planetary boundary layer (PBL) model is developed by including inertial acceleration with the Ekman momentum approximation, but a nonlinear eddy viscosity based on Blackadar’s scheme was included to improve the theoretical model proposed by Tan and Wu (1993). The model could keep the same complexity as the classical Ekman model in numerical, but extends the conventional Ekman model to include the horizontal accelerated flow with the Ekman momentum approximation. A comparison between this modified Ekman model and other simplified accelerating PBL models is made. Results show that the Ekman model overestimates (underestimates) the wind speed and pumping velocity in the cyclonic (anticyclonic) shear flow due to the neglect of the acceleration flow, however, the semi–geostrophic Ekman model overestimates the acceleration effects resulting from the underestimating (overestimating) of the wind speed and pumping velocity in the cyclonic (anticyclonic) shear flow. The Ekman momentum approximation boundary layer model could be applied to the baroclinic atmosphere. The baroclinic Ekman momentum approximation boundary layer solution has both features of classical baroclinic Ekman layer and the Ekman momentum approximate boundary lager.

A numerical study for boundary layer current and sheet flow transport induced by a skewed asymmetric wave

An analytical model with essential parameters given by a two-phase numerical model is utilized to study the net boundary layer current and sediment transport under skewed asymmetric oscillatory sheet flows. The analytical model is the first instantaneous type model that can consider phase-lag and asymmetric boundary layer development. The two-phase model supplies the essential phase-lead, instantaneous erosion depth and boundary layer development for the analytical model to enhance the understanding of velocity skewness and acceleration skewness in sediment flux and transport rate. The sediment transport difference between onshore and offshore stages caused by velocity skewness or acceleration skewness is shown to illustrate the determination of net sediment transport by the analytical model. In previous studies about sediment transport in skewed asymmetric sheet flows, the generation of net sediment transport is mainly concluded to the phase-lag effect.However, the phase-lag effect is shown important but not enough for the net sediment transport, while the skewed asymmetric boundary layer development generated net boundary layer current and mobile bed effect are key important in the transport process.

Fifth-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems (5th-CASAM-N): II. Paradigm Application to a Bernoulli Model Comprising Uncertain Parameters

This work presents the application of the recently developed “Fifth-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems (5th-CASAM-N)” to a simplified Bernoulli model. The 5th-CASAM-N builds upon and incorporates all of the lower-order (i.e., the first-, second-, third-, and fourth-order) adjoint sensitivities analysis methodologies. The Bernoulli model comprises a nonlinear model response, uncertain model parameters, uncertain model domain boundaries and uncertain model boundary conditions, admitting closed-form explicit expressions for the response sensitivities of all orders. Illustrating the specific mechanisms and advantages of applying the 5th-CASAM-N for the computation of the response sensitivities with respect to the uncertain parameters and boundaries reveals that the 5th-CASAM-N provides a fundamental step towards overcoming the curse of dimensionality in sensitivity and uncertainty analysis.

Electroosmotic flow of Eyring fluid in slit microchannel with slip boundary condition

In consideration of the electroosmotic flow in a slit microchannel, the con-stitutive relationship of the Eyring fluid model is utilized. Navier＇s slip condition is used as the boundary condition. The governing equations are solved analytically, yielding the velocity distribution. The approximate expressions of the velocity distribution are also given and discussed. Furthermore, the effects of the dimensionless parameters, the electrokinetic parameter, and the slip length on the flow are studied numerically, and appropriate conclusions are drawn.

An enhanced treatment of boundary conditions in implicit smoothed particle hydrodynamics

In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme （Implicit SPH）. Three types of virtual particles, i.e., boundary particles, image particles and mirror particles, are used to impose boundary conditions. Boundary particles are fixed on the solid boundary, and each boundary particle is associated with two fixed image particles inside the fluid domain and two fixed mirror particles outside the fluid domain. The image particles take the flow properties through fluid particles with moving least squares （MLS） interpolation and the properties of mirror particles can be obtained by the corresponding image particles. A repulsive force is also applied for boundary particles to prevent fluid particles from unphysical penetra- tion through solid boundaries. The new boundary treatment method has been validated with five numerical examples. All the numerical results show that Implicit SPH with this new boundary-treatment method can obtain accurate results for non-Newtonian fluids as well as Newtonian fluids, and this method is suitable for complex solid boundaries and can be easily extended to 3D problems.

Slicing and support structure generation for 3D printing directly on B-rep models

Traditional 3D printing is based on stereolithography or standard tessellation language models,which contain many redundant data and have low precision.This paper proposes a slicing and support structure generation algorithm for 3D printing directly on boundary representation(B-rep)models.First,surface slicing is performed by efficiently computing the intersection curves between the faces of the B-rep models and each slicing plane.Then,the normals of the B-rep models are used to detect where the support structures should be located and the support structures are generated.Experimental results show the efficiency and stability of our algorithm.

Building Morphological Characteristics and Their Effect on the Wind in Beijing

An urban boundary layer model （UBLM） is improved by incorporating the effect of buildings with a sectional drag coefficient and a height-distributed canopy drag length scale. The improved UBLM is applied to simulate the wind fields over three typical urban blocks over the Beijing area with different height-towidth ratios. For comparisons, the wind fields over the same blocks are simulated by an urban sub-domain scale model resolving the buildings explicitly. The wind fields simulated from the two different methods are in good agreement. Then, two-dimensional building morphological characteristics and urban canopy parameters for Beijing are derived from detailed building height data. Finally, experiements are conducted to investigate the effect of buildings on the wind field in Beijing using the improved UBLM.

A boundary surrogate model for micro/nano grooved surface structure applied in turbulence flow control over airfoil

The application of grooved surface structure is an emerging and effective means in turbulence flow control.However,for a realistic configuration,the global flow field described directly by simple application of massive grids makes it unfeasible to simulate.In this paper,a boundary surrogate model reproducing the effect of microscopic near-wall region is proposed to improve computational efficiency.The surrogate model trained with Lattice Boltzmann Method(LBM)considering the rarefied effect based on real micro/nanoflow field is new among literature,which accurately shows flow characteristics of the micro/nano structure.With this approach,numerical simulations via Reynolds-averaged Navier Stokes equations with modified wall boundary condition are performed in subsonic and transonic flow.The results show that micro/nano grooved surface structure has the effect of delaying transition from laminar to turbulence,thus reducing the skin friction significantly.Analysis of turbulence intensity and turbulence kinetic energy shows that the near-wall flow field of grooved airfoil is more stable compared with that of the smooth airfoil.The reducing rate of maximum turbulent intensity reaches 13.39%.The paper shows a perspective for further application of micro/nano groove structure to turbulence flow control in aircraft design by providing an accurate and efficient simulation method.

Improved Parallel Processing Function for High-Performance Large-Scale Astronomical Cross-Matching

Astronomical cross-matching is a basic method for aggregating the observational data of different wavelengths. By data aggregation, the properties of astronomical objects can be understood comprehensively. Aiming at decreasing the time consumed on I/O operations, several improved methods are introduced, including a processing flow based on the boundary growing model, which can reduce the database query operations; a concept of the biggest growing block and its determination which can improve the performance of task partition and resolve data-sparse problem; and a fast bitwise algorithm to compute the index numbers of the neighboring blocks, which is a significant efficiency guarantee. Experiments show that the methods can effectively speed up cross-matching on both sparse datasets and high-density datasets.

A Remark on Global Existence in H^4 for a Nonlinear Thermoviscoelasticity Equations with Non-monotone Pressure

We consider a one-dimensional continuous model of nutron star, described by a compressible thermoviscoelastic system with a non-monotone equation of state, due to the effective Skyrme nuclear interaction between particles. We will prove that, despite a possible destabilizing influence of the pressure, which is non-monotone and not always positive, the presence of viscosity and a sufficient thermal dissipation describe the global existence of solutions in H^4 with a mixed free boundary problem for our model.

Global Existence in H^2 for a 1-d Non-monotone Fluid

In this paper, we prove the existence in H＋^2, an incomplete metric subspace of H^2×H^2×H^2, of global solutions to the system for a one-dimensional non-monotone fluid in bounded domainΩ=（0,1）. The results in this paper have improved those previously related results.

NUMERICAL SIMULATION OF SOUND RADIATION ON GEARBOX'S HOUSING

The fluid-structure coupling finite element model and the boundary element model of a complex gearbox＇s housing are built based on the theory of fluid-structure coupling finite element method and boundary element method. At the same time, the exciting forces of the housing are analyzed and applied to the finite element models. Firstly, vibration of the housing is calculated by the fluid-structure coupling finite element model; secondly, the calculated result is verified by the experiment; finally, sound radiation of the housing is calculated by the boundary element. According to the calculated results, the housing adds some ribs not only to increase the strength, but also to reduce the sound radiation of the housing. At last, the sound radiation of the modified housing is calculated, which shows that the sound radiation of the modified housing with ribs is lower.

Second-Order Adjoint Sensitivity Analysis Methodology for Computing Exactly Response Sensitivities to Uncertain Parameters and Boundaries of Linear Systems: Mathematical Framework

This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2nd-CASAM)” for the efficient and exact computation of 1st- and 2nd-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (i.e., model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2nd-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1st-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2nd-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (N2/2 + 3 N/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2nd-LASS requires very little additional effort beyond the construction of the 1st-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1st-LASS and 2nd-LASS requires the same computational solvers as needed for solving (i.e., “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1st-LASS and the 2nd-LASS. Since neither the 1st-LASS nor the 2nd-LASS involves any differentials of the operators underlying the original system, the 1st-LASS is designated as a “first-level” (as opposed to a “first-order”) adjoint sensitivity system, while the 2nd-LASS is designated as a “second-level” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2nd-LASS that involve the imprecisely known boundary parameters. Notably, the 2nd-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.

A boundary layer model for capture of inclusions by steel-slag interface in a turbulent flow

A boundary layer model was developed to predict the capture of inclusions by steel-slag interface in a turbulent fluid flow,which is based on the detailed analysis of inclusion trajectories.The effective boundary layer for inclusion removal was proposed by a statistical method.It is noticed that the capture of inclusions by steel-slag interface is not only dependent on the diameter of inclusions but also related to the local turbulent conditions.In high turbulent flow fields,the transport of inclusions is mainly dominated by the turbulent flow,and thus,the effective boundary layer thickness is mainly affected by the level of turbulent kinetic energy and is almost independent of the inclusion diameter.The inertia of inclusions gradually takes over the stochastic effect of turbulent flow,and the effect of inclusion diameter on effective boundary layer thickness becomes more noticeable with the decrease in the level of turbulent kinetic energy.Besides,the effective boundary layer thickness is more susceptible to the inclusion diameter for larger inclusions due to its greater inertia under the same turbulent condition while it principally depends on the level of turbulent kinetic energy for smaller inclusions.As the characteristic velocity increases,the time for inclusions transport and interaction with steel-slag interface decreases,and thus,the effective boundary layer thickness decreases.Moreover,the graphical user interface was developed by using the cubic spline interpolation for ease of coupling the current boundary layer model with the macro-scale model of a turbulent fluid flow in the metallurgical vessel.

Quantifying a critical marl thickness for vertical fracture extension using field data and numerical experiments

In fractured reservoirs characterized by low matrix permeability,fracture networks control the main fluid flow paths.However,in layered reservoirs,the vertical extension of fractures is often restricted to single layers.In this study,we explored the effect of changing marl/shale thickness on fracture extension using comprehensive field data and numerical modeling.The field data were sampled from coastal exposures of Liassic limestone-marl/shale alternations in Wales and Somerset(Bristol Channel Basin,UK).The vertical fracture traces of more than 4000 fractures were mapped in detail.Six sections were selected to represent a variety of layer thicknesses.Besides the field data also thin sections were analyzed.Numerical models of fracture extension in a two-layer limestone-marl system were based on field data and laboratory measurements of Young's moduli.The modeled principal stress magnitude σ3 along the lithological contact was used as an indication for fracture extension through marls.Field data exhibit good correlation(R^2=0.76) between fracture extension and marl thickness,the thicker the marl layer the fewer fractures propagate through.The model results show that almost no tensile stress reaches the top of the marl layer when the marls are thicker than 30 cm.For marls that are less than 20 cm,the propagation of stress is more dependent on the stiffness of the marls.The higher the contrast between limestone and marl stiffness the lower the stress that is transmitted into the marl layer.In both model experiments and field data the critical marl thickness for fracture extension is ca.15-20 cm.This quantification of critical marl thicknesses can be used to improve predictions of fracture networks and permeability in layered rocks.Up-or downsampling methods often ignore spatially continuous impermeable layers with thicknesses that are under the detection limit of seismic data.However,ignoring these layers can lead to overestimates of the overall permeability.Therefore,the understanding of how fractures propagate and terminate through impermeable layers will help to improve the characterization of conventional reservoirs.

A Peridynamic Approach for the Evaluation of Metal Ablation under High Temperature

In this paper,the evaluations of metal ablation processes under high temperature,i.e.,the Al plate ablated by a laser and a heat carrier and the reactor pressure vessel ablated by a core melt,are studied by a novel peridynamic method.Above all,the peridynamic formulation for the heat conduction problem is obtained by Taylor’s expansion technique.Then,a simple and efficient moving boundary model in the peridynamic framework is proposed to handle the variable geometries,in which the ablated states of material points are described by an additional scalar field.Next,due to the automatic non-interpenetration properties of peridynamic method,a contact algorithm is established to determine the contact relationship between the ablated system and the additional heat carrier.In addition,the corresponding computational procedure is listed in detail.Finally,several numerical examples are carried out and the results verify the validity and accuracy of the present method.

Growth Ring-dependent Fracture Toughness of Sea Urchin Spines Estimated by Boundary Effect Model

Although the fracture behavior of sea urchin spines has been extensively investigated,there is as yet a lack of quantitative estimation on the effect of growth rings on the fracture properties of sea urchin spines.In sea urchin spines,much denser pores present in growth rings rather than porous layers.The tensile strength and fracture toughness of sea urchin spine samples with different numbers of growth rings are measured by the Boundary Effect Model(BEM).The experimental results of single-edge notched three-point bending tests indicate that the BEM is an appropriate method to estimate the fracture toughness of the present porous sea urchin spines,and the number of growth rings plays an important role in the fracture properties of spines.Specifically,the tensile strength and fracture toughness of sea urchin spines can be significantly improved with the increase in the number of growth rings,and their fracture toughness can even reach a relatively high value compared with some other porous materials with an identical porosity.The present research findings are expected to provide a fundamental insight into the design of high-performance bionic materials with a highly porous structure.

RECOVERY OF URANIURN FROM CARBONATE SOLUTIONS USING STRONGLY BASIC ANION EXCHANGER 3. THE MECHANISMS OF RECOVERY PROCESSES

A moving boundary model under considering the volume change of spherical resin beads during ion exchange processes was employed to recognize the mechanisms of recovering uranium from carbonate solutions using strongly basic anion exchanger Two important factors, swelling and ion exchange, which directly affect the volume of ion exchangers were taken into account. An ion exchange mechanism has been found for the forward reaction RCl/[UO2 CO3)3]4- at pH 8.5～10.0, where the rate is controlled by liquid-film diffusion at low concentration of the tetravalent uranyl complex anion [UO2 (CO3)3]4-, and is partical diffusion governing at high concentration of the complex anion. The mechanism of RCl/U(Ⅳ) at pH 5. 5～7. 5 is a chemical reaction taking place at the moving boundary of the unreacted nucleus. For the reverse reaction RnU/NaCl, the uranyl tricarbonate complex anion in the resin phase is replaced by Cl- ions with an ion exchange mechanism always determined by particle diffusion. The other forms of uranium in the solid phase loaded on the resin at pH 5. 5～7. 5 should belong to non-exchangeable uranium. The mechanism of the reverse reaction RnU/HCl is always chemical reaction which is not restricted to the moving boundary of the unreacted core.

Dynamic Control of Defective Gap Mode Through Defect Location

A one dimensional model is developed for defective gap mode（DGM）with two types of boundary conditions：conducting mesh and conducting sleeve.For a periodically modulated system without defect,the normalized width of spectral gaps equals to the modulation factor,which is consistent with previous studies.For a periodic system with local defects introduced by the boundary conditions,it shows that the conducting-mesh-induced DGM is always well confined by spectral gaps while the conducting-sleeve-induced DGM is not.The defect location can be a useful tool to dynamically control the frequency and spatial periodicity of DGM inside spectral gaps.This controllability can be potentially applied to the interaction between gap eigenmodes and energetic particles in fusion plasmas,and optical microcavities and waveguides in photonic crystals.