Numerical simulation study on the mold strength of magnetic mold casting based on a coupled electromagnetic-structural method

The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.

A coupled Legendre-Laguerre polynomial method with analytical integration for the Rayleigh waves in a quasicrystal layered half-space with an imperfect interface

The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.

Application of the N + 2 Transversal Network Method to the Study of a Coupled Resonator Filter

This paper presents a new approach to synthesize admittance function polynomials and coupling matrices for coupled resonator filters. The N + 2 transversal network method is applied to study a coupled resonator filter. This method allowed us to determine the polynomials of the reflection and transmission coefficients. A study is made for a 4 poles filter with 2 transmission zeros between the N + 2 transversal network method and the one found in the literature. A MATLAB code was designed for the numerical simulation of these coefficients for the 6, 8, and 10 pole filter with 4 transmission zeros.

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.

Numerical simulation of aluminum holding furnace with fluid-solid coupled heat transfer

To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.

Tolerance of edge cascades with coupled map lattices methods

This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating some indices including the number of failed edges, dynamic edge tolerance capacity and the perturbation of edge. In general, it assumes that the perturbation is attributed to the normal distribution in adopted simulations. By investigating the effectiveness of edge tolerance in scale-free and random networks, it finds that the larger tolerance parameter λ can more efficiently delay the cascading failure process for scale-free networks than random networks. These results indicate that the cascading failure process can be effectively controlled by increasing the tolerance parameter λ. Moreover, the simulations also show that, larger variance of perturbation can easily trigger the cascading failures than the smaller one. This study may be useful for evaluating efficiency of whole traffic systems, and for alleviating cascading failure in such systems.

Thermodynamic analysis and simulation for gas baffle entrance collimator of EAST-NBI system based on thermo-fluid coupled method

The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the beam channel for absorbing the divergence beam during the beam transmission process in the EAST-NBI system.The gas baffle entrance collimator(GBEC) is a typical high-heat-flux component located at the entrance of gas baffle. An efficient and accurate analysis of its thermodynamic performance is of great significance to explore the working limit and to ensure safe operation of the system under a high-parameter steady-state condition. Based on the thermo-fluid coupled method, thermodynamic analysis and simulation of GBEC is performed to get the working states and corresponding operating limits at different beam extraction conditions. This study provides a theoretical guidance for the next step to achieve long pulse with highpower experimental operation and has an important reference to ensure the safe operation of the system.

Transfer matrix method for determination of the natural vibration characteristics of elastically coupled launch vehicle boosters

The analysis of natural vibration characteristics has become one of important steps of the manufacture and dynamic design in the aerospace industry. This paper presents a new scenario called virtual cutting in the context of the transfer matrix method of linear multibody systems closed- loop topology for computing the free vibration characteristics of elastically coupled flexible launch vehicle boosters. In this approach, the coupled system is idealized as a triple-beam system-like structure coupled by linear translational springs, where a non-uniform free-free Euler-Bemoulli beam is used. A large thrust-to-weight ratio leads to large axial accelera- tions that result in an axial inertia load distribution from nose to tail. Consequently, it causes the development of significant compressive forces along the length of the launch vehicle. Therefore, it is important to take into account this effect in the transverse vibration model. This scenario does not need the global dynamics equations of a system, and it has high computational efficiency and low memory requirements. The validity of the presented scenario is achieved through com- parison to other approaches published in the literature.

Average vector field methods for the coupled Schrdinger KdV equations

The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.

Optimal control of attitude for coupled-rigid-body spacecraft via Chebyshev-Gauss pseudospectral method

The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a dis- crete nonlinear programming （NLP） problem using the Chebyshev-Gauss pseudospectral method （CGPM）. Solutions of the NLP problem can be obtained using the sequential quadratic programming （SQP） algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss （CG） points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform （FFT）. To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and con- trol variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.

Dynamic analysis of beam-cable coupled systems using Chebyshev spectral element method

The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a double Timoshenko beam system interconnected by discrete springs. Utilizing Chebyshev series expansion and meshing the system according to the locations of its connections, numerical results of the natural frequencies and mode shapes are obtained using only a few elements, and the results are validated by comparing them with the results of a finite-element method. Then the effects of the cable parameters and layout of connections on the natural frequencies and mode shapes of a fixed-pinned beam are studied. The results show that the modes of a beam-cable coupled system can be classified into two types, beam mode and cable mode, according to the dominant deformation. To avoid undesirable vibrations of the cable, its parameters should be controlled in a reasonable range, or the layout of the connections should be optimized.

MPS-FEM Coupled Method for Study of Wave-Structure Interaction

Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.

The determination of 52 elements in marine geological samples by an inductively coupled plasma optical emission spectrometry and an inductively coupled plasma mass spectrometry with a high-pressure closed digestion method

An improved analytical method to determine the content of 52 major, minor and trace elements in marine geological samples, using a HF-HCl-HNO_3 acid system with a high-pressure closed digestion method（HPCD）, is studied by an inductively coupled plasma optical emission spectrometry（ICP-OES） and an inductively coupled plasma mass spectrometry（ICP-MS）. The operating parameters of the instruments are optimized, and the optimal analytical parameters are determined. The influences of optical spectrum and mass spectrum interferences, digestion methods and acid systems on the analytical results are investigated. The optimal spectral lines and isotopes are chosen, and internal standard element of rhodium is selected to compensate for matrix effects and analytical signals drifting. Compared with the methods of an electric heating plate digestion and a microwave digestion, a high-pressure closed digestion method is optimized with less acid, complete digestion,less damage for digestion process. The marine geological samples are dissolved completely by a HF-HCl-HNO_3 system, the relative error（RE） for the analytical results are all less than 6.0%. The method detection limits are 2–40μg/g by the ICP-OES, and 6–80 ng/g by ICP-MS. The methods are used to determine the marine sediment reference materials（GBW07309, GBW07311, GBW07313）, rock reference materials（GBW07103, GBW07104,GBW07105）, and cobalt-rich crust reference materials（GBW07337, GBW07338, GBW07339）, the obtained analytical results are in agreement with the certified values, and both of the relative standard deviation（RSD） and the relative error（RE） are less than 6.0%. The analytical method meets the requirements for determining 52 elements contents of bulk marine geological samples.

Local discontinuous Galerkin method for solving Burgers and coupled Burgers equations

In the current work, we extend the local discontinuous Galerkin method to a more general application system. The Burgers and coupled Burgers equations are solved by the local discontinuous Galerkin method. Numerical experiments are given to verify the efficiency and accuracy of our method. Moreover the numerical results show that the method can approximate sharp fronts accurately with minimal oscillation.

Simulation Analysis of Torsion Beam Hydroforming Based on the Fluid-Solid Coupling Method

Hydroformed parts are widely used in industrial automotive parts because of their higher stiffness and fatigue strength and reduced weight relative to their corresponding cast and welded parts.This paper reports a hydraulicforming experimental platform for rectangular tube fittings that was constructed to conduct an experiment on the hydraulic forming of rectangular tube fittings.A finite element model was established on the basis of the fluid–solid coupling method and simulation analysis.The correctness of the simulation analysis and the feasibility of the fluid–solid coupling method for hydraulic forming simulation analysis were verified by comparing the experimental results with the simulation results.On the basis of the simulation analysis of the hydraulic process of the torsion beam using the fluid–solid coupling method,a sliding mold suitable for the hydroforming of torsion beams was designed for its structural characteristics.The effects of fluid characteristics,shaping pressure,axial feed rate,and friction coefficient on the wall thicknesses of torsions beams during formation were investigated.Fluid movement speed was related to tube deformation.Shaping pressure had a significant effect on rounded corners and straight edges.The axial feed speed was increased,and the uneven distribution of wall thicknesses was effectively improved.Although the friction coefficient had a nonsignificant effect on the wall thickness of the ladder-shaped region,it had a significant influence on a large deformation of wall thickness in the V-shaped area.In this paper,a method of fluid-solid coupling simulation analysis and sliding die is proposed to study the high pressure forming law in torsion beam.

Applications of F-expansion method to the coupled KdV system

An extended F-expansion method for finding periodic wave solutions of nonlinear evolution equations in mathematical physics is presented, which can be thought of as a concentration of extended Jacobi elliptic function expansion method proposed more recently. By using the homogeneous balance principle and the extended F-expansion, more periodic wave solutions expressed by Jacobi elliptic functions for the coupled KdV equations are derived. In the limit cases, the solitary wave solutions and the other type of travelling wave solutions for the system are also obtained.

Dynamic Coupling Analysis of Semisubmersible Platform Float-over Method for Docking Case

In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible platform was performed.Relevant hydrodynamic parameters were obtained using the retardation function method of three-dimensional frequency-domain potential flow theory.The results of the hydrodynamic analysis were highly consistent with the test findings,verifying the accuracy of the multifloating hydrodynamic coupling analysis,and key hydrodynamic parameters were solved for different water depths and the coupling effect.According to the obtained results,the hydrodynamic influence was the largest in shallow waters when the coupling effect was considered.Furthermore,the coupled motion equation combined with viscous damping,fender system,and mooring system was established,and the hydrodynamics,floating body motion,and dynamic response of the fender system were analyzed.Motion analysis revealed good agreement among the surge,sway,and yaw motions of the two floating bodies.However,when the wave period reached 10 s,the motion of the two floating bodies showed severe shock,and a relative motion was also observed.Therefore,excessive constraints should be added between the two floating bodies during construction to ensure construction safety.The numerical analysis and model test results of the semisubmersible platform and HYSY 229 barge at a water depth of 42 m and sea conditions of 0°,45°,and 90° were in good agreement,and the error was less than 5%.The maximum movement of the HYSY 229 barge reached 2.61 m in the sway direction,whereas that of the semisubmersible platform was 2.11 m.During construction,excessive constraints should be added between the two floating bodies to limit their relative movement and ensure construction safety.

The Hirota bilinear method for the coupled Burgers equation and the high-order Boussinesq-Burgers equation

This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink （soliton） solutions and multiple-singular-kink （soliton） solutions are derived for the two equations.

Numerical manifold method modeling of coupled processes in fractured geological media at multiple scales

The greatest challenges of rigorously modeling coupled hydro-mechanical(HM)processes in fractured geological media at different scales are associated with computational geometry.These challenges include dynamic shearing and opening of intersecting fractures at discrete fracture scales as a result of coupled processes,and contact alteration along rough fracture surfaces that triggers structural and physical changes of fractures at micro-asperity scale.In this paper,these challenges are tackled by developing a comprehensive modeling approach for coupled processes in fractured geological media based on numerical manifold method(NMM)at multiple scales.Based on their distinct geometric features,fractures are categorized into three different scales:dominant fracture,discrete fracture,and discontinuum asperity scales.Here the scale is relative,that of the fracture relative to that of the research interest or domain.Different geometric representations of fractures at different scales are used,and different governing equations and constitutive relationships are applied.For dominant fractures,a finite thickness zone model is developed to treat a fracture as a porous nonlinear domain.Nonlinear fracture mechanical behavior is accurately modeled with an implicit approach based on strain energy.For discrete fractures,a zero-dimensional model was developed for analyzing fluid flow and mechanics in fractures that are geometrically treated as boundaries of the rock matrix.With the zero-dimensional model,these fractures can be modeled with arbitrary orientations and intersections.They can be fluid conduits or seals,and can be open,bonded or sliding.For the discontinuum asperity scale,the geometry of rough fracture surfaces is explicitly represented and contacts involving dynamic alteration of contacts among asperities are rigorously calculated.Using this approach,fracture alteration caused by deformation,re-arrangement and sliding of rough surfaces can be captured.Our comprehensive model is able to handle the computational challenges with accurate representation of intersections and shearing of fractures at the discrete fracture scale and rigorously treats contacts along rough fracture surfaces at the discontinuum asperity scale.With future development of three-dimensional(3D)geometric representation of discrete fracture networks in porous rock and contacts among multi-body systems,this model is promising as a basis of 3D fully coupled analysis of fractures at multiple scales,for advancing understanding and optimizing energy recovery and storage in fractured geological media.

Numerical Simulation of the Seismic Response of a Horizontal Storage Tank Based on a SPH-FEM Coupling Method

A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a horizontal storage tank featuring a free liquid surface under seismic action was constructed using the SPH–FEM coupling method.The stored liquid was discretized using SPH particles,while the tank and supports were discretized using the FEM.The interaction between the stored liquid and the tank was simulated by using the meshless particle contact method.Then,the numerical simulation results were compared and analyzed against seismic simulation shaking table test data to validate the method.Subsequently,a series of numerical models,considering different liquid storage volumes and seismic effects,were constructed to obtain time history data of base shear and top center displacement,which revealed the seismic performance of horizontal storage tanks.Numerical simulation results and experimental data showed good agreement,with an error rate of less than 18.85%.And this conformity signifies the rationality of the SPH-FEM coupling method.The base shear and top center displacement values obtained by the coupled SPH-FEM method were only 53.3% to 69.1% of those calculated by the equivalent mass method employed in the current code.As the stored liquid volume increased,the seismic response of the horizontal storage tank exhibited a gradual upward trend,with the seismic response increasing from 73% to 388% for every 35% increase in stored liquid volume.The maximum von Mises stress of the tank and the supports remained below the steel yield strength during the earthquake.The coupled SPH-FEM method holds certain advantages in studying the seismic problems of tanks with complex structural forms,particularly due to the representation of the flow field distribution during earthquakes by involving reservoir fluid participation.