This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through...This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of arbitrary r balls of the 3-Torus. The continuity equation is satisfied together with spatially periodic boundary conditions. The yicomponent forcing terms consist of a function F as part of its expression that is arbitrarily small in an r ball where it is associated with a singular forcing expression both for inviscid and viscous cases. As a result, a significant simplification occurs with a v3(vifor all velocity components) only governing PDE resulting. The extension of three restricted subspaces in each of the principal directions in the Cartesian plane is shown as the Cartesian product ℋ=Jx,t×Jy,t×Jz,t. On each of these subspaces vi,i=1,2,3is continuous and there exists a linear independent subspace associated with the argument of the W function. Here the 3-Torus is built up from each compact segment of length 2R on each of the axes on the 3 principal directions x, y, and z. The form of the scaled velocities for non zero scaled δis related to the definition of the W function such that e−W(ξ)=W(ξ)ξwhere ξdepends on t and proportional to δ→0for infinite time t. The ratio Wξis equal to 1, making the limit δ→0finite and well defined. Considering r balls where the function F=(x−ai)2+(y−bi)2+(z−ci)2−ηset equal to −1e+rwhere r>0. is such that the forcing is singular at every distance r of centres of cubes each containing an r-ball. At the centre of the balls, the forcing is infinite. The main idea is that a system of singular initial value problems with infinite forcing is to be solved for where the velocities are shown to be locally Hölder continuous. It is proven that the limit of these singular problems shifts the finite time blowup time ti∗for first and higher derivatives to t=∞thereby indicating that there is no finite time blowup. Results in the literature can provide a systematic approach to study both large space and time behaviour for singular solutions to the Navier Stokes equations. Among the references, it has been shown that mathematical tools can be applied to study the asymptotic properties of solutions.展开更多
Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes...Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.展开更多
The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric an...The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric and Magnetic fields. Also, every moving particle has a De Broglie wavelength determined by its mass and velocity. This paper shows that all of these properties of a particle can be derived from a single wave function equation for that particle. Wave functions for the Electron and the Positron are presented and principles are provided that can be used to calculate the wave functions of all the fundamental particles in Physics. Fundamental particles such as electrons and positrons are considered to be point particles in the Standard Model of Physics and are not considered to have a structure. This paper demonstrates that they do indeed have structure and that this structure extends into the space around the particle’s center (in fact, they have infinite extent), but with rapidly diminishing energy density with the distance from that center. The particles are formed from Electromagnetic standing waves, which are stable solutions to the Schrödinger and Classical wave equations. This stable structure therefore accounts for both the wave and particle nature of these particles. In fact, all of their properties such as mass, spin and electric charge, can be accounted for from this structure. These particle properties appear to originate from a single point at the center of the wave function structure, in the same sort of way that the Shell theorem of gravity causes the gravity of a body to appear to all originate from a central point. This paper represents the first two fully characterized fundamental particles, with a complete description of their structure and properties, built up from the underlying Electromagnetic waves that comprise these and all fundamental particles.展开更多
Effective transportation systems lead to the efficient movement of goods and people, which significantly contribute to the quality of life in every society. In the heart of every economic and social development, there...Effective transportation systems lead to the efficient movement of goods and people, which significantly contribute to the quality of life in every society. In the heart of every economic and social development, there is always a transportation system. Mathematically the problem of modeling vehicle traffic flow can be solved at two main observation scales: The microscopic and the macroscopic levels. In the microscopic level, every vehicle is considered individually, and therefore, for every vehicle, we have an equation that is usually an ordinary differential equation (ODE). At a macroscopic level, we use from the dynamics models, where we have a system of partial differential equation, which involves variables such as density, speed, and flow rate of traffic stream with respect to time and space. Therefore, considering above content, this study has tried to compare solution of equation of macroscopic flow considering linear form (speed-density) and applying boundary condition that resulting to form solved is non-linear one-order partial differential equation (sharpy method) with non-linear assuming (speed and density) and consequently homographic nonlinear relation (speed-density). The recent case clearly gives more significant speeds than linear case of speed and density that can be a good scientific basis. In terms of safety for accidents and traffic signal, just as a reminder, but it is resulted of the reality that generally solutions of partial differential equations can have different forms. Therefore, the solution of partial differential equation (macroscopic flow) can have different answers and solutions so that all of these solutions apply in PDE (equation of macroscopic flow). Thus, under this condition, we can have solution of linear equation similar to greenberg or greenshield & android that are explained in logarithm and exponential function, but this article is based mostly on nonlinear solution of macroscopic equation, provided that existing nonlinear relationship between speed and density (homographic the second degree function). As mentioned above, as it gives more reliable and reasonable speeds than greenshield case, it will have more safety. This article has been provided in this field.展开更多
This paper presents an extensive survey of the most commonly used tools for diagnosing unbalanced flow in the atmosphere, namely the Lagrangian Rossby number, Psi vector, divergence equation, nonlinear balance equatio...This paper presents an extensive survey of the most commonly used tools for diagnosing unbalanced flow in the atmosphere, namely the Lagrangian Rossby number, Psi vector, divergence equation, nonlinear balance equation, generalized omega-equation, and departure from fields obtained by potential vorticity (PV) inversion. The basic thoery, assumptions as well as implementation and limitations for each of the tools are all discussed. These tools are applied to high—resolution mesoscale model data to assess the role of unbalanced dynamics in the generation of a mesoscale gravity wave event over the East Coast of the United States. Comparison of these tools in this case study shows that these various methods agree to a large extent with each other though they differ in details. Key words Unbalanced flow - Geostrophic adjustment - Gravity waves - Nonlinear balance equation - Potential vorticity inversion - Omega equations - Rossby number This research was conducted under support from NSF grant ATM-9700626 of the United States. The numerical computations described herein were performed on the Cray T90 at the North Carolina Supercomputing Center and the Cray supercomputer at the NCAR Scientific Computing Division, which also provided the initialization fields for the MM5. Thanks are extended to Mark Stoelinga at University of Washington for the RIP post-processing package.展开更多
A constitutive equation theory of Oldroyd fluid B type,i.e.the co-rotational derivative type,is developed for the anisotropic-viscoelastic fluid of liquid crystalline(LC)polymer.Analyzing the influence of the orientat...A constitutive equation theory of Oldroyd fluid B type,i.e.the co-rotational derivative type,is developed for the anisotropic-viscoelastic fluid of liquid crystalline(LC)polymer.Analyzing the influence of the orientational motion on the material behavior and neglecting the influence,the constitutive equation is applied to a simple case for the hydrodynamic motion when the orientational contribution is neglected in it and the anisotropic relaxation,retardation times and anisotropic viscosi- ties are introduced to describe the macroscopic behavior of the anisotropic LC polymer fluid.Using the equation for the shear flow of LC polymer fluid,the analytical expressions of the apparent viscosity and the normal stress differences are given which are in a good agreement with the experimental results of Baek et al.For the fiber spinning flow of the fluid,the analytical expression of the extensional viscosity is given.展开更多
A trigonometric series expansion method and two similar modified methods for the Orr-Sommerfeld equation are presented. These methods use the trigonometric series expansion with an auxiliary function added to the high...A trigonometric series expansion method and two similar modified methods for the Orr-Sommerfeld equation are presented. These methods use the trigonometric series expansion with an auxiliary function added to the highest order derivative of the unknown function and generate the lower order derivatives through successive integra- tions. The proposed methods are easy to implement because of the simplicity of the chosen basis functions. By solving the plane Poiseuille flow (PPF), plane Couette flow (PCF), and Blasius boundary layer flow with several homogeneous boundary conditions, it is shown that these methods yield results with the same accuracy as that given by the conventional Chebyshev collocation method but with better robustness, and that ob- tained by the finite difference method but with fewer modal number.展开更多
In this paper, problems of the flow over a fat plate in the large Reynolds numbercase are studied by using the method of multiple scales ̄[1,2].We have obtained N-orderuniformly valid asymptotic solutions of the Naver...In this paper, problems of the flow over a fat plate in the large Reynolds numbercase are studied by using the method of multiple scales ̄[1,2].We have obtained N-orderuniformly valid asymptotic solutions of the Naver-Stodes equations.展开更多
The authors propose a numerical algorithm for the two-dimensional Navier-Stokes equations written in stream function-vorticity formulation. The total time derivative term is treated with a first order characteristics ...The authors propose a numerical algorithm for the two-dimensional Navier-Stokes equations written in stream function-vorticity formulation. The total time derivative term is treated with a first order characteristics method. The space approximation is based on a piecewise continuous finite element method. The proposed algorithm is used to simulate the mechanical aeration process in lakes. Such process is used to combat the degradation of the water quality due to the eutrophication phenomena. For this application high computing facilities and capacities are required. In order to optimize the computing time and make possible the simulation of real applications, the authors propose a parallel implementation of the numerical algorithm. The parallelization technique is performed using the Message Passing Interface. The efficiency of the proposed numerical algorithm is illustrated by some numerical results.展开更多
By applying bromide ion as tracer, the channeling flow has been quantitatively described in saline rice soil and alkaline soil of Da'an City, Jilin Province of China. Breakthrough curves of bromide ion in the saline ...By applying bromide ion as tracer, the channeling flow has been quantitatively described in saline rice soil and alkaline soil of Da'an City, Jilin Province of China. Breakthrough curves of bromide ion in the saline rice soils after 1-year cultivation and 5-year cultivation and alkaline soil have been attained. Results show that the rice cultivation practice can improve the alkaline soil structure, however, it can accelerate the development of channeling flow pathway. Therefore, the channeling flow pathway has been developed widely in saline rice soil, but rarely in the alkaline soil. Three models of convection-dispersion equation (CDE), transfer functional model (TFM) and Back-Progation Network (BP Network) were used to simulate the transportation process of bromide ion. The peaks of probability density function of saline rice soil are higher with left skewed feature compared with that of the alkaline soil. It shows that the TIM and CDE can simulate the transportation process of the bromide ion in saline rice soil after 5-year cultivation, however, some deviation exists when it was used to simulate transportation process of bromide ion in saline rice soil after 1-year cultivation and alkaline soil; BP network can effectively simulate transportation process of bromide ion in both saline rice soil and alkaline soil.展开更多
To reduce computational costs, an improved form of the frequency domain boundary element method(BEM) is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitr...To reduce computational costs, an improved form of the frequency domain boundary element method(BEM) is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation(BIE) representation solves the two-dimensional convected Helmholtz equation(CHE) and its fundamental solution, which must satisfy a new Sommerfeld radiation condition(SRC) in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green's kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.展开更多
This work presents a new approach for simulating the random waves in viscous fluids and the associated bottom shear stresses. By generating the incident random waves in a numerical wave flume and solving the unsteady ...This work presents a new approach for simulating the random waves in viscous fluids and the associated bottom shear stresses. By generating the incident random waves in a numerical wave flume and solving the unsteady two-dimensional Navier-Stokes equations and the fully nonlinear free surface boundaiy conditions for the fluid flows in the flume, the viscous flows and laminar bottom shear stresses induced by random waves axe determined. The deterministic spectral amplitude method implemented by use of the fast Fourier transform algorithm was adopted to generate the incident random waves. The accuracy of the numerical scheme is confirmed by comparing the predicted wave spectrum with the target spectrum and by comparing the nanlerical transfer function between the shear stress and the surface elevation with the theoretical transfer function. The maximum bottom shear stress caused by random waves, computed by this wave model, is compared with that obtained by Myrhaug' s model (1995). The transfer function method is also employed to determine the maximum shear stress, and is proved accurate.展开更多
In ihis paper.an instability problem of an unsteady oscillationo flow is studied.In particultar,the phase.function of the disturbance wave.system is soived by using the charocteristic theory of partial differential eq...In ihis paper.an instability problem of an unsteady oscillationo flow is studied.In particultar,the phase.function of the disturbance wave.system is soived by using the charocteristic theory of partial differential equation and an expansion of Orysommerfeid eigenvalue problem.instead of using the disturbance model which is given previously The.flow considered is a combination of plane Poiseuille.flow with aflow oscillating periodically and its instability is found for a special initial value of a developing wave due to continuous oscillationg source.展开更多
In this paper, Beltrami vector fields in several orthogonal coordinate systems are obtained analytically and numerically. Specifically, axisymmetric incompressible inviscid steady state Beltrami (Trkalian) fluid flows...In this paper, Beltrami vector fields in several orthogonal coordinate systems are obtained analytically and numerically. Specifically, axisymmetric incompressible inviscid steady state Beltrami (Trkalian) fluid flows are obtained with the motivation to model flows that have been hypothesized to occur in tornadic flows. The studied coordinate systems include those that appear amenable to modeling such flows: the cylindrical, spherical, paraboloidal, and prolate and oblate spheroidal systems. The usual Euler equations are reformulated using the Bragg-Hawthorne equation for the stream function of the flow, which is solved analytically or numerically in each coordinate system under the assumption of separability of variables. Many of the obtained flows are visualized via contour plots of their stream functions in the <em>rz</em>-plane. Finally, the results are combined to provide a qualitative quasi-static model for a progression of tornado-like flows that develop as swirl increases. The results in this paper are equally applicable in electromagnetics, where the equivalent concept is that of a force-free magnetic field.展开更多
This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity dis...This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity distribution function obtained by the DSMC method converges to a modified form of the Boltzmann equation,which is the equation of the gas-kinetic unified algorithm to directly solve the molecular velocity distribution function.Their convergence is derived through mathematical treatment.The collision frequency is presented using various molecular models and the local equilibrium distribution function is obtained by Enskog expansion using the converged equation of the DSMC method.These two expressions agree with those used in the unified algorithm.Numerical validation of the converging consistency between these two approaches is illustrated by simulating the pressure driven Poiseuille flow in the slip transition flow regime and the two-dimensional and three-dimensional flows around a circular cylinder and spherical-cone reentry body covering the whole flow regimes from low speed micro-channel flow to high speed non-equilibrium aerothermodynamics.展开更多
The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate me...The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate method of the gas kinetic theory is studied and applied to simulate the complex multi-scale flows. Based on the uncoupling technique on molecular movement and colliding in the DSMC method, the gas-kinetic finite difference scheme is constructed to directly solve the discrete velocity distribution functions by extending and applying the unsteady time-splitting method from computational fluid dynamics. The Gauss-type discrete velocity numerical quadrature technique for different Mach number flows is developed to evaluate the macroscopic flow parameters in the physical space. As a result, the gas-kinetic numerical algorithm is established to study the three-dimensional complex flows from rarefied transition to continuum regimes. The parallel strategy adapted to the gas-kinetic numerical algorithm is investigated by analyzing the inner parallel degree of the algorithm, and then the HPF parallel processing program is developed. To test the reliability of the present gas-kinetic numerical method, the three-dimensional complex flows around sphere and spacecraft shape with various Knudsen numbers are simulated by HPF parallel computing. The computational results are found in high resolution of the flow fields and good agreement with the theoretical and experimental data. The computing practice has confirmed that the present gas-kinetic algorithm probably provides a promising approach to resolve the hypersonic aerothermodynamic problems with the complete spectrum of flow regimes from the gas-kinetic point of view of solving the Boltzmann model equation.展开更多
In this paper,we consider the asymptotic dynamics of the skew-product semiflow generated by the following time almost periodically-forced scalar reaction-diffusion equation:u_(t)=u_(xx)+f(t,u,u_(x)),t>0,0<x<L...In this paper,we consider the asymptotic dynamics of the skew-product semiflow generated by the following time almost periodically-forced scalar reaction-diffusion equation:u_(t)=u_(xx)+f(t,u,u_(x)),t>0,0<x<L(0.1)with the periodic boundary condition u(t,0)=u(t,L),u_(x)(t,0)=u_(x)(t,L),(0.2)where f is uniformly almost periodic in t.In particular,we study the topological structure of the limit sets of the skew-product semiflow.It is proved that any compact minimal invariant set(throughout this paper,we refer to it as a minimal set)can be residually embedded into an invariant set of some almost automorphically-forced flow on a circle S^(1)=R/LZ(see Definition 2.4 for“residually embedded”).Particularly,if f(t,u,p)=f(t,u,-p),then the flow on a minimal set can be embedded into an almost periodically-forced minimal flow on R(see Definition 2.4 for“embedded”).Moreover,it is proved that the ω-limit set of any bounded orbit contains at most two minimal sets that cannot be obtained from each other by phase translation.In addition,we further consider the asymptotic dynamics of the skew-product semiflow generated by(0.1)with the Neumann boundary condition u_(x)(t,0)=u_(x)(t,L)=0 or the Dirichlet boundary condition u(t,0)=u(t,L)=0.For such a system,it has been known that theω-limit set of any bounded orbit contains at most two minimal sets.By applying the new results for(0.1)+(0.2),under certain direct assumptions on f,we prove in this paper that the flow on any minimal set of(0.1)with the Neumann boundary condition or the Dirichlet boundary condition can be embedded into an almost periodically-forced minimal flow on R.Finally,a counterexample is given to show that even for quasi-periodically-forced equations,the results we obtain here cannot be further improved in general.展开更多
文摘This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of arbitrary r balls of the 3-Torus. The continuity equation is satisfied together with spatially periodic boundary conditions. The yicomponent forcing terms consist of a function F as part of its expression that is arbitrarily small in an r ball where it is associated with a singular forcing expression both for inviscid and viscous cases. As a result, a significant simplification occurs with a v3(vifor all velocity components) only governing PDE resulting. The extension of three restricted subspaces in each of the principal directions in the Cartesian plane is shown as the Cartesian product ℋ=Jx,t×Jy,t×Jz,t. On each of these subspaces vi,i=1,2,3is continuous and there exists a linear independent subspace associated with the argument of the W function. Here the 3-Torus is built up from each compact segment of length 2R on each of the axes on the 3 principal directions x, y, and z. The form of the scaled velocities for non zero scaled δis related to the definition of the W function such that e−W(ξ)=W(ξ)ξwhere ξdepends on t and proportional to δ→0for infinite time t. The ratio Wξis equal to 1, making the limit δ→0finite and well defined. Considering r balls where the function F=(x−ai)2+(y−bi)2+(z−ci)2−ηset equal to −1e+rwhere r>0. is such that the forcing is singular at every distance r of centres of cubes each containing an r-ball. At the centre of the balls, the forcing is infinite. The main idea is that a system of singular initial value problems with infinite forcing is to be solved for where the velocities are shown to be locally Hölder continuous. It is proven that the limit of these singular problems shifts the finite time blowup time ti∗for first and higher derivatives to t=∞thereby indicating that there is no finite time blowup. Results in the literature can provide a systematic approach to study both large space and time behaviour for singular solutions to the Navier Stokes equations. Among the references, it has been shown that mathematical tools can be applied to study the asymptotic properties of solutions.
文摘Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.
文摘The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric and Magnetic fields. Also, every moving particle has a De Broglie wavelength determined by its mass and velocity. This paper shows that all of these properties of a particle can be derived from a single wave function equation for that particle. Wave functions for the Electron and the Positron are presented and principles are provided that can be used to calculate the wave functions of all the fundamental particles in Physics. Fundamental particles such as electrons and positrons are considered to be point particles in the Standard Model of Physics and are not considered to have a structure. This paper demonstrates that they do indeed have structure and that this structure extends into the space around the particle’s center (in fact, they have infinite extent), but with rapidly diminishing energy density with the distance from that center. The particles are formed from Electromagnetic standing waves, which are stable solutions to the Schrödinger and Classical wave equations. This stable structure therefore accounts for both the wave and particle nature of these particles. In fact, all of their properties such as mass, spin and electric charge, can be accounted for from this structure. These particle properties appear to originate from a single point at the center of the wave function structure, in the same sort of way that the Shell theorem of gravity causes the gravity of a body to appear to all originate from a central point. This paper represents the first two fully characterized fundamental particles, with a complete description of their structure and properties, built up from the underlying Electromagnetic waves that comprise these and all fundamental particles.
文摘Effective transportation systems lead to the efficient movement of goods and people, which significantly contribute to the quality of life in every society. In the heart of every economic and social development, there is always a transportation system. Mathematically the problem of modeling vehicle traffic flow can be solved at two main observation scales: The microscopic and the macroscopic levels. In the microscopic level, every vehicle is considered individually, and therefore, for every vehicle, we have an equation that is usually an ordinary differential equation (ODE). At a macroscopic level, we use from the dynamics models, where we have a system of partial differential equation, which involves variables such as density, speed, and flow rate of traffic stream with respect to time and space. Therefore, considering above content, this study has tried to compare solution of equation of macroscopic flow considering linear form (speed-density) and applying boundary condition that resulting to form solved is non-linear one-order partial differential equation (sharpy method) with non-linear assuming (speed and density) and consequently homographic nonlinear relation (speed-density). The recent case clearly gives more significant speeds than linear case of speed and density that can be a good scientific basis. In terms of safety for accidents and traffic signal, just as a reminder, but it is resulted of the reality that generally solutions of partial differential equations can have different forms. Therefore, the solution of partial differential equation (macroscopic flow) can have different answers and solutions so that all of these solutions apply in PDE (equation of macroscopic flow). Thus, under this condition, we can have solution of linear equation similar to greenberg or greenshield & android that are explained in logarithm and exponential function, but this article is based mostly on nonlinear solution of macroscopic equation, provided that existing nonlinear relationship between speed and density (homographic the second degree function). As mentioned above, as it gives more reliable and reasonable speeds than greenshield case, it will have more safety. This article has been provided in this field.
文摘This paper presents an extensive survey of the most commonly used tools for diagnosing unbalanced flow in the atmosphere, namely the Lagrangian Rossby number, Psi vector, divergence equation, nonlinear balance equation, generalized omega-equation, and departure from fields obtained by potential vorticity (PV) inversion. The basic thoery, assumptions as well as implementation and limitations for each of the tools are all discussed. These tools are applied to high—resolution mesoscale model data to assess the role of unbalanced dynamics in the generation of a mesoscale gravity wave event over the East Coast of the United States. Comparison of these tools in this case study shows that these various methods agree to a large extent with each other though they differ in details. Key words Unbalanced flow - Geostrophic adjustment - Gravity waves - Nonlinear balance equation - Potential vorticity inversion - Omega equations - Rossby number This research was conducted under support from NSF grant ATM-9700626 of the United States. The numerical computations described herein were performed on the Cray T90 at the North Carolina Supercomputing Center and the Cray supercomputer at the NCAR Scientific Computing Division, which also provided the initialization fields for the MM5. Thanks are extended to Mark Stoelinga at University of Washington for the RIP post-processing package.
基金The project supported by the National Natural Science Foundation of China(19832050 and 10372100)
文摘A constitutive equation theory of Oldroyd fluid B type,i.e.the co-rotational derivative type,is developed for the anisotropic-viscoelastic fluid of liquid crystalline(LC)polymer.Analyzing the influence of the orientational motion on the material behavior and neglecting the influence,the constitutive equation is applied to a simple case for the hydrodynamic motion when the orientational contribution is neglected in it and the anisotropic relaxation,retardation times and anisotropic viscosi- ties are introduced to describe the macroscopic behavior of the anisotropic LC polymer fluid.Using the equation for the shear flow of LC polymer fluid,the analytical expressions of the apparent viscosity and the normal stress differences are given which are in a good agreement with the experimental results of Baek et al.For the fiber spinning flow of the fluid,the analytical expression of the extensional viscosity is given.
基金supported by the National Natural Science Foundation of China(Nos.11221062,11521091,and 91752203)
文摘A trigonometric series expansion method and two similar modified methods for the Orr-Sommerfeld equation are presented. These methods use the trigonometric series expansion with an auxiliary function added to the highest order derivative of the unknown function and generate the lower order derivatives through successive integra- tions. The proposed methods are easy to implement because of the simplicity of the chosen basis functions. By solving the plane Poiseuille flow (PPF), plane Couette flow (PCF), and Blasius boundary layer flow with several homogeneous boundary conditions, it is shown that these methods yield results with the same accuracy as that given by the conventional Chebyshev collocation method but with better robustness, and that ob- tained by the finite difference method but with fewer modal number.
文摘In this paper, problems of the flow over a fat plate in the large Reynolds numbercase are studied by using the method of multiple scales ̄[1,2].We have obtained N-orderuniformly valid asymptotic solutions of the Naver-Stodes equations.
文摘The authors propose a numerical algorithm for the two-dimensional Navier-Stokes equations written in stream function-vorticity formulation. The total time derivative term is treated with a first order characteristics method. The space approximation is based on a piecewise continuous finite element method. The proposed algorithm is used to simulate the mechanical aeration process in lakes. Such process is used to combat the degradation of the water quality due to the eutrophication phenomena. For this application high computing facilities and capacities are required. In order to optimize the computing time and make possible the simulation of real applications, the authors propose a parallel implementation of the numerical algorithm. The parallelization technique is performed using the Message Passing Interface. The efficiency of the proposed numerical algorithm is illustrated by some numerical results.
基金Under the auspices of the Key Innovation Project of Chinese Academy of Sciences (No. KZCX1-SW-19-02)
文摘By applying bromide ion as tracer, the channeling flow has been quantitatively described in saline rice soil and alkaline soil of Da'an City, Jilin Province of China. Breakthrough curves of bromide ion in the saline rice soils after 1-year cultivation and 5-year cultivation and alkaline soil have been attained. Results show that the rice cultivation practice can improve the alkaline soil structure, however, it can accelerate the development of channeling flow pathway. Therefore, the channeling flow pathway has been developed widely in saline rice soil, but rarely in the alkaline soil. Three models of convection-dispersion equation (CDE), transfer functional model (TFM) and Back-Progation Network (BP Network) were used to simulate the transportation process of bromide ion. The peaks of probability density function of saline rice soil are higher with left skewed feature compared with that of the alkaline soil. It shows that the TIM and CDE can simulate the transportation process of the bromide ion in saline rice soil after 5-year cultivation, however, some deviation exists when it was used to simulate transportation process of bromide ion in saline rice soil after 1-year cultivation and alkaline soil; BP network can effectively simulate transportation process of bromide ion in both saline rice soil and alkaline soil.
基金supported by National Engineering School of Tunis (No.13039.1)
文摘To reduce computational costs, an improved form of the frequency domain boundary element method(BEM) is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation(BIE) representation solves the two-dimensional convected Helmholtz equation(CHE) and its fundamental solution, which must satisfy a new Sommerfeld radiation condition(SRC) in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green's kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.
基金the Science Council (Grant No. NSC95-2221-E-006-474)
文摘This work presents a new approach for simulating the random waves in viscous fluids and the associated bottom shear stresses. By generating the incident random waves in a numerical wave flume and solving the unsteady two-dimensional Navier-Stokes equations and the fully nonlinear free surface boundaiy conditions for the fluid flows in the flume, the viscous flows and laminar bottom shear stresses induced by random waves axe determined. The deterministic spectral amplitude method implemented by use of the fast Fourier transform algorithm was adopted to generate the incident random waves. The accuracy of the numerical scheme is confirmed by comparing the predicted wave spectrum with the target spectrum and by comparing the nanlerical transfer function between the shear stress and the surface elevation with the theoretical transfer function. The maximum bottom shear stress caused by random waves, computed by this wave model, is compared with that obtained by Myrhaug' s model (1995). The transfer function method is also employed to determine the maximum shear stress, and is proved accurate.
文摘In ihis paper.an instability problem of an unsteady oscillationo flow is studied.In particultar,the phase.function of the disturbance wave.system is soived by using the charocteristic theory of partial differential equation and an expansion of Orysommerfeid eigenvalue problem.instead of using the disturbance model which is given previously The.flow considered is a combination of plane Poiseuille.flow with aflow oscillating periodically and its instability is found for a special initial value of a developing wave due to continuous oscillationg source.
文摘In this paper, Beltrami vector fields in several orthogonal coordinate systems are obtained analytically and numerically. Specifically, axisymmetric incompressible inviscid steady state Beltrami (Trkalian) fluid flows are obtained with the motivation to model flows that have been hypothesized to occur in tornadic flows. The studied coordinate systems include those that appear amenable to modeling such flows: the cylindrical, spherical, paraboloidal, and prolate and oblate spheroidal systems. The usual Euler equations are reformulated using the Bragg-Hawthorne equation for the stream function of the flow, which is solved analytically or numerically in each coordinate system under the assumption of separability of variables. Many of the obtained flows are visualized via contour plots of their stream functions in the <em>rz</em>-plane. Finally, the results are combined to provide a qualitative quasi-static model for a progression of tornado-like flows that develop as swirl increases. The results in this paper are equally applicable in electromagnetics, where the equivalent concept is that of a force-free magnetic field.
基金supported by the National Natural Science Foundation of China (Grant Nos. 91016027 and 91130018)
文摘This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity distribution function obtained by the DSMC method converges to a modified form of the Boltzmann equation,which is the equation of the gas-kinetic unified algorithm to directly solve the molecular velocity distribution function.Their convergence is derived through mathematical treatment.The collision frequency is presented using various molecular models and the local equilibrium distribution function is obtained by Enskog expansion using the converged equation of the DSMC method.These two expressions agree with those used in the unified algorithm.Numerical validation of the converging consistency between these two approaches is illustrated by simulating the pressure driven Poiseuille flow in the slip transition flow regime and the two-dimensional and three-dimensional flows around a circular cylinder and spherical-cone reentry body covering the whole flow regimes from low speed micro-channel flow to high speed non-equilibrium aerothermodynamics.
基金Supported by the National Natural Science Foundation of China (Grant Nos. 90205009 and 10321002)the National Parallel Computing Center
文摘The Boltzmann simplified velocity distribution function equation describing the gas transfer phenomena from various flow regimes will be explored and solved numerically in this study. The discrete velocity ordinate method of the gas kinetic theory is studied and applied to simulate the complex multi-scale flows. Based on the uncoupling technique on molecular movement and colliding in the DSMC method, the gas-kinetic finite difference scheme is constructed to directly solve the discrete velocity distribution functions by extending and applying the unsteady time-splitting method from computational fluid dynamics. The Gauss-type discrete velocity numerical quadrature technique for different Mach number flows is developed to evaluate the macroscopic flow parameters in the physical space. As a result, the gas-kinetic numerical algorithm is established to study the three-dimensional complex flows from rarefied transition to continuum regimes. The parallel strategy adapted to the gas-kinetic numerical algorithm is investigated by analyzing the inner parallel degree of the algorithm, and then the HPF parallel processing program is developed. To test the reliability of the present gas-kinetic numerical method, the three-dimensional complex flows around sphere and spacecraft shape with various Knudsen numbers are simulated by HPF parallel computing. The computational results are found in high resolution of the flow fields and good agreement with the theoretical and experimental data. The computing practice has confirmed that the present gas-kinetic algorithm probably provides a promising approach to resolve the hypersonic aerothermodynamic problems with the complete spectrum of flow regimes from the gas-kinetic point of view of solving the Boltzmann model equation.
基金supported by National Science Foundation of USA(Grant No.DMS1645673)supported by National Natural Science Foundation of China(Grant Nos.11825106,11771414 and 12090012)+2 种基金Wu Wen-Tsun Key Laboratory of Mathematics,Chinese Academy of Sciences and University of Science and Technology of Chinasupported by National Natural Science Foundation of China(Grant Nos.11971232,12071217 and 11601498)the Chinese Scholarship Council(Grant No.201906845011)for its financial support。
文摘In this paper,we consider the asymptotic dynamics of the skew-product semiflow generated by the following time almost periodically-forced scalar reaction-diffusion equation:u_(t)=u_(xx)+f(t,u,u_(x)),t>0,0<x<L(0.1)with the periodic boundary condition u(t,0)=u(t,L),u_(x)(t,0)=u_(x)(t,L),(0.2)where f is uniformly almost periodic in t.In particular,we study the topological structure of the limit sets of the skew-product semiflow.It is proved that any compact minimal invariant set(throughout this paper,we refer to it as a minimal set)can be residually embedded into an invariant set of some almost automorphically-forced flow on a circle S^(1)=R/LZ(see Definition 2.4 for“residually embedded”).Particularly,if f(t,u,p)=f(t,u,-p),then the flow on a minimal set can be embedded into an almost periodically-forced minimal flow on R(see Definition 2.4 for“embedded”).Moreover,it is proved that the ω-limit set of any bounded orbit contains at most two minimal sets that cannot be obtained from each other by phase translation.In addition,we further consider the asymptotic dynamics of the skew-product semiflow generated by(0.1)with the Neumann boundary condition u_(x)(t,0)=u_(x)(t,L)=0 or the Dirichlet boundary condition u(t,0)=u(t,L)=0.For such a system,it has been known that theω-limit set of any bounded orbit contains at most two minimal sets.By applying the new results for(0.1)+(0.2),under certain direct assumptions on f,we prove in this paper that the flow on any minimal set of(0.1)with the Neumann boundary condition or the Dirichlet boundary condition can be embedded into an almost periodically-forced minimal flow on R.Finally,a counterexample is given to show that even for quasi-periodically-forced equations,the results we obtain here cannot be further improved in general.