A Fixed-Point Fast Sweeping WENO Method with Inverse Lax-Wendroff Boundary Treatment for Steady State of Hyperbolic Conservation Laws

Fixed-point fast sweeping WENO methods are a class of efficient high-order numerical methods to solve steady-state solutions of hyperbolic partial differential equations(PDEs).The Gauss-Seidel iterations and alternating sweeping strategy are used to cover characteristics of hyperbolic PDEs in each sweeping order to achieve fast convergence rate to steady-state solutions.A nice property of fixed-point fast sweeping WENO methods which distinguishes them from other fast sweeping methods is that they are explicit and do not require inverse operation of nonlinear local systems.Hence,they are easy to be applied to a general hyperbolic system.To deal with the difficulties associated with numerical boundary treatment when high-order finite difference methods on a Cartesian mesh are used to solve hyperbolic PDEs on complex domains,inverse Lax-Wendroff(ILW)procedures were developed as a very effective approach in the literature.In this paper,we combine a fifthorder fixed-point fast sweeping WENO method with an ILW procedure to solve steadystate solution of hyperbolic conservation laws on complex computing regions.Numerical experiments are performed to test the method in solving various problems including the cases with the physical boundary not aligned with the grids.Numerical results show highorder accuracy and good performance of the method.Furthermore,the method is compared with the popular third-order total variation diminishing Runge-Kutta(TVD-RK3)time-marching method for steady-state computations.Numerical examples show that for most of examples,the fixed-point fast sweeping method saves more than half CPU time costs than TVD-RK3 to converge to steady-state solutions.

POSITIVE STEADY STATES AND DYNAMICS FOR A DIFFUSIVE PREDATOR-PREY SYSTEM WITH A DEGENERACY

In this article, we consider positive steady state solutions and dynamics for a spatially heterogeneous predator-prey system with modified Leslie-Gower and Holling-Type II schemes. The heterogeneity here is created by the degeneracy of the intra-specific pressures for the prey. By the bifurcation method, the degree theory, and a priori estimates, we discuss the existence and multiplicity of positive steady states. Moreover, by the comparison argument, we also discuss the dynamical behavior for the diffusive predator-prey system.

A Stochastic Extension of Macroscopic Stability Criterion of Nonequilibrium Steady State in Chemical Reaction Systems Governed by Master Equation

By means of both the separation of the perturbation in accordance with characteristic parnmeters and the Kramers Moyal-expansion of the master equation, it is shown that the time derivative of the partial excess quantity of stochastic entropy due to the deviation from the most probable path is related to the responsibility of a system to the external macroscopic perturbations. This evolution rate of the partial excess stochastic entropy is equivalent to the partlal excess stochastic entropy production, as well as the stochastic excess entropy production rate based on the stochastic potential npproach. It appears also as an eqivalent quantity of the Gibbs excess entropy production for the Polsson distribution. The macroscopic stability of chemical reaction systems is dominnted by this new stochastic quantity when the local equilibrium thermodynamics is broken down .

Existence for positive steady states of an eco-epidemiological model

An eco-epidemiological model with an epidemic in the predator and with a Holling type Ⅱ function is considered.A system with diffusion under the homogeneous Neumann boundary condition is studied.The existence for a positive solution of the corresponding steady state problem is mainly discussed.First,a prior estimates（positive upper and lower bounds） of the positive steady states of the reaction-diffusion system is given by the maximum principle and the Harnack inequation.Then,the non-existence of non-constant positive steady states by using the energy method is given.Finally,the existence of non-constant positive steady states is obtained by using the topological degree.

Rheological behavior of semi-solid 7075 aluminum alloy at steady state

The further application of semi-solid processing lies in the in-depth fundamental study like rheological behavior. In this research, the apparent viscosity of the semi-solid slurry of 7075 alloy was measured using a Couette type viscometer. The effects of solid fraction and shearing rate on the apparent viscosity of this alloy were investigated under different processing conditions. It can be seen that the apparent viscosity increases with an increase in the solid fraction from 10% to 50% （temperature 620 ℃ to 630 ℃） at steady state. When the solid fraction was fixed, the apparent viscosity can be decreased by altering the shearing rate from 61.235 s-1 to 489.88 s-1 at steady state. An empirical equation that shows the effects of solid fraction and shearing rate on the apparent viscosity is fitted： ηα = [0.871 - 0.00849. γ0.74924]. exp（3.7311, fs） . The microstructure of quenched samples was examined to understand the alloy＇s rheological behavior.

Experimental and Finite Element Study of Steady State Micro-cutting Characteristics of Aluminum Alloy (2A12)

This paper studies the micro-cutting characteristics of aluminum alloy (2A12) based on a series of orthogonal experiments and finite element method (FEM) simulations. An energy-based ductile failure law was proposed in the FEM simulation. The simulated cutting forces and chip morphology were compared with experimental results. The simulation result indicates that there is a close relationship between the cutting force and cutting heat. The micro-cutting force decreases as the heat flux vector increases. Both the cutting heat and the micro-cutting force need a finite time to achieve a steady state. It is observed that with the cutting speed of 169.95 m/min and uncut chip thickness of 6 μm, the heat flux vector in the workpiece increases to a stable value after 0.06 ms; meanwhile, the principal cutting force decreases to a steady state correspondingly, i.e., the micro-cutting process achieves the steady state. It is concluded that the steady state micro-cutting simulation can reflect the cutting process accurately.

40 Hz Auditory Steady State Response to Linguistic Features of Stimuli during Auditory Hallucinations

The auditory steady state response （ASSR） may reflect activity from different regions of the brain, depending on the modulation frequency used. In general, responses induced by low rates （_〈40 Hz） emanate mostly from central structures of the brain, and responses from high rates （〉80 Hz） emanate mostly from the peripheral auditory nerve or brainstem structures. Besides, it was reported that the gamma band ASSR （30-90 Hz） played an important role in working memory, speech understanding and recognition. This paper investigated the 40 Hz ASSR evoked by modulated speech and reversed speech. The speech was Chinese phrase voice, and the noise-like reversed speech was obtained by temporally reversing the speech. Both auditory stimuli were modulated with a frequency of 40 Hz. Ten healthy subjects and 5 patients with hallucination symptom participated in the experiment. Results showed re- duction in left auditory cortex response when healthy subjects listened to the reversed speech compared with the speech. In contrast, when the patients who experienced auditory hallucinations listened to the reversed speech, the auditory cortex of left hemispheric responded more actively. The ASSR results were consistent with the behavior results of patients. Therefore, the gamma band ASSR is expected to be helpful for rapid and objective diagnosis of hallucination in clinic.

A NEW APPROACH TO THE EQUATIONS FOR THE STEADY STATE CURRENTS AT MICROELECTRODES:EC'PROCESSES(SECOND ORDER REACTION)

With a new approach,the general current expressions of two typical second order catalytic reactions at microelectrodes are obtained.This allows the study of fast chemical reactions and systems where the reactants are present in similar concentrations.

Realization of minute-long steady-state H-mode discharges on EAST

In the 2016 EAST experimental campaign,a steady-state long-pulse H-mode discharge with an ITER-like tungsten divertor lasting longer than one minute has been obtained using only RF heating and current drive,through an integrated control of the wall conditioning,plasma configuration,divertor heat flux,particle exhaust,impurity management,and effective coupling of multiple RF heating and current drive sources at high injected power.The plasma current（Ip - 0.45 MA） was fully-noninductively driven（Vloop 〈 0.0 V） by a combination of-2.5 MW LHW,-0.4 MW ECH and -0.8 MW ICRF.This result demonstrates the progress of physics and technology studies on EAST,and will benefit the physics basis for steady state operation of ITER and CFETR.

Sawtooth-like oscillations and steady states caused by the m/n=2/1 double tearing mode

The sawtooth-like oscillations resulting from the m/n=2/1 double tearing mode(DTM)are numerically investigated through the three-dimensional,toroidal,nonlinear resistive-MHD code(CLT).We find that the nonlinear evolution of the m/n=2/1 DTM can lead to sawtooth-like oscillations,which are similar to those driven by the kink mode.The perpendicular thermal conductivity and the external heating rate can significantly alter the behaviors of the DTM driven sawtooth-like oscillations.With a high perpendicular thermal conductivity,the system quickly evolves into a steady state with m/n=2/1 magnetic islands and helical flow.However,with a low perpendicular thermal conductivity,the system tends to exhibit sawtooth-like oscillations.With a sufficiently high or low heating rate,the system exhibits sawtooth-like oscillations,while with an intermediate heating rate,the system quickly evolves into a steady state.At the steady state,there exist the non-axisymmetric magnetic field and strong radial flow,and both are with helicity of m/n=2/1.Like the steady state with m/n=1/1 radial flow,which is beneficial for preventing the helium ash accumulation in the core,the steady state with m/n=2/1 radial flow might also be a good candidate for the advanced steady state operations in future fusion reactors.We also find that the behaviors of the sawtooth-like oscillations are almost independent of tokamak geometry,which implies that the steady state with saturated m/n=2/1 islands might exist in different tokamaks.

STEADY STATE MOTIONS OF SHALLOW ARCH UNDER PERIODIC FORCE WITH 1:2 INTERNAL RESONANCEON THE PLANE OF PHYSICAL PARAMETERS

The bifurcation dynamics of shallow arch which possesses initial deflection under periodic excitation for the case of 1:2 internal resonance is studied in this paper. The whole parametric plane is divided into several different regions according to lire types of motions; then the distribution of steady state motions of shallow arch on the plane of physical parameters is obtained. Combining with numerical method, the dynamics of the system in different regions, especially in the Hopf bifurcation region, is studied in detail. The rule of the mode interaction and the route to chaos of the system is also analysed at the end.

Modeling and optimization of the cyclic steady state operation of adsorptive reactors

Adsorptive reactors(AR),in which an adsorptive functionality is incorporated into the catalytic reactors,offer enhanced performance over their conventional counterparts due to the effective manipulation of concentration and temperature profiles.The operation of these attractive reactors is,however,inherently unsteady state,complicating the design and operation of such sorption-enhanced processes.In order to capture,comprehend and capitalize upon the rich dynamic texture of adsorptive reactors,it is necessary to employ cyclic steady state algorithms describing the entire reaction-adsorption/desorption cycle.The stability of this cyclic steady state is of great importance for the design and operation of adsorptive reactors.In this paper,the cyclic steady state of previously proposed novel adsorptive reactor designs has been calculated and then optimized to give maximum space–time yields.The results obtained revealed unambiguously that an improvement potential of up to multifold level could be attained under the optimized cyclic steady state conditions.This additional improvement resulted from the reduction of the regeneration time well below the reaction-adsorption time,which means,in turn,more space–time yield.

Steady State Properties of One-Dimensional Non-uniform Granular Gases Subjected to Gaussian White Noise Driving

We present a model of non-uniform granular gases in one-dimensional case, whose granularity distribution has the fractal characteristic. We have studied the nonequilibrium properties of the system by means of Monte Carlo method. When the typical relaxation time T of the Brownian process is greater than the mean collision time To, the energy evolution of the system exponentially decays, with a tendency to achieve a stable asymptotic value, and the system finally reaches a nonequilibrium steady state in which the velocity distribution strongly deviates from the Gaussian one. Three other aspects have also been studied for the steady state： the visualized change of the particle density, the entropy of the system and the correlations in the velocity of particles. And the results of simulations indicate that the system has strong spatial clustering; Furthermore, the influence of the inelasticity and inhomogeneity on dynamic behaviors have also been extensively investigated, especially the dependence of the entropy and the correlations in the velocity of particles on the restitute coefficient e and the fractal dimension D.

Evaluation of The Thermal Performance of Multi-Element Doped Graphite under Steady-State High Heat Flux

Multi-element doped graphite, GBST1308 has been developed as a plasma facing material (PFM) for high heat flux components of the HT-7U device. The thermal performance of the material under steady-state (SS) high heat flux was evaluated under actively cooling conditions, the specimens were mechanically joined to copper heat sink with supercarbon sheet as a compliant layer between the interfaces. The experiments have been performed in a facility of ACT (actively cooling test stand) with a 100 kW electron gun in order to test the suitability and the loading limit of such materials. The surface temperature and bulk temperature distribution of the specimens were investigated. The experimental results are very encouraging that when heat flux is not more than 6 MW/m2, the surface temperature of GBST1308 is less than 1000℃, which is the lowest, compared with IG-430U and even with CX-2002U (CFC); The primary results indicate that the mechanically-joined material system by such a proper design as thin tile, super compliant layer, GBST as a PFM and copper-alloy heat sink, can be used as divertor plates for HT-7U in the first phase.

The Existence and Stability of Nontrivial Steady-State for a Delay Competition-Diffusion Model

This paper investigates the competition diffusion model with delay. The existence and stability of nontrivial steady state are studied.

Mass-Spring-Damper Model with Steady State Parameters for Predicting the Movement of Liquid Column and Temperature Oscillation in Loop Heat Pipe

In order to investigate the mechanism of the temperature oscillation in loop heat pipes,this paper investigated the movement of the phase interface as the changed input power by a mass-spring-damper model.The model was solved with MATLAB and was used to explain the high-frequency and low-amplitude temperature oscillation.Temperature variation with the input power from 20 W to 75 W was investigated based on a LHP prototype in a literature.The model agreed well with the experimental data in the literature.The simulation results suggested that the movement of the liquid column was caused by the fluctuation of pressure difference applied on the liquid column and the stiffness coefficients of the vapor springs increasing with the input power.According to parameter analyses,the temperature oscillation at the outlet of the condenser can be weakened by increasing the mass of the liquid column and keeping the temperature at the outlet of the condenser steady.

Steady State Creep Rate Equation of Inconel 718 Superalloy

The steady state creep rate equdtion of a nickel base superalloy Inconel 718, strengthened by coherent ordered disc-shaped bct γ^(11) phase and coherent spherical fcc γ~1 phase precipitates, has been established in the stress and temperature ranges of 620-840 MN m^(-2) and 853-943K, respecti- vely. Constant stress tensile creep tests were used to medsure the values of steady state creep rate, ε_s, and the consecutive stress reduction method was used to measure the back stress during creep deformation. The values of effective stress exponent, n_e, were detemined from the slopes of the lgε_s vs. lg(σ_a-σ_0)/G plots. The effect of grain size, d, on steady state creep rdte has been also studied in this investigation, and the grain size sensitive exponents m were detemined from the slopes of lgε_s vs. lg(b/d) plots. The creep rate equations of Inconel 718, in the above stress and temperature ranges, have been proposed to be ε_s=1.6×10^(-5)(D_1Gb/KT) (b/d )^(0.19)[(σ_a-σ_0)/G]^(1.35) in diffusional creep region, and ε_s =75(D_1Gb/KT) (b/d)^(-0.42)[(σ_a-σ_0)/G]^(5.5) in dislocation power law creep region.

A Numerical Model for Ion Charge Distribution of Plasmas in Collisional Radiative Steady State

A numerical model for the charge state distribution of plasmas in a collisional radiative steady state （CRSS） is established by averaging over the atomic process rate coefficients in universal kinetic equations. It is used to calculate the mean ion charge and ion population for a given temperature and density of the plasmas, ranging from low Z to high Z elements. The comparisons of the calculated results with those of other non-local thermodynamic equilibrium kinetics codes show that this model possesses acceptable precision. Furthermore, the NLTE effects are investigated by virtue of the model, and the differences between CRSS and LTE models for low density plasmas are quite evident.

Multiple steady state phenomenon in martensitic transformation

Based on the basic facts that the martensitic transformation is a physical phenomenon which occurs in non equilibrium conditions and there exists the feedback mechanism in the martensitic transformation, the dynamical processes of the isothermal and athermal martensitic transformations were analyzed by using nonlinear theory and a bifurcation theory model was established. It is shown that a multiple steady state phenomenon can take place as austenite is cooled, and the transitions of the steady state temperature between the branches of stable steady states can be considered the transformation from austenite to martensite. This model can estimate the starting temperature of the martensitic transformation and explain some experimental features of the martensitic transformation such as the effects of cooling rate, fluctuation and austenitic grain size on the martensitic transformation. [

Functional near-infrared spectroscopy can detect low-frequency hemodynamic oscillations in the prefrontal cortex during steady-state visual evoked potentialinducing periodic facial expression stimuli presentation

Brain oscillations are vital to cognitive functions,while disrupted oscillatory activity is linked to various brain disorders.Although high-frequency neural oscillations(>1 Hz)have been extensively studied in cognition,the neural mechanisms underlying low-frequency hemodynamic oscillations(LFHO)<1 Hz have not yet been fully explored.One way to examine oscillatory neural dynamics is to use a facial expression(FE)paradigm to induce steady-state visual evoked potentials(SSVEPs),which has been used in electroencephalography studies of high-frequency brain oscillation activity.In this study,LFHO during SSVEP-inducing periodic flickering stimuli presentation were inspected using functional near-infrared spectroscopy(fNIRS),in which hemodynamic responses in the prefrontal cortex were recorded while participants were passively viewing dynamic FEs flickering at 0.2 Hz.The fast Fourier analysis results demonstrated that the power exhibited monochronic peaks at 0.2 Hz across all channels,indicating that the periodic events successfully elicited LFHO in the prefrontal cortex.More importantly,measurement of LFHO can effectively distinguish the brain activation difference between different cognitive conditions,with happy FE presentation showing greater LFHO power than neutral FE presentation.These results demonstrate that stimuli flashing at a given frequency can induce LFHO in the prefrontal cortex,which provides new insights into the cognitive mechanisms involved in slow oscillation.