Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an...Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.展开更多
In order to improve the design results for the reconfigurable frequency response masking FRM filters an improved design method based on second-order cone programming SOCP is proposed.Unlike traditional methods that se...In order to improve the design results for the reconfigurable frequency response masking FRM filters an improved design method based on second-order cone programming SOCP is proposed.Unlike traditional methods that separately design the proposed method takes all the desired designing modes into consideration when designing all the subfilters. First an initial solution is obtained by separately designing the subfilters and then the initial solution is updated by iteratively solving a SOCP problem. The proposed method is evaluated on a design example and simulation results demonstrate that jointly designing all the subfilters can obtain significantly lower minimax approximation errors compared to the conventional design method.展开更多
In this paper, based on the second-order Taylor series expansion and the difference of convex functions algo- rithm for quadratic problems with box constraints （the DCA for QB）, a new method is proposed to solve the...In this paper, based on the second-order Taylor series expansion and the difference of convex functions algo- rithm for quadratic problems with box constraints （the DCA for QB）, a new method is proposed to solve the static response problem of structures with fairly large uncertainties in interval parameters. Although current methods are effective for solving the static response problem of structures with interval parameters with small uncertainties, these methods may fail to estimate the region of the static response of uncertain structures if the uncertainties in the parameters are fairly large. To resolve this problem, first, the general expression of the static response of structures in terms of structural parameters is derived based on the second-order Taylor series expansion. Then the problem of determining the bounds of the static response of uncertain structures is transformed into a series of quadratic problems with box constraints. These quadratic problems with box constraints can be solved using the DCA approach effectively. The numerical examples are given to illustrate the accuracy and the efficiency of the proposed method when comparing with other existing methods.展开更多
This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>...This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.展开更多
This work illustrates the application of the 1<sup>st</sup>-CASAM to a paradigm heat transport model which admits exact closed-form solutions. The closed-form expressions obtained in this work for the sens...This work illustrates the application of the 1<sup>st</sup>-CASAM to a paradigm heat transport model which admits exact closed-form solutions. The closed-form expressions obtained in this work for the sensitivities of the temperature distributions within the model to the model’s parameters, internal interfaces and external boundaries can be used to benchmark commercial and production software packages for simulating heat transport. The 1<sup>st</sup>-CASAM highlights the novel finding that response sensitivities to the imprecisely known domain boundaries and interfaces can arise both from the definition of the system’s response as well as from the equations, interfaces and boundary conditions that characterize the model and its imprecisely known domain. By enabling, in premiere, the exact computations of sensitivities to interface and boundary parameters and conditions, the 1<sup>st</sup>-CASAM enables the quantification of the effects of manufacturing tolerances on the responses of physical and engineering systems.展开更多
This work continues the illustrative application of the “Second Order Comprehensive Adjoint Sensitivity Analysis Methodology” (2<sup>nd</sup>-CASAM) to a benchmark mathematical model that can simulate th...This work continues the illustrative application of the “Second Order Comprehensive Adjoint Sensitivity Analysis Methodology” (2<sup>nd</sup>-CASAM) to a benchmark mathematical model that can simulate the evolution and/or transmission of particles in a heterogeneous medium. The model response considered in this work is a reaction-rate detector response, which provides the average interactions of particles with the respective detector or, alternatively, the time-average of the concentration of a mixture of substances in a medium. The definition of this model response includes both uncertain boundary points of the benchmark, thereby providing both direct and indirect contributions to the response sensitivities stemming from the boundaries. The exact expressions for the 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to the boundary and model parameters obtained in this work can serve as stringent benchmarks for inter-comparing the performances of all (deterministic and statistical) sensitivity analysis methods.展开更多
This work illustrates the application of the “Second Order Comprehensive Adjoint Sensitivity Analysis Methodology” (2<sup>nd</sup>-CASAM) to a mathematical model that can simulate the evolution and/or tr...This work illustrates the application of the “Second Order Comprehensive Adjoint Sensitivity Analysis Methodology” (2<sup>nd</sup>-CASAM) to a mathematical model that can simulate the evolution and/or transmission of particles in a heterogeneous medium. The model response is the value of the model’s state function (particle concentration or particle flux) at a point in phase-space, which would simulate a pointwise measurement of the respective state function. This paradigm model admits exact closed-form expressions for all of the 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to the model’s uncertain parameters and domain boundaries. These closed-form expressions can be used to verify the numerical results of production and/or commercial software, e.g., particle transport codes. Furthermore, this paradigm model comprises many uncertain parameters which have relative sensitivities of identical magnitudes. Therefore, this paradigm model could serve as a stringent benchmark for inter-comparing the performances of all deterministic and statistical sensitivity analysis methods, including the 2<sup>nd</sup>-CASAM.展开更多
A series of CdS nanoparticles with different surfaces were prepared by colloidal chemical method and reverse micelle method. Their second-order nonlinear optical (NLO) properties were experimentally studied in solutio...A series of CdS nanoparticles with different surfaces were prepared by colloidal chemical method and reverse micelle method. Their second-order nonlinear optical (NLO) properties were experimentally studied in solution by newly developed hyper-Rayleigh scattering (HRS) technique. The results show that 'per particle' first-order hyperpolarizability beta values are sensitive To the synthetic method and the surface chemical modification.展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
We develop a simple analytic calculation for the first order wave function of helium in a model in which nuclear charge screening is caused by repulsive coulomb interaction. The perturbation term, first-order correlat...We develop a simple analytic calculation for the first order wave function of helium in a model in which nuclear charge screening is caused by repulsive coulomb interaction. The perturbation term, first-order correlation energy, and first-order wave function are divided into two components, one component associated with the repulsive coulomb interaction and the other proportional to magnetic shielding. The resulting first-order wave functions are applied to calculate second-order energies within the model. We find that the second-order energies are independent of the nuclear charge screening constant in the unperturbed Hamiltonian with a central coulomb potential.展开更多
This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is...This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.展开更多
Photovoltaic(PV)inverter,as a promising voltage/var control(VVC)resource,can supply flexible reactive power to reduce microgrid power loss and regulate bus voltage.Meanwhile,active power plays a significant role in mi...Photovoltaic(PV)inverter,as a promising voltage/var control(VVC)resource,can supply flexible reactive power to reduce microgrid power loss and regulate bus voltage.Meanwhile,active power plays a significant role in microgrid voltage profile.Price-based demand response(PBDR)can shift load demand via determining time-varying prices,which can be regarded as an effective means for active power shifting.However,due to the different characteristics,PBDR and inverter-based VVC lack systematic coordination.Thus,this paper proposes a PBDR-supported three-stage hierarchically coordinated voltage control method,including day-ahead PBDR price scheduling,hour-ahead reactive power dispatch of PV inverters,and realtime local droop control of PV inverters.Considering their mutual influence,a stochastic optimization method is utilized to centrally or hierarchically coordinate adjacent two stages.To solve the bilinear constraints of droop control function,the problem is reformulated into a second-order cone programming relaxation model.Then,the concave constraints are convexified,forming a penalty convex-concave model for feasible solution recovery.Lastly,a convex-concave procedure-based solution algorithm is proposed to iteratively solve the penalty model.The proposed method is tested on 33-bus and IEEE 123-bus distribution networks and compared with other methods.The results verify the high efficiency of the proposed method to achieve power loss reduction and voltage regulation.展开更多
Based on the theory of semi-global C 1 solution and the local exact boundary controllability for first-order quasilinear hyperbolic systems,the local exact boundary controllability for a kind of second-order quasiline...Based on the theory of semi-global C 1 solution and the local exact boundary controllability for first-order quasilinear hyperbolic systems,the local exact boundary controllability for a kind of second-order quasilinear hyperbolic systems is obtained by a constructive method.展开更多
Reliability and optimization are two key elements for structural design. The reliability~ based topology optimization （RBTO） is a powerful and promising methodology for finding the optimum topologies with the uncert...Reliability and optimization are two key elements for structural design. The reliability~ based topology optimization （RBTO） is a powerful and promising methodology for finding the optimum topologies with the uncertainties being explicitly considered, typically manifested by the use of reliability constraints. Generally, a direct integration of reliability concept and topol- ogy optimization may lead to computational difficulties. In view of this fact, three methodologies have been presented in this study, including the double-loop approach （the performance measure approach, PMA） and the decoupled approaches （the so-called Hybrid method and the sequential optimization and reliability assessment, SORA）. For reliability analysis, the stochastic response surface method （SRSM） was applied, combining with the design of experiments generated by the sparse grid method, which has been proven as an effective and special discretization technique. The methodologies were investigated with three numerical examples considering the uncertainties including material properties and external loads. The optimal topologies obtained using the de- terministic, RBTOs were compared with one another; and useful conclusions regarding validity, accuracy and efficiency were drawn.展开更多
Based on the theory of semi-global classical solutions to quasilinear hyperbolic systems, the local exact boundary observability for a kind of second-order quasilinear hyperbolic systems is obtained by a constructive ...Based on the theory of semi-global classical solutions to quasilinear hyperbolic systems, the local exact boundary observability for a kind of second-order quasilinear hyperbolic systems is obtained by a constructive method.展开更多
This paper presents a reliability analysis of the pseudo-static seismic bearing capacity of a strip foundation using the limit equilibrium theory. The first-order reliability method(FORM) is employed to calculate the ...This paper presents a reliability analysis of the pseudo-static seismic bearing capacity of a strip foundation using the limit equilibrium theory. The first-order reliability method(FORM) is employed to calculate the reliability index. The response surface methodology(RSM) is used to assess the Hasofer e Lind reliability index and then it is optimized using a genetic algorithm(GA). The random variables used are the soil shear strength parameters and the seismic coefficients(khand kv). Two assumptions(normal and non-normal distribution) are used for the random variables. The assumption of uncorrelated variables was found to be conservative in comparison to that of negatively correlated soil shear strength parameters. The assumption of non-normal distribution for the random variables can induce a negative effect on the reliability index of the practical range of the seismic bearing capacity.展开更多
Ca_(3)Ti_(2)O_(7) with Ruddlesden-Popper structure exhibits the largest polarization among the known hybrid improper ferroelectrics.However,the high Curie temperature impedes the thorough study of phase transition thr...Ca_(3)Ti_(2)O_(7) with Ruddlesden-Popper structure exhibits the largest polarization among the known hybrid improper ferroelectrics.However,the high Curie temperature impedes the thorough study of phase transition through dielectric characterization.According to the previous theoretical design rule,the Curie temperature can be suppressed by increasing the tolerance factor.So,in the present work,high-quality Ca_(3-x)LaxTi_(2-x)Al_(x)O_(7)(x=0.0,0.1,0.2,0.3)ceramics with increased tolerance factors were successfully prepared.The amplitude of oxygen octahedral tilt mode indeed decreases with increasing tolerance factors,leading to a degeneration of ferroelectric polarization.However,the unexpected rigid rotation mode causes the composition-invariable coercive fields.The Curie temperatures decrease linearly with increasing x and tolerance factors.The variable-temperature dielectric constant confirms first-order improper ferroelectric transitions in Ca_(3)Ti_(2)O_(7)-based ceramics.The results of variable temperature Xray diffraction reveal the coexistence of two-phases below Curie temperature.The present work confidently confirms the first-order improper ferroelectric transition in Ca_(3)Ti_(2)O_(7)-based ceramics by combining results of variable-temperature dielectric response and in-situ X-ray powder diffraction.展开更多
In this study,the frequency characteristics of the turbulent wind and the effects of wind-wave coupling on the low-and high-frequency responses of semi-submersible floating offshore wind turbines(FOWT)are investigated...In this study,the frequency characteristics of the turbulent wind and the effects of wind-wave coupling on the low-and high-frequency responses of semi-submersible floating offshore wind turbines(FOWT)are investigated.Various wave load components,such as first-order wave loads,combined first-and second-order difference-frequency wave loads,combined first-and second-order sum-frequency wave loads,and first-and complete second-order wave loads are taken into consideration,while different turbulent environments are considered in aerodynamic loads.The com-parison is based on time histories and frequency spectra of platform motions and structural load responses and statistical values.The findings indicate that the second-order difference-frequency wave loads will significantly increase the natural frequency of low-frequency motion in the responses of the platform motion and structure load of the semi-submersible platform,which will cause structural fatigue damage.Under the action of turbulent wind,the influences of second-order wave loads on the platform motion and structural load response cannot be ignored,especially under extreme sea conditions.Therefore,in order to evaluate the dynamic responses of semi-submersible FOWT more accurately,the actual environment should be simulated more realistically.展开更多
In this paper,we discuss the numerical accuracy of asymptotic homogenization method(AHM)and multiscale finite element method(MsFEM)for periodic composite materials.Through numerical calculation of the model problems f...In this paper,we discuss the numerical accuracy of asymptotic homogenization method(AHM)and multiscale finite element method(MsFEM)for periodic composite materials.Through numerical calculation of the model problems for four kinds of typical periodic composite materials,the main factors to determine the accuracy of first-order AHM and second-order AHM are found,and the physical interpretation of these factors is given.Furthermore,the way to recover multiscale solutions of first-order AHM and MsFEM is theoretically analyzed,and it is found that first-order AHM and MsFEM provide similar multiscale solutions under some assumptions.Finally,numerical experiments verify that MsFEM is essentially a first-order multiscale method for periodic composite materials.展开更多
基金Supported by the National Natural Science Foundation of China (Grant Nos.52088102 and 51879287)National Key Research and Development Program of China (Grant No.2022YFB2602301)。
文摘Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.
基金The National Natural Science Foundation of China(No.61231002,61273266,61375028)the Ph.D.Programs Foundation of Ministry of Education of China(No.20110092130004)
文摘In order to improve the design results for the reconfigurable frequency response masking FRM filters an improved design method based on second-order cone programming SOCP is proposed.Unlike traditional methods that separately design the proposed method takes all the desired designing modes into consideration when designing all the subfilters. First an initial solution is obtained by separately designing the subfilters and then the initial solution is updated by iteratively solving a SOCP problem. The proposed method is evaluated on a design example and simulation results demonstrate that jointly designing all the subfilters can obtain significantly lower minimax approximation errors compared to the conventional design method.
基金supported by the National Natural Science Foundation of China (Grants 11002013, 11372025)the Defense Industrial Technology Development Program (Grants A0820132001, JCKY2013601B)+1 种基金the Aeronautical Science Foundation of China (Grant 2012ZA51010)111 Project (Grant B07009) for support
文摘In this paper, based on the second-order Taylor series expansion and the difference of convex functions algo- rithm for quadratic problems with box constraints （the DCA for QB）, a new method is proposed to solve the static response problem of structures with fairly large uncertainties in interval parameters. Although current methods are effective for solving the static response problem of structures with interval parameters with small uncertainties, these methods may fail to estimate the region of the static response of uncertain structures if the uncertainties in the parameters are fairly large. To resolve this problem, first, the general expression of the static response of structures in terms of structural parameters is derived based on the second-order Taylor series expansion. Then the problem of determining the bounds of the static response of uncertain structures is transformed into a series of quadratic problems with box constraints. These quadratic problems with box constraints can be solved using the DCA approach effectively. The numerical examples are given to illustrate the accuracy and the efficiency of the proposed method when comparing with other existing methods.
文摘This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.
文摘This work illustrates the application of the 1<sup>st</sup>-CASAM to a paradigm heat transport model which admits exact closed-form solutions. The closed-form expressions obtained in this work for the sensitivities of the temperature distributions within the model to the model’s parameters, internal interfaces and external boundaries can be used to benchmark commercial and production software packages for simulating heat transport. The 1<sup>st</sup>-CASAM highlights the novel finding that response sensitivities to the imprecisely known domain boundaries and interfaces can arise both from the definition of the system’s response as well as from the equations, interfaces and boundary conditions that characterize the model and its imprecisely known domain. By enabling, in premiere, the exact computations of sensitivities to interface and boundary parameters and conditions, the 1<sup>st</sup>-CASAM enables the quantification of the effects of manufacturing tolerances on the responses of physical and engineering systems.
文摘This work continues the illustrative application of the “Second Order Comprehensive Adjoint Sensitivity Analysis Methodology” (2<sup>nd</sup>-CASAM) to a benchmark mathematical model that can simulate the evolution and/or transmission of particles in a heterogeneous medium. The model response considered in this work is a reaction-rate detector response, which provides the average interactions of particles with the respective detector or, alternatively, the time-average of the concentration of a mixture of substances in a medium. The definition of this model response includes both uncertain boundary points of the benchmark, thereby providing both direct and indirect contributions to the response sensitivities stemming from the boundaries. The exact expressions for the 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to the boundary and model parameters obtained in this work can serve as stringent benchmarks for inter-comparing the performances of all (deterministic and statistical) sensitivity analysis methods.
文摘This work illustrates the application of the “Second Order Comprehensive Adjoint Sensitivity Analysis Methodology” (2<sup>nd</sup>-CASAM) to a mathematical model that can simulate the evolution and/or transmission of particles in a heterogeneous medium. The model response is the value of the model’s state function (particle concentration or particle flux) at a point in phase-space, which would simulate a pointwise measurement of the respective state function. This paradigm model admits exact closed-form expressions for all of the 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to the model’s uncertain parameters and domain boundaries. These closed-form expressions can be used to verify the numerical results of production and/or commercial software, e.g., particle transport codes. Furthermore, this paradigm model comprises many uncertain parameters which have relative sensitivities of identical magnitudes. Therefore, this paradigm model could serve as a stringent benchmark for inter-comparing the performances of all deterministic and statistical sensitivity analysis methods, including the 2<sup>nd</sup>-CASAM.
基金National Natural Science Foundation of China! (No.59582005)
文摘A series of CdS nanoparticles with different surfaces were prepared by colloidal chemical method and reverse micelle method. Their second-order nonlinear optical (NLO) properties were experimentally studied in solution by newly developed hyper-Rayleigh scattering (HRS) technique. The results show that 'per particle' first-order hyperpolarizability beta values are sensitive To the synthetic method and the surface chemical modification.
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
文摘In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
文摘We develop a simple analytic calculation for the first order wave function of helium in a model in which nuclear charge screening is caused by repulsive coulomb interaction. The perturbation term, first-order correlation energy, and first-order wave function are divided into two components, one component associated with the repulsive coulomb interaction and the other proportional to magnetic shielding. The resulting first-order wave functions are applied to calculate second-order energies within the model. We find that the second-order energies are independent of the nuclear charge screening constant in the unperturbed Hamiltonian with a central coulomb potential.
基金Project supported by the Executive Agency for Higher Education Research Development and Innovation Funding of Romania(No.PN-III-P4-PCE-2021-0993)。
文摘This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.
基金supported in part by the National Natural Science Foundation of China(No.52307091)in part by the Natural Science Foundation of Jiangsu Province(No.BK20230952)+1 种基金in part by the China Postdoctoral Science Foundation(No.2023M740976)in part by the Start Up Grant of City University of Hong Kong(No.9380163)。
文摘Photovoltaic(PV)inverter,as a promising voltage/var control(VVC)resource,can supply flexible reactive power to reduce microgrid power loss and regulate bus voltage.Meanwhile,active power plays a significant role in microgrid voltage profile.Price-based demand response(PBDR)can shift load demand via determining time-varying prices,which can be regarded as an effective means for active power shifting.However,due to the different characteristics,PBDR and inverter-based VVC lack systematic coordination.Thus,this paper proposes a PBDR-supported three-stage hierarchically coordinated voltage control method,including day-ahead PBDR price scheduling,hour-ahead reactive power dispatch of PV inverters,and realtime local droop control of PV inverters.Considering their mutual influence,a stochastic optimization method is utilized to centrally or hierarchically coordinate adjacent two stages.To solve the bilinear constraints of droop control function,the problem is reformulated into a second-order cone programming relaxation model.Then,the concave constraints are convexified,forming a penalty convex-concave model for feasible solution recovery.Lastly,a convex-concave procedure-based solution algorithm is proposed to iteratively solve the penalty model.The proposed method is tested on 33-bus and IEEE 123-bus distribution networks and compared with other methods.The results verify the high efficiency of the proposed method to achieve power loss reduction and voltage regulation.
基金supported by the Excellent Doctoral Research Foundation for Key Subject of Fudan University (No.EHH1411208)
文摘Based on the theory of semi-global C 1 solution and the local exact boundary controllability for first-order quasilinear hyperbolic systems,the local exact boundary controllability for a kind of second-order quasilinear hyperbolic systems is obtained by a constructive method.
基金Project supported by the National Natural Science Foundation of China(Nos.51275040 and 50905017)the Programme of Introducing Talents of Discipline to Universities(No.B12022)
文摘Reliability and optimization are two key elements for structural design. The reliability~ based topology optimization （RBTO） is a powerful and promising methodology for finding the optimum topologies with the uncertainties being explicitly considered, typically manifested by the use of reliability constraints. Generally, a direct integration of reliability concept and topol- ogy optimization may lead to computational difficulties. In view of this fact, three methodologies have been presented in this study, including the double-loop approach （the performance measure approach, PMA） and the decoupled approaches （the so-called Hybrid method and the sequential optimization and reliability assessment, SORA）. For reliability analysis, the stochastic response surface method （SRSM） was applied, combining with the design of experiments generated by the sparse grid method, which has been proven as an effective and special discretization technique. The methodologies were investigated with three numerical examples considering the uncertainties including material properties and external loads. The optimal topologies obtained using the de- terministic, RBTOs were compared with one another; and useful conclusions regarding validity, accuracy and efficiency were drawn.
基金supported by the National Natural Science Foundation of China(No.11526050)
文摘Based on the theory of semi-global classical solutions to quasilinear hyperbolic systems, the local exact boundary observability for a kind of second-order quasilinear hyperbolic systems is obtained by a constructive method.
基金the Ministry of Higher Education and Scientific Research of Algeria for supporting this work by offering an 11-month scholarship to the first author at the 3SR laboratory of Grenoble Alpes University,France
文摘This paper presents a reliability analysis of the pseudo-static seismic bearing capacity of a strip foundation using the limit equilibrium theory. The first-order reliability method(FORM) is employed to calculate the reliability index. The response surface methodology(RSM) is used to assess the Hasofer e Lind reliability index and then it is optimized using a genetic algorithm(GA). The random variables used are the soil shear strength parameters and the seismic coefficients(khand kv). Two assumptions(normal and non-normal distribution) are used for the random variables. The assumption of uncorrelated variables was found to be conservative in comparison to that of negatively correlated soil shear strength parameters. The assumption of non-normal distribution for the random variables can induce a negative effect on the reliability index of the practical range of the seismic bearing capacity.
基金financially supported by the National Natural Science Foundation of China under Grant Nos.51772266 and 51790493the National Key R&D Program of China under Grant No.2016YFA0300101.
文摘Ca_(3)Ti_(2)O_(7) with Ruddlesden-Popper structure exhibits the largest polarization among the known hybrid improper ferroelectrics.However,the high Curie temperature impedes the thorough study of phase transition through dielectric characterization.According to the previous theoretical design rule,the Curie temperature can be suppressed by increasing the tolerance factor.So,in the present work,high-quality Ca_(3-x)LaxTi_(2-x)Al_(x)O_(7)(x=0.0,0.1,0.2,0.3)ceramics with increased tolerance factors were successfully prepared.The amplitude of oxygen octahedral tilt mode indeed decreases with increasing tolerance factors,leading to a degeneration of ferroelectric polarization.However,the unexpected rigid rotation mode causes the composition-invariable coercive fields.The Curie temperatures decrease linearly with increasing x and tolerance factors.The variable-temperature dielectric constant confirms first-order improper ferroelectric transitions in Ca_(3)Ti_(2)O_(7)-based ceramics.The results of variable temperature Xray diffraction reveal the coexistence of two-phases below Curie temperature.The present work confidently confirms the first-order improper ferroelectric transition in Ca_(3)Ti_(2)O_(7)-based ceramics by combining results of variable-temperature dielectric response and in-situ X-ray powder diffraction.
基金supported by the Natural Science Foundation of Zhejiang Province(Grant No.LHZ21E090003)the National Nature Science Foundation of China(Grant No.52171279)+1 种基金Zhoushan Science&Technology Project(Grant No.2021C21002)supported by CNPq(Conselho Nacional de Desenvolvimento Científico e Tecnológico,Grant No.301474/2017-6).
文摘In this study,the frequency characteristics of the turbulent wind and the effects of wind-wave coupling on the low-and high-frequency responses of semi-submersible floating offshore wind turbines(FOWT)are investigated.Various wave load components,such as first-order wave loads,combined first-and second-order difference-frequency wave loads,combined first-and second-order sum-frequency wave loads,and first-and complete second-order wave loads are taken into consideration,while different turbulent environments are considered in aerodynamic loads.The com-parison is based on time histories and frequency spectra of platform motions and structural load responses and statistical values.The findings indicate that the second-order difference-frequency wave loads will significantly increase the natural frequency of low-frequency motion in the responses of the platform motion and structure load of the semi-submersible platform,which will cause structural fatigue damage.Under the action of turbulent wind,the influences of second-order wave loads on the platform motion and structural load response cannot be ignored,especially under extreme sea conditions.Therefore,in order to evaluate the dynamic responses of semi-submersible FOWT more accurately,the actual environment should be simulated more realistically.
基金the National Natural Science Foundation of China(No.11501449 and 11471262)the Center for high performance computing of Northwestern Polytechnical University.
文摘In this paper,we discuss the numerical accuracy of asymptotic homogenization method(AHM)and multiscale finite element method(MsFEM)for periodic composite materials.Through numerical calculation of the model problems for four kinds of typical periodic composite materials,the main factors to determine the accuracy of first-order AHM and second-order AHM are found,and the physical interpretation of these factors is given.Furthermore,the way to recover multiscale solutions of first-order AHM and MsFEM is theoretically analyzed,and it is found that first-order AHM and MsFEM provide similar multiscale solutions under some assumptions.Finally,numerical experiments verify that MsFEM is essentially a first-order multiscale method for periodic composite materials.
