Is it possible to intervene in the parameteric sweep of time-dependent solver?
Posted Apr 17, 2024, 3:16 a.m. EDT LiveLink for MATLAB, Studies & Solvers, Modeling Workflow 8 Replies
Please login with a confirmed email address before reporting spam
Hi there,
I want to model the evolution of structure configuration toward its steady state under various specific parameter combinations and given initial conditions.
The problem is that the period from the initial condition to the steady state is unknown before modeling. For now, I use the parameteric sweep of time-dependent solver and specify a shared pretty long time period Ttotal to each of the time-dependent study. After a trial run, I find that for some parameter combinations, the steady state has already been reached far before Ttotal, while for some other parameter combinations, the steady state has not been reached yet.
So, I was wondering if it is possible (maybe through scripts) to intervene in the parameteric sweep of time-dependent solver, so that when a customized criterion (such as energy convergence) is satisfied, comsol can stop the calculation of current parameter combination and start to the next combination, automatically.
Best, HC L
Please login with a confirmed email address before reporting spam
Check the options of the Events interface. You find it under Physics->Mathematics->ODE and DAE Interfaces.
Cheers Edgar
-------------------Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
Please login with a confirmed email address before reporting spam
Thank you @Edgar for your suggestion.
Indeed, the Events interface provides a useful tool to solve my problem.
After check the related knowledge base, presumably, I need to drop the parametric sweep node and only use a single time dependent study step, instead. During the time-dependent modeling, the parameters, such as boundary displacement, will be changed upon the trigger of events.
Following this idea, to employ the event interface, the variance of strain energy Us is denoted by Usvary. When Usvary is smaller than a given tolerance Ustol, Us is deemed to be constant. In other words, the system converges to its steady state, and the event should be triggered. Usflag is defined as Usvary-Ustol under the 'indicator states 1' node.The table below gives the values of variables near the sign change of Usvary-Ustol. It's shown that the value of Usvary-Ustol do changes its sign from >0 to <0, but Usflag doesn't change. Therefore, the event cannot be triggered. I try to assgin 1e-8 to the 'event tolerance' under 'Solver Configurations-->Time-Dependent Solver', but it doesn't work. Moreover, I change the definition of Usflag as Usvary/Ustol-1. Similary table can be given below. It seems Usvary always takes 0 under the 'indicator states 1' node when evaluate Usflag, leaving a constant Usflag. Therefore, the event cannot be triggered. Thus, I give it another shot and directly use the expression '(Usvary/Ustol-1)<0' as the condition under the 'implicit event 1' node, instead of 'Usflag<0'. It doesn't work, neither.
Besides, How to terminate the time-dependent solver when the final event is triggered? For example, the original swept parameter is the boundary displacement ub. Every time an event is triggered, ub will increase by a prescribed step, until that ub reaches its final value, ubf. After reaches steady state under ubf, the modeling should be terminated.
Thank you so much for your reading. I would appreciate it if you could give me some hints.
Best regard
HC. Liu
Check the options of the Events interface. You find it under Physics->Mathematics->ODE and DAE Interfaces.
Cheers Edgar
Please login with a confirmed email address before reporting spam
Check the options of the Events interface. You find it under Physics->Mathematics->ODE and DAE Interfaces.
Cheers Edgar
Thank you @Edgar for your advice.
I think event interface could be the solution. Some questions (posted in this thread) are encountered when using event interface.
I would appreciate it if you could give some more hints.
Best regards,
HC Liu
Please login with a confirmed email address before reporting spam
A possible simpler solution could be to add a Stop Condition node under Time-Dependent Solver. Then, you could use a stop criterion based on the time derivative of something, like
(t>1E-3)&&(abs(d(comp1.solid.Ws_tot,t))<200)
-------------------Henrik Sönnerlind
COMSOL
Please login with a confirmed email address before reporting spam
A possible simpler solution could be to add a Stop Condition node under Time-Dependent Solver. Then, you could use a stop criterion based on the time derivative of something, like
(t>1E-3)&&(abs(d(comp1.solid.Ws_tot,t))
Thank you so much @Henrik. 'Stop condition' exactly solves the problem.
Please login with a confirmed email address before reporting spam
A possible simpler solution could be to add a Stop Condition node under Time-Dependent Solver. Then, you could use a stop criterion based on the time derivative of something, like
(t>1E-3)&&(abs(d(comp1.solid.Ws_tot,t))
Hi Henrik,
Another question occurs. After the parametric sweep, I try to access the results of the last timestep under each parameter value through matlab commands like 'strain=mphevalpoint(model,{'ee'},"dataset","dset2",'outersolnum','all','selection',2,'t',Trelax)', where Trelax is the end step time at each parameter value.
'Stop Condition' makes Trelax an unknow variable. I try to use commands like 'strain=mphevalpoint(model,{'ee'},"dataset","dset2",'outersolnum','all','selection',2,'solnum','end')'. It causes error like
Invalid property value. - Property: solnum (Time)
'Time' is an array of integers. - : "1", "2", "3", "4"
Is there any way to circumvent?
Best regards,
HCL
Please login with a confirmed email address before reporting spam
Hi HCL
I can't find a model that reproduces your problem. Some of our wrapper functions do not work well with parametric sweeps where the number of returned values (e.g. times or solnum) vary accross the parametric sweep.
Try this: for i=1:number_of_outer_solnums strain(i)=mphevalpoint(model,{'ee'},"dataset","dset2",'outersolnum',i,'selection',2,'solnum','end') end
If that doesn't work for you send please send your model to support@comsol.com and mention my name.
-------------------Lars Gregersen
Comsol Denmark
Please login with a confirmed email address before reporting spam
Hi HCL
I can't find a model that reproduces your problem. Some of our wrapper functions do not work well with parametric sweeps where the number of returned values (e.g. times or solnum) vary accross the parametric sweep.
Try this: for i=1:number_of_outer_solnums strain(i)=mphevalpoint(model,{'ee'},"dataset","dset2",'outersolnum',i,'selection',2,'solnum','end') end
If that doesn't work for you send please send your model to support@comsol.com and mention my name.
Thank you @Lars for your advice! It works!
Best regards,
HCL
Reply
Please read the discussion forum rules before posting.
Discussion Forum Rules and Guidelines
The goal of COMSOL Access is to provide a forum for you to communicate effectively with COMSOL as well as your colleagues within the multiphysics simulation community. This involves providing you with access to technical support and downloads of the latest {:comsol} software releases, as well as the ability to share your comments and work with other users of the {:comsol} software through forums such as the blog, discussion forum, and Application Exchange. In order to make this an efficient and pleasant experience for you and other members of COMSOL Access, we ask that you follow a few rules and guidelines.
When registering for COMSOL Access, you agree to provide your complete and truthful information for all fields requested on your COMSOL Access account registration page. You also agree to maintain the accuracy of all information associated with you on your COMSOL Access account. You agree to maintain your COMSOL Access account for use solely by you, not to share your username and password with anyone else, and to take appropriate precautions to restrict access to your username and password from others. Furthermore, you agree not to submit any information relating to your employer through your COMSOL Access account without your employer’s authorization. Should you use a COMSOL Access account associated with an employer, you agree to immediately discontinue using that account upon termination of that employment.
The moderators of the forums will remove any generally objectionable material as quickly as possible. You acknowledge that all posts made to these forums express the views and opinions of the author and not the administrators, moderators, or webmaster (except for posts by these people). Hence, they will not be held liable.
You agree not to post or link to any material that is abusive, obscene, vulgar, slanderous, hateful, threatening, sexually oriented, or that infringes upon or violates any third-party rights or any other material that may violate any applicable laws. You agree that you will not otherwise use your COMSOL Access account to violate or to assist anyone in violating any law. Engaging in any activity in violation of these COMSOL Access rules and guidelines may lead to you being immediately and permanently banned from COMSOL Access. The IP address of all posts is recorded to aid in enforcing these conditions. You agree that the webmaster, administrator, and moderators of the forums have the right to remove, move, or close any topic at any time as they see fit. As a user of COMSOL Access, you agree to any information you have entered into any of the forums being stored in a database.
You agree that you will not use your COMSOL Access account in violation of any applicable export control laws. You represent and warrant that you are not subject to any comprehensive sanction or embargo by the U.S. or any other country, nor are you identified on any list maintained by the U.S. government that identifies persons for which the U.S. maintains restrictions. Further, you represent and warrant that you are not subject to any restriction on the receipt of technology or products under the export control laws of the U.S. or any other country.
Basic Rules
- Flaming: Do not post any messages that harass, insult, belittle, threaten, or flame another member or guest. Debates are fine, but argue with the point, not the person.
- Trolling: Do not post with the purpose of starting a dispute. Note that a person disagreeing with your opinion is not considered trolling; keep it civil, even if you are sure that the other person is wrong. Anything seen as trolling will result in you being banned.
- Spamming: Posts without content or containing nonsense waste space and everyone's time. Spam will be removed from the forum.
- Offensive Posts, Links or Images: Do not use profanity, racial, ethnic, religious, or other slurs or any other offensive material.
- Advertising: Posting of advertisements for products or services, links to auctions, affiliate links, links to promote websites, and so forth is not allowed.
- Disclosure of Personal Information: Do not disclose any other member's email, real name, address, phone number, IP address, or other personal information. This includes posting contents of emails and private messages without the sender's consent. Do not bring personal disputes into the forums. Bans and warnings issued to other users are considered personal information.
- Misleading Titles: The subject line of the post should be as informative as possible about the content of the post. This also reduces the amount of duplicate posts.
- Thread Hijacking: Taking a thread off-topic to pursue one's own agenda is not permitted.
- Piracy: Do not upload, post, email, transmit, or otherwise make available any content that infringes upon any patent, trademark, trade secret, copyright, or other proprietary rights ("Rights") of any party. This includes COMSOL model files created by a forged or illegitimate license.
- Honesty: Users must provide truthful information in creating their COMSOL Access account.
These rules are subject to change. The moderators reserve the right to remove, edit, or move posts at their discretion. The COMSOL Access administrators will reserve the right to permanently remove a user account without notice if any of the rules are not followed. Particular services accessible with your COMSOL Access account may be subject to additional rules. You agree to comply with all rules applicable to each service you access through your COMSOL Access account.
Posting Guidelines
When posting, understand that you are trying to communicate with other people. Although many COMSOL Access members are not fluent in English, the official language of this forum is English.
Here are some important guidelines of language:
- Write in English.
- If you are familiar with LaTeX, please use this to write mathematical equations.
- Always do a quick check for spelling/grammar mistakes.
- Format your post in a legible manner. Use the Preview button often.
- Be concise and articulate as much as possible.
- Use the Enter key to create paragraphs.
- Capitalize correctly. It is difficult to read posts that are written entirely in uppercase or lowercase.
- Use correct punctuation. Avoid run-on sentences.
- Try to avoid using “text speak”, “net speak”, or slang. The purpose of language is to be understood.
- Never invent acronyms and use as few acronyms as possible. For example, write "COMSOL Multiphysics" and not "CMP".
- Review your post before publishing it. Try to catch typos.
- Please check to see if a topic has already been posted. Do not post multiple threads on the same topic.
Disclaimer
By submitting content to the forums, you hereby grant COMSOL a nonexclusive, royalty-free, perpetual, worldwide, and unrestricted license to reproduce, publicly display, publicly distribute, and prepare derivative works of the content. COMSOL hereby grants you a license to copy and/or use content from the forums solely for your own internal purposes. COMSOL provides the forum service for the benefit of our users to share content with the community. All content is provided "as is" without warranty of any kind, express or implied, including without limitation, warranties of merchantability, noninfringement, design, operation, and fitness for a particular purpose, and the entire risk as to the quality and performance of the programs is with you.
Neither COMSOL, the authors, nor the copyright owners of submitted materials warrant that the programs will be error-free, uninterrupted, virus-free, secure, and suitable for your needs, produce specific results, or that errors or failures will be corrected. Comments on supplied content should be sent to the author or copyright owner through the tools provided in the forums.
Please log in to post a reply.
Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.
Suggested Content
- KNOWLEDGE BASE What hardware do you recommend for COMSOL Multiphysics®?
- FORUM Time-dependent problem: noisy temperature boundary condition. Is it possible to change the noise magnitude in a sweep analysis?
- BLOG Simulating Radiation Effects in Semiconductor Devices
- BLOG RF Signal Transformation Between the Time and Frequency Domains
- KNOWLEDGE BASE Running parametric sweeps, batch sweeps, and cluster sweeps from the command line